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); 返回错误 "Pixel Representation not found" 失败，我不明白为什么。

悲剧是如果在访问压缩数据之前，我用fileformat.saveFile("compressedPixelData.dcm"("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

好吧，也许我已经解决了问题。我已经添加了说明

dataset->removeAllButCurrentRepresentations();

在访问压缩像素数据之前。

我也可以 PixelData->removeAllButCurrentRepresentations(); 而不是之前的指令，并且以相同的方式工作....

但我真的不明白它为什么起作用....请你解释一下好吗？谢谢

Answer 2

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

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

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

DICOM 访问压缩数据 (DCMTK)

DICOM Accessing Compressed Data (DCMTK)

c++

compression

dicom

visual-studio-2013

dcmtk