Landsat栅格文件I/O
项目描述
这是一个开源的Python包,用于简单地将Landsat影像加载为NumPy数组。当从USGS Earth Explorer下载Landsat影像时,数据集包含许多波段(.tif文件)和一些元数据文件(.txt和.xml文件)。espatools被构建来解析.xml元数据文件,读取该数据集的所有波段,并在Python环境中提供方便直观的方式访问这些元数据,同时访问原始数据。espatools可在GitHub和PyPI上找到。
连接
该包大量使用properties以一致和声明的方式创建强类型对象。
该包实现了一种将这些数据集转换为PyVista数据集(vtkImageData)的方法。
PVGeo已为espatools实现了一个接口,通过XML元数据文件读取Landsat影像。有关更多详细信息,请查看PVGeo的Landsat Reader。
入门
espatools 可在 PyPI 上获取。
$ pip install espatools
使用方法
我们认为 espatools 非常易于使用;试用一下并告诉我们您的看法,因为这只是一个alpha版本!
首先,检查 这个 Jupyter Notebook,了解在Python环境中读取Landsat影像后的简单绘图演示。
然后看看 RasterSet 对象上的 .to_pyvista() 方法,在 PyVista/VTK 中查看影像的3D数据集。
然后查看 Landsat Reader,在 PVGeo 的文档中,其中 espatools 提供了在 ParaView 中直接使用的接口。
示例假彩色
import espatools
import matplotlib.pyplot as plt
# Create the reader to manage I/O
reader = espatools.RasterSetReader(filename='metadata.xml')
# Perform the read and yield a raster set
raster = reader.read()
# Get an RGB color scheme
color = raster.get_rgb('false_a')
# Now plot the false color image
plt.imshow(color)
上述代码的结果产生以下假彩色图像
您还可以使用 PyVista 在 3D 中查看数据集
mesh = raster.to_pyvista()
mesh.plot(scalars='false_a', rgb=True, cpos='xy')
项目详情
关闭
espatools-0.1.1.tar.gz 的哈希值
算法 | 哈希摘要 | |
---|---|---|
SHA256 | 22ce6957fe49d084d3bbc475b0c0d017adf3e775c8c062f116515da9ae57775d |
|
MD5 | 8f6740461aab6c2717ab3b3ff9298019 |
|
BLAKE2b-256 | fe258bac655deaaf33bfb895482a85c108f44e7cc320f479bca79de5c487f4f2 |