c++ - 使用 dcmtk 覆盖 dicom 文件中的图像/像素数据

Question

我使用 dcmtk 读取 dicom 文件并将图像提取为 .tiff 格式。进行一些图像处理后，我有一张图像，我想将其保存在源 dicom 文件中。这就是用我的新图像覆盖旧图像/像素数据，同时保留其余数据(uid、患者姓名等)一样。

我用下面的代码读取dicom

OFCondition status = src_fileformat.loadFile(src_path);

    if (status.good())
    {
        Sint32 instanceNumber = 0;
        if (src_fileformat.getDataset()->findAndGetSint32(DCM_InstanceNumber, instanceNumber).good())
        {
            cout << "instance Number N: " << instanceNumber << endl;
            sprintf(instanceNum, "%d", instanceNumber);
            printf("%s\n", instanceNum);
        }

        else
            cerr << "Error: cannot access Instance Number!" << endl;
    }
    else
        cerr << "Error: cannot read DICOM file (" << status.text() << ")" << endl;

    src_dcm = new DicomImage(src_path);
    if (src_dcm != NULL)
    {
        if (src_dcm->getStatus() == EIS_Normal)
        {
            if (src_dcm->isMonochrome())
            {
                src_dcm->setMinMaxWindow();
                Uint8 *pixelData = (Uint8 *)(src_dcm->getOutputData(16 /* bits */));
                if (pixelData != NULL)
                {
                    src_dcm->writeBMP("source.tiff", 24);  /* do something useful with the pixel data */
                }
            }
        }
        else
            cerr << "Error: cannot load DICOM image (" << DicomImage::getString(src_dcm->getStatus()) << ")" << endl;
    }

图像处理后，我有一张图像要覆盖到这个源 dicom 文件上。我查看了 image2dcm，但无法获得正确的语法/方法。任何人帮助我.. :)

编辑-1

Image2Dcm i2d; 
I2DOutputPlug *outPlug = new I2DOutputPlugSC(); 
 I2DImgSource *inputPlug = new I2DJpegSource(); 
 E_TransferSyntax writeXfer; 
 inputPlug->setImageFile(jpgFile); 
DcmDataset *dataset = NULL; 
OFCondition result = i2d.convert(inputPlug, outPlug, dataset, writeXfer); 

   // Saving output DICOM image 
  if (result.good()) 
      { 
        dataset->putAndInsertString(DCM_PhotometricInterpretation,"RGB"); 

            DcmFileFormat dcmff(dataset); 
    result = dcmff.saveFile(dcmFile, writeXfer); 
       }

我试过上面显示的语法，但不能完全理解

这是处理后的图像(上图)

这是我要覆盖的原始 dicom 图像。伙计们，有什么想法或帮助吗？？

Answer 1

基本方法应该是:

1. 从文件中加载 DICOM 数据集
2. 替换数据集中的像素数据
3. 修改各种其他元素值(例如 SOP 实例 UID)
4. 将修改后的 DICOM 数据集保存到新文件

在未压缩图像的情况下，第二步可以按照与第三步相同的方式执行，即通过对数据集适当调用 putAndInsertXXX() 方法。当然，像素数据属性的元素值应该是正确的 DICOM 格式。参见 DICOM standard第 3 部分和第 5 部分了解详细信息。