质心算法

3.1 质心检测算法

系统采用质心法进行数据处理能提高测试精度。因为质心法能使CCD 上的图像分辨率达到光敏元尺寸的1/10,那么成像亮线中心在CCD 上所对应的光敏源序号就可以是小数,而非一定是整数,这样通过计算可知,精度提高了0.1个百分点。虽然测量系统的精度有提高,但0.11%的相对误差仍不能令人满意,从误差公式可知,系统误差的改善主要取决于CCD 的像元尺寸。随着CCD 技术的不断发,像元尺寸也会不断改善,系统误差也将会有大幅度减小。

质心法图像预处理算法步骤如下[5]：（1）对图像通过灰度化和反色后阈值选择得到光斑特征区域；（2）模糊去噪（mean blur ），消除热噪声以及像素不均匀产生的噪声；（3）再次进行阈值选择，得到更清晰的光斑区域；（4）形态学处理，选择disk 中和合适的领域模板，对图像进行腐蚀和填充处理，以得到连通域的规则形状图形；（5）边缘检测得到图像边缘，反复实验证明canny 边缘检测算法最好；（6）对边缘再进行形态学strel －imerode －imclose －imfill 相关运算得到更连通的边缘曲线，调用regionprops （L ，properties ）函数，根据质心法计算质心。

下面介绍几种常用的质心算法

（1）普通质心算法

(,)ij ij ij

c c

ij

ij x I

x y I =∑∑ （3-1）

其中ij I 为二维图像上每个像素点所接收到的光强，该算法适用于没有背景噪声，背景噪声一致或信噪比较高的情况。

（2）强加权质心算法

0000000000000000,/2,/2

,/2,/2

,/2,/2

,/2,/2y w y x w x i ij j y w y i x w x c y w y x w x ij j y w y i x w x x I w x I w ++=-=-++=-=-=∑∑∑∑

0000000000000000,,22,,22

,,22,,22y

x x W y W j ij x y i x W j y W c y x x W y W ij x y i x W j y W y I w y I w ++=-=-++=-=-=∑∑∑∑

（3-2）

该算法中的加权函数包含3种形式，即a W I =，W I P =+,W I P =⨯其中a p 和均为强度值。其原理是将光斑中心较近的部分区域像素值增强，使得对光斑影响较大点的影响力进一步加大，提高质心探测精度，即使光斑形状不近似高斯分布，该方法依然能较准确的探测光斑质心。

（3）阈值质心算法

先对图像进行阈值分割，再使用COG 算法，此时的质心探测误差要小于只用COG 算法时的质心探测误差。该算法关键在阈值的选取。选取阈值的方法有很多。常用的有固定值阈值，经验阈值和直方图阈值。目前应用较广泛的是自适应阈值。自适应阈值具有阈值随图像的变化而变化以达到最好效果的优点。自适应阈值中最常用的是小波自适应阈值。该方法可以很好的去除高斯白噪声。但传统的小波变换在信号降噪中计算不足，所以可以将小波降噪算法进行改进。

（4）距离质心算法

i ij ij ij

c ij ij

ij x I W x I W =∑∑

ij ij ij ij

c ij ij

ij y I W y I W *=∑∑

11ij W s == （3-3）

其中(,)i j x y 为当前被测像素的坐标，(,)c c x y 为光斑的中心坐标，(,)c

c x y **为计算出的光斑质心坐标，ij I 为当前像素的像素值。由于距离光斑中心越远的像素对质心判断的影响越小，所以取距离的倒数作为加权函数，这样可以减弱远离光

斑的噪声，减少噪声对质心探测的干扰，提高质心的探测精度。当光斑不近似高斯分布时，该算法依然有较高的精确度，能有效的将光斑与背景分离。从而消除噪声干扰。但该算法必须在能准确求出光斑中心坐标的前提下才具有高精度探测质心的能力

[6]。

图像质心即图像灰度的重心，设图像有i ，j 两个方向，m ，n 分别为i ，j 方向像素的数量，(,)g i j 为像素点(,)i j 处的灰度值，则图像质心位置坐标表达

11

11

(,)(,)

n m j i n m j i g i j i

x g i j ====⨯=∑∑∑∑ 11

11(,)(,)

n m j i n m j i g i j j

y g i j ====⨯=∑∑∑∑ （3-4）

传统质心算法

为了便于讨论，星图采用如图3-1所示的坐标系XOY 。图中的实线方格代表图像传感器的像元，一个像元的宽度对应x 轴一个坐标单位，一个像元的高度对应y 轴一个坐标单位，第i 行j 列像元的坐标记为(,)i j x y ，则,i j x i y j ==，其中，i=0,1,2……,j=0,1,2……。

传统质心算法是在工程中经常采用的星点定位方法。它首先对星点的灰度值求面积的矩，然后在该区域内做面积平均。若第i 行j 列像元的坐标为(,)i j x y ，灰度值为(,)i j G x y ，则星点的质心可表示为

2211

2211(,)ˆ(,)i j i i j

i i j j i j i j

i i j j x G x y x G x y =====∑∑∑∑ （3-5）

传统质心算法定位精度的分析

下面以 x 方向为例分析传统质心算法的定位精度，令2

1(,)j i i j j j G G x y ==∑，因

星点像元的灰度值不可避免地存在噪声，即