如何生成可用于d3-geo的GeoJSON格式椭圆多边形？

刚好能解决你需要的类似d3.geoCircle()的椭圆生成需求，核心思路是通过计算椭圆上的经纬度采样点，再组装成标准GeoJSON Polygon结构，下面是具体实现和说明：

核心思路 #

GeoJSON的多边形本质是闭合的经纬度坐标数组，所以我们只需要：

• 基于给定的中心、半轴长度、旋转角度，生成椭圆上的一系列采样点
• 将这些点的笛卡尔坐标（米为单位）转换回经纬度
• 封装成符合GeoJSON规范的Polygon对象

基于Cesium的实现方案（推荐） #

既然你提到了Cesium的EllipseGeometry，可以直接利用Cesium的坐标转换和几何计算能力，写出一个返回GeoJSON的函数，准确性和易用性都很高：

function geoEllipse(options) {
  // 解构参数，设置默认值
  const {
    center, // [经度, 纬度]
    semiMajorAxis, // 半长轴（米）
    semiMinorAxis, // 半短轴（米）
    rotation = 0, // 旋转角度（度，默认0）
    steps = 64 // 采样点数量，越多椭圆越平滑
  } = options;

  const [centerLon, centerLat] = center;
  // 将中心经纬度转为WGS84椭球下的笛卡尔坐标
  const centerCartesian = Cesium.Cartesian3.fromDegrees(centerLon, centerLat);
  
  // 获取中心位置的东、北方向向量，并缩放至半轴长度
  const enuFrame = Cesium.Ellipsoid.WGS84.eastNorthUpToFixedFrame(centerCartesian);
  const eastAxis = Cesium.Cartesian3.normalize(
    Cesium.Cartesian3.multiplyByScalar(enuFrame.east, semiMajorAxis, new Cesium.Cartesian3()),
    new Cesium.Cartesian3()
  );
  const northAxis = Cesium.Cartesian3.normalize(
    Cesium.Cartesian3.multiplyByScalar(enuFrame.north, semiMinorAxis, new Cesium.Cartesian3()),
    new Cesium.Cartesian3()
  );

  // 对轴向量应用旋转角度
  const rotatedEast = Cesium.Cartesian3.rotateAround(
    eastAxis,
    centerCartesian,
    Cesium.Cartesian3.UNIT_Z,
    Cesium.Math.toRadians(rotation),
    new Cesium.Cartesian3()
  );
  const rotatedNorth = Cesium.Cartesian3.rotateAround(
    northAxis,
    centerCartesian,
    Cesium.Cartesian3.UNIT_Z,
    Cesium.Math.toRadians(rotation),
    new Cesium.Cartesian3()
  );

  // 生成椭圆上的采样点
  const coordinates = [];
  for (let i = 0; i <= steps; i++) {
    const angle = (i / steps) * 2 * Math.PI;
    // 椭圆参数方程计算当前点的偏移量
    const xOffset = Math.cos(angle);
    const yOffset = Math.sin(angle);

    // 计算当前点的笛卡尔坐标
    const pointCartesian = Cesium.Cartesian3.add(
      centerCartesian,
      Cesium.Cartesian3.add(
        Cesium.Cartesian3.multiplyByScalar(rotatedEast, xOffset, new Cesium.Cartesian3()),
        Cesium.Cartesian3.multiplyByScalar(rotatedNorth, yOffset, new Cesium.Cartesian3()),
        new Cesium.Cartesian3()
      ),
      new Cesium.Cartesian3()
    );

    // 转换回经纬度
    const cartographic = Cesium.Cartographic.fromCartesian(pointCartesian);
    coordinates.push([
      Cesium.Math.toDegrees(cartographic.longitude),
      Cesium.Math.toDegrees(cartographic.latitude)
    ]);
  }

  // 返回标准GeoJSON Polygon
  return {
    type: "Polygon",
    coordinates: [coordinates] // Polygon的coordinates是数组的数组，外层是环，内层是点
  };
}

和d3-geo结合使用 #

生成的GeoJSON可以直接用d3.geoPath()渲染，示例代码：

// 创建椭圆GeoJSON
const myEllipse = geoEllipse({
  center: [-75.59777, 40.03883],
  semiMajorAxis: 500000,
  semiMinorAxis: 300000,
  rotation: 60,
  steps: 80 // 增加采样点让椭圆更平滑
});

// 用d3绘制
const svg = d3.select("svg");
const pathGenerator = d3.geoPath();

svg.append("path")
  .datum(myEllipse)
  .attr("d", pathGenerator)
  .attr("fill", "rgba(0, 120, 255, 0.2)")
  .attr("stroke", "#0078ff")
  .attr("stroke-width", 2);

无Cesium依赖的替代思路 #

如果不想引入Cesium，可以自己实现球面坐标转换：

1. 将中心经纬度转为球面坐标（弧度）
2. 利用球面三角学计算每个采样点相对于中心的偏移，得到经纬度
3. 组装成GeoJSON

不过这种方法需要处理更多球面几何细节，小范围椭圆可以用平面近似（将经纬度转成墨卡托投影坐标，生成椭圆后再转回经纬度），但大范围会有明显误差，所以更推荐用Cesium的方案来保证准确性。

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