Opened 13 years ago

Closed 12 years ago

#1581 closed enhancement (fixed)

[raster] ST_Clip should not impose a nodata value when the source raster do not have one

Reported by: pracine Owned by: pracine
Priority: medium Milestone: PostGIS 2.0.2
Component: raster Version: master
Keywords: history Cc:

Description

In the case where:

1) the input raster do not have a nodata value defined

2) no alternative nodata values are provided

The behavior right now is to set a nodata value to the resulting raster using ST_MinPossibleValue(ST_BandPixelType(rast)). We should not do that...

Most users would expect to get a raster with no nodata value defined since this is what they provided. However in the case the geometry is not a rectangle aligned with the raster (e.g a circle), there must be a way to set some pixels to nodata.

When the provided raster have a nodata defined, this nodata value is used. No problem.

When an array of nodata value is provided, those nodata values are used. No problem.

There is no such need in the case where the geometry is a rectangular fraction of the raster extent or when the geometry encompasses the raster and hence there is no need to set a nodata value in the resulting raster but we have no way to know that.

When both of those conditions are not met (the provided raster do not have a nodata defined AND no array of nodata value is provided) and hence absolutely no nodata value is known, instead of imposing ST_MinPossibleValue() we should just NOTICE the user that we could not set a nodata value and that if he wants one to be defined he must use one of the variant accepting nodata values (or set a nodata value to its raster before clipping).

Change History (5)

comment:1 by pracine, 13 years ago

Note that this is not the case for ST_Intersection(raster, raster) since in this case we produce a new raster and we do not expect the resulting raster to have exact same characteristics as the input rasters. This is the case for ST_Clip: the user expect to get a raster similar to the one provided except for its extent (including nodata areas).

comment:2 by pracine, 13 years ago

Milestone: PostGIS Raster 2.0.1PostGIS 2.0.1

comment:3 by Bborie Park, 12 years ago

Milestone: PostGIS 2.0.1PostGIS 2.0.2

comment:4 by Bborie Park, 12 years ago

Fixed in -trunk (2.1) as of r10568.

comment:5 by Bborie Park, 12 years ago

Keywords: history added
Resolution: fixed
Status: newclosed

Fixed in 2.0 as of r10570

Note: See TracTickets for help on using tickets.