画像処理についてあれこれ

Friday, February 17, 2012

Luaで線の幅を指定する

Luaで線の幅を指定するには、以下のコードを実行します。

require"imlua"
require"cdlua"
require"cdluaim"

-- キャンバスを作成
image = im.ImageCreate(200, 200, im.RGB, im.BYTE)
canvas = image:cdCreateCanvas()

canvas:Activate()  
-- 背景色
canvas:Background(cd.EncodeColor(255, 255, 255))
canvas:Clear()
-- 線の色
canvas:Foreground(cd.EncodeColor(0, 0, 255))
-- 太線を描画
canvas:LineWidth(10)
canvas:Line(0, 0, 200, 200)
canvas:Kill()

-- ファイルを保存
image:Save("linewidth.png", "PNG")

出力画像

動作環境
Lua for Windows 5.1.4

Tuesday, February 14, 2012

Luaでラインスタイルを設定する

Luaでラインスタイルを設定するには、以下のコードを実行します。

require"imlua"
require"cdlua"
require"cdluaim"

-- キャンバスを作成
image = im.ImageCreate(200, 200, im.RGB, im.BYTE)
canvas = image:cdCreateCanvas()

canvas:Activate()  
-- 背景色
canvas:Background(cd.EncodeColor(255, 255, 255))
canvas:Clear()
-- 線の色
canvas:Foreground(cd.EncodeColor(0, 0, 255))
-- 破線を描画
canvas:LineStyle(cd.DASHED)
canvas:Line(40, 0, 40, 200)

-- 点線を描画
canvas:LineStyle(cd.DOTTED)
canvas:Line(80, 0, 80, 200)

-- 線を描画(DASH_DOT)
canvas:LineStyle(cd.DASH_DOT)
canvas:Line(120, 0, 120, 200)

-- 線を描画(DASH_DOT_DOT)
canvas:LineStyle(cd.DASH_DOT_DOT)
canvas:Line(160, 0, 160, 200)

canvas:Kill()

-- ファイルを保存
image:Save("linestyle.png", "PNG")

出力画像

動作環境
Lua for Windows 5.1.4

Saturday, February 11, 2012

Luaで直線を描画する

Luaで直線を描画するには、以下のコードを実行します。

require"imlua"
require"cdlua"
require"cdluaim"

-- キャンバスを作成
image = im.ImageCreate(200, 200, im.RGB, im.BYTE)
canvas = image:cdCreateCanvas()

canvas:Activate()  
-- 背景色
canvas:Background(cd.EncodeColor(255, 255, 255))
canvas:Clear()
-- 線の色
canvas:Foreground(cd.EncodeColor(0, 0, 255))
-- 線を描画
canvas:Line(0, 0, 200, 200)
canvas:Kill()

-- ファイルを保存
image:Save("line.png", "PNG")

出力画像

動作環境
Lua for Windows 5.1.4

Friday, January 27, 2012

GeoToolsでラインを描画する

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

import java.awt.*;
import java.awt.image.*;
import java.io.*;
import java.net.*;
import java.util.*;
import javax.imageio.*;
import org.geotools.data.*;
import org.geotools.data.shapefile.*;
import org.geotools.data.simple.*;
import org.opengis.filter.*;
import org.geotools.factory.*;
import org.geotools.feature.*;
import org.geotools.feature.simple.*;
import org.geotools.geometry.*;
import org.geotools.geometry.jts.*;
import org.geotools.map.*;
import org.geotools.renderer.lite.*;
import org.geotools.styling.*;
import org.opengis.feature.simple.*;
import org.opengis.feature.type.*;
import org.opengis.geometry.coordinate.*;
import org.opengis.referencing.crs.*;
import com.vividsolutions.jts.geom.*;

