IDL编程技巧

1 Absolute Distance

; win is radius of moving window
D_D_all=((win-indgen(win*2+1)#(intarr(1,win*2+1)+1))^2+(win-(intarr(win*2+1)+1)#indgen(1,win*2+1))^2)^0.5

2 Mask

;quality control using QA
index=where(QA gt 3.0, count)
if count gt 0 then tmpLST[index] = !VALUES.F_NAN

3 时空插补像元

• 索引转坐标
• 确定窗口
• 相似像元查找
  前后两个影像，若第一个找0、1，若最后找n-1、n-2
• 空间权重
• 时间权重
• 依权重插值