使用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…