Examples of org.geotools.coverage.grid.GridCoverage2D

org.geotools.coverage.grid.GridCoverage2D
Basic access to grid data values backed by a two-dimensional {@linkplain RenderedImage rendered image}. Each band in an image is represented as a {@linkplain GridSampleDimension sample dimension}.
Grid coverages are usually two-dimensional. However, {@linkplain #getEnvelope their envelope}may have more than two dimensions. For example, a remote sensing image may be valid only over some time range (the time of satellite pass over the observed area). Envelopes for such grid coverage can have three dimensions: the two usual ones (horizontal extent along x and y), and a third one for start time and end time (time extent along t). However, the {@linkplain GeneralGridRange grid range} for all extra-dimension musthave a {@linkplain GeneralGridRange#getLength size} not greater than 1. In other words, a{@code GridCoverage2D} can be a slice in a 3 dimensional grid coverage. Each slice can have anarbitrary width and height (like any two-dimensional images), but only 1 voxel depth (a "voxel" is a three-dimensional pixel).
Serialization note:
Because it is serializable, {@code GridCoverage2D} can be included as method argument or asreturn type in Remote Method Invocation (RMI). However, the pixel data are not sent during serialization. Instead, the image data are transmitted "on-demand" using socket communications. This mechanism is implemented using JAI {@link SerializableRenderedImage}class. While serialization (usually on server side) should work on J2SE 1.4 and above, deserialization (usually on client side) of {@code GridCoverage2D} instances requires J2SE 1.5. @since 2.1 @source $URL$ @version $Id$ @author Martin Desruisseaux (IRD)


        // getting a reader
        GeoTiffReader reader = new GeoTiffReader(input);


        // reading the coverage
        GridCoverage2D coverage = (GridCoverage2D) reader.read(null);


        // check coverage and crs
        assertNotNull(coverage);
        assertNotNull(coverage.getCoordinateReferenceSystem());
        reader.dispose();


        // get a writer
        final File writedir = new File("./target", "testWriter");
        writedir.mkdirs();
        final File output = new File(writedir, "target.tif");
        GeoTiffWriter writer = new GeoTiffWriter(output);




        ParameterValue<Boolean> value = GeoTiffFormat.RETAIN_AXES_ORDER.createValue();
        value.setValue(true);
        
        // switching axes
        final CoordinateReferenceSystem latLon4267=CRS.decode("EPSG:4267");
        assertEquals(CRS.getAxisOrder(latLon4267), AxisOrder.NORTH_EAST);
        final GeneralEnvelope envelope = (GeneralEnvelope) CRS.transform(coverage.getEnvelope(), latLon4267);
        envelope.setCoordinateReferenceSystem(latLon4267);
        
        coverage=CoverageFactoryFinder.getGridCoverageFactory(GeoTools.getDefaultHints()).create(
                coverage.getName(),
                coverage.getRenderedImage(),
                envelope
                );
        
        
        
        writer.write(coverage, new GeneralParameterValue[] { value });
        writer.dispose();
        coverage.dispose(true);


        // getting a reader
        reader = new GeoTiffReader(output,null);// this way I do not impose the lonlat ordering
        final GridCoverage2D gc= reader.read(null);
        final MathTransform g2w_ = gc.getGridGeometry().getGridToCRS();
        assertTrue(g2w_ instanceof AffineTransform2D);
        AffineTransform2D g2w = (AffineTransform2D) g2w_;
        assertTrue(XAffineTransform.getSwapXY(g2w)==-1);
        assertEquals(AxisOrder.NORTH_EAST,CRS.getAxisOrder(gc.getCoordinateReferenceSystem()));
        RenderedImage ri = gc.getRenderedImage();
        assertEquals(ri.getWidth(), 120 );
        assertEquals(ri.getHeight(), 121 );
        assertTrue(((GeneralEnvelope)gc.getEnvelope()).equals(coverage.getEnvelope(),XAffineTransform.getScaleX0(g2w)*1E-1,false));
        reader.dispose();


       
    }

