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

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

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

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


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