ccsm_refer.pro

来自「idl环境下运行」· PRO 代码 · 共 19 行

pro ccsm_refer,check=check,refer=refer,rms_refer=rms_refer,cns=cns,cnl=cnl,cnb=cnb,p=p
 ;    ;计算交叉相关系数
    for n=0l,cns*cnl-1 do begin
;     ;计算第一项积的和m_sum,平方和r_square,c_square
     m_sum=0l
     r_sum=0l
     c_sum=0l
        FOR i =n*cnb,(n+1)*cnb-1 DO begin
        m_sum=float(m_sum+refer[i]*check[i])
        r_sum=float(r_sum+refer[i]^2)
        c_sum=float(c_sum+check[i]^2)
        endfor
     a=float(cnb*m_sum)-float(total(check[n*cnb:(n+1)*cnb-1])*total(refer[n*cnb:(n+1)*cnb-1]))
     b=float(cnb*r_sum)-float(total(refer[n*cnb:(n+1)*cnb-1])^2)     
     c=float(cnb*c_sum)-float(total(check[n*cnb:(n+1)*cnb-1])^2)
     rms_refer[p*cns*cnl+n]=float(a)/float(sqrt(b*c+0.00001))
    
     endfor
end