【python】sarscape输出的哨兵一号雷达影像干涉图展示

1.参考链接

envi sarscape 帮助文档
np.arctan() math.tan() math.atan()
【python】读取envi输出的复数二进制影像文件

2.实现过程

1. 读取复数二进制影像文件参考 参考链接3
2. 根据sarscape的文档介绍，干涉图相位的计算公式如下图
   
   3.使用np.arctan()来计算干涉图复数影像对应的相位，得到相位的弧度值

import gdal
import numpy as np
def img2array(path):
    dataset = gdal.Open(path)
    width = dataset.RasterXSize #栅格矩阵的列数
    height = dataset.RasterYSize #栅格矩阵的行数
    data = dataset.ReadAsArray(0,0,width,height)#单波段
    return data

path = r'H:\1.cd\8.winter_sbas_SBAS_processing\work\work_interferogram_stacking\IS_20181106_m_0_20190105_s_2_reflat_fint'
dataset = img2array(path)
ds = dataset[200:202,200:202]
print('complex array is:\n',ds,'\n')
print('real value of complex array is:\n',ds.real,'\n')
print('imag value of complex array is:\n',ds.imag,'\n')
print('interferometric phase array is:\n',np.arctan(ds.imag/ds.real),'\n')#计算得到是弧度值