将图像的rgb数据转成DICOM医学图像格式

dcmtk官方文档：https://support.dcmtk.org/docs/
dcmtk最新源码下载：https://www.dcmtk.org/en/dcmtk/dcmtk-software-development/
dcmtk旧版本源码下载：https://dicom.offis.de/download/dcmtk/

用DCMTK库实现将图像转成dcm格式

dcmtk库的编译这里就不叙述了，网上有很多这方面的详细内容。直接上代码：

#include <iostream>
#include "opencv2/opencv.hpp"
#include "dcmtk/config/osconfig.h"
#include "dcmtk/dcmdata/dctk.h"
#include "dcmtk/dcmimgle/dcmimage.h"
#include "dcmtk/dcmdata/libi2d/i2djpgs.h"
#include "dcmtk/dcmjpeg/djencode.h"
#include "dcmtk/dcmjpeg/djrplol.h"

using namespace std;
bool ToDicom()
{
	DcmFileFormat dicomfile;
	DcmDataset *dataset = dicomfile.getDataset();

	Uint32 all_len = 0;
	E_TransferSyntax ts = EXS_LittleEndianExplicit;

	cv::Mat mt = cv::imread("test.jpg");
	cv::cvtColor(mt, mt, cv::COLOR_BGR2RGB);

	/*	添加患者信息	*/
	dataset->putAndInsertUint16(DCM_AccessionNumber, 0);
	dataset->putAndInsertString(DCM_PatientName, "李与白", true);
	dataset->putAndInsertString(DCM_PatientID, "201210409217");
	dataset->putAndInsertString(DCM_PatientBirthDate, "19900612");
	dataset->putAndInsertString(DCM_PatientSex, "男");
	dataset->putAndInsertString(DCM_PatientAge, "32");

	dataset->putAndInsertString(DCM_InstitutionName, "四川大学华西医院");
	dataset->putAndInsertString(DCM_InstitutionalDepartmentName, "InstitutionalDepartmentName");

	/*	添加Study信息	*/
	dataset->putAndInsertString(DCM_StudyDate, "StudyDate");
	dataset->putAndInsertString(DCM_StudyTime, "StudyTime");
	dataset->putAndInsertString(DCM_StudyDescription, "StudyDescription");
	char uid[100];
	dcmGenerateUniqueIdentifier(uid, SITE_STUDY_UID_ROOT);
	dataset->putAndInsertString(DCM_StudyInstanceUID, uid);
	dataset->putAndInsertString(DCM_StudyID, "SF11240921450001");

	/*	添加Series信息	*/
	dataset->putAndInsertString(DCM_SeriesDate, "SeriesDate");
	dataset->putAndInsertString(DCM_SeriesTime, "SeriesTime");
	dataset->putAndInsertString(DCM_SeriesDescription, "SeriesDescription");
	memset(uid, 0, sizeof(char) * 100);
	dcmGenerateUniqueIdentifier(uid, SITE_SERIES_UID_ROOT);
	dataset->putAndInsertString(DCM_SeriesInstanceUID, uid);

	/*	添加Image信息	*/
	dataset->putAndInsertString(DCM_PlanarConfiguration, "0");
	dataset->putAndInsertString(DCM_ImageType, "ORIGINAL\\PRIMARY\\AXIAL");
	dataset->putAndInsertString(DCM_ContentDate, "ContentDate");
	dataset->putAndInsertString(DCM_ContentTime, "ContentTime");
	dataset->putAndInsertString(DCM_InstanceCreationDate, "20231125");//年月日
	dataset->putAndInsertString(DCM_InstanceCreationTime, "090750");//09:07:50
	dataset->putAndInsertString(DCM_AcquisitionDate, "20231126");
	dataset->putAndInsertString(DCM_AcquisitionTime, "090751");
	dataset->putAndInsertString(DCM_InstanceNumber, "1");
	dataset->putAndInsertString(DCM_PixelSpacing, "0.3\\0.3");
	dataset->putAndInsertString(DCM_WindowCenter, "600"); //源图像的窗位
	dataset->putAndInsertString(DCM_WindowWidth, "800");//源图像的窗宽
	dataset->putAndInsertString(DCM_RescaleIntercept, "0");//灰度偏移
	dataset->putAndInsertString(DCM_RescaleSlope, "1");

	dataset->putAndInsertString(DCM_NumberOfFrames, "1");
	dataset->putAndInsertUint16(DCM_SamplesPerPixel, 3);
	dataset->putAndInsertUint16(DCM_Rows, mt.rows);
	dataset->putAndInsertUint16(DCM_Columns, mt.cols);
	dataset->putAndInsertUint16(DCM_BitsAllocated, 8);
	dataset->putAndInsertUint16(DCM_BitsStored, 8);
	dataset->putAndInsertUint16(DCM_HighBit, 7);
	dataset->putAndInsertUint8Array(DCM_PixelData, (Uint8*)mt.data, mt.channels()*mt.total());
	dataset->putAndInsertOFStringArray(DCM_PhotometricInterpretation, "RGB");
	if (dataset->canWriteXfer(ts))
	{
		OFCondition status = dicomfile.saveFile("test.dcm", ts);
		if (status.good())
		{
			std::cout << "save ok" << std::endl;
			return true;
		}
		else
		{
			std::cerr << "Error saving DICOM file: " << status.text() << std::endl;
			return false;
		}
	}
	else
	{
		cout << "can not write" << endl;
		return false;
	}
}

