shpファイルから座標を読み取り、距離を計算する

Question

地点からshpファイル（ports.shp）までの距離を自然地理データから計算したい。

例えば、私は、ファイルの特徴をロードしています：

... 
String filename = "10m_cultural/ne_10m_ports.shp"; 
... 


public static void Calcs(String filename) 
    throws IOException, NoSuchAuthorityCodeException, FactoryException, TransformException { 

    HashMap<String, Object> params = new HashMap<>(); 
    params.put("url", DataUtilities.fileToURL(new File(filename))); 
    DataStore ds = DataStoreFinder.getDataStore(params); 

    String name = ds.getTypeNames()[0]; 
    SimpleFeatureSource source = ds.getFeatureSource(name); 
    SimpleFeatureCollection features = source.getFeatures(); 

}

今、私は距離を計算したいから、例えば点は次のとおりです。

GeometryFactory gf = JTSFactoryFinder.getGeometryFactory(); 
Point p = gf.createPoint(new Coordinate(43, 18)); 

私は距離を計算することを知っています：

 CoordinateReferenceSystem crs = CRS.decode("EPSG:4326");   


     Point start = gf.createPoint(new Coordinate(43, 18)); 
     Point dest = gf.createPoint(new Coordinate(?????)); 

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


     double distance = gc.getOrthodromicDistance(); 

ファイルをロードするから

Point dest = gf.createPoint(new Coordinate(?????));

私はfeatures持っているが、それはどんなgetCoordinates()を持っていない：

が、私は先の点の座標（ports.shpファイル）を見つける方法がわかりません方法。

また、ports.shpは多くのPOINTジオメトリで構成されています。私は何らかの理由で基準点のすべての点を計算してから、最も近い点を選択する必要がありますか？

Answer 1

フィーチャーには、必要なポイントを与えるgetDefaultGeometryメソッドがあります。次にポイントから座標を取得できます。

EDITあなたの問題は、あなたがユニットのミスマッチ（、度でそう約360）バウンディングボックスの幅にMinDistを設定するが、メートル単位での距離と比較することしてた

（その程度780万）ので、あなたは保存するのに十分近いところにポイントを見つけることはありませんでした。

最初の検索範囲を制限して検索を効率的にすることでプレイを開始しましたが、それが役立つかどうかは実際には分かりません。

final double MAX_SEARCH_DISTANCE = Math.max(index.getBounds().getWidth(), index.getBounds().getHeight()); 
    double searchDist = 0.01; 

    while (searchDist < MAX_SEARCH_DISTANCE) { 
     // start point (user input) 
     Coordinate coordinate = p.getCoordinate(); 

     ReferencedEnvelope search = new ReferencedEnvelope(new Envelope(coordinate), 
       index.getSchema().getCoordinateReferenceSystem()); 

     search.expandBy(searchDist); 
     BBOX bbox = ff.bbox(ff.property(index.getSchema().getGeometryDescriptor().getName()), (BoundingBox) search); 
     SimpleFeatureCollection candidates = index.subCollection(bbox); 

     double minDist = Double.POSITIVE_INFINITY; // can't use 
                // MAX_Search_dist here 
                // as it is degrees and 
                // dists are meters 
     Coordinate minDistPoint = null; 
     double dist = 0; 
     Point dest = null; 
     SimpleFeatureIterator itr = candidates.features(); 
     CoordinateReferenceSystem crs = DefaultGeographicCRS.WGS84; 
     try { 
      SimpleFeature feature = null; 
      while (itr.hasNext()) { 
       feature = itr.next(); 

       // destination point 
       dest = (Point) feature.getDefaultGeometry(); 
       GeodeticCalculator gc = new GeodeticCalculator(crs); 
       gc.setStartingPosition(JTS.toDirectPosition(p.getCoordinate(), crs)); 
       gc.setDestinationPosition(JTS.toDirectPosition(dest.getCoordinate(), crs)); 
       // Calculate distance between points 
       dist = gc.getOrthodromicDistance(); 
       // System.out.println(feature.getID()+": "+dist); 
       if (dist < minDist) { 
        minDist = dist; 
        minDistPoint = dest.getCoordinate(); 
        lastMatched = feature; 
       } 
      } 

     } finally { 
      itr.close(); 
     } 
     Point ret = null; 

     if (minDistPoint == null) { 
      searchDist *= 2.0; 
      System.out.println("repeat search"); 
     } else { 
      ret = gf.createPoint(minDistPoint); 
      return ret; 
     } 
    } 
    return gf.createPoint(new Coordinate()); 
}