Opened 9 years ago

Last modified 7 years ago

#3318 assigned defect

ST_MapAlgebra with mask seems to hard-code distance 1 instead of inferring distance from mask matrix

Reported by: robe Owned by: dustymugs
Priority: medium Milestone: PostGIS Fund Me
Component: raster Version: master
Keywords: Cc:

Description (last modified by robe)

I've moved this from the now closed mask documentation #2603 to here, I'm still having issue with this even as of r14198 even though dusty has made some changes.

DROP TABLE IF EXISTS shapes;
CREATE TABLE shapes(rid integer, rast raster);  

INSERT INTO shapes(rid,rast)
VALUES ( 1, ST_AsRaster(
		ST_Buffer(
			ST_GeomFromText('LINESTRING(50 50,150 150,150 50)'), 10,'join=bevel'), 
			200,200,ARRAY['8BUI'], ARRAY[118], ARRAY[0]) );

 SELECT st_mapalgebra(rast,1,'st_mean4ma(double precision[], int[], text[])'::regprocedure,'{{1,0,1,1,1}, {1,0,1,1,1}, {1,0,1,1,1}, {1,0,1,1,1}, {1,0,1,1,1}}'::double precision[],false)
FROM shapes;

Gives error:

ERROR:  rt_pixel_set_array: mask dimensions 5 x 5 do not match given dims 3 x 3

It's my understanding based on dusty's comments that the neighborhood size should be inferred from the mask size, so there seems to me something wrong here.

Bborie's comments copied from other ticket


Distance is two pixels left, right, up, down from the pixel of interest…

2 2 2 2 2
2 1 1 1 2
2 1 0 1 2
2 1 1 1 2
2 2 2 2 2

So the distance x/y are correct. As for that error message, I'll need to dig.


Change History (13)

comment:1 by robe, 9 years ago

Description: modified (diff)

comment:2 by robe, 9 years ago

Component: postgisraster
Owner: changed from pramsey to dustymugs

comment:3 by robe, 9 years ago

Milestone: PostGIS 2.3.0PostGIS 2.2.1

comment:4 by dustymugs, 9 years ago

Status: newassigned

comment:5 by dustymugs, 9 years ago

Priority: mediumcritical

comment:6 by robe, 9 years ago

Milestone: PostGIS 2.2.1PostGIS 2.3.0

Bborie,

I'm guessing you won't have time to look at this within the next day or so, so pushing to 2.3. If you do feel free to push back. I'll push to 2.2.2 after release.

comment:7 by dustymugs, 9 years ago

Wow I forgot about this… I'll take a look tonight.

comment:8 by robe, 9 years ago

Milestone: PostGIS 2.3.0PostGIS 2.2.1

okay I'll switch back then. Hurry Hurry PostgreSQL 9.5.0 is already been stamped.

comment:9 by robe, 9 years ago

Milestone: PostGIS 2.2.1PostGIS 2.2.2

comment:10 by pramsey, 9 years ago

Milestone: PostGIS 2.2.2PostGIS 2.2.3

comment:11 by robe, 8 years ago

Milestone: PostGIS 2.2.3PostGIS 2.4.0

comment:12 by robe, 7 years ago

Milestone: PostGIS 2.4.0PostGIS Future
Priority: criticalmedium

comment:13 by robe, 7 years ago

Milestone: PostGIS FuturePostGIS Fund Me

Milestone renamed

Note: See TracTickets for help on using tickets.