Values in raster change after CookieCutter function

[rickemery]

When the CookieCutter function is used to cut out a feature from a shapefile and raster with the same projection, the values in the raster are not maintained.

An example can be found in ftp://volga.ags.io/files/zs_large_raster.zip using zs_large_raster.tif and zs_large_raster.shp:

1. Download zs_large_raster.zip and extract it's contents
2. Open zs_large_raster.tif, examine values for pixels in FID 0
3. Run the following:

inimg = GeoImage('/path/to/zs_large_raster.tif')
feature = GeoVector('/path/to/zs_large_raster.shp')[0]
CookieCutter(
    GeoImages(['/path/to/zs_large_raster.tif']),
    feature,
    './cut_raster.tif',
    abs(inimg.Affine()[1]),
    abs(inimg.Affine()[5]),
    True
)

4. Open cut_raster.tif, examine values for pixels in FID 0