1

我有一个具有PALETTE_COLOR光度解释的 RLE 压缩 DICOM 文件。使用 GDCM，我可以使用以下代码获取片段：

gdcm.ImageReader imagereader = new gdcm.ImageReader();
imagereader.SetFileName(fileName);
if (imagereader.Read())
  {
  DataElement compressedPixelData = imagereader.GetFile().GetDataSet().GetDataElement(Helper.PIXEL_DATA);
  SequenceOfFragments sf = compressedPixelData.GetSequenceOfFragments();
  if (sf == null) throw new Exception("Cannot get fragments!");
  Fragment frag = sf.GetFragment((uint)frameNumber);
  uint bufLen = Convert.ToUInt32(frag.GetByteValue().GetLength().toString());
  byte[] buffer = new byte[bufLen];
  frag.GetByteValue().GetBuffer(buffer, bufLen);
  }

现在我正在尝试从缓冲区创建图像。既然是PALETTE_COLOR，我必须将 LUT 应用到缓冲区。我使用这段代码：

gdcm.LookupTable lut = imagereader.GetImage().GetLUT();

int size = ImageWidth * ImageHeight * 3;
byte[] decodedBuffer = new byte[size];

bool worked = lut.Decode(decodedBuffer, (uint)size, buffer, (uint)bufLen);

if (worked)
   return decodedBuffer;
else
   return buffer;

decodedBuffer 大小是Width * Height * 3因为我在应用 LUT 后期望 RGB 像素。但生成的图像不正确。

如果我使用ImageApplyLookupTable类，我可以正确显示图像。

我不想使用ImageApplyLookupTable类，因为它会将整个图像（所有片段！）解码为原始 RGB 像素并消耗大量内存。我想逐帧解码以最小化内存使用量。

您能否为我指出如何gdcm.LookupTable正确使用类一次解码一帧的正确方向？我使用的示例文件在这里。

更新：使用ImageRegionReader适用于 8 位调色板颜色，但不适用于 16 位，您能检查一下原因吗？我这里有 16 位的示例。

4

1 回答 1

0

GDCM 附带了一堆示例，在您的情况下，我建议您开始阅读：ExtractImageRegionWithLUT.cs。

该代码非常愚蠢，因为它将在前一个帧上一遍又一遍地写入提取的帧。但无论如何，这应该给你基本的想法。这是我在这里所做的：

首先确保您的 C# 示例已构建（我使用 make 在我的情况下）：

$ make GDCMCSharpExample

然后：

$ ExtractImageRegionWithLUT.exe PAL-16_RLE.dcm

在 UNIX 上，这将创建一个原始文件，您可以使用以下命令将其转换回 DICOM gdcmimg：

$ gdcmimg --depth 16 --spp 3 --size 800,600 /tmp/frame_rgb.raw /tmp/frame_rgb.raw.dcm

然后您可以简单地查看它：

$ gdcmviewer /tmp/frame_rgb.raw.dcm

您甚至可以使用以下方法来说服自己：

$ gdcmimg -C 1.2.840.10008.5.1.4.1.1.3.1 --depth 16 --spp 3 --size 800,600 /tmp/frame_rgb.raw /tmp/frame_rgb.raw.dcm
于 2016-10-27T10:18:35.860 回答