矩阵运算,局部标准差!
IDL这个语言中常规计算局部标准差的方法都是采用双循环的方式,本人提出一种矩阵计算得到不同尺寸的“滑块”的标准差(3×3、5×5、7×7等都可以实现),话不多说直接上代码,代码中有详细的说明注释。
下面的函数是用来计算局部标准差:
FUNCTION uniformity_convol_method, raw_data, filter_h, filter_w, invalid_value
;raw_data 是原始的二维数组
;fileter_h 是滑块的高度
;filter_w 是滑块的宽度
; invalid_value 是数据中的无效值、空值等,根据需要其是否参与计算
; 函数返回一个与初始矩阵大小相同的,标准差矩阵
w = N_ELEMENTS(raw_data[