图像中点扩散函数的获取

在图像恢复技术中, 点扩展函数( PSF) 是影响图像恢复结果的关键因素, 所以常常利用先验知识和后验判断方法估计PSF函数来恢复图像。

#include 
#include 
#include 
//#include “cdjpeg.h“		/* Common decls for cjpeg/djpeg applications */
#include “basic_op.h“

//#define NULL 0
/********************************
 * roger: some basic operations *
 ********************************/

/* 确定运动模糊角度的两个函数 */
void matrix_orient（float angledouble *motion_matrix）
{
	

	if（angle>=-90.0&&angle<-60.0）
	{
		motion_matrix[0]=-1-2*sin（angle*PI/180.0）+2*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[1]=-2*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[2]= 0.0;
		motion_matrix[3]= 2+2*sin（angle*PI/180.0）-4*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[4]= 4*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[5]= 0.0;
		motion_matrix[6]=-1.0;
		motion_matrix[7]= 0.0;
		motion_matrix[8]= 0.0;
	}
	else if（angle>=-60.0&&angle<-30.0）
	{
		motion_matrix[0]= 0.0;
		motion_matrix[1]=-2-4*sin（angle*PI/180.0）+2*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[2]= 1+2*sin（angle*PI/180.0）-2*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
        motion_matrix[3]= 0.0;
		motion_matrix[4]= 4+4*sin（angle*PI/180.0）-4*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[5]=-2-2*sin（angle*PI/180.0）+4*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[6]=-1.0;
		motion_matrix[7]= 0.0;
		motion_matrix[8]= 0.0;
	}
	else if（angle>=-30.0&&angle<0.0）
	{
		motion_matrix[0]= 0.0;
		motion_matrix[1]= 0.0;
		motion_matrix[2]= 0.0;
		motion_matrix[3]= 0.0;
		motion_matrix[4]=-4*sin（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[5]= 2*sin（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[6]=-1.0;
		motion_matrix[7]= 2+4*sin（angle*PI/180.0）-2*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[8]=-1-2*sin（angle*PI/180.0）+2*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
	}
	else if（angle>=0.0&&angle<30.0）
	{
		motion_matrix[0]=-1.0;
		motion_matrix[1]= 2-4*sin（angle*PI/180.0）-2*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[2]=-1+2*sin（angle*PI/180.0）+2*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[3]= 0.0;
		motion_matrix[4]= 4*sin（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[5]= -2*sin（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[6]= 0.0;
		motion_matrix[7]= 0.0;
		motion_matrix[8]= 0.0;
	}
	else if（angle>=30.0&&angle<60.0）
	{
		motion_matrix[0]=-1.0
		motion_matrix[1]= 0.0;
		motion_matrix[2]= 0.0;
		motion_matrix[3]= 0.0;
		motion_matrix[4]= 4-4*sin（angle*PI/180.0）-4*cos（angle*PI/180.0）+4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[5]=-2+2*sin（angle*PI/180.0）+4*cos（angle*PI/180.0）-4*sin（angle*PI/180.0）*cos（angle*PI/180.0）;
		motion_matrix[6]= 0

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        784  2008-07-24 18:02  psf\basic_op.h

     文件       7854  2008-07-29 21:58  psf\bmpio.c

     文件        855  2008-07-18 17:49  psf\bmpio.h

     文件       5529  2008-07-09 16:15  psf\bmpReadWrite.c

     文件       2688  2008-07-31 19:29  psf\deblur.c

     文件       3026  2008-07-30 14:05  psf\edgedetect.c

     文件       4154  2008-07-29 11:47  psf\fft.c

     文件       4211  2008-08-06 10:29  psf\fft2.c

     文件       4083  2008-07-24 18:05  psf\img.c

     文件       3091  2008-07-18 17:02  psf\img.h

     文件      18418  2008-07-25 09:17  psf\img_base_lib.c

     文件      24413  2008-08-06 10:37  psf\img_base_lib1.c

     文件        350  2008-07-15 15:22  psf\img_template_data.h

     文件       5354  2008-08-06 17:23  psf\liu.dsp

     文件        514  2008-07-16 16:47  psf\liu.dsw

     文件     181248  2010-04-07 19:25  psf\liu.ncb

     文件       6049  2010-04-07 19:20  psf\liu.plg

     文件      17853  2008-08-06 17:22  psf\psf.c

     文件        159  2008-07-23 10:12  psf\psf.h

     文件       2842  2008-08-05 15:08  psf\psf2otf.c

     文件       9517  2008-08-07 10:54  psf\wiener .c

     文件       5056  2008-07-29 14:07  psf\win.c

     文件       7025  2008-08-06 17:28  psf\win1.c

     文件       6632  2008-08-06 17:41  psf\win2.c

     文件          0  2008-07-30 18:33  psf\~VC17C.tmp

     文件      95232  2008-08-07 15:20  psf\~VC196.tmp

     文件      95232  2008-08-07 11:54  psf\~VC197.tmp

     文件        512  2008-08-06 14:12  psf\~VC222.tmp

     文件          0  2008-08-06 14:12  psf\~VC223.tmp

     文件          0  2008-09-10 19:14  psf\~VC438.tmp

............此处省略30个文件信息