Show direction of linestring on map - auto zoom on map

冷暖自知 提交于 2019-12-10 22:11:14

问题


I have this code, which plots on a map a linestring which is the track of 2 coordinate points that the user supplies.

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);


}

I have 2 questions:

1) How can I show the direction of the line?From start point to end?

2) When you run the program and see the map , you must manually search for the linestring(red line) and then zoom to map in order to find it.Is there a way to automatically zoom to the line (the coordinates) when the map appears?


回答1:


For the style you need something like the SLD described here, in code that becomes:

    // 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);

Zooming into the line is just a case of setting the map viewport to the bounds of the line:

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

If you want you might want to grow those bounds by a little (10%?) so that you can see the surroundings too.

EDIT

To avoid the deprecated methods in StyleBuilder you can use:

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

Expanding the bounding box is just a case of adding some distance to the envelope:

    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 );


来源:https://stackoverflow.com/questions/45235244/show-direction-of-linestring-on-map-auto-zoom-on-map

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!