gdal.Buildvrt和gdal.Warp工具实现TIF影像拼接、矢量裁剪

最新推荐文章于 2024-09-18 10:17:37 发布

见贤思齐547

最新推荐文章于 2024-09-18 10:17:37 发布

阅读量6k

点赞数 4

分类专栏： Python学习过程 Python地理数据处理文章标签： python

本文链接：https://blog.csdn.net/weixin_40501429/article/details/114838275

版权

本文记录了使用GDAL库的gdalwarp和gdal.Buildvrt+gdal.Warp工具进行TIF影像拼接和裁剪的过程，通过Python实现。对比了不同操作结果，并探讨了数据类型对图像显示的影响。

摘要由CSDN通过智能技术生成

看到GDAL库有很多工具可以直接使用，避免写长串的代码。这里算是记录一下自己的学习过程吧~

1.比如gdalwarp，李老师的博客中写到的：（https://liminlu.blog.csdn.net/article/details/8531979）

以上是用命令行的方式实现，自己敲了几行代码在pycharm中实现（毕竟不太习惯使用命令行的方式）：

from osgeo import gdal,gdalnumeric
from osgeo.gdalconst import *

# 查看DEM图的属性信息(无投影信息)
# imageInformation = gdal.Info('D:\ProfessionalProfile\DEMdata\OutputDataPycharm2\ASTGTM2_N34E114_dem.tif')
# 以下数据有投影信息UTM-47N
# imageInformation2 = gdal.Info('D:\ProfessionalProfile\MODISrelevant\MOD13Q1mrtProcess\MOD13Q1.A2020161.250m_16_days_NDVI.tif')

# 打开TXT文件路径，此文件存储了整个山东省的DEM图
f = open(r'D:\ProfessionalProfile\DEMdata\OutputDataPython\StatisticsDEMfile.txt', 'r')