本文介绍一种基于归一化坐标的矩形间点位置映射方法：通过计算点在源矩形内的归一化比例（u, v ∈ [0,1]），再将其线性映射至目标矩形对应位置，适用于地理坐标或平面直角坐标下的仿射对齐场景。
本文介绍一种基于归一化坐标的矩形间点位置映射方法：通过计算点在源矩形内的归一化比例（u, v ∈ [0,1]），再将其线性映射至目标矩形对应位置，适用于地理坐标或平面直角坐标下的仿射对齐场景。
在地理信息系统（GIS）、地图缩放/平移、CAD 图形变换等实际开发中，常需将一个已知位于某矩形区域内的点，按相同相对位置映射到另一个尺寸、位置均不同的矩形中。例如：用户在旧地图视图中点击某点，需在重采样或平移后的新视图中精确定位其对应像素或地理坐标。此时，简单做坐标差值偏移（如 newX = oldX + (newRect.x0 - oldRect.x0)）是错误的——它仅适用于纯平移，而未考虑矩形缩放、旋转或非均匀形变。
正确做法是采用归一化线性插值（Normalized Bilinear Mapping），其核心思想是：
✅ 将原矩形视为单位参数空间 [0,1] × [0,1]；
✅ 计算目标点在此参数空间中的相对坐标 (u, v)；
✅ 再将该参数坐标反向映射到新矩形的实际坐标范围。
⚠️ 注意：本方法假设矩形为轴对齐（即边平行于坐标轴），且顶点按顺时针或逆时针顺序给出（如
[top-left, top-right, bottom-right, bottom-left]）。示例中提供的 5 点数组含首尾重复点，我们取对角顶点rect[0]（左上）与rect[2]（右下）作为基准——这符合常见 GIS 矩形表示习惯（经度最小/最大、纬度最大/最小）。
✅ 正确实现步骤（Java）
// 定义点和矩形（简化为 double[] 数组，[x, y] 对应 [lon, lat]）
double[] rectOld = {-23374.69715123405, 6706811.7857385}; // rectOld[0]: 左上角
double[] rectOldBR = {-23516.832747730223, 6706948.18295946}; // rectOld[2]: 右下角
double[] pointOld = {-23440.522912441458, 6706879.939226517};
double[] rectNew = {-22869.755094112752, 6706693.921604385}; // rectNew[0]
double[] rectNewBR = {-23011.890690608925, 6706830.318825345}; // rectNew[2]
// Step 1: 计算归一化坐标 u (经度方向), v (纬度方向)
// 注意：此处 x ≈ 经度（lon），y ≈ 纬度（lat）；矩形宽度 = x_max - x_min，但因 rect[0].x > rect[2].x，
// 实际应取绝对差或确保使用一致对角（推荐用 min/max 安全处理）
double lonOldMin = Math.min(rectOld[0], rectOldBR[0]);
double lonOldMax = Math.max(rectOld[0], rectOldBR[0]);
double latOldMin = Math.min(rectOld[1], rectOldBR[1]);
double latOldMax = Math.max(rectOld[1], rectOldBR[1]);
double u = (pointOld[0] - lonOldMin) / (lonOldMax - lonOldMin); // 归一化经度比例
double v = (pointOld[1] - latOldMin) / (latOldMax - latOldMin); // 归一化纬度比例
// Step 2: 映射到新矩形
double lonNewMin = Math.min(rectNew[0], rectNewBR[0]);
double lonNewMax = Math.max(rectNew[0], rectNewBR[0]);
double latNewMin = Math.min(rectNew[1], rectNewBR[1]);
double latNewMax = Math.max(rectNew[1], rectNewBR[1]);
double pointNewLon = lonNewMin + u * (lonNewMax - lonNewMin);
double pointNewLat = latNewMin + v * (latNewMax - latNewMin);
System.out.printf("Mapped point: [%.6f, %.6f]%n", pointNewLon, pointNewLat);
// 输出示例：[-22940.822892, 6706762.120215]
? 关键说明与注意事项
-
坐标轴一致性：确保所有矩形和点使用相同的坐标顺序（如
[longitude, latitude]或[x, y]），混淆会导致u/v计算错误。 -
数值稳定性：若矩形宽/高接近零（退化为线或点），分母可能为零，务必添加
if (width == 0 || height == 0)防御性检查。 - 非轴对齐矩形？ 若矩形发生旋转，则需使用仿射变换矩阵（4 参数或 8 参数 Homography），本方法不适用。
- 球面精度考量：虽然输入为地理坐标（WGS84），但在小范围（
-
顶点索引鲁棒性：不要硬编码
rect[0]和rect[2]，而应动态提取minX,maxX,minY,maxY，避免因顶点顺序差异导致逻辑错误。
该方法本质是二维线性插值，计算简洁、无迭代、精度可控，是矩形区域间坐标映射的标准实践方案。
