图像空间频率和清晰度

/********************************************************************************
*函数描述： SpaceFreq 计算并返回一幅图像的空间频率     
*函数参数： IplImage *img 单通道8位图像
*函数返回值：double   
*********************************************************************************/
double SpaceFreq(IplImage *img)
{
double RF = 0;
double CF = 0;
double SF = 0;

int i,j;//循环变量
int height=img->height;
int width=img->width;
int step=img->widthStep/sizeof(uchar);
uchar *data=(uchar*)img->imageData;
double num = width*height;

//行频计算
for(i=0;i<height;i++)
{
for(j=0;j<width;j++)
{
RF += (data[i*step+j+1]-data[i*step+j])*(data[i*step+j+1]-data[i*step+j]);
}
}
    RF = sqrt(1.0*RF/num);

//列频计算
for(i=0;i<height;i++)