public class Geotools16
{
  public static void main(String args[])
    throws Exception
  {
    URL url = new URL("file://C:/share/geotools/world.shp");
    ShapefileDataStore shapefile = new ShapefileDataStore(url);

    SimpleFeatureSource fs = shapefile.getFeatureSource();
    // 国の名前でフィルタリング
    FilterFactory ff = CommonFactoryFinder.getFilterFactory2(null);
    Filter filter = ff.equals(
      ff.property("NAME"),
      ff.literal("JAPAN")
    );
    SimpleFeatureCollection sfc = fs.getFeatures(filter);

    SimpleFeatureType schema = sfc.getSchema();
    CoordinateReferenceSystem crs = 
      schema.getGeometryDescriptor().getCoordinateReferenceSystem();
    MapLayer layers[] = {};
    DefaultMapContext map = new DefaultMapContext(layers, crs);

    //スタイルを作成
    StyleFactory sf = CommonFactoryFinder.getStyleFactory(null);
    // ポリゴンの線の色
    org.geotools.styling.Stroke stroke1 = sf.createStroke(
      ff.literal(new Color(0x59, 0x58, 0x55)),
      ff.literal(1)
    );
    // ポリゴンの塗りつぶし色
    org.geotools.styling.Fill fill1 = sf.createFill(
      ff.literal(new Color(0xD2, 0xDB, 0xD5))
    );
    PolygonSymbolizer sym1 = sf.createPolygonSymbolizer(stroke1, fill1, null);

    // フィルタに一致する場合のルールを作成
    Rule rule1 = sf.createRule();
    rule1.symbolizers().add(sym1);
    rule1.setFilter(filter);
    Rule rules[] = {rule1};
    FeatureTypeStyle fts = sf.createFeatureTypeStyle(rules);
    Style style1 = sf.createStyle();
    style1.featureTypeStyles().add(fts);

    map.addLayer(new FeatureLayer(sfc, style1));

    // --------------------------------------------------
    // 別のレイヤーを作成
    // ポイント用のタイプを作成
    SimpleFeatureType pointtype = DataUtilities.createType(
      "Location", 
      "the_geom:LineString," + 
      "name:String"
    );
    SimpleFeatureBuilder sfb = new SimpleFeatureBuilder(pointtype);

    // 東京を示すポイントを作成
    // 経度
    double longitude1 = 139.6917064d;
    // 緯度
    double latitude1 = 35.6894875d;

    // 那覇を示すポイントを作成
    // 経度
    double longitude2 = 127.67898559570312d;
    // 緯度
    double latitude2 = 26.19364441820758d;

    com.vividsolutions.jts.geom.GeometryFactory gf = JTSFactoryFinder.getGeometryFactory(null);
    Coordinate cords[] = {
      new Coordinate(longitude1, latitude1), 
      new Coordinate(longitude2, latitude2)
    };
    com.vividsolutions.jts.geom.LineString lineseg = gf.createLineString(cords);
    sfb.add(lineseg);
    sfb.add("tokyo and naha");

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

    // 線の色
    org.geotools.styling.Stroke stroke2 = sf.createStroke(
      ff.literal(new Color(0xC8, 0x46, 0x63)),
      ff.literal(3)
    );

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

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

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

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

    BufferedImage image = new BufferedImage((int)rect.width, 
      (int)rect.height, 
      BufferedImage.TYPE_INT_RGB);
    Graphics2D 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"));
  }
}

出力画像

動作環境
JDK7 Update1, geotools 2.7.4

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

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

Wednesday, January 25, 2012

GeoToolsで指定した地点にクロスマークを描画する

GeoToolsで指定した地点にクロスマークを描画するには、以下のコードを実行します。

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

