在多帧图像中插入修改过的像素数据

Inserting modified pixel data in a multiframe dicom image?

本文关键字:像素 像素数 数据 修改 插入 图像      更新时间:2023-10-16

首先我要说明的是,我使用的是Dcmtk库。

我已经设法学习如何修改单帧dicom图像的像素数据。现在,我尝试在多帧图像中做同样的事情。我可以提取所有必要的信息,甚至可以为每帧单独提取像素数据,并可以修改它们。但是,当我必须插入修改后的像素数据时,问题就出现了。在单帧的情况下,我使用DcmDataset中的方法:

putAndInsertUint8Array()

但是对于多帧图像,我看不到任何这样的选项。我在dcelement:

中使用此方法获取每帧的像素数据
getUncompressedFrame()

,我只需要把帧索引,以获得相应的像素数据。但是在插入的时候,我找不到任何这样的选项。我的编程代码如下:

int main()
{
   MdfDatasetManager file;
   if(EC_Normal==file.loadFile("test.dcm",ERM_autoDetect,EXS_Unknown))
   {
      DcmDataset *dataset = file.getDataset();
      E_TransferSyntax xfer= dataset->getOriginalXfer();
      bool OriginallyCompressed=false;
      if(xfer!=0 && xfer !=1 && xfer!=2 && xfer!=3)
      {
         OriginallyCompressed=true;
         DJDecoderRegistration::registerCodecs();
         if(EC_Normal==dataset->chooseRepresentation(EXS_LittleEndianExplicit, NULL))
         {
            if(dataset->canWriteXfer(EXS_LittleEndianExplicit))
            {
               cout<<"Originally it's a compressed image, but now decompressed!n";
            }
         }
      }
      DcmElement* element=NULL;
      Uint16 rows = 0;
      Uint16 cols = 0;
      Uint16 samplePerPixel = 0;
      Uint16 planarConfiguration = 0;
      int index=0;
      // I've fixed these values but later I will change them to dinaymic and make it work as per user's wish.
      int ymin=50;//minimum rows
      int ymax=500;//maximum rows
      int xmin=100;//Minimum columns
      int xmax=600;//Maximum columns
      if(EC_Normal==dataset->findAndGetUint16(DCM_Rows, rows))
      {
         if(EC_Normal==dataset->findAndGetUint16(DCM_Columns, cols))
         {
            if(EC_Normal==dataset->findAndGetUint16(DCM_SamplesPerPixel,samplePerPixel))
            {   
               if(EC_Normal==dataset->findAndGetUint16(DCM_PlanarConfiguration,planarConfiguration))
               {
                  if(EC_Normal==dataset->findAndGetElement(DCM_PixelData,element))
                  {
                     Uint32 startFragment=0;
                     Uint32 sizeF=0;
                     element->getUncompressedFrameSize(dataset,sizeF);
                     long int numOfFrames=0;
                     dataset->findAndGetLongInt(DCM_NumberOfFrames,numOfFrames);
                     for(int i=0;i<int(numOfFrames);i++)
                     {
                        Uint8 * buffer = new Uint8[int(sizeF)];
                        OFString decompressedColorModel=NULL;
                        DcmFileCache * cache=NULL;
                        if(EC_Normal==element->getUncompressedFrame(dataset,i,startFragment,buffer,sizeF,decompressedColorModel,cache))
                        {
                           Uint8 * newBuffer = new Uint8[int(sizeF)];
                           if(buffer != NULL)
                           {
                              for(unsigned long y = 0; y < rows; y++)
                              {
                                 for(unsigned long x = 0; x < cols; x++)
                                 {
                                    if(planarConfiguration==0)
                                    {
                                       if(x>xmin && x<xmax && y>ymin && y<ymax)
                                       {
                                          index=(x + y +  y*(cols-1))*samplePerPixel;
                                          newBuffer[index]  = 0;
                                          newBuffer[index + 1]  = 0;
                                          newBuffer[index +2]  = 0;
                                       }
                                       else
                                       {
                                          index=(x + y +  y*(cols-1))*samplePerPixel;
                                          newBuffer[index]  = buffer[index];
                                          newBuffer[index + 1]  = buffer[index + 1];
                                          newBuffer[index + 2]  = buffer[index + 2];
                                       }
                                    }
                                 }
                              }
                           }
                           delete newBuffer;
                        }
                        delete buffer;
                     }
                  }
               }
            }
         }
      }
   }
   return 0;
}

如果我只是设法找到一种方法来插入修改像素数据的每一帧,这个程序将是完整的。请建议我该做什么。或者请告诉我,如果你知道,多帧dicom图像中所有帧的整个像素数据是如何存储在一起的。然后我可以从所有帧中取出所有的像素数据并修改它们然后尝试将所有修改过的像素数据一起插入

case#1,未压缩像素数据:

///. get PixelData element in DCM dataset
pDcmDataSet->findAndGetPixelData(...);
///. get pixels in PixelData element
pDcmPixelDataSet->findAndGetOW(...);

您将获得所有帧的完整像素数据。

case#2,压缩像素数据:

///. get PixelData element in DCM dataset
pDcmDataSet->findAndGetPixelData(...);
///. get PixelSequence in PixelData element
pPixelData->getEncapsulatedRepresentation(...)
///. get PixelItem in PixelSequence
pDcmPixelSequence->getItem(...);
///. get frame in Pixel Item
pPixelItem->getUint8Arrary(...);

你将得到一帧压缩图像