View Full Code Here

        
                // getting a reader
                GeoTiffReader reader = new GeoTiffReader(input);
        
                // reading the coverage
                GridCoverage2D coverage = (GridCoverage2D) reader.read(null);
        
                // check coverage and crs
                assertNotNull(coverage);
                assertNotNull(coverage.getCoordinateReferenceSystem());
                reader.dispose();
        
                // get a writer
                final File output = new File(TestData.file(GeoTiffReaderTest.class, "."), "bigtiff" + i + ".tif");
                GeoTiffWriter writer = new GeoTiffWriter(output);
                
                GeoTiffWriteParams params = new GeoTiffWriteParams();
                params.setForceToBigTIFF(true);
                ParameterValue<GeoToolsWriteParams> value = GeoTiffFormat.GEOTOOLS_WRITE_PARAMS.createValue();
                value.setValue(params);
                
                writer.write(coverage, new GeneralParameterValue[]{value});
                writer.dispose();
                coverage.dispose(true);
                
                // getting a reader
                reader = new GeoTiffReader(output);
                RenderedImage ri = reader.read(null).getRenderedImage();
                assertEquals(ri.getWidth(), i == 0 ? 120 : 12);

View Full Code Here

        
        GeneralEnvelope coverageEnvelope2D = new GeneralEnvelope(gg2D.getEnvelope());
        GridSampleDimension[] sampleDimensions = new GridSampleDimension[]{new GridSampleDimension("testDim")};
        boolean getGeophysics = false;
        // Creation of the Coverage
        GridCoverage2D coverage1 = (GridCoverage2D) Utilities.createCoverageFromImage(coverageFactory , coverageName , imageIndex, image ,
                raster2Model , spatialReferenceSystem2D, null, sampleDimensions,
                getGeophysics);
        
        GridCoverage2D coverage2 = (GridCoverage2D) Utilities.createCoverageFromImage(coverageFactory , coverageName , imageIndex, image ,
                null , null, coverageEnvelope2D, sampleDimensions,
                getGeophysics);
        
        // Ensure the two coverages have the same envelope and raster dimensions
        assertTrue(coverage1.getEnvelope2D().boundsEquals(coverage2.getEnvelope2D(), 0, 1, 0.01d));
        assertSame(coverage1.getRenderedImage(), coverage2.getRenderedImage());
    }

View Full Code Here

          CoordinateReferenceSystem crs=reader.getCoordinateReferenceSystem();
          
          final File prj= TestData.file(GeoTiffReaderTest.class, "override/sample.prj");
          final CoordinateReferenceSystem crs_=new PrjFileReader(new FileInputStream(prj).getChannel()).getCoordinateReferenceSystem();
          assertTrue(CRS.equalsIgnoreMetadata(crs, crs_));
          GridCoverage2D coverage=reader.read(null);
          assertTrue(CRS.equalsIgnoreMetadata(coverage.getCoordinateReferenceSystem(), crs_));


          coverage.dispose(true);
      }

View Full Code Here

                                            ((GeoTiffReader) reader).getMetadata()
                                                            .getRootNode());
                            System.out.println(iIOMetadataDumper.getMetadata());            
            }
            // reading the coverage
            GridCoverage2D coverage1 = (GridCoverage2D) reader.read(null);


            // check coverage and crs
            assertNotNull(coverage1);
            assertNotNull(coverage1.getCoordinateReferenceSystem());
            assertNotSame(coverage1.getCoordinateReferenceSystem(),crs_);
            reader.dispose();


            coverage1.dispose(true);
        }

