QGIS中tif切片导出为空白或缺失瓦片怎么办？

Question

在QGIS中使用“导出为瓦片（XYZ）”功能导出TIFF影像切片时，常出现部分或全部瓦片为空白（纯黑/白/灰）、缺失、分辨率异常或404错误。常见原因包括：①源TIFF未正确设置CRS或存在坐标系不匹配；②未启用“渲染图层”或图层被意外关闭/透明度设为100%；③金字塔（Overviews）缺失导致缩放层级瓦片生成失败；④输出路径含中文或特殊字符，引发GDAL写入中断；⑤内存不足或瓦片范围超出原始影像有效区域（如误选全球范围）。此外，QGIS 3.28+版本中“导出为瓦片”工具默认跳过空瓦片（优化选项），若未勾选“保留空瓦片”，会导致逻辑上“缺失”。建议操作前验证图层可见性、执行“构建金字塔”、统一使用WGS84/Pseudo-Mercator投影，并优先选用处理大图更稳定的插件（如QTiles或GDAL命令行`gdal_translate + gdal_tile`).

Answer 1

```html

一、现象层：瓦片导出异常的典型表现

• 部分层级（如z=10–14）出现大面积纯黑/纯白/50%灰瓦片（255,255,255 或 0,0,0）
• 浏览器加载时返回 404 Not Found（尤其在非整数缩放级别或边缘区域）
• 瓦片尺寸不一致（如本应为256×256，却生成128×128或全零字节空文件）
• QGIS日志中反复出现 GDAL: Unhandled exception in GDALDataset::RasterIO 或 Failed to compute overview

二、配置层：核心参数与隐式依赖关系

QGIS“导出为瓦片（XYZ）”并非原子操作，其行为受以下5个关键配置链式影响：

配置项	默认值（QGIS 3.28+）	风险触发条件
CRS of output tiles	Pseudo-Mercator (EPSG:3857)	源TIFF为WGS84（EPSG:4326）但未重投影→坐标拉伸/裁剪失真
Extent source	Canvas extent	误用“Full layer extent”导致超出影像有效范围→生成大量空瓦片
Empty tiles	✅ Skip empty tiles（默认启用）	未勾选“Preserve empty tiles”→逻辑缺失（非物理丢失）

三、数据层：TIFF元数据完整性诊断流程

# 使用GDAL验证TIFF健康度（终端执行）
gdalinfo -stats -proj4 your_image.tif | grep -E "(PROJCS|GEOGCS|OVERVIEWS|Band|Size)"
# 关键输出应包含：
#   Projection: PROJCS["WGS 84 / Pseudo-Mercator",...]
#   Overviews: 2x, 4x, 8x, 16x, 32x
#   Band 1 Block=256x256 Type=UInt16, ColorInterp=Gray

四、执行层：QGIS内置工具失效时的降级路径

graph TD
  A[QGIS XYZ Tiles Export] -->|失败| B{检查金字塔}
  B -->|缺失| C[右键图层 → Properties → Pyramids → Build overviews]
  B -->|存在| D[切换至QTiles插件]
  D --> E[支持自定义压缩/空值填充/并发线程]
  A -->|仍失败| F[GDAL命令行方案]
  F --> G[gdal_translate -of GTiff -co TILED=YES -co COMPRESS=LZW]
  F --> H[gdal_tile -zoom 0-15 -webviewer all -tileIndex index.shp]

五、工程层：生产环境推荐工作流（含容错设计）

1. 【预处理】用 gdalwarp -t_srs EPSG:3857 -r bilinear -co TILED=YES -co COMPRESS=DEFLATE 统一重投影并分块
2. 【增强】执行 gdaladdo -r average --config COMPRESS_OVERVIEW DEFLATE input_3857.tif 2 4 8 16 32
3. 【导出】在QGIS中加载重投影后TIFF → 确认图层可见性/透明度=0% → 勾选“Preserve empty tiles”
4. 【验证】用Python脚本批量校验瓦片：
   import os, numpy as np; for f in os.listdir('tiles/z12'): arr = np.fromfile(f, dtype=np.uint8); assert arr.mean() > 5 and arr.std() > 10

六、进阶层：GDAL 3.8+内存映射与并发优化策略

针对>10GB TIFF，需绕过QGIS GUI内存瓶颈：

• 启用GDAL虚拟内存：设置环境变量 GDAL_CACHEMAX=2048（单位MB）
• 强制内存映射：在gdal_translate中添加 -oo NUM_THREADS=ALL_CPUS
• 瓦片空值智能填充：使用-srcnodata 0 -a_nodata 0 -dstalpha避免黑边传播

七、监控层：自动化诊断清单（Checklist）

检查项	验证命令/操作	合格标准
CRS一致性	QGIS图层属性→Source→Coordinate Reference System	显示EPSG:3857且无警告图标
金字塔状态	图层右键→Properties→Pyramids→Overview levels	显示≥5级且状态为“Built”
输出路径合规性	检查路径是否含中文/空格/emoji	仅含ASCII字母、数字、下划线、连字符

```