0

我必须将 GPS 坐标（纬度和经度）转换为 Lambert Conformal Conic 投影。地图指的是意大利。由于我不是这个主题的专家，我发现（谷歌搜索）Java 的最佳库之一是 GeoToolkit。不幸的是，我没有找到任何示例或文档（只有 Javadoc）。这里的投影参数（proj=Lcc Lat_2=45 Lat_1=40 Lat_0=42.5 Lon_0=9.2 x_0=600000 y_0=2200000）：

PARAM_MT["Lambert_Conformal_Conic_2SP", 
PARAMETER["semi_major", 6378137.0], 
PARAMETER["semi_minor", 6356752.314], 
PARAMETER["central_meridian", 9.2], 
PARAMETER["latitude_of_origin", 42.5], 
PARAMETER["false_easting", 600000.0], 
PARAMETER["false_northing", 2200000.0], 
PARAMETER["standard_parallel_2", 45.0], 
PARAMETER["standard_parallel_1", 40.0]]

有谁知道比 GeoToolkit 更好的库？

我编写了以下代码来转换坐标：

// 1. Get the MathTransform Factory
MathTransformFactory factory;

FactoryRegistry registry = new FactoryRegistry(
        MathTransformFactory.class);

factory = registry.getServiceProvider(MathTransformFactory.class, null,
        null, Hints.MATH_TRANSFORM_FACTORY);

// 2. Define source and target coord reference systems
GeographicCRS sourceCRS = DefaultGeographicCRS.WGS84;

ParameterValueGroup parameters = factory
        .getDefaultParameters("Lambert_Conformal_Conic_2SP");
parameters.parameter("semi_major").setValue(WGS84_AXIS_MAJOR);
parameters.parameter("semi_minor").setValue(WGS84_AXIS_MINOR);
parameters.parameter("central_meridian").setValue(9.2);
parameters.parameter("latitude_of_origin").setValue(42.5);
parameters.parameter("standard_parallel_1").setValue(40.0);
parameters.parameter("standard_parallel_2").setValue(45.0);
parameters.parameter("false_easting").setValue(600000);
parameters.parameter("false_northing").setValue(2200000);

// 3. Create the MathTransform with the given parameters
Conversion conversion = new DefiningConversion("GPS to Lambert",
        parameters);

CRSFactory crsFactory = FactoryFinder.getCRSFactory(null);
Map<String, ?> properties = Collections.singletonMap(
        ProjectedCRS.NAME_KEY, "LatLon2Lambert");
ProjectedCRS targetCRS = crsFactory.createProjectedCRS(properties,
        sourceCRS, conversion, DefaultCartesianCS.GENERIC_2D);

MathTransform tr = CRS.findMathTransform(sourceCRS, targetCRS);

DirectPosition sourcePt = new DefaultDirectPosition(DefaultGeographicCRS.WGS84);
sourcePt.setOrdinate(0, 45.0);
sourcePt.setOrdinate(1, 8.0);

DirectPosition targetPt = tr.transform(p.getValue(), null);

但是，它不起作用。我哪里错了？提前致谢。

4

1 回答 1

0

提供的代码是正确的。唯一的问题是输入：GeoToolkit 使用位置符号来识别纬度和经度。它要求第一个参数是经度。例如，在上述代码的末尾，您可以执行转换。

下面是一个例子：

DirectPosition sourcePt = new GeneralDirectPosition(DefaultGeographicCRS.WGS84);
sourcePt.setOrdinate(0, Double.parseDouble(8.751426)); // longitude
sourcePt.setOrdinate(1, Double.parseDouble(45.12478)); // latitude

DirectPosition targetPt = tr.transform(p.getValue(), null);

在 targetPt 变量中，您在 Lambert 中有正确的坐标。

于 2012-05-15T06:49:56.260 回答