栅格数据创建与保存

栅格数据创建与保存

作者：阿振

邮箱：tanzhenyugis@163.com

博客：https://blog.csdn.net/theonegis/article/details/80089375

修改时间：2018-05-24

声明：本文为博主原创文章，转载请注明原文出处

---

思路与方法

使用Python进行栅格数据处理，很多时候，我们会将GDAL的Dataset对象转化为NumPy的ndarray对象，这样我们可以使用很多通用的Python库对数据进行处理，然后再借助GDAL库将数据写回到文件。

不同于普通的二进制文件，空间栅格数据的写需要注意两点：

1. 数据的投影信息（确定了平面坐标系）
2. 数据的地理坐标信息（确定了图像在给定坐标系下的位置）

在GDAL中，我们首先需要创建Dataset对象，然后给Dataset对象填充数据以及元数据信息就OK了。

Driver或者说GDALDriver（Python版本的API中对象名称好像都去掉了前缀GDAL，而C/C++版本的API很多对象前面都是有GDAL前缀的，如GDALDataset对象在Python中对应的是Dataset对象）有两个方法：Create()和CreateCopy()

所以，相应地，我们也有两种思路去创建一个Dataset对象：

1. 如果我们有一个原型数据，比如我们对原始数据进行了处理，处理之后，空间信息，波段等都没有变化，则可以将原始数据作为原型数据，使用CreateCopy()方法创建一个和原始数据一样的Dataset对象，然后在创建好的对象中填充一个ndarray