在地图上显示线串的方向 - 在地图上自动缩放

时间:2017-07-21 10:25:32

标签: java gis geotools

我有这个代码,它在地图上绘制一个线串,它是用户提供的2个坐标点的轨迹。

public class Quickstart {

public static void main(String[] args) throws Exception {
    // display a data store file chooser dialog for shapefiles
    File file = JFileDataStoreChooser.showOpenFile("shp", null);
    if (file == null) {
        return;
    }

    FileDataStore store = FileDataStoreFinder.getDataStore(file);
    SimpleFeatureSource featureSource = store.getFeatureSource();

  GeometryFactory gf = JTSFactoryFinder.getGeometryFactory();

  double latitude, longitude, latitudeDest, longitudeDest;
  Scanner reader = new Scanner(System.in);
  reader.useLocale(Locale.US);
  System.out.println("Enter reference longitude and latitude:\n");
  longitude = reader.nextDouble();
  latitude = reader.nextDouble();
  System.out.println("Enter destination longitude and latitude:\n");
  longitudeDest = reader.nextDouble();
  latitudeDest = reader.nextDouble();
  reader.close();

    final String EPSG4326 = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\"," +
            "\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\", " + 
            "0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]]";
    CoordinateReferenceSystem crs = CRS.parseWKT(EPSG4326);        


    Point start = gf.createPoint(new Coordinate(longitude, latitude));
    Point end = gf.createPoint(new Coordinate(longitudeDest, latitudeDest));

    GeodeticCalculator gc = new GeodeticCalculator(crs);
    gc.setStartingPosition(JTS.toDirectPosition(start.getCoordinate(), crs));
    gc.setDestinationPosition(JTS.toDirectPosition(end.getCoordinate(), crs));

    // Calculate distance between points
    double distance = gc.getOrthodromicDistance();

    int totalmeters = (int) distance;
    int km = totalmeters / 1000;
    int meters = totalmeters - (km * 1000);
    float remaining_cm = (float) (distance - totalmeters) * 10000;
    remaining_cm = Math.round(remaining_cm);
    float cm = remaining_cm / 100;

    System.out.println("Distance = " + km + "km " + meters + "m " + cm + "cm");

    Coordinate[] coordinates = {start.getCoordinate(), end.getCoordinate()};       
    LineString line = gf.createLineString(coordinates);

    SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
    builder.setName("TwoDistancesType");
    builder.setCRS(DefaultGeographicCRS.WGS84);
    //builder.add("start", Point.class);
    //builder.add("end", Point.class);
    builder.add("line", LineString.class);
    // build the type
    final SimpleFeatureType TYPE = builder.buildFeatureType();

    SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(TYPE);
    featureBuilder.add(line);

    SimpleFeature feature = featureBuilder.buildFeature(null);
    DefaultFeatureCollection featureCollection = new DefaultFeatureCollection("internal", TYPE);
    featureCollection.add(feature);

    // Create style for the line
    //Style style = SLD.createSimpleStyle(TYPE, Color.red);
    Style style = SLD.createLineStyle(Color.red, 2.0f);
    Layer layer = new FeatureLayer(featureCollection, style);

    // Create style for the file
    Style shpStyle = SLD.createSimpleStyle(TYPE, Color.blue);
    Layer shpLayer = new FeatureLayer(featureSource, shpStyle);

    // Create a map content and add our shapefile to it
    MapContent map = new MapContent();
    map.setTitle("TEST");
    map.addLayer(layer);
    map.addLayer(shpLayer);

    // Now display the map
    JMapFrame.showMap(map);


}

我有两个问题:

1)如何显示线的方向?从起点到终点?

2)当你运行程序并看到地图时,你必须手动搜索线串(红线),然后缩放到地图才能找到它。有没有办法自动缩放到线条(坐标) )当地图出现时?

1 个答案:

答案 0 :(得分:1)

对于样式,您需要类似SLD描述的here,代码变为:

    // Create style for the line
    // Style style = SLD.createSimpleStyle(TYPE, Color.red);
    org.geotools.styling.Style style = SLD.createLineStyle(Color.red, 2.0f);
    StyleBuilder sb = new StyleBuilder();
    FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2();
    PointSymbolizer point = sb.createPointSymbolizer();
    Mark mark = sb.createMark("shape://oarrow");
    mark.setFill(sb.createFill(Color.RED));
    mark.setStroke(sb.createStroke(Color.red));

    Graphic graphic = sb.createGraphic(null, mark, null);
    graphic.setRotation(ff.function("endAngle", ff.property("line")));
    point.setGraphic(graphic);

    point.setGeometry(ff.function("endpoint", ff.property("line")));

    Rule rule = sb.createRule(point);
    style.getFeatureTypeStyles()[0].addRule(rule );
    Layer layer = new FeatureLayer(featureCollection, style);

放大直线只是将地图视口设置为直线边界的情况:

    MapViewport viewport = new MapViewport(featureCollection.getBounds());
    map.setViewport(viewport );

如果你想要,你可能想要将这些边界增长一点(10%?),这样你就可以看到周围环境。

enter image description here

修改

要避免StyleBuilder中已弃用的方法,您可以使用:

style.featureTypeStyles().get(0).rules().add(rule);

扩展边界框只是在信封上添加一些距离的情况:

    ReferencedEnvelope bounds = featureCollection.getBounds();
    double delta = bounds.getWidth()/20.0; //5% on each side
    bounds.expandBy(delta );
    MapViewport viewport = new MapViewport(bounds);
    map.setViewport(viewport );