如何将仅含NAXIS头的FITS图像中DAOStarFinder输出坐标转ra/dec

完全可以做xcentroid、ycentroid到RA/Dec的转换，但受限于当前掌握的已知信息，转换结果是带系统误差的近似值，无法达到专业天文定标的精度。

你需要先补全的必要参数 #

你目前已经掌握的锚点信息（中心像素位置、对应中心天球坐标、图像尺寸）不足以支撑精确转换，还需要确定两类核心参数：

• 图像的像素角比例尺：即单个像素在天球上对应的角跨度，通常需要分别明确赤经、赤纬两个方向的单像素角度值，单位常用角秒/像素
• 图像的方位标定参数：包括图像坐标轴相对天球北方向的旋转角度、是否存在光学切向畸变、像素轴正方向与赤经/赤纬方向的对应关系（比如x轴正方向指向东还是西、y轴正方向指向北还是南）

小视场场景下的近似转换流程 #

你拍摄的是M52星团的2000×2000像素图像，这类常规观测属于小视场范围，天球的球面投影畸变可以忽略，可直接用平面近似完成简化转换，步骤为：

1. 以给定的中心锚点为基准：像素(1000,1000)对应天球坐标RA=351.2°, Dec=61.5933°
2. 对DAOStarFinder输出的每个源，计算它和中心像素的相对偏移：
   • 水平偏移量：dx = xcentroid - 1000
   • 垂直偏移量：dy = ycentroid - 1000
3. 代入选定的角比例尺，结合中心赤纬做投影修正（赤经方向的实际角距离需要乘以cos(Dec_center)，抵消高赤纬区域的天球投影压缩效应），再根据图像的实际朝向调整偏移量的正负号，得到每个源相对中心的天球坐标差
4. 把坐标差叠加到中心锚点的RA/Dec上，即可得到对应源的近似天球坐标

实际使用时的注意事项 #

• 因为当前FITS头除了NAXIS没有任何WCS相关字段，所有转换需要的比例尺、方位参数都需要从外部补充（比如同设备的定标结果、从已知M52的角直径反推的估算比例尺），这套近似方法的误差完全取决于你补入的参数准确度，如果没有经过实地校准，转换得到的坐标可能存在几十角秒甚至更大的系统偏差
• 如果需要达到和外部星表交叉匹配的精度，建议从DAOStarFinder识别的源中挑选几颗亮星，和已知的亮星位置做交叉匹配，拟合出完整的世界坐标系（WCS）参数后再做坐标转换，不要直接使用无校准的近似结果
• 计算前一定要确认DAOStarFinder输出的像素坐标的索引规则：部分天文处理软件的像素坐标以图像最左上角为(0,0)（0索引），部分以最左上角为(1,1)（1索引），如果和你已知的中心像素的坐标系定义不一致，会引入整像素级的固定偏移误差

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