GDAL学习笔记03：获取遥感影像基本信息

全力以赴做好自己该做的事情。

---

前言

1. 版本

   1.1 2022年7月24日，Version 1

2. 摘要

GDAL，栅格数据基本信息

3. 说明

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

4. 微信公众号GISRSGeography

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

---

一、获取栅格数据基本信息

• 基于GDAL获取栅格数据基本信息的代码如下：

# 1. 相关模块的导入
import gdal

# 2. 路径处理和基本变量定义
rootdir = r'E:\03_学习盘\04_数据分析实例\01_Python\06_GDAL\gdata\\'
#print(rootdir)
tiffile = rootdir + 'geotiff_file.tif' # 栅格文件

# 3. 读取栅格文件
tifsource = gdal.Open(tiffile)
#dir(tifsource) # 通过dir函数查看tifsource可以进行的相关操作
  # 3.1 查看影像元数据
print('影像元数据信息如下：')
print(tifsource.GetMetadata())
  # 3.2 获取栅格数据的信息
print('影像数据存储信息如下：')
print(tifsource.GetDescription())
  # 3.3 获取栅格数据的波段
print('栅格数据波段信息如下：')
print(tifsource.RasterCount)
  # 3.4 获取影像的行列数目
print('影像的行列数目如下：')
img_width,img_height = tifsource.RasterXSize,tifsource.RasterYSize
print(img_width,img_height)
  # 3.5 获取影像的空间参考
print('影像的空间参考如下：')
print(tifsource.GetGeoTransform())
  # 3.6 获取影像的投影信息
print('影像的投影信息如下：')    
print(tifsource.GetProjection())

• 运行结果如下：

影像元数据信息如下：
{'AREA_OR_POINT': 'Area', 'PyramidResamplingType': 'NEAREST'}
影像数据存储信息如下：
E:\03_学习盘\04_数据分析实例\01_Python\06_GDAL\gdata\\geotiff_file.tif
栅格数据波段信息如下：
3
影像的行列数目如下：
1500 900
影像的空间参考如下：
(1868454.913, 30.0, 0.0, 5353126.266, 0.0, -30.0)
影像的投影信息如下：
PROJCS["Albers_Beijing54",GEOGCS["Unknown datum based upon the Krassowsky 1940 ellipsoid",DATUM["Not_specified_based_on_Krassowsky_1940_ellipsoid",SPHEROID["Krassowsky 1940",6378245,298.3,AUTHORITY["EPSG","7024"]],AUTHORITY["EPSG","6024"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["latitude_of_center",0],PARAMETER["longitude_of_center",105],PARAMETER["standard_parallel_1",25],PARAMETER["standard_parallel_2",47],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]

二、参考资料

1. 《Python与开源GIS》