Вставить измененные данные пикселей в многокадровое изображение dicom?

Прежде всего я должен упомянуть, что для этой цели я использую библиотеку Dcmtk.

Мне уже удалось узнать, как я могу изменить пиксельные данные в одном кадре изображения. Теперь я пытаюсь сделать то же самое в случае многокадровых изображений. Я могу извлечь всю необходимую информацию и даже может извлечь данные пикселей отдельно для каждого кадра и может изменить их. Но проблема возникает, когда мне нужно вставить измененные данные пикселей. В случае одного кадра я использую метод в 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;
}

Если мне удастся найти способ вставить измененные данные пикселей для каждого кадра, эта программа будет завершена. Пожалуйста, предложите мне, что я должен сделать. Или, пожалуйста, скажите мне, если вы знаете, как все пиксельные данные всех кадров в многокадровом двухцветном изображении хранятся вместе. Тогда, может быть, я смогу собрать все данные о пикселях вместе из всех кадров и изменить их, а затем попытаться вставить все данные об измененных пикселях вместе.

1

Решение

случай № 1, несжатые данные пикселей:

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

Вы получите целые пиксельные данные всех кадров в одном куске.

случай № 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(...);

Вы получите один кадр сжатого изображения.

0

Другие решения

Других решений пока нет …

По вопросам рекламы [email protected]