Python GDAL工具使用及使用VRT格式数据处理

最新推荐文章于 2025-01-15 12:56:43 发布
最新推荐文章于 2025-01-15 12:56:43 发布
阅读量3.3k 收藏 8
点赞数
分类专栏： python 遥感 地理信息 文章标签： python 大数据
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/tk20190411/article/details/120525665
版权
本文介绍了如何利用GDAL的VRT(Virtual Raster Tile)功能来管理和处理地理空间数据。通过构建虚拟栅格文件，结合gdalbuildvrt和gdalwarp命令，实现了数据列表的拼接和重投影，显著减少了内存占用和处理时间。这种方法适用于多幅DEM影像的整合，并且可以直接在QGIS中打开或用GDAL读取，提高了工作效率。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

前言

VRT即虚拟栅格文件。在这里插入图片描述

使用方法

先找到gadl工具exe文件的路径

D:\Anaconda\envs\geopandas\Library\bin

对数据列表构建vrt：
首先在工具路径下打开cmd，或者打开cmd输入完整工具路径
构建虚拟栅格

gdalbuildvrt -input_file_list E:\GeoLearn\dem_vrt_test\data_list.txt E:\GeoLearn\dem_vrt_test\data_vrt.vrt

在这里插入图片描述

对虚拟栅格进行重投影, 并指定输出分辨率和插值方法，可以覆盖已存在文件

gdalwarp -t_srs EPSG:32651 -r bilinear -tr 90 90 -overwrite E:\GeoLearn\dem_vrt_test\data_vrt.vrt E:\GeoLearn\dem_vrt_test\data_vrt_utm.vrt

在这里插入图片描述
原数据是两幅DEM影像，下图为拼接并重投影后的vrt文件，只有5 kb
在这里插入图片描述
可以直接在QGIS中打开，或者使用GDAL读取，非常节省数据处理时间和内存空间
在这里插入图片描述
更多使用方法可以查看GDAL的官方文档https://gdal.org/programs/index.html#raster-programs

