1.参考链接
envi sarscape 帮助文档
np.arctan() math.tan() math.atan()
【python】读取envi输出的复数二进制影像文件
2.实现过程
- 读取复数二进制影像文件参考 参考链接3
- 根据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')#计算得到是弧度值