1.一种基于非对称光学系统点扩散函数模型的图像复原方法,其特征在于,所述方法包括以下步骤:步骤一:预处理:对由于光学系统自身的不适应性导致退化的图像 进行预处理,其中M和N为整数,分别表示图像矩阵的行数和列数,预处理包括滤除噪声的干扰以及几何畸变的校正,其中对模糊图像的噪声滤除选择最佳陷波滤波的效果较好,同时根据仿射变换来校正畸变的图像;
步骤二:退化PSF的估计模型:根据构造公式建立可变脉宽VPW的近似PSF估计模型f(x,y),VPW脉冲的函数表达式为:f(t)=fs(t)+fa(t),t∈[0,T) (1)其中
这里fs(t)和fa(t)分别表示VPW函数f(t)的对称部分和非对称部分,T表示函数的时域长度,VPW函数由a、d、t′、r四个参数确定退化过程,而高斯函数仅由一个参数σ来决定图像模糊的平滑程度,相比较而言,VPW函数模型具有更高的灵活性和更广的适用范围以及更能考虑到与实际模糊过程的近似误差;
步骤三:参数选取:根据光学系统的成像具有中心像素值高,向四周扩散过程中像素值逐渐降低的特点选择脉宽r、脉冲时延t′;根据图像退化模型要在形态上近似实际点成像的特点选择合适的幅值a和b;
步骤四:建立VPW估计模型:利用上述步骤的构造公式以及选取的参数来建立VPW的近似PSF估计模型;
步骤五:利用上述步骤建立的VPW模糊函数模型,进而构造相应的3×3或者5×5的PSF卷积核h,构造方法为:在生成一个3×3或者5×5的矩阵后将每个网格的坐标值作为VPW函数模型f(x,y)的坐标值,计算每点的权重并进行元素值的归一化处理后生成最终的3×3或者5×5的PSF卷积核h;
步骤六:选取相应的复原算法,这里选取传统的Lucy-Richardson图像复原算法来恢复图像;
步骤七:去模糊的初始估计:利用上述步骤的卷积核h和初始估计图像g∈RM×N得到去模糊的初始图像F1,变换方法如下:F1=h*g (4)
其中F1表示第1次迭代的去模糊图像的估计,h为PSF卷积核,g∈RM×N为经过预处理得到的模糊图像,符号*表示卷积操作;
步骤八:求取期望因子:利用上一步骤的估计图像Fk和输入的模糊图像g∈RM×N以及PSF卷积核h变换得到估计图像的期望最大化因子,Fk为第k次迭代的去模糊图像的估计,变换方式如公式(5)所示:其中符号·表示相关运算,符号*表示空间卷积;
步骤九:输出复原图像:不断更新迭代的估计图像,倘若达到迭代次数k,其一般为10-
20次迭代,输出最终的去模糊图像,更新方法如公式(6)所示:Fk+1Fk×μ (6)
其中Fk+1为输出的复原图像。