Matlab读取复数（SLC）数据

就让往事随风去...

已于 2024-12-17 14:01:17 修改

阅读量2.7k

点赞数 2

分类专栏： Matlab 相关文章标签： matlab

于 2022-12-26 12:30:18 首次发布

本文链接：https://blog.csdn.net/weixin_41652574/article/details/128442314

版权

Matlab 相关专栏收录该内容

7 篇文章

订阅专栏

本文介绍了一种使用Matlab进行InSAR数据处理的方法，包括复杂短数据读取、多视函数应用、高程数据读取及拼接等步骤，并提供了具体代码示例。此外还分享了荷兰Delft大学开发的InSAR图像处理源代码及其工具箱。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

complex_short 读取

        [data,count]=freadbk(infile,lines,bkformat,ro,rn,co,cn);
        master=freadbk(‘...\\master.raw',5000,'cpxint16');
        master_m=abs(master);         # 幅度
        master_p=angle(master)；        # 相位
        Figure; imagesc(abs(master_m))

多视函数

m_master=cpxmultilook(master,1,5)；

hgt读取
[phase,ampli]=freadhgt(infile,numlines);
[phase,ampli]=freadhgt('...\\master_slave.uint',500);

SRTM DEM读取
        SRTM下载地址：
        http://rmw.recordist.com/index.html

       % 打开文件
        fid=fopen('*.hgt');
       % 读取文件
        A=fread(fid,size,precision,skip,machineformat)
        N29E057=fread(fid,[1201,1201],'int16','b')；

拼接
dem=[[N29E057',N29E058'];[HN28E057',N28E058']]; % 矩阵拼接

写数据
        fid=fopen('...\bam.dem','w')；
        fwrite(fid,dem','int16')；
        fclose(fid)；

验证是否正确
        fid=fopen('...\bam.dem')；
        Bam_ dem=fread(fid,[2402 2402],'int16')；
        fclose(fid)；
        figure,imagesc(Bam_dem),colorbar

注：上面部分m函数文件为 “荷兰Delft大学开发的insar图像处理源代码”，见下面分享