赤道坐标系下星座边界点数组插值合并实现方法咨询

核心思路说明 #

你之前使用Array.prototype.map达不到预期，本质是因为map是一对一映射的操作，每个输入元素只能对应一个输出元素，而你的场景需要每个边界点（除最后一个）对应输出「当前点+N个中间插值点」多个元素，属于一对多映射场景，优先使用flatMap或reduce实现。
实现时要注意避免相邻点的端点重复，只在每对处理时保留前一个端点，最后单独补最后一个端点即可，不会出现重复或遗漏。
我们先统一定义用到的类型和函数签名方便演示：

// 赤道坐标点类型
type EquatorialCoordinate = {
  ra: number; // 赤经
  dec: number; // 赤纬
}
// 你的插值函数签名，返回两点之间的中间插值点（不含两个输入端点）
declare function getInterpolatedEquatorialCoordinates(
  p1: EquatorialCoordinate,
  p2: EquatorialCoordinate,
  precision: number
): EquatorialCoordinate[]

实现1：使用Array.prototype.flatMap（最简洁，推荐） #

flatMap的逻辑是先执行map操作，再自动把结果数组展平一层，刚好适配一对多映射的场景：

// 入参定义
const originalBoundary: EquatorialCoordinate[] = // 你的IAU官方原始边界坐标数组
const interpolationPrecision = 0.1 // 替换为你实际需要的插值精度参数
// 插值拼接逻辑
const interpolatedBoundary = originalBoundary.flatMap((currentPoint, index, arr) => {
  // 最后一个点直接返回自身，不需要再做插值
  if (index === arr.length - 1) return [currentPoint]
  const nextPoint = arr[index + 1]
  // 返回当前点 + 两点之间的所有中间插值点
  return [currentPoint, ...getInterpolatedEquatorialCoordinates(currentPoint, nextPoint, interpolationPrecision)]
})

实现2：使用Array.prototype.reduce（兼容性更好） #

如果运行环境不支持ES2019的flatMap，可以用reduce实现，逻辑同样清晰，性能和flatMap几乎一致：

const interpolatedBoundary = originalBoundary.reduce((acc, currentPoint, index, arr) => {
  // 先存入当前原始点
  acc.push(currentPoint)
  // 非最后一个点时，追加和下一个点之间的插值结果
  if (index !== arr.length - 1) {
    const nextPoint = arr[index + 1]
    const interpolatedPoints = getInterpolatedEquatorialCoordinates(currentPoint, nextPoint, interpolationPrecision)
    acc.push(...interpolatedPoints)
  }
  return acc
}, [] as EquatorialCoordinate[])

注意事项 #

• 确保你的getInterpolatedEquatorialCoordinates返回的插值点不包含输入的两个端点，否则会出现坐标点重复的问题。如果你的插值函数本身已经包含了端点，需要调整逻辑，只取中间的点或者做去重处理。
• 如果需要处理闭合的星座边界，即最后一个点和第一个点也需要做插值的话，可以先在原始数组末尾追加第一个点，再执行上述插值逻辑，完成后再去掉最后一个重复的点即可。
• 性能层面：如果原始边界点数量很大（比如超过10000个），上述两种实现的性能差异可以忽略，都远优于多次遍历数组、反复拼接的实现方案。

内容的提问来源于stack exchange，提问作者Micheal J. Roberts