python读取InSAR复数数据

InSAR复数数据的读取，如：fint(滤波后相位)，dint(滤波前相位)等

注：注意行列号输入正确！

转载请标明出处，谢谢！

注：有任何问题请私信或评论！

import numpy as np
import math
import os
import matplotlib.pyplot as plt
import struct
import matplotlibdef xshow(filename, nx, nz):f = open(filename, 'rb')pic = np.zeros((nx, nz))for i in range(nx):for j in range(nz):data = f.read(4)elem = struct.unpack("f", data)[0]pic[i][j] = elemf.close()return picdef ReadFloatComplex(file_path, col: int, row: int):""":param file_path:  float_complex:param col: samples:param row: lines:return: phase，intensity data"""# 读取数据row = row*2Data_RealImag = xshow(filename=file_path, nx=col, nz=row)coldiv2 = math.floor(row / 2.0)# 实部、虚部数据在复数数据中的索引Real_Index = list(range(0, row - 1, 2))Imag_Index = list(range(1, row, 2))# 提取实部数据real = []for i in range(coldiv2):Col_Ctl_Real = Real_Index[i]Data_Real = Data_RealImag[:, Col_Ctl_Real]real.append(Data_Real)# list转numpyreal = np.array(real)# print('real is {}'.format(real))# 提取虚部数据imag = []for j in range(coldiv2):Col_Ctl_Imag = Imag_Index[j]Data_Imag = Data_RealImag[:, Col_Ctl_Imag]imag.append(Data_Imag)# list转numpyimag = np.array(imag)# print('imag is {}'.format(imag))# 合成复数数据Data = real + 1j * imag# 相位数据phase = np.angle(Data)phase = np.transpose(phase)# 零值转为nanphase[phase == 0] = np.nan# 强度数据intensity = abs(Data)intensity = np.transpose(intensity)# 零值转为nanintensity[intensity == 0] = np.nanreturn phase, intensityif __name__ == '__main__':# -------------------------------------------------参数设置---------------------------------------------------------file = r''  # 复数数据路径col =   # 行列号row =   # 行列号phase, intensity = ReadFloatComplex(file_path=file, col=col, row=row)norm = matplotlib.colors.Normalize(vmin=-1, vmax=1)plt.imshow(phase, cmap='jet')plt.colorbar()plt.show()plt.imshow(intensity, cmap='gray', norm=norm)plt.colorbar()plt.show()

python读取InSAR复数数据相关推荐

如何用python读取文本中指定行的内容
如何用python读取文本中指定行的内容 搜索资料 我来答分享新浪微博 QQ空间浏览 5284 次查看全文 http://www.taodudu.cc/news/show-64036.ht ...
Window10 64位，通过Python读取.mif, shp文件
1. 读取依赖库及安装依赖模块osgeo ,osgeo依赖于GDAL .mif 是Mapinfo支持的信息 .shp是 postgis存储地理空间信息的文件格式: 尝试了很多安装方法:包括安装Ana ...
python 读取excel文件效率时间格式_python读取Excel文件中的时间数据
在使用python读取Excel文件中的时间格式,碰到的时间格式转换问题: 读取这样的表格: 输出这样的数据结果: 然而这样的结果却不是我们想要的,我们需要的是这样的结果: 1.安装python官方库 ...
技术图文：C# VS. Python 读取CSV文件指南
背景 CSV 是一种以逗号进行特征分隔的文本文件类型,在数据库或电子表格中是一种非常常见的导入导出格式.本篇图文就以泰坦尼克号船员获救预测( Kaggle)中使用的数据集为例来说明 C#.Python ...
python读取一个文件夹/子文件夹下的所有文件名字
python读取一个文件夹/子文件夹下的所有文件名字示例代码: import osfile_path = './images/' all_file_name = os.listdir(file_pa ...
python读取配置文件获取所有键值对_python笔记14-读取yaml配置文件(pyyaml)
yaml简介 1.yaml [ˈjæməl]: Yet Another Markup Language :另一种标记语言.yaml 是专门用来写配置文件的语言,非常简洁和强大,之前用ini也能写配置文 ...
python读取大文件目录_Python读取系统文件夹内所有文件并统计数量的方法
大家先看一下Python os模块中的部分函数 python 路径相关的函数 os.listdir(dirname):列出dirname下的目录和文件 os.getcwd():获得当前工作目录 os. ...
Python读取多个excel文件（删除字段、数据格式转换、dataframe多表合并）并写入ElasticSearch实战（自动创建索引、写入ElasticSearch、探索性数据分析）
Python读取多个excel文件(删除字段.数据格式转换.dataframe多表合并)并写入ElasticSearch实战(自动创建索引.写入ElasticSearch.探索性数据分析) 目录
python删除重复值所在的行数_使用python读取txt文件的内容,并删除重复的行数方法...
注意,本文代码是使用在txt文档上,同时txt文档中的内容每一行代表的是图片的名字. #coding:utf-8 import shutil readDir = "原文件绝对路经" ...

python读取InSAR复数数据

python读取InSAR复数数据

python读取InSAR复数数据相关推荐

最新文章

热门文章