DAOPHOT中Sharpness在Photutils与Stetson定义差异及计算偏差咨询

1. 对齐H_00的计算逻辑 #

photutils的daofinder中，sharpness依赖归一化的PSF加权二阶矩H_00，而非简单的区域像素均值。你需要：

• 对比自定义代码与photutils的H_00计算步骤：photutils会基于PSF模型对目标区域像素做加权处理，再计算归一化矩，若你误用原始像素均值替代，必然出现结果偏差。
• 可以直接调用photutils内部函数辅助验证：

  from photutils.detection.daofinder import _compute_sharpness
  # 传入你的PSF数据、目标位置等参数，获取标准计算结果
  sharpness_photutils = _compute_sharpness(...)
  

2. 确认峰值定义一致 #

• Stetson原版DAOPHOT用的是目标区域内的原始像素最大值，但photutils的daofinder可能采用PSF拟合后的中心峰值而非原始像素峰值。检查你自定义代码的峰值提取逻辑是否与photutils匹配。

3. 匹配周边区域选取规则 #

• 原版DAOPHOT的“周边像素”是固定半径的环形区域，而photutils可能使用PSF支持范围内的加权区域。核对两者的区域划分（包括半径、是否排除中心区域）是否完全一致。

4. 检查归一化步骤 #

• photutils的sharpness计算可能包含额外归一化操作（比如除以PSF半高宽相关参数），若你的自定义代码遗漏这一步，结果会出现数量级差异。查看photutils源码中sharpness的最终计算式，确认是否存在归一化系数。

5. 统一PSF数据预处理 #

• 确保测试用PSF数据在传入photutils和自定义代码前的预处理完全一致：
  • 是否做了背景扣除？
  • 是否将PSF总通量归一化为1？
  • 像素单位（电子数/ADU）是否统一？

快速定位技巧：打印photutils计算过程中的中间变量（H_00、峰值、加权均值），与自定义代码的中间结果逐一对比，差异点会直接显现。

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