package mil.nga.giat.geowave.adapter.vector.plugin; import org.geotools.data.Query; import org.geotools.data.simple.SimpleFeatureCollection; import org.geotools.factory.Hints; import org.geotools.geometry.jts.ReferencedEnvelope; import org.geotools.process.ProcessException; import org.geotools.process.factory.DescribeParameter; import org.geotools.process.factory.DescribeProcess; import org.geotools.process.factory.DescribeResult; import org.geotools.process.gs.GSProcess; import org.opengis.coverage.grid.GridGeometry; /** * This class can be used as a GeoTools Render Transform ('geowave:Subsample') * within an SLD on any layer that uses the GeoWave Data Store. An example SLD * is provided (example-slds/SubsamplePoints.sld). The pixel-size allows you to * skip more than a single pixel. For example, a pixel size of 3 would skip an * estimated 3x3 pixel cell in GeoWave's row IDs. Note that rows are only * skipped when a feature successfully passes filters. * */ @DescribeProcess(title = "SubsampleAtScreenResolution", description = "This process will enable GeoWave to subsample WMS requests based on pixel resolution to not oversample data. This will efficiently render overlapping point geometry that would otherwise be hidden but it assumes an opaque style and does not take transparency into account. It will use the centroid for other geometry types than point which can produce visual artifacts - distributed rendering is an alternative approach to efficiently render lines and polygons") public class SubsampleProcess { public static final Hints.Key PIXEL_SIZE = new Hints.Key( Double.class); public static final Hints.Key OUTPUT_BBOX = new Hints.Key( ReferencedEnvelope.class); public static final Hints.Key OUTPUT_WIDTH = new Hints.Key( Integer.class); public static final Hints.Key OUTPUT_HEIGHT = new Hints.Key( Integer.class); @DescribeResult(name = "result", description = "This is just a pass-through, the key is to provide enough information within invertQuery to perform a map to screen transform") public SimpleFeatureCollection execute( @DescribeParameter(name = "data", description = "Feature collection containing the data") final SimpleFeatureCollection features, @DescribeParameter(name = "outputBBOX", description = "Georeferenced bounding box of the output") final ReferencedEnvelope argOutputEnv, @DescribeParameter(name = "outputWidth", description = "Width of the output raster") final Integer argOutputWidth, @DescribeParameter(name = "outputHeight", description = "Height of the output raster") final Integer argOutputHeight, @DescribeParameter(name = "pixelSize", description = "The pixel size to base subsampling on") final Double pixelSize ) throws ProcessException { // vector-to-vector render transform that is just a pass through - the // key is to add map to screen transform within invertQuery return features; } public Query invertQuery( @DescribeParameter(name = "outputBBOX", description = "Georeferenced bounding box of the output") final ReferencedEnvelope argOutputEnv, @DescribeParameter(name = "outputWidth", description = "Width of the output raster") final Integer argOutputWidth, @DescribeParameter(name = "outputHeight", description = "Height of the output raster") final Integer argOutputHeight, @DescribeParameter(name = "pixelSize", description = "The pixel size to base subsampling on") final Double pixelSize, final Query targetQuery, final GridGeometry targetGridGeometry ) throws ProcessException { // add to the query hints targetQuery.getHints().put( OUTPUT_WIDTH, argOutputWidth); targetQuery.getHints().put( OUTPUT_HEIGHT, argOutputHeight); targetQuery.getHints().put( OUTPUT_BBOX, argOutputEnv); if (pixelSize != null) { targetQuery.getHints().put( PIXEL_SIZE, pixelSize); } return targetQuery; } }