Scripting GeoServer

Download as pdf or txt
Download as pdf or txt
You are on page 1of 28

Scripting GeoServer with

GeoScript
Justin Deoliveira, OpenGeo

GeoScript
Spatial capabilities for scripting languages.

GeoTools
Library providing tools for geospatial data.

GeoServer
Application server for publishing, editing, and
processing spatial data.

GeoServer Scripting Extension


Bringing the convenience of GeoScript together
with the power of GeoServer.

JSR-223

Scripting Hooks

Scripting Hooks

WPS Process

WFS Transaction Callback

HTTP Endpoint

Filter Function

Output Format *

Data Format *

* In progress

Demos

Read PostGIS Features - Java


...
import org.geotools.data.DataStore;
import org.geotools.data.postgis.PostgisNGDataStoreFactory;
...
Map<String,Serializable> params = new HashMap<String, Serializable>();
params.put(JDBCDataStoreFactory.HOST.key, "localhost");
params.put(JDBCDataStoreFactory.PORT.key, 5432);
params.put(JDBCDataStoreFactory.DATABASE.key, "geoscript");
params.put(JDBCDataStoreFactory.DBTYPE.key, "postgis");
params.put(JDBCDataStoreFactory.USER.key, "jdeolive");
PostgisNGDataStoreFactory factory = new PostgisNGDataStoreFactory();
DataStore pg = factory.createDataStore(params);
FeatureCollection features = pg.getFeatureSource("states").getFeatures();
FeatureIterator it = features.features();
try {
while(it.hasNext()) {
Feature f = it.next();
System.out.println(f.getProperty("STATE_NAME").getValue());
}
}
finally {
it.close();
}

Read PostGIS Features - Python

from geoscript.workspace import PostGIS


pg = PostGIS('geoscript')
for f in pg['states'].features():
print f['STATE_NAME']

Read PostGIS Features - JavaScript

var ws = require("geoscript/workspace");
var pg = ws.PostGIS("geoscript");
for (var f in pg.get("states").features) {
print(f.get("STATE_NAME"));
}

Geometry

>>> from geoscript import geom


>>> geom.Point(30, 10)
POINT(30 10)
>>> poly = geom.Point(30,10).buffer(5)
>>> poly.area
78.03612880645133

Geometry - I/O

>>> from geoscript import geom


>>> point = geom.Point(30, 10)
>>> geom.writeKML(point)
<kml:Point
xmlns:kml="http://earth.google.com/kml/2.1">
<kml:coordinates>0.0,0.0</kml:coordinates>
</kml:Point>
>>> geom.writeJSON(point)
{"type":"Point","coordinates":[30,10]}

Geometry - Visualization
js> var geom = require("geoscript/geom")
js> require("geoscript/viewer").bind()
js> var poly1 = geom.Point([0, 0]).buffer(1)
js> var poly2 = poly1.transform({dx: 0.5, dy: 0.5})
js> poly1.difference(poly2)
<Polygon [[[0.9095298326166407, -0.409529...>

Projection

js> var proj = require("geoscript/proj");


js> var p = proj.Projection("epsg:4326");
js> p.wkt
GEOGCS["WGS 84",
DATUM["World Geodetic System 1984",
...

Projection
>>> from geoscript import geom
>>> from geoscript.proj import Projection
>>> p = Projection('epsg:4326')
>>> p.transform((-111, 45.7), 'epsg:26912')
(500000.0, 5060716.313515949)
>>> g = geom.Point(0, 0).buffer(4)
>>> g = ...
>>> p.transform(g, 'epsg:26912')
>>> p.transform(g, 'epsg:3005')

WGS
84

UTM

Albers

Data Access - Workspace


>> from geoscript.workspace import PostGIS
>> pg = PostGIS('spearfish')
>> pg.layers()
['archsites', 'bugsites', ..., 'streams']
>> l = pg['archsites']

Data Access - Layers

>>> from geoscript.layer import Shapefile


>>> states = Shapefile('states.shp')
>>> states = states.reproject('epsg:3005')

Data Access - Features


js> var ws = require("geoscript/workspace")
js> var dir = ws.Directory("Data")
js> var states = dir.get("states")
js> for (var f in states.query("STATE_ABBR LIKE 'M%'")) {
>
print(f.get("STATE_NAME"));
> }
Maryland
Missouri
Mississippi
Michigan
...

Data Access - Rasters

>>> from geoscript.layer import GeoTIFF


>>> dem = GeoTIFF('sfdem.tif')
>>> dem.bands
[GRAY_INDEX]
>>> dem.extrema()
([-9.999999933815811e+36], [1840.0])
>>> h = dem.histogram(0, 1850)

Data Access - Plotting


>>> from geoscript.layer import GeoTIFF
>>> dem = GeoTIFF('sfdem.tif')
>>> h = dem.histogram(0, 1850)
>>>
>>>
>>>
>>>

from geoscript.plot import bar


y = h.counts()
x = map(lambda (x,y): y, h.bins())
bar.xy(zip(x, y)).show()

Styling - Stroke
>>> from geoscript.style import Stroke
>>> Stroke('#000000', width=2)
>>> Stroke('black', width=2, dash=[5,5])
>>> Stroke((0,0,0),width=2).hatch('vertline')

Styling - Shape and Icon


js> var style = require("geoscript/style");
js> style.Shape({name: "star", fill: "yellow"})
<Shape name: 'star', size: 6>
js> style.Icon("rainy.svg")
<Icon url: 'rainy.svg'>

Styling - Labels
>>> from geoscript.style import *
>>> font = 'bold 16px Arial'
>>> Shape() + Label('name',font)
.point(displace=(20,0))
>>> Stroke() + Label('name',font)
.linear(offset=10)
>>> Fill() + Label('name',font).halo('white',2)

Styling - Theming

>>> from geoscript.style *


>>> style = Stroke() + Label('STATE_ABBR', 14, 'Serif')
>>> style += Fill('#4DFF4D', 0.7).where('PERSONS < 2E6')
>>> style += Fill('#FF4D4D', 0.7).where('PERSONS BETWEEN 2E6 AND 4E6)
>>> style += Fill('#4D4DFF', 0.7).where('PERSONS > 4E6)

Process
>>> from geoscript.layer import Shapefile
>>> from geoscript.process import Process
>>> cities = Shapefile('world_cities.shp')
>>> p = Process.lookup('vec:Heatmap')
>>> p.run(cities, weightAttr='population')

Process

>>> p = Process.lookup('ras:PolygonExtraction')
>>> p.run(hm, range=(0.25, 0.85))

Thanks!

You might also like