groovyとgeotoolsでラインを描画する

groovyとgeotoolsでラインを描画するには、以下のコードを実行します。

import java.awt.*
import java.awt.image.*
import javax.imageio.*
import org.geotools.data.shapefile.*
import org.geotools.factory.*
import org.geotools.feature.*
import org.geotools.feature.simple.*
import org.geotools.filter.*
import org.geotools.map.*
import org.geotools.renderer.lite.*
import org.geotools.styling.*
import org.geotools.geometry.jts.*
import com.vividsolutions.jts.geom.*
import org.geotools.data.*
import org.opengis.feature.type.*

url = new URL("file://C:/share/geotools/world.shp")
shapefile = new ShapefileDataStore(url)

fs = shapefile.getFeatureSource()
schema = fs.getSchema()
crs = schema.getGeometryDescriptor().getCoordinateReferenceSystem()
map = new DefaultMapContext([] as MapLayer[], crs)

ff = FilterFactoryFinder.createFilterFactory()
sf = CommonFactoryFinder.getStyleFactory()

// フィルタを作成
filter = ff.createCompareFilter(FilterType.COMPARE_EQUALS)
filter.addLeftValue(ff.createAttributeExpression("NAME"))
filter.addRightValue(ff.createLiteralExpression("JAPAN"))
fs = fs.getFeatures(filter)

// ポリゴンの線の色1
stroke1 = sf.createStroke(
  ff.literal(new Color(0x59, 0x58, 0x55)),
  ff.literal(1)
)

// ポリゴンの塗りつぶし色1
fill1 = sf.createFill(
  ff.literal(new Color(0xD2, 0xDB, 0xD5))
)
sym1 = sf.createPolygonSymbolizer(stroke1, fill1, null)

// フィルタに一致する場合のルールを作成
rule1 = sf.createRule()
rule1.symbolizers().add(sym1)

fts = sf.createFeatureTypeStyle([rule1] as Rule[])
style1 = sf.createStyle()
style1.featureTypeStyles().add(fts)

// レイヤーとして追加
map.addLayer(new FeatureLayer(fs, style1))


// ------------------------------------------
// レイヤーを作成
// タイプを作成
pointtype = DataUtilities.createType(
  "Location", 
  "the_geom:LineString," + 
  "name:String"
)
sfb = new SimpleFeatureBuilder(pointtype)

// 東京を示すポイントを作成
// 経度
longitude1 = 139.6917064d
// 緯度
latitude1 = 35.6894875d
// 那覇を示すポイントを作成
// 経度
longitude2 = 127.67898559570312d
// 緯度
latitude2 = 26.19364441820758d
gf = JTSFactoryFinder.getGeometryFactory(null)
lineseg = gf.createLineString([
  new Coordinate(longitude1, latitude1),
  new Coordinate(longitude2, latitude2)] as Coordinate[]
)
sfb.add(lineseg)
sfb.add("tokyo and naha")

col = FeatureCollections.newCollection()
feature1 = sfb.buildFeature(null)
col.add(feature1)

// 線の色
stroke2 = sf.createStroke(
  ff.literal(new Color(0xC8, 0x46, 0x63)),
  ff.literal(4)
)

GeometryDescriptor geomDesc = fs.getSchema().getGeometryDescriptor()
geometryAttributeName = geomDesc.getLocalName()
sym2 = sf.createLineSymbolizer(stroke2, geometryAttributeName)

rule2 = sf.createRule()
rule2.symbolizers().add(sym2)
fts2 = sf.createFeatureTypeStyle([rule2] as Rule[])
style2 = sf.createStyle()
style2.featureTypeStyles().add(fts2)
// レイヤーとして追加
map.addLayer(new FeatureLayer(col, style2))

// レンダリング
renderer = new StreamingRenderer()
renderer.setContext(map)

width = 400
bounds = map.getLayerBounds()
rect = new Rectangle(0, 0, width, 
  (int)(width * bounds.getHeight() / bounds.getWidth()))

image = new BufferedImage((int)rect.width, (int)rect.height, 
  BufferedImage.TYPE_INT_RGB)
gr = image.createGraphics()
gr.setPaint(Color.WHITE)
gr.fill(rect)
renderer.paint(gr, rect, bounds)

ImageIO.write(image, "jpeg", new File("tokyo_and_naha.jpg"))

出力画像

※世界地図のシェイプファイルは以下からダウンロード
World map for APRS
http://aprsworld.net/gisdata/world/

※日本地図は国土地理院のウェブサイトからダウンロードできます。
ダウンロード・利用規約などは以下を参照。
http://www1.gsi.go.jp/geowww/globalmap-gsi/download/index.html

※.groovy/libからGROOVY_HOME/libにgeotools-2.7.3-bin.zipに
含まれるjarをコピーする

動作環境
groovy 1.8.2, JDK6 Update27, GeoTools 2.7.3

関連情報
・GeoToolsのウェブサイト
http://geotools.org/