你的习惯决定了你会成为什么样的人。

---

前言

1. 版本

   1.1 山东青岛，2022年7月17日，Version 1

2. 摘要

GDAL基础知识

3. 说明

  此系列文章是GDAL的学习笔记，主要是对 《Python与开源GIS》 和 《Python地理数据处理》 的摘录和学习，版权归上述书籍作者所有，如有不妥，请告知，必第一时间删除，谢谢。

4. 微信公众号GISRSGeography

• 公众号 GISRSGeography的内容涉及GIS，遥感和作物模型等的内容，会坚持更新，
  欢迎大家关注，谢谢！。
  

一、GDAL简介

  GDAL (Geospatial Data Abstraction Library)是一个在X/MIT协议下的开源栅格数据库，可以进行空间数据的读写^[1]，严格定义可以参考GDAL官网。

二、导入GDAL

1、在Python中导入GDAL通常可以基于以下三种方法

# 1. 导入方法1
import gdal# 2. 导入方法2
from osgeo import gdal# 3. 导入方法3
try:import gdal
except:from osgeo import gdal

2、除gdal包之外，gdalconst包亦建议导入

# 导入gdalconst
from osgeo.gdalconst import *

2、gdal中的栅格数据驱动

• 要先读取某种类型的数据，需要先注册数据驱动，即初始化一个对象，让其知道某种数据结构。可以使用以下语句注册所有数据驱动，但是只能读不能写：

gdal.AllRegister()

• 单独注册某一类型的数据驱动，即可以读，也可以写，甚至可以创建数据集，下面语句注册了TIFF数据类型

# 注册'GTiff'数据驱动
driver = gdal.GetDriverByName('GTiff')

3、查看GDAL支持的数据驱动

# 查看GDAL支持的数据驱动
drv_count = gdal.GetDriverCount()
print(drv_count)

三、读取遥感影像的信息

1、打开GeoTIFF文件

• 对于GIS栅格数据而言，一个数据集通常由若干个文件组成，并且使用额外的信息组织他们之间的关系。
• 对于GDAL而言，栅格数据集是由栅格的波段数据以及所有波段的共同属性构成的。
• 打开GeoTIFF栅格数据集的操作如下：

from osgeo import gdal
# 打开GeoTIFF数据集
rds = gdal.Open('/gdata/geotiff_file.tif')
# 通过dir函数查看当前对象可用操作
dir(rds)