• 大小: 2KB
    文件类型: .m
    金币: 1
    下载: 0 次
    发布日期: 2021-05-17
  • 语言: Matlab
  • 标签: Goldstein  matlab  

资源简介

完全可以使用,滤波效果很好。用时记得改改我程序中的加载的图片的名字,改成自己的要处理的图片就行。程序会覆盖原文件,请注意备份。

资源截图

代码片段和文件信息

%-------------------------------
% InSAR 干涉条纹图Goldstein滤波  
% Author: 龙学军xuejun.long@gmail.com
% Date  : 2008.9
%-------------------------------
% 本程序用于实现 Goldstein and Werner (1998)提出的局域频率加权滤波方法
%
% 使用低通滤波器w对局域窗口进行滤波
% 加权系数a决定滤波的深度
clear all;
close all;
%-------------------------------------------------------------------------
% 读取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(‘test1.bmp‘);
%image=double(image(500:1024500:1024));%取部分图像
[widthheight]=size(image);%获取图像宽度和高度

image=double(image);
%转换为正余弦图
image=image/255*2*pi;
image=cos(image)+i*sin(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+1he

评论

共有 条评论