Question

我从雷达获得了很多坐标，这些坐标格式与谷歌地图坐标非常相似，我不太确定，但我问公司，他们告诉我为了获得阿尔伯投影的坐标，我需要做以下事情：

然后，您将能够创建2个CoordinateReferenceSystem（即坐标系），设置一个默认值（即lon / lat），并设置另一个是WKT字符串（将投影x / y）。那么你可以轻松创建2个MathTransform进行转换方向。

这是Alber预测的OGC WKT：

PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Albers_Conic_Equal_Area"],
    PARAMETER["standard_parallel_1",31.996308],
    PARAMETER["standard_parallel_2",33.996308],
    PARAMETER["latitude_of_center",32.996308],
    PARAMETER["longitude_of_center",35.415901],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,AUTHORITY["EPSG","9001"]]]

所以从我得到的是，我需要从长/ lat转换到WKT投影，以显示在Alber的投影地图图像上。

所以在GeoTools我使用了以下代码：

    CoordinateReferenceSystem source = CRS.decode("EPSG:4326");
    CoordinateReferenceSystem target = CRS.parseWKT("PROJCS[\"unnamed\", GEOGCS[\"WGS 84\", DATUM[\"WGS_1984\", SPHEROID[\"WGS 84\",6378137,298.2572235629972, AUTHORITY[\"EPSG\",\"7030\"]], AUTHORITY[\"EPSG\",\"6326\"]], PRIMEM[\"Greenwich\",0], UNIT[\"degree\",0.0174532925199433], AUTHORITY[\"EPSG\",\"4326\"]], PROJECTION[\"Albers_Conic_Equal_Area\"], PARAMETER[\"standard_parallel_1\",31.996308], PARAMETER[\"standard_parallel_2\",33.996308], PARAMETER[\"latitude_of_center\",32.996308], PARAMETER[\"longitude_of_center\",35.415901], PARAMETER[\"false_easting\",0], PARAMETER[\"false_northing\",0], UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]]]");
    MathTransform transform = CRS.findMathTransform(source, target, true);

    Coordinate c = JTS.transform(new Coordinate(34, 35), new Coordinate(), transform);
    System.out.println(c.toString());

这就是我得到的输出：

(-38422.86847540497, 111410.0483012808, NaN)

现在，可能是因为错误的源坐标系，但他的默认长/拉系统是什么意思？

即使我解决了这个问题，我怎样才能让它在我的地图图像上显示这些点？我的意思是它必须知道图像的宽度/高度不是吗？

Answer 1

这个表示WKT投影从EPSG 4326投影到距离零点的距离（参见参数）。我们怎么知道它的米？因为WKT有一个设置的UNIT标签，显示为米单位。

那我是怎么用的？

我公司给了我两个文件，地图JPEG文件和地图常量的XML文件。

该地图常量XML文件包含从该地图的零点到地图角落的距离。因此，如果您在地图上至少有一个点，那么您可以找到所有内容。

将其转换为地图X / Y需要知道的事情：

每个像素多少公里（在我的情况下是1.6正好1英里）
地图上距离零点的距离上的至少一个已知点

我就这样做了：

        MathTransform transform = CRS.findMathTransform(epsg4326, targetWKT, true);

        DirectPosition2D srcDirectPosition2D
                = new DirectPosition2D(epsg4326, latitude, longitude);
        DirectPosition2D destDirectPosition2D
                = new DirectPosition2D();
        transform.transform(srcDirectPosition2D, destDirectPosition2D);

        double transX = destDirectPosition2D.x;
        double transY = destDirectPosition2D.y;

        int kmPerPixel = mapImage.getWidth / 1024; // It is known to me that my map is 1024x1024km ...

        double x = zeroPointX + ( (transX * 0.001) * kmPerPixel);
        double y = zeroPointY + ( ( (transX * -1) * 0.001) * kmPerPixel);

你可以通过相同的计算得到零点X和Y，而不是添加，就像我在地图上的距离一样，就像我在角落里做的那样。

可以帮助一些人，所以我发布了我想出来的东西。

将经度和纬度转换为Alber投影上的X，Y像素坐标图像

1 个答案: