单个dcm文件含有多帧数据，如何拆分成多个dcm文件

最新推荐文章于 2024-08-12 16:55:46 发布
最新推荐文章于 2024-08-12 16:55:46 发布
阅读量1.7k 收藏 6
点赞数
分类专栏： dcmtk 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/laziji/article/details/123071673
版权

最近研究了题目中的问题：

bool Parse2Dcms(std::string infilePath, std::string output_path)
{
	 
	if (!this->filepath.length())
	{
		cout << "No No parent dcm filepath! error" << endl;
		return false;
	}
	DcmFileFormat fileformat;
	
	OFCondition oc = fileformat.loadFile(this->filepath.c_str());              //读取Dicom图像  
	if (!oc.good())     //判断Dicom文件是否读取成功  
	{
		std::cout << "file Load error" << std::endl;
		return false;
	}
	/*               dicom  meta infomation            */
	DcmMetaInfo* metainfo =  fileformat.getMetaInfo();

	OFString FMIG;
	metainfo->findAndGetOFString(DCM_FileMetaInformationGroupLength, FMIG);
	OFString strVersion;
	metainfo->findAndGetOFString(DCM_FileMetaInformationVersion, strVersion);
	OFString strMSSOPClassUID;
	metainfo->findAndGetOFString(DCM_MediaStorageSOPClassUID, strMSSOPClassUID);
	OFString strMMSOPIUID;
	metainfo->findAndGetOFString(DCM_MediaStorageSOPInstanceUID, strMMSOPIUID);
	OFString strTransferSyntax;
	metainfo->findAndGetOFString(DCM_TransferSyntaxUID, strTransferSyntax);
	OFString strImpUID;
	metainfo->findAndGetOFString(DCM_ImplementationClassUID, strImpUID);
	OFString strImpVerName;
	metainfo->findAndGetOFString(DCM_ImplementationVersionName, strImpVerName);
 
	/*      dicom dataset infomation */
	 
	DcmDataset *dataset = fileformat.getDataset();                             //得到Dicom的数据集  
	E_TransferSyntax xfer = dataset->getOriginalXfer();                        //得到传输语法 
	unsigned short img_bits(0);
	dataset->findAndGetUint16(DCM_SamplesPerPixel, img_bits);                    //单个像素占用多少byte 
	unsigned short bit_count(0);
	dataset->findAndGetUint16(DCM_BitsStored, bit_count);                        //获取像素的位数 bit 
	OFString isRGB;
	dataset->findAndGetOFString(DCM_PhotometricInterpretation, isRGB);           //DCM图片的图像模式  
	OFString strPrimaryAngle;
	dataset->findAndGetOFStringArray(DCM_PositionerPrimaryAngleIncrement, strPrimaryAngle);
	std::string tagValue = strPrimaryAngle.c_str();
	while (true)
	{
		size_t nPos = tagValue.find('\\');
		if (std::string::npos != nPos) {
			vAngles.emplace_back(atof(tagValue.substr(0, nPos).c_str()));
			tagValue = tagValue.substr(nPos + 1, tagValue.length());
		}
		else {
			vAngles.emplace_back(atof(tagValue.c_str()));
			break;
		}
	}  //end while
	 
	//0008
	OFString strSChSet;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0005), strSChSet);
	OFString strImageType;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0008), strImageType);
	OFString strUniqueID;
	dataset->findAndGetOFString(DcmTagKey(0x0008,0x0016), strUniqueID);
	OFString strSOPInstanceUID;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0018), strSOPInstanceUID);
	OFString strStudyDate;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0020), strStudyDate);
	OFString strSeriesDate;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0021), strSeriesDate);
	OFString strContentDate;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0023), strContentDate);
	OFString strAcquisitionDateTime;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x002a), strAcquisitionDateTime);
	OFString strStudyTime;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0030), strStudyTime);
	OFString strSeriesTime;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0031), strSeriesTime);
	OFString strContentTime;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0033), strContentTime);
	OFString strAccessionNumber;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0050), strAccessionNumber);
	OFString strModality;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0060), strModality);
	OFString strManufacturer;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0070), strManufacturer);
	OFString strInstitutionName;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0080), strInstitutionName);
	OFString strReferringPhysicianName;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x0090), strReferringPhysicianName);
	OFString strStationName;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1010), strStationName);

	OFString strStudyDescription;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1030), strStudyDescription);
	OFString strSeriesDescription;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x103e), strSeriesDescription);
	OFString strPerformingPhysicianName;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1050), strPerformingPhysicianName);
	OFString strManufacturerModelName;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1090), strManufacturerModelName);
	OFString strReferencedStudySequence;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1110), strReferencedStudySequence);
	OFString strReferencedSOPClassUID;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1150), strReferencedSOPClassUID);
	OFString strReferencedSOPInstanceUID;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1155), strReferencedSOPInstanceUID);
	OFString strReferencedPerformedProcedureStepSequence;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x1111), strReferencedPerformedProcedureStepSequence);
	OFString strStartTrim;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x2142), strStartTrim);
	OFString strStopTrim;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x2143), strStopTrim);
	OFString strPatientIdentityRemoved;
	dataset->findAndGetOFString(DcmTagKey(0x0012, 0x0062), strPatientIdentityRemoved);
	OFString strIrradiationEventUID;
	dataset->findAndGetOFString(DcmTagKey(0x0008, 0x3010), strIrradiationEventUID);
	//0010
	OFString strPatientName;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x0010), strPatientName);
	OFString strPatientID;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x0020), strPatientID);
	OFString strIssuerOfPatientID;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x0021), strIssuerOfPatientID);
	OFString strPatientBirthDate;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x0030), strPatientBirthDate);
	OFString strPatientSex;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x0040), strPatientSex);
	OFString strPatientSize;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x1020), strPatientSize);
	OFString strPatientWeight;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x1030), strPatientWeight);
	OFString strPatientComments;
	dataset->findAndGetOFString(DcmTagKey(0x0010, 0x4000), strPatientComments);
		 

	//0018
	OFString strImageAndFluoros;
	dataset->findAndGetOFString(DcmTagKey(0x0018,0x115e), strImageAndFluoros);
	OFString strImagePixelSpace;
	dataset->findAndGetOFStringArray(DcmTagKey(0x0018,0x1164), strImagePixelSpace);
	OFString strPositionMotion;
	dataset->findAndGetOFString(DcmTagKey(0x0018,0x1500), strPositionMotion);
	OFString strPositionPrimaryAngle;
	OFString strPositionSecondAngle;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1510), strPositionPrimaryAngle);
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1511), strPositionSecondAngle);
	OFString strContrastBolus;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x0010),strContrastBolus);
	OFString strPartExamed;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x0015), strPartExamed);
	OFString strCineRate;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x0040), strCineRate);
	OFString strKVP;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x0060), strKVP);
	OFString strDeviceSerialNumber;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1000), strDeviceSerialNumber);
	OFString strSoftwareVer;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1020), strSoftwareVer);
	OFString strProtocolName;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1030), strProtocolName);
	OFString strDistanceToDetector;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1110), strDistanceToDetector);
	OFString strDistanceToPatient;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1111), strDistanceToPatient);
	OFString strExposureTime;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1150), strExposureTime);
	OFString strXRayTube;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1151), strXRayTube);
	OFString strRadiation;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x1155), strRadiation);
	OFString strRadiationMode;
	dataset->findAndGetOFString(DcmTagKey(0x0018, 0x115a), strRadiationMode);
	
	//0020
	OFString strStudyInstanceUID;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x000d), strStudyInstanceUID);
	OFString strSeriesInstanceUID;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x000e), strSeriesInstanceUID);
	OFString strStudyID;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x0010), strStudyID);
	OFString strSeriesNumber;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x0011), strSeriesNumber);
	OFString strAcquisitionNumber;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x0012), strAcquisitionNumber);
	OFString strInstanceNumber;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x0013), strInstanceNumber);
	OFString strPatientOrientation;
	dataset->findAndGetOFStringArray(DcmTagKey(0x0020, 0x0020), strPatientOrientation);
	OFString strImagePositionPatient;
	dataset->findAndGetOFStringArray(DcmTagKey(0x0020, 0x0032), strImagePositionPatient);
	OFString strLaterality;
	dataset->findAndGetOFString(DcmTagKey(0x0020, 0x0060), strLaterality);
    //0028
	OFString strPixelSpacing;
	dataset->findAndGetOFStringArray(DcmTagKey(0x0028, 0x0030), strPixelSpacing);
	unsigned short bitsAllocated(0);
	dataset->findAndGetUint16(DcmTagKey(0x0028, 0x0100), bitsAllocated);
	unsigned short bitsStored(0);
	dataset->findAndGetUint16(DcmTagKey(0x0028, 0x0101), bitsStored);
	unsigned short hightBit(0);
	dataset->findAndGetUint16(DcmTagKey(0x0028, 0x0102), hightBit);
	unsigned short pixelRep(0);
	dataset->findAndGetUint16(DcmTagKey(0x0028, 0x0103), pixelRep);
	OFString strPixelRelationShip;
	dataset->findAndGetOFString(DcmTagKey(0x0028, 0x1040), strPixelRelationShip);
	OFString strWC;
	dataset->findAndGetOFString(DcmTagKey(0x0028, 0x1050), strWC);
	OFString strWW;
	dataset->findAndGetOFString(DcmTagKey(0x0028, 0x1051), strWW);
	OFString strFrameIncrementPointer;
	dataset->findAndGetOFString(DcmTagKey(0x0028, 0x0009), strFrameIncrementPointer);
 
	//0040
	OFString strPerformedProcedureStepStartDate;
	dataset->findAndGetOFString(DcmTagKey(0x0040, 0x0244), strPerformedProcedureStepStartDate);
	OFString strPerformedProcedureStepStartTime;
	dataset->findAndGetOFString(DcmTagKey(0x0040, 0x0245), strPerformedProcedureStepStartTime);
	OFString strPerformedProcedureStepID;
	dataset->findAndGetOFString(DcmTagKey(0x0040, 0x0253), strPerformedProcedureStepID);
	OFString strPerformedProcedureStepDescription;
	dataset->findAndGetOFString(DcmTagKey(0x0040, 0x0254), strPerformedProcedureStepDescription);
	  

	unsigned short m_width;                                                     //获取图像的窗宽高  
	unsigned short m_height;
	OFString framecount;
  //DCM_MinimumStoredValueMapped  
 // DCM_MaximumStoredValueMapped  
	dataset->findAndGetOFString(DCM_NumberOfFrames, framecount);				//DCM图片的帧数  
	unsigned short number = atoi(framecount.c_str());
	dataset->findAndGetUint16(DCM_Rows, m_height);
	dataset->findAndGetUint16(DCM_Columns, m_width);
	 
	int rows(m_height);
	int cols(m_width);
	const int length(rows*cols*2);
	Uint16* img_pixel = new Uint16[length];
	
	DicomImage* m_pDicomImage = new DicomImage("文件名", CIF_UsePartialAccessToPixelData, 0, 1);
	if (m_pDicomImage == NULL)
		throw "Open DICOM File failed";
	double min = 0, max = 0;
	m_pDicomImage->getMinMaxValues(min, max);
	unsigned long FrameCnt = m_pDicomImage->getFrameCount();
	unsigned long indexFirstFrame = m_pDicomImage->getFirstFrame();
	//
	DcmElement* element = NULL;
	OFCondition result = dataset->findAndGetElement(DCM_PixelData, element); //获取图像数据  

	if (result.bad() || element == NULL)
		return false;
	Uint16* image16_data = nullptr;
	result = element->getUint16Array(image16_data);      //获取>8位的图像数据  

	do
	{
		DcmFileFormat *fileFormat = new DcmFileFormat();
		DcmDataset *dataSet = fileFormat->getDataset();
		DcmMetaInfo* metaInfo = fileformat.getMetaInfo();
		memset(img_pixel, 0, length);

		
		//m_pDicomImage->getOutputData(img_pixel, length, Img_bitCount, indexFirstFrame);
		//const DiPixel* pixel = m_pDicomImage->getInterData();
		memcpy_s(img_pixel, length, image16_data + rows*cols*indexFirstFrame, length);
		
		//0x0002
		metaInfo->putAndInsertString(DCM_FileMetaInformationGroupLength, FMIG.c_str());
		metaInfo->putAndInsertString(DCM_FileMetaInformationVersion, strVersion.c_str());
		metaInfo->putAndInsertString(DCM_MediaStorageSOPClassUID, strMSSOPClassUID.c_str());
		 
		metaInfo->putAndInsertString(DCM_MediaStorageSOPInstanceUID, strMMSOPIUID.c_str());
		metaInfo->putAndInsertString(DCM_TransferSyntaxUID, strTransferSyntax.c_str());
		metaInfo->putAndInsertString(DCM_ImplementationClassUID, strImpUID.c_str());
		metaInfo->putAndInsertString(DCM_ImplementationVersionName, strImpVerName.c_str());
		metaInfo->computeGroupLengthAndPadding(EGL_withGL, EPD_noChange, META_HEADER_DEFAULT_TRANSFERSYNTAX, EET_UndefinedLength);
		//0x0008
		dataSet->putAndInsertString(DCM_SpecificCharacterSet, strSChSet.c_str());
		dataSet->putAndInsertString(DCM_ImageType, strImageType.c_str());
	 	dataSet->putAndInsertString(DCM_SOPClassUID, strUniqueID.c_str());
		dataSet->putAndInsertString(DCM_SOPInstanceUID, strSOPInstanceUID.c_str());
		dataSet->putAndInsertString(DCM_StudyDate, strStudyDate.c_str());
		dataSet->putAndInsertString(DCM_SeriesDate, strSeriesDate.c_str());
		dataSet->putAndInsertString(DCM_ContentDate, strContentDate.c_str());
		dataSet->putAndInsertString(DCM_AcquisitionDateTime, strAcquisitionDateTime.c_str());
		dataSet->putAndInsertString(DCM_StudyTime, strStudyTime.c_str());
		dataSet->putAndInsertString(DCM_SeriesTime, strSeriesTime.c_str());
		dataSet->putAndInsertString(DCM_ContentTime, strContentTime.c_str());
		dataSet->putAndInsertString(DCM_AccessionNumber, strAccessionNumber.c_str());
		dataSet->putAndInsertString(DCM_Modality, strModality.c_str());
		dataSet->putAndInsertString(DCM_Manufacturer, strManufacturer.c_str());
		dataSet->putAndInsertString(DCM_InstitutionName, strInstitutionName.c_str());
		dataSet->putAndInsertString(DCM_ReferringPhysicianName, strReferringPhysicianName.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1010), strStationName.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1030), strStudyDescription.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x103e), strSeriesDescription.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1050), strPerformingPhysicianName.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1090), strManufacturerModelName.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1110), strReferencedStudySequence.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1150), strReferencedSOPClassUID.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1155), strReferencedSOPInstanceUID.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x1111), strReferencedPerformedProcedureStepSequence.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x2142), strStartTrim.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x2143), strStopTrim.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0012, 0x0062), strPatientIdentityRemoved.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0008, 0x3010), strIrradiationEventUID.c_str());
		
		//0010
		 
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x0010), strPatientName.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x0020), strPatientID.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x0021), strIssuerOfPatientID.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x0030), strPatientBirthDate.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x0040), strPatientSex.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x1020), strPatientSize.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x1030), strPatientWeight.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0010, 0x4000), strPatientComments.c_str());

		//0x0020
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x000d), strStudyInstanceUID.c_str(),true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x000e), strSeriesInstanceUID.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0010), strStudyID.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0011), strSeriesNumber.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0012), strAcquisitionNumber.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0013), strInstanceNumber.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0020), strPatientOrientation.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0032), strImagePositionPatient.c_str(), true);
		dataSet->putAndInsertString(DcmTagKey(0x0020, 0x0060), strLaterality.c_str(), true);
		//0x0018
		dataSet->putAndInsertString(DCM_ImageAndFluoroscopyAreaDoseProduct, strImageAndFluoros.c_str(), true);
	 	dataSet->putAndInsertString(DCM_ImagerPixelSpacing, strImagePixelSpace.c_str(), true);
	 	dataSet->putAndInsertString(DCM_PositionerMotion, strPositionMotion.c_str(), true);
		dataSet->putAndInsertString(DCM_PositionerPrimaryAngle, strPositionPrimaryAngle.c_str(), true);
		dataSet->putAndInsertString(DCM_PositionerSecondaryAngle, strPositionSecondAngle.c_str(), true);
		dataSet->putAndInsertString(DCM_PositionerPrimaryAngleIncrement, std::to_string(vAngles[indexFirstFrame]).c_str());
	 
		dataSet->putAndInsertString(DCM_ContrastBolusAgent, strContrastBolus.c_str(), true);
		dataSet->putAndInsertString(DCM_BodyPartExamined, strPartExamed.c_str(), true);
		dataSet->putAndInsertString(DCM_CineRate, strCineRate.c_str(), true);
		dataSet->putAndInsertString(DCM_KVP, strKVP.c_str(), true);
		dataSet->putAndInsertString(DCM_DeviceSerialNumber, strDeviceSerialNumber.c_str(), true);
		dataSet->putAndInsertString(DCM_SoftwareVersions, strSoftwareVer.c_str(), true);
		dataSet->putAndInsertString(DCM_ProtocolName, strProtocolName.c_str(), true);
		dataSet->putAndInsertString(DCM_DistanceSourceToDetector, strDistanceToDetector.c_str(), true);
		dataSet->putAndInsertString(DCM_DistanceSourceToPatient, strDistanceToPatient.c_str(), true);
		dataSet->putAndInsertString(DCM_ExposureTime, strExposureTime.c_str(), true);
		dataSet->putAndInsertString(DCM_XRayTubeCurrent, strXRayTube.c_str(), true);
		dataSet->putAndInsertString(DCM_RadiationSetting, strRadiation.c_str(), true);
		dataSet->putAndInsertString(DCM_RadiationMode, strRadiationMode.c_str(), true);
	
		//0x0028
		dataSet->putAndInsertUint16Array(DCM_PixelData, OFreinterpret_cast(Uint16*, img_pixel), length);
		dataSet->putAndInsertString(DCM_NumberOfFrames, "1");
		dataSet->putAndInsertUint16(DCM_Rows, rows);
		dataSet->putAndInsertUint16(DCM_Columns, cols);
		dataSet->putAndInsertUint16(DCM_SamplesPerPixel, img_bits);
		dataSet->putAndInsertUint16(DCM_BitsStored, bit_count);
		dataSet->putAndInsertUint16(DCM_BitsAllocated,bitsAllocated);
		dataSet->putAndInsertString(DcmTagKey(0x0028, 0x0009), strFrameIncrementPointer.c_str());
		dataSet->putAndInsertString(DCM_PixelSpacing, strPixelSpacing.c_str(), true);
		dataSet->putAndInsertUint16(DCM_HighBit, hightBit);
		dataSet->putAndInsertUint16(DCM_PixelRepresentation, pixelRep);
		dataSet->putAndInsertString(DCM_PixelIntensityRelationship, strPixelRelationShip.c_str(), true);
		dataSet->putAndInsertString(DCM_WindowCenter, strWC.c_str(), true);
		dataSet->putAndInsertString(DCM_WindowWidth, strWW.c_str(), true);
		//0x0040
		 
		dataSet->putAndInsertString(DcmTagKey(0x0040, 0x0244), strPerformedProcedureStepStartDate.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0040, 0x0245), strPerformedProcedureStepStartTime.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0040, 0x0253), strPerformedProcedureStepID.c_str());
		dataSet->putAndInsertString(DcmTagKey(0x0040, 0x0254), strPerformedProcedureStepDescription.c_str());

		dataSet->putAndInsertString(DCM_PhotometricInterpretation, isRGB.c_str());
		
	 
		
		//保存文件
		char t[10];
		string filename;
		sprintf_s(t, "%d", indexFirstFrame);
		filename = t;
		filename = output_path + filename + ".dcm";
		OFCondition status = fileFormat->saveFile(filename.c_str(), xfer);
		if (status.bad())
			cout << "Error: cannot write DICOM file (" << status.text() << ")" << endl;

		delete fileFormat;
		indexFirstFrame = m_pDicomImage->getFirstFrame();

		//this->ShowImage(filename);
	}while (m_pDicomImage->processNextFrames());
	 
	delete[] img_pixel;
	return true;
}

