在上传到我的云存储之前,我需要在工作站中压缩dicom图像。我正在计划获取像素阵列并将其保存为jp2格式。然后,我从中创建PixelData并将其附加回dicom ds对象。我正在使用scipy保存和读取图像。但是,我看到正在创建的图像不一致。其中一些往往与原始图像不同。
import pydicom
import scipy.misc
ds = pydicom.dcmread('CT2_J2KR')
pixData1 = ds.pixel_array.tostring()
scipy.misc.imsave('xout.jp2', ds.pixel_array)
outarr = scipy.misc.imread('xout.jp2')
newPixelData = outarr.tostring()
这是压缩dicom图像的正确方法吗?如果没有,还有什么更好的方法?