View Full Code Here

        final AbstractGridFormat format = new GeoTiffFormat();
        assertTrue(format.accepts(noCrs));
        GeoTiffReader reader = (GeoTiffReader) format.getReader(noCrs);
        CoordinateReferenceSystem crs=reader.getCoordinateReferenceSystem();
        assertTrue(CRS.equalsIgnoreMetadata(crs, DefaultEngineeringCRS.GENERIC_2D));
        GridCoverage2D coverage=reader.read(null);
        assertTrue(CRS.equalsIgnoreMetadata(coverage.getCoordinateReferenceSystem(), DefaultEngineeringCRS.GENERIC_2D));
        


        // hint for CRS
        crs = CRS.decode("EPSG:32632", true);
        final Hints hint = new Hints();
        hint.put(Hints.DEFAULT_COORDINATE_REFERENCE_SYSTEM, crs);


        // getting a reader
        reader = new GeoTiffReader(noCrs, hint);


        // reading the coverage
        GridCoverage2D coverage1 = (GridCoverage2D) reader.read(null);


        // check coverage and crs
        assertNotNull(coverage1);
        assertNotNull(coverage1.getCoordinateReferenceSystem());
        assertEquals(CRS.lookupIdentifier(coverage1.getCoordinateReferenceSystem(), true),
                "EPSG:32632");
        reader.dispose();


        //
        // use prj and wld
        //
        final File wldprjFile = TestData.file(GeoTiffReaderTest.class, "no_crs_no_envelope.tif");
        assertTrue(format.accepts(wldprjFile));


        // getting a reader
        reader = new GeoTiffReader(wldprjFile);


        // reading the coverage
        GridCoverage2D coverage2 = (GridCoverage2D) reader.read(null);


        // check coverage and crs
        assertNotNull(coverage2);
        assertNotNull(coverage2.getCoordinateReferenceSystem());
        assertEquals(CRS.lookupIdentifier(coverage2.getCoordinateReferenceSystem(), true),
                "EPSG:32632");
        reader.dispose();


        //
        // use prj and hint
        //
        final File wldFile = TestData.file(GeoTiffReaderTest.class, "no_crs_no_envelope2.tif");
        assertTrue(format.accepts(wldFile));


        // getting a reader
        reader = new GeoTiffReader(wldFile, hint);


        // reading the coverage
        GridCoverage2D coverage3 = (GridCoverage2D) reader.read(null);


        // check coverage and crs
        assertNotNull(coverage3);
        assertNotNull(coverage3.getCoordinateReferenceSystem());
        assertEquals(CRS.lookupIdentifier(coverage3.getCoordinateReferenceSystem(), true),
                "EPSG:32632");
        reader.dispose();
        
        coverage1.dispose(true);
        coverage2.dispose(true);
        coverage3.dispose(true);
    }

View Full Code Here

  }


  static void checkCoverage(final ImageMosaicReader reader,
      GeneralParameterValue[] values, String title, Rectangle rect) throws IOException {
    // Test the coverage
    final GridCoverage2D coverage = getCoverage(reader, values, true);
    testCoverage(reader, values, title, coverage, rect);
  }

View Full Code Here

    testCoverage(reader, values, title, coverage, rect);
  }


  static GridCoverage2D getCoverage(final ImageMosaicReader reader,
      GeneralParameterValue[] values, final boolean checkForNull) throws IOException {
    final GridCoverage2D coverage = (GridCoverage2D) reader.read(values);
    if (checkForNull) {
        Assert.assertNotNull(coverage);
    }
    return coverage;
  }

