【python】读取envi输出的复数二进制影像文件-程序员宅基地

技术标签： python envi insar InSAR 图像处理算法 sarscape 二进制影像文件

1.参考链接

2.复数二进制影像文件

此处读取envi sarscape时序干涉处理过程中生成的复数二进制文件。

读取二进制影像文件，使用gdal.Open(path)即可读取，path不需要加入后缀名
读取得到的数组的 dtype为complex64
提取复数的实部和虚部进行后续计算

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

path = r'..\0.data\0.intensity\burst_IW3_8_slc'
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')

在这里插入图片描述

3. 需要有对应的头文件

有头文件
没有头文件

4. 两个实数矩阵合并得到一个复数矩阵

【python】将两个实数矩阵合并为一个复数矩阵

# 输入：复数三维向量
import numpy as np
real = np.random.random(size=(2,3))
imag = np.random.random(size=(2,3))
inp = real+imag*1j
print(inp)

输出
在这里插入图片描述

5.求复数的幅度和幅角

np.abs()求复数的幅度
原博客
使用np.abs求复数的模

import numpy as np
a = 1-2j
print(abs(a)) #2.23606797749979
print(np.sqrt(5)) #2.23606797749979

np.angle()求复数的角度
原博客

import numpy as np
print('结果用弧度制表示：{}'.format(np.angle([1.0, 1.0j, 1+1j])))
print('结果用角度制表示：{}'.format(np.angle(1+1j, deg=True)))
# 结果用弧度制表示：[0. 1.57079633 0.78539816]
# 结果用角度制表示：45.0

6. 几何形式(实部和虚部)和指数形式(幅度和角度)的相互转换

b = 1 + 1j
amp = np.abs(b)
ang = np.angle(b)
print("幅度为：\n",amp,"\n角度为:\n",ang)
real = amp*np.cos(ang)
imag = amp*np.sin(ang)
print("计算得到的实部为：\n",real,"\n计算得到的虚部为:\n",imag)

输出为：
在这里插入图片描述

本文链接：https://blog.csdn.net/qq_45110581/article/details/109891201

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

数仓实战作业_select campseg_id, s_start_date,count(1) from i_h_-程序员宅基地

文章浏览阅读464次。作业题：在会员分析中计算最近七天连续三天活跃会员数。项目的数据采集过程中，有哪些地方能够优化，如何实现？连续值求解问题对表 dws_member_start_day 中的数据使用 row_number()函数进行排序，并且根据设备号 device_id分组，登录日期dt排序，用日期减去排名得到groupId 并限定如期最近7天，根据设备号和groupId分组，求出统计值cnt 大于等于 3的记录。WITH tmp as( SELECT ..._select campseg_id, s_start_date,count(1) from i_h_iop_2_93001_day_dm group b

nagios用NsClient自定义windows监控-程序员宅基地

文章浏览阅读76次。Nagios用NsClient自定义windows监控NsClient++来监控windows主机有三种方式：check_nt,check_nrpe,nsca.相对而言，check_nrpe已经可以满足大部分的监控任务要求了。这三种方式在nsclient的PDF中说的很详细，但全是Ｅ文，根据我自己的部署过程，边做边写下自己理解：check_nt的使用方法：这是NsClient++默认的使用方法，但..._nagios自定义监控windows

docker- 构建 oracle2c-r2(12.2.0.1) 的镜像-程序员宅基地

文章浏览阅读423次。需求由于公司的数据库需要使用新的oracle版本（12c-r2 -->12.2.0.1），从之前的oracle11g迁移到12c。所以，便有了我们今天的内容。首先，我们就先来介绍一下如何构建oracle12c的镜像（docker image）。如果大家有使用的需求而又不是正式的项目，可以直接到docker hub 上面 pull 一个别人家的。在这里附上链接：https://hu..._/opt/database/install/unzip -qqqo ../stage/components/oracle.jdk/1.8.0.91.0/

寻（光阴的故事）-程序员宅基地

文章浏览阅读101次。寻寻一束芬芳的鲜花插于窗前在我疲劳的时候看它一眼寻一把眷恋的故土带在身边在我思乡的时候捧于胸前寻一段纯真的恋情藏在心间在我寂寞的时候把她思念 ...

ZStack实践汇 | ZStack平台的使用心得-程序员宅基地

文章浏览阅读2.1k次。作者：ZStack 社区王彬Iaas云服务的普及,让我们在使用服务器的时候享受了飞一般的感觉，新兴企业在构建自己的系统时，往往都会选择购买云厂商的云服务器(虚拟机)进行使用，使用这样的虚拟机企业不需要购置任何硬件，不需要考虑机房、网络、服务器维护相关的工作便可以获取到一个低成本、安全免维护、伸缩性强、可灵活迁移的云服务器。在这个云服务器上我们可以快速的构建企业的业务系统。随着企业的不断发展，...

Spec2006使用说明-程序员宅基地

文章浏览阅读1.7k次。Spec2006使用说明五 10 十月 2014Bypenglee工具介绍SPEC CPU 2006 benchmark是SPEC新一代的行业标准化的CPU测试基准套件。重点测试系统的处理器，内存子系统和编译器。这个基准测试套件包括的SPECint基准和SPECfp基准。其中SPECint2006基准包含12个不同的基准测试和SPECfp2006年基准包含19个不..._spec2006如何手动停止