Matlab读取复数（SLC）数据

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图像处理源代码”，见下面分享

InSAR_Matlab工具箱

链接：https://pan.baidu.com/s/17NIo-49nneZS2faqhiCaPQ

提取码：ddjx