goldstein 滤波

%-------------------------------
% 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