c++ - DICOM 访问压缩数据 (DCMTK)

Question

我正在努力使用 DCMTK 3.6.1 库，现在我遇到了一种常见情况:访问 DICOM 图像的压缩数据像素。

正如图书馆的作者在这里建议的那样 http://support.dcmtk.org/redmine/projects/dcmtk/wiki/Howto_AccessingCompressedData ，这是获取压缩数据的正确方法。

此代码有效，但它们从文件开始创建数据集。在我的代码中，我以这种方式填充了我的数据集

status = data->putAndInsertUint8Array(DCM_PixelData, pSource, sizeBuffer);

其中 pSource 包含我未压缩的像素数据。在此之后，我添加图像 DICOM 标签(第 28 组)并使用

进行压缩

status = data->chooseRepresentation(EXS_JPEGProcess14SV1, &param);

从这个数据集开始，我想访问压缩数据

status = data->findAndGetElement(DCM_PixelData, element);
    DcmPixelData *dpix = NULL;
    dpix = OFstatic_cast(DcmPixelData*, element);
    /* Since we have compressed data, we must utilize DcmPixelSequence
    in order to access it in raw format, e. g. for decompressing it
    with an external library.
    */
    DcmPixelSequence *dseq = NULL;
    E_TransferSyntax xferSyntax = EXS_Unknown;
    const DcmRepresentationParameter *rep = NULL;
    // Find the key that is needed to access the right representation of the data within DCMTK
    dpix->getOriginalRepresentationKey(xferSyntax, rep);

    // Access original data representation and get result within pixel sequence
    status = dpix->getEncapsulatedRepresentation(xferSyntax, rep, dseq);
    Uint32 length;
    if (status.good())
    {
        DcmPixelItem* pixitem = NULL;
        // Access first frame (skipping offset table)
        dseq->getItem(pixitem, 1);
        if (pixitem == NULL)
            return 1;
        Uint8* pixData = NULL;
        // Get the length of this pixel item (i.e. fragment, i.e. most of the time, the lenght of the frame)
        length = pixitem->getLength();
        if (length == 0)
            return 1;
        // Finally, get the compressed data for this pixel item
        status = pixitem->getUint8Array(pixData);
        // Metto i Pixel Data compressi su pSorgCompr
        pSorgCompr = (LPBYTE)pixData;
    }
    ////////////////////////////

    DJEncoderRegistration::cleanup();
    DJDecoderRegistration::cleanup();

但是 status = dpix->getEncapsulatedRepresentation(xferSyntax, rep, dseq); 行返回失败并显示错误“未找到像素表示”，我不明白为什么。

奇怪的事情是，如果在访问压缩数据之前，我用fileformat.saveFile("compressedPixelData.dcm", EXS_JPEGProcess14SV1); 保存压缩文件，然后下一步我使用 result = fileformat.loadFile("compressedPixelData.dcm"); 加载文件，一切正常

好像loadfile函数解决了这个问题，不知道怎么办，也许是填充了一些标签？

我在调用 chooseRepresentation 函数之前填写的标签是:

• DCM_Rows
• DCM_Columns
• DCM_BitsStored
• DCM_SamplesPerPixel
• DCM_PlanarConfiguration
• DCM_HighBit
• DCM_BitsAllocated
• DCM_PixelRepresentation
• DCM_RescaleIntercept
• DCM_RescaleSlope
• DCM_PhotometricInterpretation
• DCM_PixelAspectRatio
• DCM_ImagerPixelSpacing
• DCM_PixelSpacing

Answer 1

在 DICOM 中，压缩帧的存储方式不同于像素数据元素 (7FE0, 0010) 下的未压缩(原始)帧。这是 DICOM 中压缩图像编码的基本概念。

封装后的像素流(压缩图像数据)在顶层数据集的像素数据元素(7FE0，0010)下被分割成一个或多个片段(Items)。在封装编码中，像素数据元素是一个序列，它包含两个或多个项目元素(Fragments)，每个片段(Item Element)都有自己明确的长度。此外，封装格式支持单帧和多帧图像编码。一个框架可以完全包含在一个片段(Item)中，也可以跨越多个片段。封装像素流的 Fragment 序列以定界符结束。
编码像素数据流之前的项目序列中的第一个项目(FFFE，E000)通常为空，或者可能包含一个基本偏移表项目，该项目保存每个帧的第一个片段的项目标签的第一个字节的字节偏移量在项目序列中。

因此，当您想从 DICOM Data Set 中提取压缩像素流时，您需要跳过 Pixel Data 序列下的第一个 Item 元素。我希望这有助于理解您引用的文档。请参阅 DICOM 标准 PS 3.5 附件 A 了解更多信息。