其实DicomImage可以不用。这个怎么和DcmDataSet或DcmFileFormat联系起来，欢迎大家讨论。

peanut_wu
关注
  • 0
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
专栏目录
dicom多帧转换_如何提取/查看多帧DICOM文件的所有帧？
weixin_39622901的博客
12-18 251
I know I can view a dicom file using following code:import dicomfrom dicom.contrib.pydicom_PIL import show_PILf = "CT-MONO2-8-abdo.dcm"ds = dicom.read_file(f, force=True)show_PIL(ds)However, how can I...
排序算法相关的时间复杂度和进程稳定性(下)
yanhaijunyan的博客
10-10 446
待继续分析
多帧Dcm图片的读取与拆分存储
m_buddy的博客
11-28 9159
1. 概述 在项目中为了保留对之前代码的兼容和加上新功能，需要将多帧的Dcm图像转换为单帧的Dcm图像，近期对其进行了实践。本篇文章中将实现的代码予以分享，不妥之处请见谅。多帧图像 2. 实现代码 之前有两篇博客分别讲了DCMTK在vs上配置过程和读取Dcm图的实现，感兴趣的读者可以看看之前的文章。 //*****************************************
python dicom图像分割_如何使用python将pydicom文件（图像）分成两半？
weixin_40007016的博客
11-21 682
I have a lot of images (pydicom files). I would like to divide in half. From 1 image, I would like 2 images: part left and part right.Input: 1000x1000Output: 500x1000 (width x height).Currently, I can...
多帧Dcm图像
11-28
里面包含30张彩色的多帧Dcm图像，这里仅供用作开发测试用途
ITK笔记——读取多帧DICOM图像
Andy的专栏
10-26 2758
上篇笔记介绍了单个DICOM图像的读取，接下来看如何读取DICOM序列以及使用。
DCM4CHE处理多帧图像
lemin_zhao的博客
02-09 2819
概要 dicom医疗影像图像有单帧和多帧的分别，一般来说普通CT CR图像都是单帧的(也有多帧的情况)，一般DXA和US图形是多帧的，今天我们就主要介绍一下多帧图像的处理 准备工作 准备一个多帧图像的dicom 查看这个dicom的压缩协议，如果是带压缩协议的dicom需要解压缩，具体的参见我的博客DCM4CHE解压缩的DICOM文件 大概思路 单帧dicom是把一幅...
ITK操作合并dcm图像为VTK等三维数据
10-20
将多个DCM图像合并成一个VTK（ Visualization Toolkit）数据，可以用于构建三维模型，便于后续的可视化和分析。 首先，我们需要理解DCM图像文件。DCM文件不仅包含图像数据，还包含元数据，如患者信息、扫描参数等。...
ITK 实现多张图像转成单个nii.gz或mha文件案例
12-17
这个案例展示了ITK在医学图像处理中的强大功能，能够方便地处理图像序列，将其整合为单个文件，这对于数据分析、可视化和后续的处理任务非常有用。通过这样的方法，用户可以更高效地管理和操作大量的医学影像数据。
读DCM数据，对DCM图按照序列号进行分类
11-20
DCM文件由多个数据元素（Data Elements）组成，每个元素都有一个特定的标签（Tag），包含了元数据（如患者信息、设备信息）和图像数据。元数据通常以ASCII编码的文本形式存在，而图像数据则可能以压缩或未压缩的二...
DICOM.zip_dicom_dicom matlab_dicom 文件
07-15
在实际应用中，如果你有一个包含多个DICOM文件的压缩包"DICOM.zip"，你需要先解压，然后才能逐一处理这些文件。MATLAB的`unzip`函数可以完成解压操作： ```matlab unzip('DICOM.zip'); ``` 解压后，你可以遍历所有...
DCM反激输出滤波电容的发热诱因分析
08-14
此时，可以通过并联多个电容来分散ESR损耗，从而降低单个电容的发热。 针对这些问题，一些常用的解决方案包括：一是增加输出电容的数量，通过并联降低每个电容承受的纹波电流；二是采用二级LC滤波，进一步减少纹波...
Python解析多帧dicom数据
技术随笔
02-20 2792
目录概述DICOM多帧数据存储pydicom应用 概述 pydicom是一个常用python DICOM parser。但是，没有提供解析多帧图的示例。本文结合相关函数和DICOM知识做一个简单说明。 DICOM多帧数据存储 DICOM标准中关于多帧数据存储的最重要一部分说明是PS3.5 Annex A.4 A.4 Transfer Syntaxes For Encapsulation of En...
DCMTK:带有多张图片的dcm文件转成多张QImage与多张图片添加到一个dcm文件
syz9011的博客
11-18 3278
随着对DICOM标准慢慢的理解，发现其实带有单张图片的dcm文件与QImage文件的互转并不是一件很难的事情，一般情况下，结果出不来的原因都是对DICOM标准不理解。然而往往一点点理解不够，可能浪费大量的事情，造成很大困扰。带有多张图的dcm文件转QImage这个最简单，在原有但张图的dcm文件转QImage文章（http://blog.csdn.net/syz9011/article/detail
DICOM医学图像处理：DICOM存储操作之“多幅BMP图像数据存入DCM文件”
热门推荐
只要踏出一步，路就在前方——zssure
12-24 3万+
背景： 本专栏“DICOM医学图像处理”受众较窄，起初只想作为自己学习积累和工作经验的简单整理。前几天无聊浏览了一下，发现阅读量两极化严重，主要集中在“关于BMP（JPG）与DCM格式转换”和“DICOM 通讯协议”，尤其是许久前的第一篇博文DCMTK开源库的学习笔记1：将DCM文件保存成BMP文件或数据流（即数组）。因此在2014年底前打算写几篇关于DCM格式转换的文章，此次主要聚焦“如何将BMP、JPG等常规图像保存成DCM文件”，以DCMTK库为基础，给出简单的实例。
vtk读取dcm文件
Jeremy写字的地方
11-21 802
import vtk reader = vtk.vtkDICOMImageReader() reader.SetFileName('...\\.dcm') reader.Update() viewer = vtk.vtkImageViewer2() viewer.SetInputConnection(reader.GetOutputPort()) iren = vtk.vtkRenderWindowInteractor() viewer.SetupInteractor(iren) viewer.Re
VTK系列14_VTK读取DCM医学图像(CT图)并显示
qq_32809093的博客
09-16 2669
实例14：VTK读取DCM医学图像并显示 #include "vtkAutoInit.h" VTK_MODULE_INIT(vtkRenderingOpenGL2); VTK_MODULE_INIT(vtkInteractionStyle); #include <vtkDICOMImageReader.h>//DCM医学文...
vtk读取.dcm文件（医学影像文件）
xi__q的专栏
12-23 2081
VTK读取dcm文件（医学影像文件）的一个核心的知识点就是，窗宽窗位。 知识点： 窗宽/窗位的概念 窗宽是图像显示的灰度范围。一般显示器的灰度范围为256级，而医学图像的灰度范围则远远大于该范围，因此通过显示器显示时不能显示所有灰度级，需要使用窗宽来定义欲显示的灰度范围。灰度值高于该范围的最大值时：白影。低于：黑影。比如，窗宽200，窗位100，则可视范围：0-200；窗位500，可视范围：400-600 QDialog *dlg = new QDialog(); dlg->resize(600
C++笔记5•内存管理•
最新发布
qq_64446190的博客
08-12 1044
内存管理
dcmtk3.6.3 dcm2pnm.exe把多帧dicom图像转换成单张png图像
07-12
要使用DCMTK 3.6.3的dcm2pnm.exe将多帧DICOM图像转换为单个PNG图像，您可以按照以下步骤操作： 1. 首先，确保已经正确安装DCMTK 3.6.3库。您可以从DCMTK的官方网站（https://dicom.offis.de/dcmtk.php.en）下载并按照说明进行安装。 2. 打开命令行提示符或终端窗口，并导航到包含DICOM图像文件的目录。 3. 运行以下命令以将多帧DICOM图像转换为单个PNG图像： ``` dcm2pnm +Fpng input.dcm output.png ``` 其中，`input.dcm`是您要转换的多帧DICOM图像文件的路径，`output.png`是生成的单个PNG图像文件的路径。 请确保在命令行中正确指定了dcm2pnm.exe的路径，或者将其添加到系统环境变量中，以便可以直接在任何目录下运行该命令。 4. 执行命令后，dcm2pnm将会解析输入的多帧DICOM图像，并将其转换为单个PNG图像。生成的PNG图像文件将保存在指定的输出路径中。 请注意，dcm2pnm是DCMTK提供的一个实用工具，可用于将DICOM图像转换为其他常见的图像格式。根据您的需求，您可以使用其他选项和参数来进一步控制转换过程。您可以通过运行 `dcm2pnm --help` 命令获取更多用法和选项的详细信息。 此外，如果您希望在自己的程序中实现DICOM图像转换，您可以使用DCMTK库提供的API来读取和处理DICOM图像，并将其保存为PNG格式。具体的实现方式将涉及更多的编程细节，可以参考DCMTK的文档和示例代码以获取更多信息。
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值