处理 150GB .jp2 图像

Question

我下载了一个 150GB 的卫星 .jp2 图像，我一次只想要其中的一小部分。 我将如何将图像平铺为可管理的块？只需提取图像的一部分就足够了。

由于我对 Python 只是有点熟悉，所以我查看了 Pillow 和 OpenCV 库，但由于图像分辨率超出了它们的限制而没有成功。 我还研究了 Openslide for Python，但无法消除错误（无法找到模块“libopenslide-0.dll”）。

Answer 1

libvips 可以高效处理大图片。

例如，对于这个 2.8gb 测试图像：

$ vipsheader 9235.jp2
9235.jp2: 107568x79650 uchar, 4 bands, srgb, jp2kload
$ ls -l 9235.jp2
-rw-r--r-- 1 john john 2881486848 Mar  1 22:37 9235.jp2

我明白了：

$ /usr/bin/time -f %M:%e \
    vips crop 9235.jp2 x.jpg 10000 10000 1000 1000
190848:0.45

因此，它在 0.5 秒内从 110,000 x 80,000 像素 jp2 图像中取出 1,000 x 1,000 像素块，并且需要不到 200mb 的内存。

有 python、ruby、node 等的绑定，因此您不必使用 CLI。

在 python 中你可以这样写：

import pyvips

image = pyvips.Image.new_from_file("9235.jp2")
tile = image.crop(10000, 10000, 1000, 1000)
tile.write_to_file("x.jpg")

它确实有点依赖于你正在读取的 jp2 图像。有些是直到（!!）并且读出一个部分可能非常慢。

也有 Windows 二进制文件，请查看“下载”页面。

如果您使用的是 Linux，vipsdisp 可以非常快速地查看像这样的巨大图像。

Answer 2

Grok JPEG 2000 工具包
能够解压缩极大图像的区域，例如
您链接到的 150 GB 图像。

示例命令：

grk_decompress -i FOO.jp2 -o FOO.tif -d 10000,10000,15000,15000 -v

解压缩区域(10000,10000,15000,15000)代码> 为 TIFF 格式。