public class Geotools15
{
  public static void main(String args[])
    throws Exception
  {
    URL url = new URL("file://C:/share/geotools/world.shp");
    ShapefileDataStore shapefile = new ShapefileDataStore(url);

    SimpleFeatureSource fs = shapefile.getFeatureSource();
    // 国の名前でフィルタリング
    FilterFactory ff = CommonFactoryFinder.getFilterFactory2(null);
    Filter filter = ff.equals(
      ff.property("NAME"),
      ff.literal("JAPAN")
    );
    SimpleFeatureCollection sfc = fs.getFeatures(filter);

    SimpleFeatureType schema = sfc.getSchema();
    CoordinateReferenceSystem crs = 
      schema.getGeometryDescriptor().getCoordinateReferenceSystem();
    MapLayer layers[] = {};
    DefaultMapContext map = new DefaultMapContext(layers, crs);

    //スタイルを作成
    StyleFactory sf = CommonFactoryFinder.getStyleFactory(null);
    // ポリゴンの線の色
    org.geotools.styling.Stroke stroke1 = sf.createStroke(
      ff.literal(new Color(0x59, 0x58, 0x55)),
      ff.literal(1)
    );
    // ポリゴンの塗りつぶし色
    org.geotools.styling.Fill fill1 = sf.createFill(
      ff.literal(new Color(0xD2, 0xDB, 0xD5))
    );
    PolygonSymbolizer sym1 = sf.createPolygonSymbolizer(stroke1, fill1, null);

    // フィルタに一致する場合のルールを作成
    Rule rule1 = sf.createRule();
    rule1.symbolizers().add(sym1);
    rule1.setFilter(filter);
    Rule rules[] = {rule1};
    FeatureTypeStyle fts = sf.createFeatureTypeStyle(rules);
    Style style1 = sf.createStyle();
    style1.featureTypeStyles().add(fts);

    map.addLayer(new FeatureLayer(sfc, style1));

    // --------------------------------------------------
    // 別のレイヤーを作成
    // ポイント用のタイプを作成
    SimpleFeatureType pointtype = DataUtilities.createType(
      "Location", 
      "the_geom:Point," + 
      "name:String"
    );
    SimpleFeatureBuilder sfb = new SimpleFeatureBuilder(pointtype);

    // 大阪を示すポイントを作成
    // 経度
    double longitude = 135.5021651d;
    // 緯度
    double latitude = 34.6937378d;
    GeometryFactory gf = JTSFactoryFinder.getGeometryFactory(null);
    com.vividsolutions.jts.geom.Point point = gf.createPoint(new Coordinate(longitude, latitude));
    sfb.add(point);
    sfb.add("osaka");

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

    // マーカーの線
    org.geotools.styling.Stroke stroke2 = sf.createStroke(
      ff.literal(new Color(0xC8, 0x46, 0x63)),
      ff.literal(3)
    );
    // マーカーの塗りつぶし色
    org.geotools.styling.Fill fill2 = sf.createFill(ff.literal(new Color(0xff, 0xC8, 0x61)));

    // マーカーの形
    Mark mark = sf.getCrossMark();
    mark.setFill(fill2);
    mark.setStroke(stroke2);

    Graphic graphic = sf.createDefaultGraphic();
    graphic.graphicalSymbols().clear();
    graphic.graphicalSymbols().add(mark);
    graphic.setSize(ff.literal(18));

    GeometryDescriptor geomDesc = fs.getSchema().getGeometryDescriptor();
    String geometryAttributeName = geomDesc.getLocalName();
    PointSymbolizer sym2 = sf.createPointSymbolizer(graphic, geometryAttributeName);

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

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

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

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

    ImageIO.write(image, "jpeg", new File("osaka.jpg"));
  }
}

出力画像

Monday, January 23, 2012

GeoToolsで表示範囲を設定する

GeoToolsで表示範囲を設定するには、以下のコードを実行します。

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

