使用 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 Instance UID）
4. 将修改后的 DICOM 数据集保存到新文件中

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