使用Rasterio旋转栅格数据
旋转栅格数据使用Rasterio
使用Python旋转卫星图像并保持地理定位准确性
栅格数据类似于普通图像数据。每个像素都与地球表面上的一个位置相连。这使得问题复杂化。如果我们要旋转数据,我们还必须考虑底层的坐标参考系统(CRS)。在不调整其地理定位的情况下扭曲栅格会导致空间分析不准确。
调整地理定位不是一件简单的事情。幸运的是,Rasterio可以帮助我们。它是一个流行的Python库,用于地理空间数据分析。我们将使用该库来:
- 旋转一个栅格
- 并将图像投影到正确的CRS中。
在此过程中,我们将讨论Python代码,并且您可以在GitHub上找到完整的项目。
本文假设读者具备栅格数据和使用其CRS的基本知识。如果您想回顾,请查阅下面的文章。它更详细地介绍了栅格数据的投影。
如何使用Python在Landsat卫星图像上绘制坐标
使用Landsat元数据和Rasterio将像素位置映射到地理坐标
towardsdatascience.com
下载Landsat场景
对于我们的栅格数据,我们将使用卫星图像。具体来说,是Landsat场景。您可以使用EarthExplorer门户网站下载其中之一。或者,如果您想使用Python,下面的文章将引导您完成整个过程:
使用Python下载Landsat卫星图像
使用landsatxplore Python包简化Landsat场景下载
towardsdatascience.com
最终,您应该拥有一个包含所有Landsat级别2科学产品可用文件的文件夹。我们将使用红色可见光波段进行操作。对于Landsat 8或9场景,这由B4波段给出。
打开栅格文件
我们使用以下代码来打开和显示这个波段。ID为这个特定场景的Landsat场景ID…