DCMTK 学习笔记一 之解析tag标签

Posted 2021-06-23 kevinlq

阅读本文大概需要 3.3 分钟

我们要解析显示一张 dcm图像，首先必须要解析其部分标签信息。

标签信息主要有三类：

• Meta info
• 普通 tag
• 像素 tag

根据需要进行读取

使用 DCMKT读取 dicom标签是很容易的一件事情

meata info tag读取

meta 信息在 0002组

有三种方式读取，这里展示最简单的一种

    DcmMetaInfo metainfo;
    OFCondition status;
    status = metainfo.loadFile("G:/dicomFile/3.dcm");

    if (status.good())
    {
        OFString sopClassUID, xferUID;
        if (metainfo.findAndGetOFString(DCM_MediaStorageSOPClassUID, sopClassUID).good())
            COUT << "SOP Class UID: " << sopClassUID << OFendl;
        if (metainfo.findAndGetOFString(DCM_TransferSyntaxUID, xferUID).good())
            COUT << "Transfer Syntax UID: " << xferUID << OFendl;

        metainfo.print(COUT);
    }

输出结果信息：

SOP Class UID: 1.2.840.10008.5.1.4.1.1.2
Transfer Syntax UID: 1.2.840.10008.1.2.4.51
# Dicom-Meta-Information-Header
# Used TransferSyntax: Little Endian Explicit
(0002,0000) UL 202                                      #   4, 1 FileMetaInformationGroupLength
(0002,0001) OB 00\\01                                    #   2, 1 FileMetaInformationVersion
(0002,0002) UI =CTImageStorage                          #  26, 1 MediaStorageSOPClassUID
(0002,0003) UI [1.2.276.0.7230010.3.1.4.2012715.17328.1511764633.841627] #  56, 1 MediaStorageSOPInstanceUID
(0002,0010) UI =JPEGExtended:Process2+4                 #  22, 1 TransferSyntaxUID
(0002,0012) UI [1.2.276.0.7230010.3.0.3.6.1]            #  28, 1 ImplementationClassUID
(0002,0013) SH [OFFIS_DCMTK_361]                        #  16, 1 ImplementationVersionName

官方参考资料

普通 tag读取

下面以读取病人姓名为例

    DcmFileFormat dcmFormat;

    OFCondition ofResult = dcmFormat.loadFile("G:/dicomFile/3.dcm");

    if (!ofResult.good())
    {
        qDebug() << "dcmFormat read fail..." << ofResult.text();
        return false;
    }

    DcmDataset *pDataSet = dcmFormat.getDataset();

    if (nullptr == pDataSet)
    {
        return false;
    }

    OFString strPatientName;
    ofResult = pDataSet->findAndGetOFString(DCM_PatientName, strPatientName);

    if (!ofResult.good())
    {
        return false;
    }

    qDebug() << "patientName:" << strPatientName.c_str();

结果输出

patientName: Zhang xxx

实际开发当中解析时需要分别对应不同类型，单独封装对应的解析方法

比如解析字符串tag

bool getStringTagValue(DcmDataset* pDataSet,const DcmTag& pTag, QString &value)
{
    if (nullptr == pDataSet)
    {
        return false;
    }

    const char* strValue = nullptr;
    Uint32 lLength = 0;
    OFCondition ofResult = pDataSet->findAndGetString(pTag, strValue, lLength);

    if (!ofResult.good())
    {
        return false;
    }

    value = QLatin1String(strValue, lLength);

    return true;
}

比如解析浮点数tag

bool getU16TagValue(DcmDataset* pDataSet,const DcmTag& pTag, quint16 &nValue)
{
    Uint16	uint16 = 0;

    if ( pDataSet->findAndGetUint16( pTag, uint16).good() )
    {
        nValue = uint16;

        return true;
    }

    return false;
}

其实 DcmFileFormat对象内部就包含了DcmMetaInfo和DcmDataset对象指针，在开始的构造函数中内部自动创建的，内部采用了双链表结构进行存储

DcmFileFormat::DcmFileFormat()
  : DcmSequenceOfItems(InternalUseTag),
    FileReadMode(ERM_autoDetect)
{
    DcmMetaInfo *MetaInfo = new DcmMetaInfo();
    DcmSequenceOfItems::itemList->insert(MetaInfo);
    MetaInfo->setParent(this);

    DcmDataset *Dataset = new DcmDataset();
    DcmSequenceOfItems::itemList->insert(Dataset);
    Dataset->setParent(this);
}

总结

当然了，上述方法均是采用 DCMTK开源库进行读取，自己造轮子也能一步一步读取，就是比较麻烦一点

GitHub上面有一个开源很久的项目 qdcm，也是基于Qt实现的，可以参考下

qdcm

比如读取文件，判断是否包含头信息代码

bool DcmFile::hasHeader(const QString &path)
{
    QFile file(path);
    if (!file.open(QFile::ReadOnly)) {
        return false;
    }

    QDataStream ds(&file);
    int s = ds.skipRawData(128);
    if (s != 128) {
        file.close();
        return false;
    }

    char signature[4];
    s = ds.readRawData(signature, 4);
    if (s != 4) {
        file.close();
        return false;
    }

    file.close();

    if (signature[0] == \'D\' && signature[1] == \'I\' && signature[2] == \'C\' && signature[3] == \'M\') {
        return true;
    }

    return false;
}

下一篇来学习下 Dicom传输语法