Как я могу использовать GeoTools / ProJ.4 (или другие API), чтобы преобразовать координаты из HK80 GRID в широту / долготу? - PullRequest
Новая
       53

Как я могу использовать GeoTools / ProJ.4 (или другие API), чтобы преобразовать координаты из HK80 GRID в широту / долготу?

0 голосов
/ 27 декабря 2011

как конвертировать сетку hk80 в лат / долг?

HK 1980 Grid равен EPSG 2326, в то время как Lat / Long WGS84 - EPSG 4326

.

Или есть что-нибудь более подходящее для этого

1 Ответ

4 голосов
/ 17 августа 2012

Я надеюсь, что следующий код может вам помочь.Мне нужно использовать GeoTools для достижения вашей цели.Кроме того, вам нужно вставить соответствующий север и восток для epsg 2326. Однако я все еще не могу решить проблему исключений.Чтобы правильно использовать этот код, вам необходимо загрузить и импортировать все файлы GeoTools, кроме gt-epsg-postgresql .

package org.geotools.tutorial;

import java.awt.geom.Point2D;
import java.io.File;

import org.geotools.data.FileDataStore;
import org.geotools.data.FileDataStoreFinder;
import org.geotools.data.simple.SimpleFeatureSource;
import org.geotools.factory.Hints;
import org.geotools.map.FeatureLayer;
import org.geotools.map.Layer;
import org.geotools.map.MapContent;
import org.geotools.referencing.CRS;
import org.geotools.referencing.ReferencingFactoryFinder;
import org.geotools.styling.SLD;
import org.geotools.styling.Style;
import org.geotools.swing.JMapFrame;
import org.geotools.swing.data.JFileDataStoreChooser;
import org.opengis.referencing.crs.CRSAuthorityFactory;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;
import org.geotools.referencing.*;
/**
 * Prompts the user for a shapefile and displays the contents on the screen in a map frame.
 * <p>
 * This is the GeoTools Quickstart application used in documentationa and tutorials. *
 */
public class quickstart {

/**
 * GeoTools Quickstart demo application. Prompts the user for a shapefile and displays its
 * contents on the screen in a map frame
 */
    public static void main(String[] args) throws Exception {
        //Point2D.Double srcProjec = new Point2D.Double (836694.050, 819069.800);
        //Point2D.Double dstProjec = new Point2D.Double (132,37);
        // display a data store file chooser dialog for shapefiles
        //Hints hints = new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE);
        //CRSAuthorityFactory factory = ReferencingFactoryFinder.getCRSAuthorityFactory("EPSG", hints);
        //CoordinateReferenceSystem crs = factory.createCoordinateReferenceSystem("EPSG:4326");
        CRSAuthorityFactory   factory = CRS.getAuthorityFactory(true);
        CoordinateReferenceSystem srcCRS = factory.createCoordinateReferenceSystem("EPSG:2326");
        CoordinateReferenceSystem dstCRS = factory.createCoordinateReferenceSystem("EPSG:4326");
        boolean lenient = true; // allow for some error due to different datums
        MathTransform transform = CRS.findMathTransform(srcCRS, dstCRS, lenient);
        double[] srcProjec = {818039, 836361};// easting, northing, 
        double[] dstProjec = {0, 0};
        transform.transform(srcProjec, 0, dstProjec, 0, 1);
        System.out.println("longitude: " + dstProjec[0] + ", latitude: " + dstProjec[1]);

    }

}
...