int main()
{
	ToDicom();
	system("pause");
	return 0;
}

代码中用到了opencv图像处理库，主要是方便读取图像数据，而且可以通过opencv将其支持的任何格式转成dcm格式。

如果转换后的dcm图像打不开则可以用dcmdump.exe进行调试：

1. 进入到dcmdump.exe所在文件夹
2. 在地址栏输入cmd回车
3. 在弹出的提示符中输入:

dcmdump.exe d:/test.dcm

然后回车就能看到文件信息

对未编码的dcm文件进行编码或改变其编码方式

int main()
{
	DJEncoderRegistration::registerCodecs(); // 注册JPEG编解码器
	DcmFileFormat fileformat;
	if (fileformat.loadFile("test.dcm").good())
	{
		DcmDataset *dataset = fileformat.getDataset();
		DcmItem *metaInfo = fileformat.getMetaInfo();

		// 创建数据集的无损JPEG版本
		if (dataset->chooseRepresentation(EXS_JPEGProcess14SV1, NULL).good() && dataset->canWriteXfer(EXS_JPEGProcess14SV1))
		{
			// store in lossless JPEG format
			OFCondition status = fileformat.saveFile("test_jpeg.dcm", EXS_JPEGProcess14SV1);
			if (status.good())
				std::cout << "encode jpeg ok" << std::endl;
			else
				std::cout << "encode jpeg err" << std::endl;
		}
		else
			std::cout << "chooseRepresentation err" << std::endl;
	}
	DJEncoderRegistration::cleanup();
	system("pause");
	return 0;
}

直接将rgb数据转DICOM并Jpeg编码

#include <iostream>
#include "opencv2/opencv.hpp"
#include "dcmtk/config/osconfig.h"
#include "dcmtk/dcmdata/dctk.h"
#include "dcmtk/dcmimgle/dcmimage.h"
#include "dcmtk/dcmdata/libi2d/i2djpgs.h"
#include "dcmtk/dcmjpeg/djencode.h"
#include "dcmtk/dcmjpeg/djrplol.h"

