[论文解读] Weighted Least Squares Fit of an Ellipse to Describe Complete or Spotty Diffraction Rings on a Planar 2D Detector

Michael Hart, Michael Drakopoulos|arXiv (Cornell University)|Nov 21, 2013

X-ray Diffraction in Crystallography参考文献 1被引用 3

一句话总结

本文提出了一种针对二维探测器上X射线粉末衍射环的加权最小二乘椭圆拟合方法，利用峰强度作为权重，能够精确建模完整环和斑点环。该方法在完整环与斑点环的拟合椭圆之间实现了亚10 μm的定位一致性，从而在ΔQ/Q水平达到10⁻⁵的高精度校准X射线束能量与探测器几何参数。

ABSTRACT

We present a procedure for fitting an ellipse to powder diffraction patterns recorded on a planar 2D detector, which uses the peak intensities as weights. This procedure works for complete and spotty diffraction rings. We outline all the steps required: the interpolation of Cartesian pixel data into polar co-ordinates based upon an approximate ellipse centre; the identification of peak locations along each radial direction of a diffraction ring; and the weighted least squares algebraic fit of an ellipse to that diffraction ring. The performance of the procedure is assessed by characterising a complete diffraction ring from a powder standard, and a spotty ring generated from a complete ring of the standard. The fit of an ellipse to a spotty ring from an alkane complex is also shown. Agreement in the calculated parameters for the two rings is good. This procedure is implemented as part of an X-ray beam energy and 2D detector calibration routine which is currently in use at Beamline I12 (Diamond Light Source, UK).

研究动机与目标

• 开发一种稳健的方法，用于对平面二维探测器上的完整和斑点X射线粉末衍射环进行椭圆拟合。
• 通过根据峰强度加权拟合过程，提高斑点环的拟合精度，从而降低低强度噪声的影响。
• 利用椭圆参数实现对X射线束能量、样品-探测器距离和探测器方向的精确校准。
• 通过拟合至峰位置而非原始像素强度，解决衍射环重叠的挑战。
• 提供一种适用于同步辐射光束线数据处理流程中常规使用的几何与解析方法。

提出的方法

• 通过衍射图像自卷积的局部极大值寻找近似环中心，采用降采样以提高计算效率。
• 使用双线性插值将笛卡尔坐标系下的探测器像素数据插值到极坐标系，其中方位角方向的射线均匀分布（例如，完整环R = 720，斑点环R最高可达3000）。
• 通过检测径向强度分布中的局部极大值识别峰位置，随后对所有射线进行直方图分析，以分离出主导的峰位置。
• 沿每条射线对强度分布拟合高斯函数，以精确确定径向和方位角方向的峰中心。
• 对峰中心坐标应用加权最小二乘代数拟合，权重来源于拟合得到的峰强度，以突出高置信度数据点。
• 该方法通过拟合孤立峰位置而非环形区域内的原始像素强度，避免了重叠环带来的偏移影响。

实验结果

研究问题

• RQ1 加权最小二乘椭圆拟合能否准确描述二维探测器上完整和斑点衍射环？
• RQ2 强度加权在存在噪声或缺失数据的情况下，如何提升椭圆拟合的鲁棒性与准确性？
• RQ3 当环为斑点状时，拟合椭圆在多大程度上能保持原始环的几何参数（如长半轴、中心）？
• RQ4 该方法能否在真实同步辐射实验中可靠地用于束流能量与探测器校准？
• RQ5 与直接对环形区域像素强度进行椭圆拟合的替代方法相比，该方法有何优势？

主要发现

• 完整衍射环（10阶CeO₂环）及其对应斑点环的拟合椭圆长半轴分别为1014.05像素和1014.06像素，表明结果高度一致。
• 完整环与斑点环的椭圆中心坐标分别为(1434.21, 1441.93)像素和(1434.27, 1441.95)像素，定位一致性达到10 μm量级。
• 尽管斑点环的总强度仅为完整环的17%，其拟合椭圆参数仍保持高度一致，表明对低信号水平具有强鲁棒性。
• 该方法成功地对来自NBD复合物的斑点环拟合出单一椭圆，即使斑点不对应单一Q值，椭圆仍能代表峰位置的集体分布。
• 采用该椭圆拟合方法的校准流程可实现ΔQ/Q ≈ 10⁻⁵的精度，支持高精度的二维到一维数据还原。
• 该方法优于直接对环形区域像素强度进行椭圆拟合的替代方法，后者可能因环重叠而产生偏移。