Skip to main content

I have an area of interest (a polygon) and I want to read only the raster data that is present in it. I am trying to use a FeatureReader (Reader=GEODATABASE_FILE_RASTER_DATASET) and set the spatial filter to "INTERSECTS". I know this option works well with vector data. But here it doesn't work and all the raster data is read.

It will not clip an area out of a raster.

Intersect is a relational query wich will result in reading all rasters crossing initiator boundary.

When selecting a raster by using intersect trough a feature reader, you probably still wil want to reduce it to the BB of your initiator. Wich means clipping. (creating a smaller rectangular object..)

Wich is wastefull; select raster by extracting its rasterextent (RasterExtenCoercer), then intersect wiht your area of interest. Then select raster by attribute (name) and clipp it. Much faster.

Selecting by feature reader is slow (for ecw's at least)


It will not clip an area out of a raster.

Intersect is a relational query wich will result in reading all rasters crossing initiator boundary.

When selecting a raster by using intersect trough a feature reader, you probably still wil want to reduce it to the BB of your initiator. Wich means clipping. (creating a smaller rectangular object..)

Wich is wastefull; select raster by extracting its rasterextent (RasterExtenCoercer), then intersect wiht your area of interest. Then select raster by attribute (name) and clipp it. Much faster.

Selecting by feature reader is slow (for ecw's at least)

Thanks for you answer!

Reply