using namespace std;
void JPGToDicom()
{
	DJEncoderRegistration::registerCodecs();

	DcmFileFormat dicomfile;
	DcmDataset *dataset = dicomfile.getDataset();
	E_TransferSyntax ts = EXS_JPEGProcess14SV1;

	cv::Mat mt = cv::imread("src.jpg");
	cv::cvtColor(mt, mt, cv::COLOR_BGR2RGB);

	dataset->putAndInsertUint16(DCM_AccessionNumber, 0);
	dataset->putAndInsertString(DCM_PatientName, "xxx", true);
	dataset->putAndInsertString(DCM_PatientID, "201210409217");
	dataset->putAndInsertString(DCM_PatientBirthDate, "19900612");
	dataset->putAndInsertString(DCM_PatientSex, "男");
	dataset->putAndInsertString(DCM_PatientAge, "32");
	dataset->putAndInsertString(DCM_InstitutionName, "四川大学华西医院");
	dataset->putAndInsertString(DCM_InstitutionalDepartmentName, "InstitutionalDepartmentName");
	dataset->putAndInsertString(DCM_StudyDate, "StudyDate");
	dataset->putAndInsertString(DCM_StudyTime, "StudyTime");
	dataset->putAndInsertString(DCM_StudyDescription, "StudyDescription");
	char uid[100];
	dcmGenerateUniqueIdentifier(uid, SITE_STUDY_UID_ROOT);
	dataset->putAndInsertString(DCM_StudyInstanceUID, uid);
	dataset->putAndInsertString(DCM_StudyID, "SF11240921450001");
	dataset->putAndInsertString(DCM_SeriesDate, "SeriesDate");
	dataset->putAndInsertString(DCM_SeriesTime, "SeriesTime");
	dataset->putAndInsertString(DCM_SeriesDescription, "SeriesDescription");
	memset(uid, 0, sizeof(char) * 100);
	dcmGenerateUniqueIdentifier(uid, SITE_SERIES_UID_ROOT);
	dataset->putAndInsertString(DCM_SeriesInstanceUID, uid);
	dataset->putAndInsertString(DCM_PlanarConfiguration, "0");
	dataset->putAndInsertString(DCM_ImageType, "ORIGINAL\\PRIMARY\\AXIAL");
	dataset->putAndInsertString(DCM_ContentDate, "ContentDate");
	dataset->putAndInsertString(DCM_ContentTime, "ContentTime");
	dataset->putAndInsertString(DCM_InstanceCreationDate, "20231125");//年月日
	dataset->putAndInsertString(DCM_InstanceCreationTime, "090750");//09:07:58
	dataset->putAndInsertString(DCM_AcquisitionDate, "20231126");
	dataset->putAndInsertString(DCM_AcquisitionTime, "090751");
	dataset->putAndInsertString(DCM_InstanceNumber, "1");
	dataset->putAndInsertString(DCM_PixelSpacing, "0.3\\0.3");
	dataset->putAndInsertString(DCM_WindowCenter, "600"); //源图像的窗位
	dataset->putAndInsertString(DCM_WindowWidth, "800");//源图像的窗宽
	dataset->putAndInsertString(DCM_RescaleIntercept, "0");//灰度偏移
	dataset->putAndInsertString(DCM_RescaleSlope, "1");
	dataset->putAndInsertString(DCM_NumberOfFrames, "1");
	dataset->putAndInsertUint16(DCM_SamplesPerPixel, 3);
	dataset->putAndInsertUint16(DCM_Rows, mt.rows);
	dataset->putAndInsertUint16(DCM_Columns, mt.cols);
	dataset->putAndInsertUint16(DCM_BitsAllocated, 8);
	dataset->putAndInsertUint16(DCM_BitsStored, 8);
	dataset->putAndInsertUint16(DCM_HighBit, 7);
	dataset->putAndInsertUint8Array(DCM_PixelData, (Uint8*)mt.data, mt.channels()*mt.total());
	dataset->putAndInsertOFStringArray(DCM_PhotometricInterpretation, "RGB");
	if (dataset->chooseRepresentation(ts, NULL).good() && dataset->canWriteXfer(ts))
	{
		OFCondition status = dicomfile.saveFile("test_jpg.dcm", ts);
		if (status.good())
			std::cout << "encode jpeg ok" << std::endl;
	}
	DJEncoderRegistration::cleanup();
}

int main()
{
	ToDicom();
	system("pause");
	return 0;
}

dcm转Tiff

TIFF文件需要提前安装libtiff，并且DCMTK编译时加上tiff选项。
png,jpeg格式同理

#include <iostream>
#include "dcmtk/config/osconfig.h"
#include "dcmtk/dcmdata/dctk.h"
#include "dcmtk/dcmimgle/dcmimage.h"
#include "dcmtk/dcmdata/libi2d/i2djpgs.h"
#include "dcmtk/dcmjpeg/djencode.h"
#include "dcmtk/dcmjpeg/djrplol.h"
#include "dcmtk/dcmimage/dipitiff.h"

using namespace std;

int dcmToTiff(DcmFileFormat* fileformat)
{
	DicomImage* pDcmImage = new DicomImage(fileformat, fileformat->getDataset()->getOriginalXfer());
	if (pDcmImage && pDcmImage->getStatus() == EIS_Normal)
		return pDcmImage->writePluginFormat(new DiTIFFPlugin(), "test.tiff");
	else
		return 0;
}

int main()
{
	DcmFileFormat fileformat;
	if (fileformat.loadFile("test.dcm").good())
	{
		cout << "write to tiff return: " << dcmToTiff(&fileformat) << endl;
	}
	else
		cerr << "load file err!" << endl;

	system("pause");
	return 0;
}