2026年5月5日
GDAL—栅格数据按瓦片读取、虚拟文件创建、VRT创建、进程和线程处理
摘要： 1.按瓦片读取 from osgeo import gdal import numpy as np ds = gdal.Open(r"E:\DroughtAnalysis\data\SPEI_raster\gdal_spei12\large_raster.tif") band = ds.GetRas 阅读全文
posted @ 2026-05-05 17:19 MyEngine 阅读(32) 评论(0) 推荐(0)
GDAL—瓦片格式栅格数据创建和修改
摘要： 瓦片栅格数据 要点：瓦片是物理分区存储方式；实际的写入操作是逻辑分区方式；通过新建瓦片栅格的方式修改瓦片设置 import numpy as np from osgeo import gdal, osr # 参数设置 cols = 25000 # 列数 rows = 20000 # 行数 bands 阅读全文
posted @ 2026-05-05 16:13 MyEngine 阅读(21) 评论(0) 推荐(0)
GDAL-栅格数据和矢量数据的创建和修改
摘要： 1.栅格数据的创建 顺序：选择驱动gdal.GetDriverByName("GTiff")→创建空数据ds = driver.Create(output_raster, cols, rows, bands, dtype)→写入数据ds.GetRasterBand(1).WriteArray(dat 阅读全文
posted @ 2026-05-05 15:13 MyEngine 阅读(19) 评论(0) 推荐(0)
GDAL-栅格数据重采样，波段运算
摘要： 1. 栅格数据重采样 from osgeo import gdal gdal.UseExceptions() input_raster = r"" output_raster = r"" # 目标分辨率（度） target_xres = 0.1 target_yres = 0.1 gdal.Warp 阅读全文
posted @ 2026-05-05 14:55 MyEngine 阅读(19) 评论(0) 推荐(0)
GDAL——栅格数据的裁剪，矢量数据的裁剪和按属性导出
摘要： 1. 栅格数据的裁剪 需要注意的小细节： 核心函数是gdal.Warp() Warp()的属性中output_raster和input_raster可以使路径也可以是Dataset,而cutlineDSName只能是路径。 srcBand属性：表示只裁剪第几个波段，不设置则裁剪全部波段。 cutli 阅读全文
posted @ 2026-05-05 14:21 MyEngine 阅读(49) 评论(0) 推荐(0)