public class Geotools14
{
  public static void main(String args[])
    throws Exception
  {
    URL url = new URL("file://C:/share/geotools/world.shp");
    ShapefileDataStore shapefile = new ShapefileDataStore(url);

    SimpleFeatureSource fs = shapefile.getFeatureSource();
    // 国の名前でフィルタリング
    FilterFactory ff = CommonFactoryFinder.getFilterFactory2(null);
    Filter filter = ff.equals(
      ff.property("NAME"),
      ff.literal("JAPAN")
    );
    SimpleFeatureCollection sfc = fs.getFeatures(filter);

    SimpleFeatureType schema = sfc.getSchema();
    CoordinateReferenceSystem crs = 
      schema.getGeometryDescriptor().getCoordinateReferenceSystem();
    MapLayer layers[] = {};
    DefaultMapContext map = new DefaultMapContext(layers, crs);

    //スタイルを作成
    StyleFactory sf = CommonFactoryFinder.getStyleFactory(null);
    // ポリゴンの線の色
    org.geotools.styling.Stroke stroke1 = sf.createStroke(
      ff.literal(new Color(0x59, 0x58, 0x55)),
      ff.literal(1)
    );
    // ポリゴンの塗りつぶし色
    org.geotools.styling.Fill fill1 = sf.createFill(
      ff.literal(new Color(0xD2, 0xDB, 0xD5))
    );
    PolygonSymbolizer sym1 = sf.createPolygonSymbolizer(stroke1, fill1, null);

    // フィルタに一致する場合のルールを作成
    Rule rule1 = sf.createRule();
    rule1.symbolizers().add(sym1);
    rule1.setFilter(filter);
    Rule rules[] = {rule1};
    FeatureTypeStyle fts = sf.createFeatureTypeStyle(rules);
    Style style1 = sf.createStyle();
    style1.featureTypeStyles().add(fts);

    map.addLayer(new FeatureLayer(sfc, style1));

    // --------------------------------------------------
    // 別のレイヤーを作成
    // ポイント用のタイプを作成
    SimpleFeatureType pointtype = DataUtilities.createType(
      "Location", 
      "the_geom:Point," + 
      "name:String"
    );
    SimpleFeatureBuilder sfb = new SimpleFeatureBuilder(pointtype);

    // 東京を示すポイントを作成
    // 経度
    double longitude = 139.6917064d;
    // 緯度
    double latitude = 35.6894875d;
    GeometryFactory gf = JTSFactoryFinder.getGeometryFactory(null);
    com.vividsolutions.jts.geom.Point point = gf.createPoint(new Coordinate(longitude, latitude));
    sfb.add(point);
    sfb.add("tokyo");

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

    // マーカーの線
    org.geotools.styling.Stroke stroke2 = sf.createStroke(
      ff.literal(new Color(0xC8, 0x46, 0x63)),
      ff.literal(3)
    );
    // マーカーの塗りつぶし色
    org.geotools.styling.Fill fill2 = sf.createFill(ff.literal(new Color(0xff, 0xC8, 0x61)));

    // マーカーの形
    Mark mark = sf.getCircleMark();
    mark.setFill(fill2);
    mark.setStroke(stroke2);

    Graphic graphic = sf.createDefaultGraphic();
    graphic.graphicalSymbols().clear();
    graphic.graphicalSymbols().add(mark);
    graphic.setSize(ff.literal(10));

    GeometryDescriptor geomDesc = fs.getSchema().getGeometryDescriptor();
    String geometryAttributeName = geomDesc.getLocalName();
    PointSymbolizer sym2 = sf.createPointSymbolizer(graphic, geometryAttributeName);

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

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

    int width = 400;
    // 表示範囲を設定
    ReferencedEnvelope bounds = new ReferencedEnvelope(
      longitude-1, longitude+1, latitude-1, latitude+1, 
      map.getViewport().getBounds().getCoordinateReferenceSystem()
    );
    Rectangle rect = new Rectangle(0, 0, width, 
      (int)(width * bounds.getHeight() / bounds.getWidth()));

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

    ImageIO.write(image, "jpeg", new File("around_tokyo.jpg"));
  }
}

出力画像

画像処理についてあれこれ

Friday, February 17, 2012

Luaで線の幅を指定する

Tuesday, February 14, 2012

Luaでラインスタイルを設定する

Saturday, February 11, 2012

Luaで直線を描画する

Friday, January 27, 2012

GeoToolsでラインを描画する

Wednesday, January 25, 2012

GeoToolsで指定した地点にクロスマークを描画する

Monday, January 23, 2012

GeoToolsで表示範囲を設定する

About Me

Labels