Search code examples
javagisgeotools

Show line on map which connects points - receiving an empty map


I want to plot a line which connects 2 points on a map.

The code I am using:

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

      // ask for current and destination positions
      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");

        SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
        builder.setName("TwoDistancesType");
        builder.setCRS(DefaultGeographicCRS.WGS84);
        builder.add("location", Point.class);

        // build the type
        final SimpleFeatureType TYPE = builder.buildFeatureType();

        SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(TYPE);
        featureBuilder.add(start);
        //featureBuilder.add(end);

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

        Style style = SLD.createSimpleStyle(TYPE, Color.red);
        Layer layer = new FeatureLayer(featureCollection, style);

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

        map.addLayer(layer);

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


    }

}

I have 2 problems:

1) I can't add a second feature to featureBuilder.It doesn't allow it.It shows Can handle 1 attributes only, index is 1.

So, how can I plot a line then?

2) With the above code, I am receiving:

org.geotools.renderer.lite.StreamingRenderer fireErrorEvent SEVERE: The scale denominator must be positive
java.lang.IllegalArgumentException: The scale denominator must be positive

------- UPDATE ------------------------

After the solution that @Michael gave for the first question , now I no longer receive the error regarding the denominator, but I am receiving an empty map (white space).

----- UPDATE according to @iant suggestion ----------------

So, I tried this.Created a coordinates which holds the coordinates of the points (start and end) ,then created a linestring with those coordinates and added it to featurebuilder.

        SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
        builder.setName("TwoDistancesType");
        builder.setCRS(DefaultGeographicCRS.WGS84);

        builder.add("line", LineString.class); //added a linestring class


        final SimpleFeatureType TYPE = builder.buildFeatureType();

        SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(TYPE);

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

and even though I am loading a map (countries.shp) it shows me an empty white space with a red line.

------ SOLUTION -------------

Ok, so the solution is (thanks to @iants comments):

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

and now you have a red line on a blue map!


Solution

  • You need to create a LineString from your points and then store that in your feature.

    You should then get a correct scale but you might want to add some other data such as a coast line to the map first. The quick start tutorial can show you how to do that.