gdal应用-gdalbuildvrt
弗里曼的小伙伴
03-31 5927
gdal应用-gdalbuildvrt sf2gis@163.com 2015年1月16日 1 gdalbuildvrt: 创建VRT（virtualgdal dataset）。 gdalbuildvrt [-tileindex field_name]              [-resolution {highest|lowest|average|user}]         
gdal.Buildvrt和gdal.Warp工具实现TIF影像拼接、矢量裁剪
weixin_40501429的博客
03-15 6323
看到GDAL库有很多工具可以直接使用，避免写长串的代码。这里算是记录一下自己的学习过程吧~ 比如gdalwarp，李老师的博客中写到的：（https://liminlu.blog.csdn.net/article/details/8531979） 以上是用命令行的方式实现，自己敲了几行代码在pycharm中实现（毕竟不太习惯使用命令行的方式）： from osgeo import gdal,gdalnumeric from osgeo.gdalconst import * # 查看DEM图的属
【待细看】关于GDAL的VRT格式
Theropod的博客
03-13 4849
最近才知道Gdal有这个功能，简单来说就是用一个xml描述对一个或一组栅格的处理（投影、mosiac、各种计算），生成一个.vrt文件，这样并没有真的进行计算，也没有存储文件，只有当要使用vrt的时候这些操作才会进行。比较省心。 Lazy raster processing with GDAL VRTs 官方文档 也可以导入到postgis里面...
利用GDAL vrt实现单波段浮点数dem包装为三波段rgb影像
CPWWHSU的博客
04-30 1016
本文介绍了利用GDAL vrt实现单波段浮点数dem包装为三波段rgb影像，具体为用vrt表达dem、将dem拉伸为0-255 区间的灰度图、灰度着色三个过程，所有中间结果可用vrt格式来表达。借助GDALTranslate函数，设置Options参数将高程值拉伸到0-255区间。借助GDALBuildVRT构建tif的vrt描述文件。
gdal影像镶嵌并裁剪
陨星落云的博客
04-05 1609
# -*- coding: utf-8 -*- """ Created on Tue Apr 5 20:18:49 2022 @author: 陨星落云 """ import os from osgeo import gdal def mergeAndClip(img_dir,clip_shp,out_dir): files = [os.path.join(img_dir,s)for s in sorted(os.listdir(img_dir))] merged_vrt = os.
Python地学分析 — GDAL对遥感影像重投影
一叶知秋的博客
12-20 9772
欢迎关注博主的微信公众号：“智能遥感”。 该公众号将为您奉上Python地学分析、爬虫、数据分析、Web开发、机器学习、深度学习等热门源代码。 本人的GitHub代码资料主页（持续更新中，多给Star，多Fork）： https://github.com/xbr2017 CSDN也在同步更新： https://blog.csdn.net/XBR_2014 “ 对遥感影像重投影是遥感...
Python GDAL VRT构建技术与STAC响应集成示例
在给定的文件信息中，涉及的知识点涵盖了遥感影像处理、地理信息系统（GIS）、Python编程语言以及遥感数据格式转换。具体知识点如下： 1. 静态（Static） 在描述中，“静态”一词可能用于形容静态数据或静态图像，...
Python工具gdal2mbt：快速转换GDAL到MBTiles格式
资源摘要信息:"gdal2mbt是一个命令行工具，它允许用户从GDAL兼容的文件格式创建MBTiles格式的瓦片地图包。...以上就是对gdal2mbt工具的详细介绍和知识点解析，希望能对使用该工具进行地理空间数据处理的用户提供帮助。
【GDAL-Python】6-在Python中使用GDAL合并栅格数据
weixin_44702962的博客
04-18 1636
（1）教程内容：介绍如何使用 Python 中的 GDAL 将同一目录中的一堆栅格切片合并到单个 geotiff 文件中； （2）技术路线：使用GDAL合并栅格数据有两个选项：1） 在Python脚本中使用子进程调用gdal_merge.py （终端也可以）2） 创建一个虚拟栅格vrt，然后使用translate转换为geotiff PS:在合并影像文件较多且数据量大时推荐2），1）方法会加载所有文件，对硬件的内存占用较高。
python && gdal开发文档
05-29
- **强制GDAL使用指定的驱动程序**: 指定使用特定的驱动程序来处理文件。 **2. 数据操作** - **获取shapefile特征计数**: 计算shapefile文件中特征的数量。 - **从PostgreSQL数据库获取所有PostGIS图层**: ...
GADL源代码（值得一读）
07-19
gdal源码，便于了解gdal底层实现
GDAL对图像文件格式的转换
m0_48011081的博客
03-21 5008
我们常常在图像处理过程中遇到不同软件或程序要求输入的图像格式不同（有些程序或软件支持的数据格式不是常用的Tiff，Img等数据格式），因此需要对不同的数据格式相互进行转换。 我这里以GTiff(.tif)数据转换为PCRaster(.map)数据为例。 首先需要安装GDAL，我这里是在Anaconda上直接安装了基于Python的GDAL，可以在下面网站自行下载，https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 例如下面对应的就是Pyth...
使用VRT裁剪栅格影像
qq_43290189的博客
03-28 2248
代码同样来自《Python地理数据处理》，这里做一个学习记录。书籍原外文名“Geoprocessing with Python”，在曼宁出版社可以搜索到，可以线上免费预览全文，下载源码和配套数据。附一个链接：Geoprocessing with Python VRT（GDAL虚拟栅格格式），允许使用XML文件来定义数据集，本身不存储数据，只描述了如何从其他文件中获得数据。我自己的理解觉得有点像数据库里视图与数据表中视图的概念。 1.实现代码和输出结果 # -*- coding: utf-8 -*-
【GDAL】如何利用GDAL批量进行等同于Arcmap中Mosaic To New Raster的操作
qq_41925877的博客
03-27 439
随后利用 lxml文件读写VRT文件，并加入自定义的Pixel function以解决重叠栅格取值的问题，下面代码以获取最大值为例。
GDAL栅格数据操作拼接
GIS探索之路
04-10 936
使用python GDAL 对下载好的瓦片（一定处理后）进行拼接。
使用GADL对高程数据进行填洼
weixin_58851039的博客
07-30 1098
以上代码演示了如何使用`gdal.FillSinks()`函数对DEM数据进行填洼处理，并保存为新的文件。请将`input_dem_file`和`output_filled_dem_file`分别替换为实际的输入和输出DEM文件的路径。4. 保存结果：填洼完成后，可以使用GADL将处理后的DEM数据保存为新的文件，可使用函数如GDALCreate()和GDALWriteRaster()。2. 检测洼地：使用GADL中的算法，例如流域填洼算法（Fill sinks），来检测DEM中的洼地。
GADL 读取栅格数据(tif)，遍历数组
qq_35662333的博客
03-22 477
简单记录下用c++读取栅格数据，遍历数组。
GDAL/OGR：使用虚拟图层 .vrt
最新发布
weixin_43416590的博客
01-15 124
与GDAL类似，OGR也支持虚拟文件格式(VRT)。OGR虚拟驱动程序功能更强大，可用于实时数据转换。使用基于XML的格式文件，可以将多个数据层合并成一个虚拟层.vrt。VRT文件还用于配置将表格数据读取为空间数据格式。
GDAL栅格数据处理教程
weixin_44785184的博客
02-20 4843
本文结合实例详细讲解了如何使用GDAL操作栅格数据，包括栅格数据的读写，以及一些常用的栅格数据处理操作，如坐标变换、裁剪、镶嵌、插值等。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值