%-------------------------------
%?InSAR?干涉條紋圖Goldstein濾波??
%?Author:?龍學(xué)軍xuejun.long@gmail.com
%?Date??:?2008.9
%-------------------------------
%?本程序用于實(shí)現(xiàn)?Goldstein?and?Werner?（1998）提出的局域頻率加權(quán)濾波方法
%
%?使用低通濾波器w對(duì)局域窗口進(jìn)行濾波
%?加權(quán)系數(shù)a決定濾波的深度
clear?all;
close?all;
clc;
%-------------------------------------------------------------------------
%?讀取InSAR復(fù)數(shù)據(jù)
%master=read_InSAR_date（‘D:\工作文件夾\matlab\insar復(fù)數(shù)據(jù)\kashi-master_multilook_256_1024.dat‘10242561）;
%slave=read_InSAR_date（‘D:\工作文件夾\matlab\insar復(fù)數(shù)據(jù)\kashi-Flattened-slave.dat‘10242561）;
%image=master.*conj（slave）;%共軛相乘
%--------------------------------------------------------------------------
%?直接讀取干涉圖
[imagemap]=imread（‘1.bmp‘）;
%image=double（image（500:1024500:1024））;%取部分圖像
image=double（image）;
%轉(zhuǎn)換為正余弦圖
image=image*2*pi/255;
image=cos（image）+i*sin（image）;

[widthheight]=size（image）;%獲取圖像寬度和高度

%--------------------------------------------------------------------------
%?設(shè)定參數(shù)并初始化變量
%?FFT局域窗口設(shè)定
m=15;
n=15;

armw=（m-1）/2;
armh=（n-1）/2;


%?設(shè)定低通濾波窗口
w=zeros（mn）;
w（armw:armw+2armh:armh+2）=1;
%?初始化變量
result=ones（width-m+1height-n+1）;
i1=ones（mn）;
i2