jericks

File file = new File("earth.tif")
GeoTIFF geotiff = new GeoTIFF(file)
Raster raster = geotiff.read("earth")

File arcGridFile = new File("earth.asc")
ArcGrid arcGrid = new ArcGrid(arcGridFile)
arcGrid.write(raster)
Raster arcGridRaster = arcGrid.read("earth")

arcGridRaster.eachCell {double value, double x, double y ->

jericks

package geoscript.layer

import geoscript.feature.Feature
import geoscript.workspace.Directory
import geoscript.workspace.Memory
import org.geotools.data.DataUtilities
import org.geotools.data.DefaultTransaction
import org.geotools.data.FeatureStore
import org.geotools.data.Transaction
import org.geotools.data.collection.ListFeatureCollection

jericks

import geoscript.layer.*
import geoscript.feature.*
import geoscript.workspace.*

import org.geotools.data.Transaction
import org.geotools.data.DefaultTransaction
import org.geotools.data.FeatureStore
import org.geotools.data.DataUtilities

/**

jericks

import geoscript.layer.*
import geoscript.feature.*
import geoscript.workspace.*

import org.geotools.data.Transaction
import org.geotools.data.DefaultTransaction
import org.geotools.data.FeatureStore
import org.geotools.data.DataUtilities
import org.geotools.data.FeatureWriter as GtFeatureWriter

jericks

import geoscript.layer.Shapefile
import geoscript.proj.Projection
import geoscript.feature.Field

def shp = new Shapefile("test_uk.shp")
def out = shp.workspace.create(shp.schema.reproject("EPSG:27700", "test_uk_valid"))
out.add(shp.cursor.collect{f ->
    f.geom = Projection.transform(f.geom, "EPSG:4326", "EPSG:27700")
    if (!f.geom.isValid()) {
        f.geom = f.geom.buffer(0)

jericks

import geoscript.geom.Bounds
import geoscript.map.Map
import geoscript.workspace.Directory
import geoscript.style.*
import geoscript.filter.Function

def dir = new Directory("/Users/jericks/Downloads/OpenGeo Workshop/data/shp-gs")
def layer = dir.get("nyc_buildings")
layer.style = new Style([
    new PolygonSymbolizer(

jericks

import geoscript.layer.*
import geoscript.proj.*
import geoscript.style.*
import geoscript.map.*
import geoscript.filter.Function

def statesShp = new Shapefile("states.shp")

statesShp.style = new Style([
    new PolygonSymbolizer(

jericks

import geoscript.feature.Field
import org.geotools.factory.CommonFactoryFinder
import org.opengis.filter.expression.Function as GtFunction
import org.geotools.filter.FunctionImpl
import org.opengis.filter.expression.Expression

/**
 * A GeoScript Function either wraps an existing GeoTools Function or
 * can easily create a new Function by passing in a Closure.
 * @author Jared Erickson

jericks

elif options.offsetcalc:
        import exifread
        from datetime import datetime
        imgDate = datetime.strptime(args[0], "%Y-%m-%d %H:%M:%S")
        f = open(args[1], 'rb')
        tags = exifread.process_file(f)
        f.close()
        exifDateStr = tags['EXIF DateTimeDigitized'].values
        exifDate = datetime.strptime(exifDateStr, "%Y:%m:%d %H:%M:%S")
        diff = exifDate - imgDate

jericks

import geoscript.layer.*
import geoscript.geom.Bounds
import geoscript.proj.Projection
import javax.imageio.ImageIO

MBTiles mbtiles = new MBTiles(new File("geography-class.mbtiles"))
Raster raster = mbtiles.getRaster(mbtiles.tiles(4)).reproject(new Projection("EPSG:4326")).resample(bbox: new Bounds(-180,-90,180,90,"EPSG:4326"), size: [800,400])
ImageIO.write(raster.image, "png", new File("geography-class.png"))