使用D3地理投影时，如何追踪自适应重采样的投影前经纬度值？

嘿，这个问题我之前折腾过好一阵子！D3的geoPath内部的重采样逻辑确实有点“藏着掖着”，直接拦截投影的point方法只能拿到线段端点，因为中间的采样点是在路径生成阶段悄悄生成的，根本不会走投影的point方法提前暴露原始经纬度。不过有两个靠谱的办法能解决这个问题：

方法一：手动生成可控的大圆弧经纬度采样点 #

最简单直接的方式是用d3.geoInterpolate自己生成大圆弧上的所有原始经纬度点，完全绕开D3内部的黑箱逻辑。你可以自由控制采样密度，还能直接拿到所有投影前的坐标：

// 定义两个城市的经纬度 [经度, 纬度]
const beijing = [116.4074, 39.9042];
const newYork = [-74.0060, 40.7128];

// 创建大圆弧插值函数
const arcInterpolator = d3.geoInterpolate(beijing, newYork);

// 生成50个分段的采样点（最终得到51个坐标）
const sampleCount = 50;
const rawArcPoints = Array.from({ length: sampleCount + 1 }, (_, index) => {
  const t = index / sampleCount; // 从0到1的插值参数
  return arcInterpolator(t);
});

// rawArcPoints就是完整的投影前经纬度序列，直接用它画路径或者做其他处理都没问题

这种方法的好处是完全可控，你可以根据地图缩放级别动态调整采样数量，比依赖D3内部的默认精度要灵活得多。

方法二：Hook D3内部的重采样逻辑（进阶玩法） #

如果你想复用D3根据当前视图自动调整采样密度的逻辑（比如缩放时自动增减采样点），可以通过扩展geoPath来捕获内部生成的原始经纬度点。核心思路是在geoStream处理地理数据时，手动触发D3的重采样逻辑：

const orthoProjection = d3.geoOrthographic();
const customGeoPath = d3.geoPath().projection(orthoProjection);

// 重写geoPath的路径生成逻辑，捕获重采样后的原始经纬度
customGeoPath._path = function(geoObject) {
  const resampledPoints = [];
  const precision = this.precision(); // 获取D3默认的采样精度（默认0.5度）

  // 定义重采样单个线段的函数
  const resampleSegment = (start, end) => {
    const segmentLength = d3.geoLength({ type: "LineString", coordinates: [start, end] });
    const sampleNum = Math.ceil(segmentLength / precision); // 计算需要的采样点数量
    return Array.from({ length: sampleNum + 1 }, (_, i) => 
      d3.geoInterpolate(start, end)(i / sampleNum)
    );
  };

  // 处理LineString类型的地理数据（大圆弧一般用这个类型）
  if (geoObject.type === "LineString") {
    for (let i = 0; i < geoObject.coordinates.length - 1; i++) {
      const segmentPoints = resampleSegment(geoObject.coordinates[i], geoObject.coordinates[i+1]);
      resampledPoints.push(...segmentPoints);
    }
    // 现在resampledPoints就是D3内部逻辑生成的原始经纬度采样点
    console.log("重采样后的原始经纬度：", resampledPoints);
  }

  // 最后调用原始的路径生成方法，保证正常绘图
  return d3.geoPath.prototype._path.call(this, geoObject);
};

这个方法需要对D3的geoStream和路径生成机制有一定了解，适合需要和D3内置采样逻辑对齐的场景。

为什么拦截投影的point方法没用？ #

你之前遇到的问题，本质是D3的执行顺序：geoPath会先对原始线段进行重采样（生成中间的经纬度点），然后才会把每个采样点传给投影的point方法转换为屏幕坐标。但这个重采样过程是在geoPath内部完成的，并没有把中间的原始经纬度暴露给投影的point方法，所以你只能拿到最终的屏幕坐标，看不到投影前的原始点。

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