GIS遥感数据工程【免费下载链接】gdalGDAL is an open source MIT licensed translator library for raster and vector geospatial data formats.项目地址https://gitcode.com/gh_mirrors/gd/gdal点击查看免费下载gdal raster update是 GDAL 自 3.12 版本引入的栅格更新子命令用于将输入栅格的像素值写入一个已经存在的目标栅格数据集并在输入/输出坐标系不一致时自动完成重投影目标数据集的 extent、尺寸、分辨率与坐标系在整个过程中保持不变。本文基于 doc/source/programs/gdal_raster_update.rst 官方文档并结合 apps/gdalalg_raster_update.cpp 源码实现与 autotest/utilities/test_gdalalg_raster_update.py 测试用例完整讲解该命令的全部选项、重采样方法、概览overview同步刷新机制、gdal raster pipeline中的组合用法以及常见错误与规避方式帮助读者在影像更新、镶嵌修补、局部替换等场景中直接上手。命令概览做什么、不做什么gdal raster update的官方用途一句话即可概括Update the destination raster with the content of the input one.用输入栅格的内容更新目标栅格这一描述同样作为GDALRasterUpdateAlgorithm的DESCRIPTION常量定义在 gdalalg_raster_update.h 中。其核心行为特征如下就地更新目标数据集必须已存在且以**更新模式GA_Update**打开。源码 gdalalg_raster_update.cpp 中通过CPLAssert(poDstDS-GetAccess() GA_Update)强制校验这一点。自动重投影当输入与输出数据集的坐标系CRS不同时像素值会在写入前被重投影到目标 CRS 下的正确位置。不改变目标栅格结构目标数据集的 extent、尺寸、分辨率、CRS 均不被修改——这是与gdalwarp生成新文件最大的区别update只负责把新值填进已有的像素格子。概览自动刷新命令默认会同步刷新目标数据集已存在的 overviews内部机制复用gdal raster overview refresh除非显式指定--no-update-overviews。版本与能力演进该命令自 GDAL 3.12 起提供从 GDAL 3.13 开始它可以作为gdal raster pipeline中的一个步骤使用此时其输入数据集来自上一环节的输出。从命令整体架构看gdal raster update归属于gdal raster这一栅格子命令家族gdal raster入口自 3.11 引入并列于convert、clip、reproject、mosaic等子命令详见 doc/source/programs/gdal_raster.rst。底层实现基于 GDALWarp 的轻量封装尽管文档以命令行角度描述理解其底层实现有助于预判行为。在 gdalalg_raster_update.cpp 的RunStep()中算法做了以下事情获取输入数据集poSrcDS与目标数据集poDstDS校验二者不同详见下文常见错误。若指定了--geometry通过GetClipGeometry()继承自GDALClipCommon见 gdalalg_clip_common.h解析裁剪几何并导出为 WKT 作为-cutline选项。将--resampling-r、--warp-option-wo、--transform-option-to、--error-threshold-et等参数统一组装成CPLStringList交给GDALWarpAppOptionsNew()构建 warp 选项。调用核心函数GDALWarp(nullptr, hDstDS, 1, hSrcDS, psOptions, nullptr)即以目标数据集为输出句柄执行 GDALWarp——这正是写入已有数据集的实现关键。warp 完成后若目标带 1 有 overview 且未禁用更新则调用GDALRasterOverviewAlgorithmRefresh对受影响区域做部分刷新。可见gdal raster update本质上是GDALWarp面向原地更新场景的高层封装重投影、重采样、裁剪等几何能力全部来自 warp 框架。程序专属选项详解以下选项为gdal raster update专属均可同时用于命令行与 pipeline 步骤中。--et, --error-threshold ERROR-THRESHOLD变换近似的误差阈值以输入像素数为单位。默认值为0.125像素即允许近似变换引入不超过 1/8 像素的误差。例外情况当指定了RPC_DEM变换选项时默认切换为精确变换即--error-threshold0。传0表示使用精确变换不做近似逼近精度最高但计算更慢。从源码看该参数经由GDALRasterReprojectUtils::AddWarpOptTransformOptErrorThresholdArg()注册在 gdalalg_raster_reproject.cpp 中限定取值必须 ≥ 0。--geometry WKT_or_GeoJSON以 WKT 或 GeoJSON 字符串形式指定一个多边形或多多边形将更新操作限制在该几何范围内几何之外的像素保持不变。CRS 判定规则官方文档明确规定若输入几何为GeoJSON默认假定其坐标为 WGS84经纬度除非 GeoJSON 几何体内自带 CRS 定义或显式给出--geometry-crs。若输入几何为WKT默认假定其坐标系为输入数据集的 CRS除非显式给出--geometry-crs。坐标轴遵循GIS 友好约定X 为经度或东向eastingY 为纬度或北向northing。测试用例 test_gdalalg_raster_update.py 演示了geometryPOLYGON ((1 -1,2 -1,2 -2,1 -2,1 -1))的用法在 3×3 全零目标上仅在多边形覆盖的中央像素写入 1。若传入无法解析的字符串则报错Clipping geometry is neither a valid WKT or GeoJSON geometry见测试test_gdalalg_raster_update_geometry_invalid。--geometry-crs CRS指定--geometry坐标值所在的 CRS。若不指定按上一节规则推断。几何会从 geometry-crs 重投影到被更新数据集的 CRS后再执行裁剪。--no-update-overviews禁止更新目标数据集已存在的概览。默认行为是自动刷新 overviews这在目标数据文件较大、且更新范围很小时会额外耗时若希望只写全分辨率像素、概览稍后统一重建可加此开关。--to, --transform-option NAMEVALUE设置变换器选项适用于GDALCreateGenImgProjTransformer2()接受的参数RPC 相关的变换选项参考GDALCreateRPCTransformerV2()的文档。可多次指定。源码中该参数同样由AddWarpOptTransformOptErrorThresholdArg()注册gdalalg_raster_reproject.cpp并按NAMEVALUE键值对解析。测试用例给出了实际用法transform_option{SRC_METHOD: GEOTRANSFORM}。-r, --resampling RESAMPLING重采样方法默认nearest。可选值及其说明来自共享文档 doc/source/programs/gdal_options/warp_resampling.rst在源码 gdalalg_raster_reproject.cpp 中注册了完整的候选集合方法说明nearest最近邻默认最快插值质量最差bilinear双线性插值cubic三次卷积插值cubicsplineB 样条插值lanczosLanczos 窗函数 sinc 插值average对所有非 NODATA 贡献像素做加权平均rms所有非 NODATA 贡献像素的均方根 / 二次均值mode众数选取出现次数最多的值并列时取首个识别到的众数max所有非 NODATA 贡献像素的最大值min所有非 NODATA 贡献像素的最小值med所有非 NODATA 贡献像素的中位数q1所有非 NODATA 贡献像素的第一四分位数q3所有非 NODATA 贡献像素的第三四分位数sum所有非 NODATA 贡献像素的加权和其中average、rms、mode、max、min、med、q1、q3、sum属于聚合类重采样适用于降采样场景例如把高分辨率输入写入低分辨率目标格网时。官方示例即使用-r cubic进行三次卷积插值。--wo, --warp-option NAMEVALUE设置 warp 选项可多次指定所有可用键由GDALWarpOptions::papszWarpOptions文档定义。源码中该参数在 gdalalg_raster_reproject.cpp 注册并利用GDALWarpGetOptionList()提供命令行自动补全。典型用途例如--wo SKIP_NOSOURCEYES源数据未覆盖的区域跳过、不写入这正是测试用例中演示的用法。标准选项gdal raster update还支持两类适用于所有gdal系列命令的标准选项--if, --input-format format指定尝试打开输入文件的驱动/格式名。通常无需指定但在自动驱动探测失败或选错驱动时可用于跳过自动探测。可以重复指定多个候选驱动注意它并不强制使用这些驱动打开数据部分驱动对文件扩展名有要求。详见 doc/source/programs/gdal_options/if.rst。--oo, --open-option NAMEVALUE数据集打开选项格式相关可重复。详见 doc/source/programs/gdal_options/oo.rst。Overview 同步刷新机制这是gdal raster update区别于裸 GDALWarp的重要增值行为。默认情况下命令会在像素更新完成后刷新目标数据集已存在的概览且只刷新受影响区域而非全量重建。源码中的实现逻辑gdalalg_raster_update.cpp检查目标数据集第 1 波段是否有 overview 且未禁用刷新。读取输入数据集的仿射变换geotransform计算其四个角点的地理范围。若输入与输出都有 CRS 或都无 CRS则创建OGRCoordinateTransformation调用TransformBounds()将输入范围变换到目标 CRS 下得到待刷新包围盒overviewRefreshBBox。若上述条件不满足例如输入无 CRS 而输出有 CRS则发出警告Overviews can not be updatedoverviews 保持原样。刷新动作委托给GDALRasterOverviewAlgorithmRefresh并把bbox与resampling参数传入实现按包围盒部分刷新。bbox部分的刷新逻辑PartialRefreshFromBBOX定义在 gdalalg_raster_overview_refresh.cpp将地理包围盒反算为像素窗口后调用GDALRegenerateOverviewsMultiBand()仅重建该窗口内的概览像素。进度回调上源码将总进度拆为两段warp 阶段占 0.00.75overview 刷新阶段占 0.751.0见 gdalalg_raster_update.cpp。测试test_gdalalg_raster_update_refresh_overviewstest_gdalalg_raster_update.py验证了该机制目标 1024×1024 含 2 倍概览输入 2×2 像素写入中心区域后概览对应位置也被更新为 76由 255 与 127 混合的平均结果而概览外围仍保留旧值 127。这说明概览刷新确实是局部且与主分辨率像素更新一致的。在gdal raster pipeline中组合使用GDAL 3.13从 GDAL 3.13 起update可作为 raster pipeline 的一个步骤其输入数据集来自上一环节的输出。在 gdalalg_raster_update.h 中该类同时声明了CanBeLastStep()返回true——update 可以作为管线的最后一步此时直接写目标数据集CanBeMiddleStep()返回true——update 也可以作为中间步骤其输出更新后的数据集继续流入后续环节IsNativelyStreamingCompatible()返回false——该步骤需要物化materialize数据不适合流式按块求值OutputDatasetAllowedBeforeRunningStep()返回true——目标数据集可在步骤运行前由外部如--output提供。管线用法以!分隔各步骤规则详见 doc/source/programs/gdal_raster_pipeline.rst其中第 322 行起即为 update 步骤的说明入口。测试用例演示了两种典型形态作为中间步骤gdal raster pipeline read src.tif ! update dst.tif ! info --stats对应测试test_gdalalg_raster_update_pipeline_intermediate_steptest_gdalalg_raster_update.pyupdate 把 1×1 的 src 写入 3×3 的 dst.tif 后info --stats读取的波段最大值为 1。作为最后一步gdal raster pipeline read src.tif ! update --output dst.tif对应测试test_gdalalg_raster_update_pipeline_last_steptest_gdalalg_raster_update.py。注意由于update不是原生流式兼容步骤在管线中使用时建议配合 doc/source/programs/gdal_raster_materialize.rst 文档中提到的materialize步骤按需物化中间结果以获得更好的性能。完整使用示例基本用法三次插值更新官方文档示例gdal_raster_update.rst$ gdal raster update -r cubic in.tif out.tif将in.tif的像素内容更新进已存在的out.tif使用 cubic三次卷积插值完成重采样/重投影。局部范围更新WKT 多边形裁剪$ gdal raster update \ --geometry POLYGON ((1 -1,2 -1,2 -2,1 -2,1 -1)) \ in.tif out.tif仅把几何覆盖范围内的像素更新几何外像素原样保留。局部范围更新GeoJSON指定 CRS$ gdal raster update \ --geometry {type:Polygon,coordinates:[[[...]]]} \ --geometry-crs EPSG:32631 \ in.tif out.tifGeoJSON 几何默认按 WGS84 解释若坐标实际为投影坐标务必用--geometry-crs显式声明。更新并跳过无源区域$ gdal raster update -r cubic \ --wo SKIP_NOSOURCEYES \ in.tif out.tifSKIP_NOSOURCEYES使输入未覆盖的像素保持目标原有值不被覆盖在更新已有数据场景中很常用测试test_gdalalg_raster_update_gtiff即使用了该组合。更新但不刷新概览$ gdal raster update --no-update-overviews in.tif out.tif作为管线步骤$ gdal raster pipeline read in.tif ! update dst.tif ! info --stats返回状态码依据共享文档 doc/source/programs/return_code.rst 的约定程序成功执行返回状态码0出错时返回非零值。以警告形式发出的非阻塞错误non-blocking errors emitted as warnings被视为成功执行——例如Overviews can not be updated这类警告不会导致命令失败测试test_gdalalg_raster_update_cannot_update_overviews中同时断言了告警与失败两种状态的分野。常见错误与规避建议结合源码校验逻辑gdalalg_raster_update.cpp与测试用例以下几点需要特别留意源与目标不能是同一数据集。无论是以相同数据集指针传入test_gdalalg_raster_update_src_is_same_as_dst_pointer还是源/目标文件名相同test_gdalalg_raster_update_src_is_same_as_dst_filename都会触发Source and destination datasets must be different错误。源码中不仅比较数据集指针还会比较非 MEM 驱动的文件名杜绝原地自更新。目标必须已存在且可更新。update不会创建目标文件目标不存在时命令失败test_gdalalg_raster_update_dst_does_not_exist。目标文件必须支持以更新模式打开。裁剪几何必须合法。无法解析为 WKT 或 GeoJSON 的字符串会直接报错见test_gdalalg_raster_update_geometry_invalid。概览刷新有前提条件。当输入/输出 CRS 情况无法建立坐标变换例如一个有 CRS、一个没有overviews 不会被刷新并给出警告。此时可先统一 CRS或接受警告后另行使用 doc/source/programs/gdal_raster_overview_refresh.rst 中gdal raster overview refresh手动全量刷新。小结gdal raster update为向既有栅格数据集写入新像素这一高频需求提供了官方一等公民命令底层复用成熟的 GDALWarp 重投影/重采样框架上层自动处理裁剪几何、概览局部刷新并可通过--no-update-overviews、--wo、--to、--et等选项精细控制行为自 GDAL 3.13 起还能无缝嵌入gdal raster pipeline组合复杂处理链。对影像拼接修补、瓦片级局部更新、传感器数据回填等场景它是比手工组合gdalwarpgdal addo更简洁可靠的替代方案。若需深入了解各 warp 选项的完整清单可进一步查阅 gdalalg_raster_update.cpp 中选项到GDALWarpAppOptions的映射逻辑以及 gdalalg_raster_reproject.cpp 中共享的参数注册实现。赞分享GIS遥感数据工程【免费下载链接】gdalGDAL is an open source MIT licensed translator library for raster and vector geospatial data formats.项目地址https://gitcode.com/gh_mirrors/gd/gdal点击查看免费下载相关推荐GDAL gdal raster overview add 命令完全指南为栅格数据集构建、刷新与复用概览金字塔GDAL gdal raster overview add 命令完全指南为栅格数据集构建、刷新与复用概览金字塔 gdal raster overview adGIS遥感数据工程GDAL gdal raster reproject 子命令完全指南栅格重投影原理、参数详解与实战GDAL gdal raster reproject 子命令完全指南栅格重投影原理、参数详解与实战 本文以 GDAL 3.11 起引入的新一代子命令 gdalGIS遥感数据工程GDAL gdal raster proximity 子命令完全指南生成栅格邻近度地图GDAL gdal raster proximity 子命令完全指南生成栅格邻近度地图 gdal raster proximity 是 GDAL3.12 起GIS遥感数据工程上一篇清华PPT模板终极指南5分钟打造专业学术演示文稿下一篇如何快速掌握League Akari英雄联盟玩家的效率提升完整指南创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考
GIS遥感数据工程【免费下载链接】gdalGDAL is an open source MIT licensed translator library for raster and vector geospatial data formats.项目地址https://gitcode.com/gh_mirrors/gd/gdal点击查看免费下载gdal raster update是 GDAL 自 3.12 版本引入的栅格更新子命令用于将输入栅格的像素值写入一个已经存在的目标栅格数据集并在输入/输出坐标系不一致时自动完成重投影目标数据集的 extent、尺寸、分辨率与坐标系在整个过程中保持不变。本文基于 doc/source/programs/gdal_raster_update.rst 官方文档并结合 apps/gdalalg_raster_update.cpp 源码实现与 autotest/utilities/test_gdalalg_raster_update.py 测试用例完整讲解该命令的全部选项、重采样方法、概览overview同步刷新机制、gdal raster pipeline中的组合用法以及常见错误与规避方式帮助读者在影像更新、镶嵌修补、局部替换等场景中直接上手。命令概览做什么、不做什么gdal raster update的官方用途一句话即可概括Update the destination raster with the content of the input one.用输入栅格的内容更新目标栅格这一描述同样作为GDALRasterUpdateAlgorithm的DESCRIPTION常量定义在 gdalalg_raster_update.h 中。其核心行为特征如下就地更新目标数据集必须已存在且以**更新模式GA_Update**打开。源码 gdalalg_raster_update.cpp 中通过CPLAssert(poDstDS-GetAccess() GA_Update)强制校验这一点。自动重投影当输入与输出数据集的坐标系CRS不同时像素值会在写入前被重投影到目标 CRS 下的正确位置。不改变目标栅格结构目标数据集的 extent、尺寸、分辨率、CRS 均不被修改——这是与gdalwarp生成新文件最大的区别update只负责把新值填进已有的像素格子。概览自动刷新命令默认会同步刷新目标数据集已存在的 overviews内部机制复用gdal raster overview refresh除非显式指定--no-update-overviews。版本与能力演进该命令自 GDAL 3.12 起提供从 GDAL 3.13 开始它可以作为gdal raster pipeline中的一个步骤使用此时其输入数据集来自上一环节的输出。从命令整体架构看gdal raster update归属于gdal raster这一栅格子命令家族gdal raster入口自 3.11 引入并列于convert、clip、reproject、mosaic等子命令详见 doc/source/programs/gdal_raster.rst。底层实现基于 GDALWarp 的轻量封装尽管文档以命令行角度描述理解其底层实现有助于预判行为。在 gdalalg_raster_update.cpp 的RunStep()中算法做了以下事情获取输入数据集poSrcDS与目标数据集poDstDS校验二者不同详见下文常见错误。若指定了--geometry通过GetClipGeometry()继承自GDALClipCommon见 gdalalg_clip_common.h解析裁剪几何并导出为 WKT 作为-cutline选项。将--resampling-r、--warp-option-wo、--transform-option-to、--error-threshold-et等参数统一组装成CPLStringList交给GDALWarpAppOptionsNew()构建 warp 选项。调用核心函数GDALWarp(nullptr, hDstDS, 1, hSrcDS, psOptions, nullptr)即以目标数据集为输出句柄执行 GDALWarp——这正是写入已有数据集的实现关键。warp 完成后若目标带 1 有 overview 且未禁用更新则调用GDALRasterOverviewAlgorithmRefresh对受影响区域做部分刷新。可见gdal raster update本质上是GDALWarp面向原地更新场景的高层封装重投影、重采样、裁剪等几何能力全部来自 warp 框架。程序专属选项详解以下选项为gdal raster update专属均可同时用于命令行与 pipeline 步骤中。--et, --error-threshold ERROR-THRESHOLD变换近似的误差阈值以输入像素数为单位。默认值为0.125像素即允许近似变换引入不超过 1/8 像素的误差。例外情况当指定了RPC_DEM变换选项时默认切换为精确变换即--error-threshold0。传0表示使用精确变换不做近似逼近精度最高但计算更慢。从源码看该参数经由GDALRasterReprojectUtils::AddWarpOptTransformOptErrorThresholdArg()注册在 gdalalg_raster_reproject.cpp 中限定取值必须 ≥ 0。--geometry WKT_or_GeoJSON以 WKT 或 GeoJSON 字符串形式指定一个多边形或多多边形将更新操作限制在该几何范围内几何之外的像素保持不变。CRS 判定规则官方文档明确规定若输入几何为GeoJSON默认假定其坐标为 WGS84经纬度除非 GeoJSON 几何体内自带 CRS 定义或显式给出--geometry-crs。若输入几何为WKT默认假定其坐标系为输入数据集的 CRS除非显式给出--geometry-crs。坐标轴遵循GIS 友好约定X 为经度或东向eastingY 为纬度或北向northing。测试用例 test_gdalalg_raster_update.py 演示了geometryPOLYGON ((1 -1,2 -1,2 -2,1 -2,1 -1))的用法在 3×3 全零目标上仅在多边形覆盖的中央像素写入 1。若传入无法解析的字符串则报错Clipping geometry is neither a valid WKT or GeoJSON geometry见测试test_gdalalg_raster_update_geometry_invalid。--geometry-crs CRS指定--geometry坐标值所在的 CRS。若不指定按上一节规则推断。几何会从 geometry-crs 重投影到被更新数据集的 CRS后再执行裁剪。--no-update-overviews禁止更新目标数据集已存在的概览。默认行为是自动刷新 overviews这在目标数据文件较大、且更新范围很小时会额外耗时若希望只写全分辨率像素、概览稍后统一重建可加此开关。--to, --transform-option NAMEVALUE设置变换器选项适用于GDALCreateGenImgProjTransformer2()接受的参数RPC 相关的变换选项参考GDALCreateRPCTransformerV2()的文档。可多次指定。源码中该参数同样由AddWarpOptTransformOptErrorThresholdArg()注册gdalalg_raster_reproject.cpp并按NAMEVALUE键值对解析。测试用例给出了实际用法transform_option{SRC_METHOD: GEOTRANSFORM}。-r, --resampling RESAMPLING重采样方法默认nearest。可选值及其说明来自共享文档 doc/source/programs/gdal_options/warp_resampling.rst在源码 gdalalg_raster_reproject.cpp 中注册了完整的候选集合方法说明nearest最近邻默认最快插值质量最差bilinear双线性插值cubic三次卷积插值cubicsplineB 样条插值lanczosLanczos 窗函数 sinc 插值average对所有非 NODATA 贡献像素做加权平均rms所有非 NODATA 贡献像素的均方根 / 二次均值mode众数选取出现次数最多的值并列时取首个识别到的众数max所有非 NODATA 贡献像素的最大值min所有非 NODATA 贡献像素的最小值med所有非 NODATA 贡献像素的中位数q1所有非 NODATA 贡献像素的第一四分位数q3所有非 NODATA 贡献像素的第三四分位数sum所有非 NODATA 贡献像素的加权和其中average、rms、mode、max、min、med、q1、q3、sum属于聚合类重采样适用于降采样场景例如把高分辨率输入写入低分辨率目标格网时。官方示例即使用-r cubic进行三次卷积插值。--wo, --warp-option NAMEVALUE设置 warp 选项可多次指定所有可用键由GDALWarpOptions::papszWarpOptions文档定义。源码中该参数在 gdalalg_raster_reproject.cpp 注册并利用GDALWarpGetOptionList()提供命令行自动补全。典型用途例如--wo SKIP_NOSOURCEYES源数据未覆盖的区域跳过、不写入这正是测试用例中演示的用法。标准选项gdal raster update还支持两类适用于所有gdal系列命令的标准选项--if, --input-format format指定尝试打开输入文件的驱动/格式名。通常无需指定但在自动驱动探测失败或选错驱动时可用于跳过自动探测。可以重复指定多个候选驱动注意它并不强制使用这些驱动打开数据部分驱动对文件扩展名有要求。详见 doc/source/programs/gdal_options/if.rst。--oo, --open-option NAMEVALUE数据集打开选项格式相关可重复。详见 doc/source/programs/gdal_options/oo.rst。Overview 同步刷新机制这是gdal raster update区别于裸 GDALWarp的重要增值行为。默认情况下命令会在像素更新完成后刷新目标数据集已存在的概览且只刷新受影响区域而非全量重建。源码中的实现逻辑gdalalg_raster_update.cpp检查目标数据集第 1 波段是否有 overview 且未禁用刷新。读取输入数据集的仿射变换geotransform计算其四个角点的地理范围。若输入与输出都有 CRS 或都无 CRS则创建OGRCoordinateTransformation调用TransformBounds()将输入范围变换到目标 CRS 下得到待刷新包围盒overviewRefreshBBox。若上述条件不满足例如输入无 CRS 而输出有 CRS则发出警告Overviews can not be updatedoverviews 保持原样。刷新动作委托给GDALRasterOverviewAlgorithmRefresh并把bbox与resampling参数传入实现按包围盒部分刷新。bbox部分的刷新逻辑PartialRefreshFromBBOX定义在 gdalalg_raster_overview_refresh.cpp将地理包围盒反算为像素窗口后调用GDALRegenerateOverviewsMultiBand()仅重建该窗口内的概览像素。进度回调上源码将总进度拆为两段warp 阶段占 0.00.75overview 刷新阶段占 0.751.0见 gdalalg_raster_update.cpp。测试test_gdalalg_raster_update_refresh_overviewstest_gdalalg_raster_update.py验证了该机制目标 1024×1024 含 2 倍概览输入 2×2 像素写入中心区域后概览对应位置也被更新为 76由 255 与 127 混合的平均结果而概览外围仍保留旧值 127。这说明概览刷新确实是局部且与主分辨率像素更新一致的。在gdal raster pipeline中组合使用GDAL 3.13从 GDAL 3.13 起update可作为 raster pipeline 的一个步骤其输入数据集来自上一环节的输出。在 gdalalg_raster_update.h 中该类同时声明了CanBeLastStep()返回true——update 可以作为管线的最后一步此时直接写目标数据集CanBeMiddleStep()返回true——update 也可以作为中间步骤其输出更新后的数据集继续流入后续环节IsNativelyStreamingCompatible()返回false——该步骤需要物化materialize数据不适合流式按块求值OutputDatasetAllowedBeforeRunningStep()返回true——目标数据集可在步骤运行前由外部如--output提供。管线用法以!分隔各步骤规则详见 doc/source/programs/gdal_raster_pipeline.rst其中第 322 行起即为 update 步骤的说明入口。测试用例演示了两种典型形态作为中间步骤gdal raster pipeline read src.tif ! update dst.tif ! info --stats对应测试test_gdalalg_raster_update_pipeline_intermediate_steptest_gdalalg_raster_update.pyupdate 把 1×1 的 src 写入 3×3 的 dst.tif 后info --stats读取的波段最大值为 1。作为最后一步gdal raster pipeline read src.tif ! update --output dst.tif对应测试test_gdalalg_raster_update_pipeline_last_steptest_gdalalg_raster_update.py。注意由于update不是原生流式兼容步骤在管线中使用时建议配合 doc/source/programs/gdal_raster_materialize.rst 文档中提到的materialize步骤按需物化中间结果以获得更好的性能。完整使用示例基本用法三次插值更新官方文档示例gdal_raster_update.rst$ gdal raster update -r cubic in.tif out.tif将in.tif的像素内容更新进已存在的out.tif使用 cubic三次卷积插值完成重采样/重投影。局部范围更新WKT 多边形裁剪$ gdal raster update \ --geometry POLYGON ((1 -1,2 -1,2 -2,1 -2,1 -1)) \ in.tif out.tif仅把几何覆盖范围内的像素更新几何外像素原样保留。局部范围更新GeoJSON指定 CRS$ gdal raster update \ --geometry {type:Polygon,coordinates:[[[...]]]} \ --geometry-crs EPSG:32631 \ in.tif out.tifGeoJSON 几何默认按 WGS84 解释若坐标实际为投影坐标务必用--geometry-crs显式声明。更新并跳过无源区域$ gdal raster update -r cubic \ --wo SKIP_NOSOURCEYES \ in.tif out.tifSKIP_NOSOURCEYES使输入未覆盖的像素保持目标原有值不被覆盖在更新已有数据场景中很常用测试test_gdalalg_raster_update_gtiff即使用了该组合。更新但不刷新概览$ gdal raster update --no-update-overviews in.tif out.tif作为管线步骤$ gdal raster pipeline read in.tif ! update dst.tif ! info --stats返回状态码依据共享文档 doc/source/programs/return_code.rst 的约定程序成功执行返回状态码0出错时返回非零值。以警告形式发出的非阻塞错误non-blocking errors emitted as warnings被视为成功执行——例如Overviews can not be updated这类警告不会导致命令失败测试test_gdalalg_raster_update_cannot_update_overviews中同时断言了告警与失败两种状态的分野。常见错误与规避建议结合源码校验逻辑gdalalg_raster_update.cpp与测试用例以下几点需要特别留意源与目标不能是同一数据集。无论是以相同数据集指针传入test_gdalalg_raster_update_src_is_same_as_dst_pointer还是源/目标文件名相同test_gdalalg_raster_update_src_is_same_as_dst_filename都会触发Source and destination datasets must be different错误。源码中不仅比较数据集指针还会比较非 MEM 驱动的文件名杜绝原地自更新。目标必须已存在且可更新。update不会创建目标文件目标不存在时命令失败test_gdalalg_raster_update_dst_does_not_exist。目标文件必须支持以更新模式打开。裁剪几何必须合法。无法解析为 WKT 或 GeoJSON 的字符串会直接报错见test_gdalalg_raster_update_geometry_invalid。概览刷新有前提条件。当输入/输出 CRS 情况无法建立坐标变换例如一个有 CRS、一个没有overviews 不会被刷新并给出警告。此时可先统一 CRS或接受警告后另行使用 doc/source/programs/gdal_raster_overview_refresh.rst 中gdal raster overview refresh手动全量刷新。小结gdal raster update为向既有栅格数据集写入新像素这一高频需求提供了官方一等公民命令底层复用成熟的 GDALWarp 重投影/重采样框架上层自动处理裁剪几何、概览局部刷新并可通过--no-update-overviews、--wo、--to、--et等选项精细控制行为自 GDAL 3.13 起还能无缝嵌入gdal raster pipeline组合复杂处理链。对影像拼接修补、瓦片级局部更新、传感器数据回填等场景它是比手工组合gdalwarpgdal addo更简洁可靠的替代方案。若需深入了解各 warp 选项的完整清单可进一步查阅 gdalalg_raster_update.cpp 中选项到GDALWarpAppOptions的映射逻辑以及 gdalalg_raster_reproject.cpp 中共享的参数注册实现。赞分享GIS遥感数据工程【免费下载链接】gdalGDAL is an open source MIT licensed translator library for raster and vector geospatial data formats.项目地址https://gitcode.com/gh_mirrors/gd/gdal点击查看免费下载相关推荐GDAL gdal raster overview add 命令完全指南为栅格数据集构建、刷新与复用概览金字塔GDAL gdal raster overview add 命令完全指南为栅格数据集构建、刷新与复用概览金字塔 gdal raster overview adGIS遥感数据工程GDAL gdal raster reproject 子命令完全指南栅格重投影原理、参数详解与实战GDAL gdal raster reproject 子命令完全指南栅格重投影原理、参数详解与实战 本文以 GDAL 3.11 起引入的新一代子命令 gdalGIS遥感数据工程GDAL gdal raster proximity 子命令完全指南生成栅格邻近度地图GDAL gdal raster proximity 子命令完全指南生成栅格邻近度地图 gdal raster proximity 是 GDAL3.12 起GIS遥感数据工程上一篇清华PPT模板终极指南5分钟打造专业学术演示文稿下一篇如何快速掌握League Akari英雄联盟玩家的效率提升完整指南创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考
