Mitre_sfr代码注解（二） sfrproc函数

   笔者学识有限。 本博客旨在对mitre_sfr核心代码进行简单注解， 详细内容读者仁者见仁智者见智。

   bsfrproc函数是ISO12233附录A代码中的同名函数， 也是mitre_sfr中最核心的函数。因此对sfrproc函数的理解至关重要， 如果读者想对mitre_sfr功能进行自定义的移植， 最重要的就是对这个函数的移植与理解。

    sfrproc函数涵盖了ISO12233 Sect 6.3.2 中绝大部分的内容， 整个SFR计算的核心内容都在这个函数里面， 包含了下列步骤的主干(具体函数的注解见后续博文， 本文以核心函数注解为主， 其他main函数涉及的其他非主要部分函数， 且有必要注解的函数见后续再另行注解)：

    1. 寻找sfr刀口的质心（locate_centroids函数）

    2. 拟合刀口划过的直线（fit函数）

    3. 刀口的超采样（ESF）（bin_to_regular_xgrid函数）

    4. ESF转换为LSF（calculate_derivative函数）

    5. 对LSF加汉明窗（apply_hamming_window函数）

    6. 计算SFR（ftwos函数）

short sfrProc (double **freq, double **sfr, 
	       int *len,
	       double *farea,
	       unsigned short size_x, int *nrows,
	       double *slope, int *numcycles, int *pcnt2, double *off, double *R2,
	       int version, int iterate, int user_angle)
{
  unsigned short i, j, col, err = 0;
  long pcnt;
  double dt, dt1, sfrc, tmp, tmp2;
  double *temp=NULL, *shifts=NULL, *edgex=NULL, *Signal=NULL;
  double *AveEdge=NULL, *AveTmp=NULL;
  long *counts=NULL;
  int nzero;
  unsigned short size_y;
  unsigned int bin_len;
  double avar, bvar, offset1, offset2, offset;
  double centroid;
  int start_row, center_row;
  double *farea_old;
  double cyclelimit;
  FILE *f