[raster] Respect imprecise real valued nodata
|Reported by:||mlt||Owned by:||Bborie Park|
PostGIS should handle well real valued, e.g., fractional NODATA values, e.g. with ST_DumpAsPolygons. I suspect rounding takes place somewhere in the code as the following query from my testing database
select (gv).val, ST_Centroid((gv).geom) geom, (gv).val = nd from ( select ST_DumpAsPolygons(ST_Union(ST_Clip(rast, geom)), 1, True) gv, ST_BandNoDataValue(rast,1) nd from dem, depressions where ST_Intersects(rast, geom) and gid=1 group by gid, nd ) foo --where (gv).val!=(select ST_BandNoDataValue(rast,1) from dem where rid=1) order by (gv).val limit 1
whereas I'd either expect False for direct comparison with NODATA due to imprecision/rounding, or different value in the first column if PostGIS properly eliminated NODATA from the output.
on Windows 7 64 bit.