2

最初に、この目的のために Dcmtk ライブラリを使用していることに言及する必要があります。

単一フレームの dicom 画像のピクセル データを変更する方法については、既に学習しました。今、マルチフレーム画像の場合も同じことをしようとしています。必要な情報をすべて抽出でき、フレームごとに個別にピクセル データを抽出して変更することもできます。しかし、変更されたピクセル データを挿入する必要があるときに問題が発生します。単一フレームの場合、DcmDataset のメソッドを使用します。

putAndInsertUint8Array()

しかし、マルチフレーム画像にはそのようなオプションはありません。DcmElement で次のメソッドを使用して、各フレームのピクセル データを取得します。

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 画像内のすべてのフレームの全ピクセル データがどのように一緒に格納されるかを知っている場合は、教えてください。次に、すべてのフレームからピクセル データ全体を取得して変更し、変更されたピクセル データ全体を一緒に挿入することができます。

4

1 に答える 1

0

ケース#1、非圧縮ピクセル データ:

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

すべてのフレームの全ピクセル データを 1 つにまとめます。

ケース#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(...);

圧縮された画像の 1 フレームが得られます。

于 2013-11-27T19:42:07.087 に答える