JPEG2000 的 PythonMagick 数据提取

Question

我正在开发一个应用程序，该应用程序(除其他外)需要读取卫星图像(每个图像只有一个波段)并处理像素数据。

格式是 JPEG-2000，因此我不能使用 PIL 库(它简化了一切)。我找到了 PythonMagick 库，我可以完美地读取图像并提取像素值。但只有一个像素!

im=PythonMagick.Image(dirimage)      # (This is very slow....)
a=im.pixelColor(j-1,i-1).intensity() # the value intensity is extracted for one pixel
a=a/int(XML_var[37][2])              # the reflectance values are normalised to the range [0,1]

因此，我需要一个 for 循环来获取所有像素值(图像非常大)。我尝试使用 Blob 函数来获取数据，但它崩溃了。

有没有更好的选择？如何快速获取JPEG2000图像的像素数据并保存到数组中？

Answer 1

使用 Blob 应该可行:

import numpy
from PythonMagick import Image, Blob

i = Image('http://www.microimages.com/gallery/jp2/potholes2.jp2')
b = Blob()
i.write(b, 'GRAY')
a = numpy.fromstring(b.data, 'uint8').reshape((i.rows(), i.columns()))