View Full Code Here

          // getting a reader
          GeoTiffReader reader = new GeoTiffReader(o, new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE));
          if (reader != null) {
    
            // reading the coverage
                  final GridCoverage2D coverage = (GridCoverage2D) reader.read(null);
    
            // Crs and envelope
            if (TestData.isInteractiveTest()) {
              buffer.append("CRS: ").append(
                  coverage.getCoordinateReferenceSystem2D()
                      .toWKT()).append("\n");
              buffer.append("GG: ").append(
                  coverage.getGridGeometry().toString()).append("\n");
            }
            // display metadata
            if (org.geotools.TestData.isExtensiveTest()) {
              IIOMetadataDumper iIOMetadataDumper = new IIOMetadataDumper(
                  ((GeoTiffReader) reader).getMetadata()
                      .getRootNode());
              buffer.append("TIFF metadata: ").append(
                  iIOMetadataDumper.getMetadata()).append("\n");
            }
            
            // layout checks
            final ImageLayout layout = reader.getImageLayout();
            assertNotNull(layout);
            assertNotNull(layout.getColorModel(null));
            assertNotNull(layout.getSampleModel(null));
            assertEquals(0,layout.getMinX(null));
            assertEquals(0,layout.getMinY(null));
            assertTrue(layout.getWidth(null) > 0);
            assertTrue(layout.getHeight(null) > 0);
            assertEquals(0,layout.getTileGridXOffset(null));
            assertEquals(0,layout.getTileGridYOffset(null));
            assertTrue(layout.getTileHeight(null) > 0);
            assertTrue(layout.getTileWidth(null) > 0);
            
            
            // showing it
            if (TestData.isInteractiveTest()){
                coverage.show();
            }
            else {
                PlanarImage.wrapRenderedImage(coverage.getRenderedImage()).getTiles();
            }
            
            // write and read back
            final File destFile = File.createTempFile("test", ".tif",writeDirectory);        
            final GeoTiffWriter writer= new GeoTiffWriter(destFile);
            writer.write(coverage, null);
            writer.dispose();
            
            // read back
            assertTrue(format.accepts(destFile));
            reader = new GeoTiffReader(destFile, new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE));
            final GridCoverage2D destCoverage = (GridCoverage2D) reader.read(null);
            reader.dispose();
            
            final double eps=XAffineTransform.getScaleX0((AffineTransform)coverage.getGridGeometry().getGridToCRS())*1E-2;
            assertTrue("CRS comparison failed:" +o.toString(),CRS.findMathTransform(coverage.getCoordinateReferenceSystem(), destCoverage.getCoordinateReferenceSystem(), true).isIdentity());
            assertTrue("CRS comparison failed:" +o.toString(),CRS.equalsIgnoreMetadata(coverage.getCoordinateReferenceSystem(), destCoverage.getCoordinateReferenceSystem()));
            assertTrue("GridRange comparison failed:" +o.toString(),coverage.getGridGeometry().getGridRange().equals(destCoverage.getGridGeometry().getGridRange()));
            assertTrue("Envelope comparison failed:" +o.toString(),((GeneralEnvelope)coverage.getGridGeometry().getEnvelope()).equals(destCoverage.getGridGeometry().getEnvelope(),eps,false));
            coverage.dispose(true);
            destCoverage.dispose(true);
          }
    
        } else
          buffer.append("NOT ACCEPTED").append("\n");
        if (TestData.isInteractiveTest())

View Full Code Here

        if (headless) {
            System.out.println("      Skipping test in headless build");
            return;
        }


        GridCoverage2D coverage = createCoverage();
        Style style = createCoverageStyle("1");


        MapContext context = new MapContext();
        context.addLayer(coverage, style);
        StreamingRenderer renderer = new StreamingRenderer();

View Full Code Here

0 1 2 3 4 5 6 7 8 9

TOP

Related Classes of org.geotools.coverage.grid.GridCoverage2D

mil.nga.giat.geowave.raster.plugin.GeoWaveRasterReader

org.geoserver.catalog.rest.CoverageTest

org.geoserver.data.test.MockData

org.geoserver.wcs.GetCoverageTest

org.geoserver.wcs.ResourceAccessManagerWCSTest

org.geoserver.wcs.response.Wcs10GetCoverageResponse

org.geoserver.wcs2_0.GetCoverage

org.geoserver.wcs2_0.kvp.CRSExtentionKVPTest

org.geoserver.wcs2_0.kvp.GeoTiffKvpTest

org.geoserver.wcs2_0.kvp.RangeSubsetKvpTest

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.