如何准确裁剪并保存 NumPy 图像数组为 JPEG 格式？

本文共计963个文字，预计阅读时间需要4分钟。

相关专题内容如下：

本文详解因数组索引顺序错误导致 pil 无法保存空图像的常见问题，重点纠正 numpy 图像切片中行（高度）与列（宽度）维度的误用，并提供可直接运行的修复代码与最佳实践。

在使用 NumPy 处理图像（如从 .svs 文件读取的多维数组）并进行分块裁剪时，一个极易被忽视却致命的错误是数组切片维度顺序与图像坐标系不匹配。你的原始代码中：

croppedImage = img[XStart:YStart, XEnd:YEnd]

这一行存在双重逻辑错误：

1. 维度混淆：NumPy 数组 img.shape 返回 (height, width, depth)（即 (H, W, C)），对应索引应为 [row_start:row_end, col_start:col_end, ...]；而 XStart/XEnd 是基于宽度（列方向）计算的，YStart/YEnd 是基于高度（行方向）计算的。但你却将 XStart:YStart 用作第一个切片维度（本该是 YStart:YEnd），把 XEnd:YEnd 用作第二个维度（本该是 XStart:XEnd）；
2. 切片语法错误：a[i:j, k:l] 表示「第 i 到 j−1 行、第 k 到 l−1 列」，而非 a[i, j, k, l]。你写成 XStart:YStart 实际上是用列偏移去切行范围——当 XStart > YStart（例如 12 > 107 不成立，但后续循环中 XStart = 12 + 256*x 快速增大）时，会生成空切片（如 img[300:50, ...] → 空数组），导致 croppedImage.size == 0，PIL 因此抛出 ValueError: cannot write empty image as JPEG。

✅ 正确写法应严格遵循 (height, width, channels) 的索引顺序：

croppedImage = img[YStart:YEnd, XStart:XEnd] # ✅ 行范围在前，列范围在后

完整修复后的核心循环如下（含边界校验与格式兼容性处理）：

import numpy as np from PIL import Image as im import tifffile as tif # ...（加载图像、计算 excess pixels 等前置代码保持不变） for y in range(YNumberOfTiles): for x in range(XNumberOfTiles): XStart = leftExcessPixels + (tilesize * x) YStart = topExcessPixels + (tilesize * y) XEnd = min(XStart + tilesize, width) # 防止越界（虽此处理论不需，但强推荐） YEnd = min(YStart + tilesize, height) # ✅ 关键修正：先高度（Y），后宽度（X） croppedImage = img[YStart:YEnd, XStart:XEnd] # 可选：验证裁剪结果有效性 if croppedImage.size == 0: print(f"Warning: Empty tile at x={x}, y={y} -> skipped") continue # 确保数据类型兼容 PIL（JPEG 需 uint8） if croppedImage.dtype != np.uint8: # 假设输入是 uint16 或 float，做安全归一化（按需调整） if croppedImage.dtype == np.uint16: croppedImage = (croppedImage // 256).astype(np.uint8) elif croppedImage.dtype in (np.float32, np.float64): croppedImage = np.clip(croppedImage * 255, 0, 255).astype(np.uint8) try: pil_img = im.fromarray(croppedImage) pil_img.save(f'Sliced/sample_x{x}_y{y}.jpg', quality=95) except Exception as e: print(f"Failed to save tile x{x}_y{y}: {e}")

? 重要注意事项：

• 通道顺序：若原始图像是 BGR（如 OpenCV 加载）或带 alpha 通道，需显式转换（如 croppedImage = croppedImage[..., :3] 或 cv2.cvtColor(..., cv2.COLOR_BGR2RGB)）；
• 内存与性能：对超大图像（如 13K×20K），避免一次性加载全部 tiles 到内存，可考虑流式处理或使用 tifffile 的 page 参数分页读取；
• 文件格式选择：JPEG 不支持透明度且会压缩，若需无损保存，改用 'Sliced/...png' 并调用 save(..., format='PNG')；
• 调试技巧：在循环内添加 print(f"Tile {x},{y}: shape={croppedImage.shape}") 可快速定位空切片来源。

遵循上述修正与规范，即可稳定生成尺寸统一、内容完整的 JPEG 图像块，彻底规避“empty image”错误。