基于给定起始经纬度与比例尺，将二维数组转换为经纬度点

前提假设 #

首先得明确数组索引和地理边界的对应关系（这是转换的核心，要是实际情况不一样，调整参数就行）：

• 数组的data[row][col]中：
  • row=719, col=0 → 左下角（LL）：纬度LL_lat=28.33346126086149，经度LL_lon=30.757889519059812
  • row=719, col=719 → 右下角（LR）：纬度LR_lat=28.633821527352005，经度LR_lon=39.242110331928565
  • row=0, col=0 → 左上角（UL）：需要你补充这个位置的纬度值UL_lat；如果是规则矩形网格，左上角和右上角的纬度应该相等
  • row=0, col=719 → 右上角（UR）：经度和右下角一致，纬度和左上角一致（规则矩形假设）

情况1：规则矩形网格（南北边缘平行） #

如果你的网格是标准矩形（同一行的纬度相同，同一列的经度相同），可以用简化的线性计算：

1.1 计算经度 #

经度从最西端（LL_lon）到最东端（LR_lon）均匀分布在720列上：

private static double getLongitude(int col) {
    int totalCols = 720;
    double LL_lon = 30.757889519059812;
    double LR_lon = 39.242110331928565;
    // 720个点对应719个间隔，计算每列的经度步长
    double lonStep = (LR_lon - LL_lon) / (totalCols - 1.0);
    return LL_lon + col * lonStep;
}

1.2 计算纬度 #

假设最北端（row=0）的纬度是UL_lat，最南端（row=719）的纬度是LL_lat，则：

private static double getLatitude(int row, double UL_lat) {
    int totalRows = 720;
    double LL_lat = 28.33346126086149;
    // 计算每行的纬度步长
    double latStep = (UL_lat - LL_lat) / (totalRows - 1.0);
    // row=0对应最高纬度，row=719对应最低纬度
    return UL_lat - row * latStep;
}

情况2：梯形网格（边缘不平行，双线性插值） #

如果右下角和左下角的纬度确实不同（网格是梯形），就得用双线性插值来计算任意索引的坐标：

核心思路 #

先在水平方向（同一行）插值得到南北边缘对应列的纬度，再在垂直方向（同一列）插值得到当前行的最终纬度，经度则直接按水平线性插值计算。

代码实现 #

public class GridToGeo {
    // 已知的边界坐标
    private static final double LL_LAT = 28.33346126086149;
    private static final double LL_LON = 30.757889519059812;
    private static final double LR_LAT = 28.633821527352005;
    private static final double LR_LON = 39.242110331928565;
    // 请替换成实际的左上角、右上角纬度值
    private static final double UL_LAT = 29.33346126086149; 
    private static final double UR_LAT = 29.633821527352005; 
    private static final int GRID_SIZE = 720;

    public static double[] getCoordinates(int row, int col) {
        if (row < 0 || row >= GRID_SIZE || col < 0 || col >= GRID_SIZE) {
            throw new IllegalArgumentException("索引超出网格范围");
        }

        // 计算归一化参数（范围0到1）
        double rowNorm = (double) row / (GRID_SIZE - 1); // row=0对应北端（0），row=719对应南端（1）
        double colNorm = (double) col / (GRID_SIZE - 1); // col=0对应西端（0），col=719对应东端（1）

        // 1. 计算经度（水平线性插值）
        double longitude = LL_LON + colNorm * (LR_LON - LL_LON);

        // 2. 计算南北边缘对应列的纬度
        double southLat = LL_LAT + colNorm * (LR_LAT - LL_LAT); // 南边缘当前列的纬度
        double northLat = UL_LAT + colNorm * (UR_LAT - UL_LAT); // 北边缘当前列的纬度

        // 3. 垂直插值得到当前行的纬度
        double latitude = northLat + rowNorm * (southLat - northLat);

        return new double[]{latitude, longitude};
    }

    // 测试示例
    public static void main(String[] args) {
        double[] llCoords = getCoordinates(719, 0);
        System.out.println("左下角坐标：lat=" + llCoords[0] + ", lon=" + llCoords[1]);

        double[] lrCoords = getCoordinates(719, 719);
        System.out.println("右下角坐标：lat=" + lrCoords[0] + ", lon=" + lrCoords[1]);
    }
}

关键提示 #

• 如果你的数组索引对应关系相反（比如row=0是最南端），只需要调整rowNorm的计算方式，改成rowNorm = 1 - (double)row/(GRID_SIZE-1)就行。
• 补充左上角的坐标后，梯形插值的结果会更精准；如果是严格的矩形网格，直接用情况1的简化计算效率更高。

内容的提问来源于stack exchange，提问作者Ben Beri