0

我正在尝试为整个世界生成矩形网格，单元格大小为 350 米（宽度和高度相等）。

我试图做的事情：

1) 在主循环中迭代 [-180 范围内的槽经度；180]

2) 在嵌套循环中迭代 [-85.06 范围内的低谷纬度；85.06]（Web 墨卡托边界）

3) 使用 Java GeographicLib Geodesic.WGS84.Direct()在东（在主循环中）和北（在嵌套循环中）增加绑定到 350 米的单元格。

Java代码：

double lon_min = -180;
double lon_max = 180;
double lat_min = -85.06;
double lat_max = 85.06;

double lon_max_c = lon_min;
double lon_min_c = lon_min;

int grid_x = 0;
int grid_y = 0;
int cell_id = 0;

while(lon_max_c < lon_max)
{
    double lat_min_c = lat_min;
    double lat_max_c = lat_min;

    lon_min_c = lon_max_c;
    GeodesicData g = Geodesic.WGS84.Direct(lat_min_c, lon_min_c, 90, 350);

    boolean isXlast = g.lon2 > lon_max || (lon_min_c > 0 && g.lon2 < 0);
    lon_max_c = isXlast? lon_max : g.lon2;

    grid_y = 0;

    while(lat_max_c < lat_max)
    {
        lat_min_c = lat_max_c;

        GeodesicData g1 = Geodesic.WGS84.Direct(lat_min_c, lon_min_c, 0, 350);

        boolean isYlast = g1.lat2 > lat_max;
        lat_max_c = isYlast? lat_max : g1.lat2;

        System.out.print("\rid: " +  cell_id + " lon: " + lon_max_c + " lat: " + lat_max_c);

        grid_y = grid_y + 1;
        cell_id = cell_id + 1;
    }
    grid_x = grid_x + 1;
}

结果，产生了约 5.3 亿个细胞。然而，整个星球的面积约为 5.1 亿平方公里。所以，细胞总数应该是40亿左右。我对测地线数学很陌生，所以出了什么问题？我应该怎么做才能正确生成这个网格？

4

1 回答 1

0

整个世界的矩形网格，单元格大小固定为 350 米（宽度和高度相等）

你想要做的事情是不可能的，因为地球不是一个矩形。有关了解更多信息的起点，请参阅https://en.wikipedia.org/wiki/Map_projection 。

虽然对于小区域可以忽略这种影响，但在“范围 [-180; 180] 范围内的经度”、“范围 [-85.06; 85.06] 范围内的纬度”等尺度上忽略它是不可行的。

于 2018-11-27T11:00:49.660 回答