Add a few Raster examples.

Author: jericks

examples/raster_crop.groovy

@@ -0,0 +1,16 @@
+import geoscript.layer.*
+import geoscript.geom.Bounds
+import geoscript.proj.Projection
+
+// Read
+def geoTiffFormat = new GeoTIFF()
+def raster = geoTiffFormat.read(new File("raster.tif"), new Projection("urn:ogc:def:crs:EPSG::EPSG:4326"))
+
+// Crop
+def croppedRaster = raster.crop(new Bounds(-10, -10, 10, 10, raster.proj))
+println croppedRaster.bounds
+println croppedRaster.proj.id
+println croppedRaster.size
+
+// Write
+geoTiffFormat.write(croppedRaster, new File("raster_cropped.tif"))

examples/raster_getvalue.groovy

@@ -0,0 +1,14 @@
+import geoscript.geom.Geometry
+import geoscript.layer.*
+import geoscript.geom.Bounds
+import geoscript.proj.Projection
+
+// Read
+def geoTiffFormat = new GeoTIFF()
+def raster = geoTiffFormat.read(new File("raster.tif"), new Projection("urn:ogc:def:crs:EPSG::EPSG:4326"))
+
+def multiPoint = Geometry.createRandomPoints(new Bounds(-180,-90,180,80).geometry, 10)
+multiPoint.points.eachWithIndex{pt,i ->
+    def value = raster.getValue(pt)
+    println "${i}). ${value} at ${pt}"
+}

examples/raster_info.groovy

@@ -0,0 +1,21 @@
+import geoscript.layer.*
+
+// Read a GeoTIFF
+def format = new GeoTIFF()
+def raster = format.read(new File("raster.tif"))
+
+// Print some information
+println "Format = ${raster.format}"
+println "Proj EPSG = ${raster.proj.id}"
+println "Proj WKT = ${raster.proj.wkt}"
+println "Bounds = ${raster.bounds.geometry.wkt}"
+println "Size = ${raster.size}"
+println "Block Size = ${raster.blockSize}"
+println "Pixel Size = ${raster.pixelSize}"
+println "Band:"
+raster.bands.eachWithIndex{b,i ->
+    println "   ${i}). ${b}"
+}
+
+
+

examples/raster_read_write.groovy

@@ -0,0 +1,13 @@
+import geoscript.layer.*
+
+// Read
+def geoTiffFormat = new GeoTIFF()
+def raster = geoTiffFormat.read(new File("raster.tif"))
+
+// Write
+def pngFormat = new WorldImage()
+pngFormat.write(raster, new File("raster.png"))
+
+// Read
+def pngRaster = pngFormat.read(new File("raster.png"))
+println pngRaster.format