转自:https://blog.csdn.net/zssureqh/article/details/8804736
dcm文件是医学领域DICOM3.0标准所对应的主要的文件格式。前两篇学习笔记中,学习了读取dcm文件的相关信息,如信息头MetaInformation元素、像素数据元素,只停留在了读取和显示dcm文件的阶段,随着学习的深入,自己开发的平台的功能逐渐增加,“修改dcm文件的相关数据元”在所难免。下面将简单的介绍一下对dcm文件的修改,主要的参考资料是:dcmtk开源库的官方文件,另外也很感谢diqiucun666在CSDN的博文http://blog.csdn.net/diqiucun666/article/details/2866908,里面详细介绍了dcmtk的主要构成,是dcmtk开源库很好的中文是说明手册。
参照博文中的:
dcmtk开源库中对dcm文件定义的最顶层的类是:DcmFileFormat,每次需要利用DcmFileFormat的对象来进行dcm文件的导入(loadFile)和导出(saveFile)。dcm文件其本质与常见的DIB文件格式相似,可总体分为“文件头”和“文件体”两部分。dcmtk开源库分别为“文件头”和“文件体“封装了相应的类:DcmMetaInfo和DcmDataset。两个类都继承自DcmItem类,其含有主要的成员变量为:elementList,一个存储相应dcm文件数据元的双向列表。(注意:DcmFileFormat类继承自DcmSequenceOfItems类,其含有的数据成员变量时itemList,即含有DcmMetaInfo类和DcmDataset类的基类的双向列表——这说明了dcm文件的数据元可以进行嵌套存储),DcmFileFormat变量的成员参见下图:
此处是一个空的DcmFileFormat类的对象mDcmFileFormat,其itemList双向链表中只有两个item成员,firstNode就是DcmMetaInfo类的实例,lastNode就是DcmDataset类的实例。
那么找到修改itemList数据成员的方法,就能够修改相应的dcm数据元,达到自由修改dcm文件的目的,但是在修改的时候需要注意,各个数据元的元素类型的匹配。
查看dcitem.h文件,可以发现dcmtk给我们定义了几种操作itemList双向链表的方法,如下图:
下面给出修改dcm文件中的患者姓名和dcm像素数据的代码:
DcmFileFormat mDcmImage;mDcmImage.loadFile("c:\\test.dcm");short* pPixelData=static_cast<short*>(mDcmImage.getPixelData());int size=mDcmImage.getWidth()*mDcmImage.getHeight();for(int i=0;i<size;++i)if(pPixelData[i]<800)pPixelData[i]=0;elsepPixelData[i]=2000;UINT16 *pData=new UINT16[size];memcpy(pData,pPixelData,size*sizeof(UINT16));DcmDataset *dataset=mDcmImage.getDataset();DcmMetaInfo *metainfo=mDcmImage.getMetaInfo();metainfo->remove(DCM_TransferSyntaxUID);delete dataset->remove(DCM_PatientName);dataset->putAndInsertString(DCM_PatientName,"TEST-THU");delete dataset->remove(DCM_PixelData);dataset->putAndInsertUint16Array(DCM_PixelData,pData,size);mDcmImage.saveFile("c:\\self.dcm",EXS_LittleEndianExplicit);
其中getPixelData()函数是对:
DcmDataset *dataset = this->getDataset();// decompress data set if compresseddataset->chooseRepresentation(EXS_LittleEndianExplicit, NULL);DcmElement* element=NULL;dataset->findAndGetElement(DCM_PixelData,element);unsigned char* pImage=NULL;element->getUint8Array(pImage);return pImage;
代码的封装。利用像素指针dcm文件的像素数据进行二值化。
然后我们利用DcmItem类中定义的putAndInsertString()修改DCM_PatinetName标签,即患者的姓名。
对别结果如下: