DICOM医学图像格式转换的C++实现 DIOCM转 BMP

最新推荐文章于 2021-04-25 15:54:10 发布
最新推荐文章于 2021-04-25 15:54:10 发布
阅读量2.3k 收藏 19
点赞数 2
分类专栏： VTK  ITK
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/sinat_33514750/article/details/104923531
版权 
//zoukeke 1670386358@qq.com 2020年3月17日
#include "stdafx.h"
#include <opencv2/opencv.hpp>
#include <windows.h>
using namespace cv;
struct TagValue
{
	unsigned short tag1;
	unsigned short tag2;
};

// 原始象素类型
enum E_SourcePixelType
{ 
	ESourcePixelType_U16,		// USHORT
	ESourcePixelType_I16,		// SHORT
};

int _tmain(int argc, _TCHAR* argv[])
{
	bool isVR=true;
	bool isLitteEndian=true;
	int file_length=0;
	char VR[3];

	unsigned int pixDataLen=0;
	unsigned int pixDataOffset=0;
	unsigned short channle=0;
	unsigned short rows=0;
	unsigned short cols=0;
	unsigned short dataLen=0;  
	unsigned short validLen=0;

	E_SourcePixelType pixelType;

	int windowsWidth=0;
	int windowsCenter=0;

	bool ZeroIsBlack=true;
	float RescaleSlope =0.06;
	float RescaleIntercept=0;

	FILE *fp;
	fp=fopen("11.dcm","rb");
	if(fp==NULL)
	{
		printf("can not open file!");
		return 0;
	} 
	fseek(fp,0,SEEK_END);
	file_length=ftell(fp);
	fseek(fp,0,SEEK_SET);

	fseek(fp,128,SEEK_SET);

	char headchar[5];
	memset(headchar,0,5);
	int read_num = fread(headchar,1,4,fp);
	if(read_num!=4)
	{
		fclose(fp);
		return 0;
	}

	if(strcmp(headchar,"DICM"))   
	{
		fclose(fp);
		printf("File is not DICM");
		return 0;
	}

	while(ftell(fp)+6<file_length)
	{
		TagValue tag;
		unsigned int len;
		memset(VR,0,3);
		fread(&tag,sizeof(TagValue),1,fp);
		int index=ftell(fp);

		if(tag.tag1==0x02)
		{
			fread(VR,1,2,fp);
			if(!strcmp(VR,"OB")||!strcmp(VR,"OW")||!strcmp(VR,"SQ"))
			{
				fseek(fp,2,SEEK_CUR);
				fread(&len,sizeof(unsigned int),1,fp);
			}
			else
			{
				unsigned short l;
				int ss = fread(&l,sizeof(unsigned short),1,fp);
				int a=ftell(fp);
				len =(unsigned int)l ;
			}
		}
		else if(tag.tag1==0xfffe)
		{
			if(tag.tag2==0xe000||tag.tag2==0xe00d||tag.tag2==0xe0dd)
			{
				fread(&len,sizeof(unsigned int),1,fp);
			}
		}
		else if(isVR==true)
		{
			fread(VR,1,2,fp);
			int a= ftell(fp);
			if(!strcmp(VR,"OB")||!strcmp(VR,"OW")||!strcmp(VR,"SQ"))
			{
				fseek(fp,2,SEEK_CUR);
				fread(&len,sizeof(unsigned int),1,fp);
			}
			else
			{
				unsigned short l;
				l=sizeof(unsigned short);
				fread(&l,sizeof(unsigned short),1,fp);
				a= ftell(fp);
				len =(unsigned int)l ;
			}
		}
		else if(isVR==false)
		{
			fread(&len,sizeof(unsigned int),1,fp);
		}

		if(tag.tag1==0x02&&tag.tag2==0x10)
		{
			char msg[124];
			memset(msg,0,124);
			fread(msg,1,len,fp);

			if(!strcmp(msg,"1.2.840.10008.1.2.1"))
			{
				isLitteEndian=true;
				isVR=true;
			}
			else if(!strcmp(msg,"1.2.840.10008.1.2.2"))
			{
				isLitteEndian=false;
				isVR=true;
			}
			else if(!strcmp(msg,"1.2.840.10008.1.2"))
			{
				isLitteEndian=true;
				isVR=false;
			}
		}
		else if(tag.tag1 ==0x28 && tag.tag2==0x103)
		{
			unsigned short m;
			fread(&m,sizeof(unsigned short),1,fp);
			if(m==0)
			{
				pixelType = ESourcePixelType_U16;
			}
			else if(m ==1)
			{
				pixelType = ESourcePixelType_U16;
			}
		}
		else if(tag.tag1==0x7fe0&&tag.tag2==0x10)
		{
			pixDataLen=len;
			pixDataOffset=ftell(fp);
			fseek(fp,len,SEEK_CUR);

		}
		else if(tag.tag1==0x28&&tag.tag2==0x10)
		{
			fread(&rows,sizeof(unsigned short),1,fp);
		}
		else if(tag.tag1==0x28&&tag.tag2==0x11)
		{
			fread(&cols,sizeof(unsigned short),1,fp);
		}
		else if(tag.tag1==0x28&&tag.tag2==0x02)
		{
			fread(&channle,sizeof(unsigned short),1,fp);
		}
		else if(tag.tag1==0x28&&tag.tag2==0x101)
		{
			fread(&validLen,sizeof(unsigned short),1,fp);
		}
		else if(tag.tag1==0x28&&tag.tag2==0x100)
		{
			fread(&dataLen,sizeof(unsigned short),1,fp);
		}
		else if(tag.tag1==0x28&&tag.tag2==0x1050)
		{
			char msg[11];
			memset(msg,0,11);
			fread(msg,1,len,fp);
			windowsCenter=atoi(msg);

		}
		else if(tag.tag1==0x28,tag.tag2==0x1051)
		{
			//fseek(fp,len,SEEK_CUR);
			char msg[40];
			memset(msg,0,40);
			fread(msg,1,len,fp);
			windowsWidth=atoi(msg);
		}
		else if(tag.tag1==0x0028&&tag.tag2==0x0004)
		{
			char msg[40];
			memset(msg,0,40);
			fread(msg,1,len,fp);
			if(!strcmp(msg,"MONOCHROME1 "))
			{
				ZeroIsBlack=false;
			}
			else if(!strcmp(msg,"MONOCHROME2 "))
			{
				ZeroIsBlack=true;
			}
		}
		else if(tag.tag1==0x0028&&tag.tag2==0x1052)
		{
			char msg[40];
			memset(msg,0,40);
			fread(msg,1,len,fp);
			RescaleIntercept=atof(msg);
		}
		else if(tag.tag1==0x0028&&tag.tag2==0x1053)
		{
			char msg[40];
			memset(msg,0,40);
			fread(msg,1,len,fp);
			RescaleSlope =atof(msg);
		}
		else
		{
			char msg[1024];
			memset(msg,0,1024);
			fread(msg,1,len,fp);


		}
	}

	fseek(fp,pixDataOffset,SEEK_SET);
	if(windowsCenter==0&&windowsWidth==0)
	{
		windowsWidth = 1 << validLen;
		windowsCenter = windowsWidth / 2;
	}
	
	int min_value,max_value;
	min_value=windowsCenter-windowsWidth/2.0+0.5;
	max_value=windowsCenter+windowsWidth/2.0+0.5;
	double pers = 255.0/(max_value-min_value);  

	Mat src;    
	Mat src2;

	int nPixel= 0;

	double fCtA = 0; 
	double fCtB = 0;
	fCtA = (double)256 /windowsWidth;
	fCtB = 128 - 256 * (double)windowsCenter / windowsWidth;
	if (fCtB < 0)
	{
		fCtB = 0;
	}
	if (fCtB > 255)
	{
		fCtB = 255;
	}

	if(channle==1)
	{
		src.create((int)rows,(int)cols,CV_8UC1);
	
		src2.create((int)rows,(int)cols, CV_16SC1);
		for (int i = 0;i<rows;i++)
		{
			for(int j=0;j<cols;j++)
			{
				unsigned short   gray=0;
				short gray2 = 0;
				unsigned char pix[2];
				fread(pix,1,2,fp);
				if(pixelType == ESourcePixelType_U16)
				{
					if(validLen<=8)
					{
						if(isLitteEndian)
						{
							gray=pix[0];
						}
						else   
						{
							gray=pix[1];
						}

					}
					else 
					{
						long temp = 0;
						if(isLitteEndian)
						{
							gray=*(unsigned short*)pix;

							if(gray > 32767)
							{
								gray = 32767;
							}
							temp =gray*RescaleSlope+RescaleIntercept;
							temp = temp * fCtA + fCtB;
						}
						else
						{
							gray=pix[1]+pix[0]*256;
						}

						int nValue = (int)temp;
						if(nValue>0xff)
						{
							nValue=0xff;
						}
						else if(nValue<0)
						{
							nValue=0;

						}
						nPixel = nValue;
					}
				}
				else if(pixelType ==ESourcePixelType_I16)
				{
					if(validLen<=8)
					{
						if(isLitteEndian)   
						{
							gray2=pix[0];
						}
						else   
						{
							gray2=pix[1];
						}

					}
					else 
					{
						long temp = 0;
						if(isLitteEndian)
						{
							gray2=*(short*)pix;

							if(gray2 > 32767)
							{
								gray2 = 32767;
							}
							if(gray2 < -32767)
							{
								gray2 = -32767;
							}
							temp =gray2*RescaleSlope+RescaleIntercept;
							temp = temp * fCtA + fCtB;
						}
						else
						{
							gray2=pix[1]+pix[0]*256;
						}

						int nValue = (int)temp;
						if(nValue>0xff)
						{
							nValue=0xff;
						}
						else if(nValue<0)
						{
							nValue=0;

						}
						nPixel = nValue;
					}
				}		
				if(!ZeroIsBlack)
				{
					nPixel=255-nPixel;
				}
				src.at<uchar>(i,j)=nPixel;
			}
			std::cout<<std::endl;
		}
	}
	else if(channle==3)
	{
		src.create((int)rows,(int)cols,CV_8UC3);
		for (int i = 0;i<rows;i++)
		{
			for(int j=0;j<cols;j++)
			{
				unsigned char pix[3];
				fread(pix,1,3,fp);
				src.at<Vec3b>(i,j)[0]=pix[2];
				src.at<Vec3b>(i,j)[1]=pix[1];
				src.at<Vec3b>(i,j)[2]=pix[0];
			}
		}
	}
	imwrite("D:\\src.bmp", src);  
	fclose(fp);
	imshow("src",src);
	waitKey(0);
	return 0;      
}
**

## *参考资料*

**

在这里插入图片描述
在这里插入图片描述

sinat_33514750
关注
专栏目录
DICOM与BMP相互转换 c++代码 VS2010
04-16
DICOM与BMP相互转换 c++代码 VS2010
C++将DICOM文件中图像转换为JPEG文件
Ericohe的博客
07-23 1013
void SaveJpg(vtkSmartPointer<vtkImageData> imageData) { int width = imageData->GetDimensions()[0]; int height = imageData->GetDimensions()[1]; QImage image(width, height, QImage::Format_RGB32); short *colorsPtr = reinterpret_c.
DICOM医学图像格式转换的C++实现
03-23
多数医疗影像无法实现不同格式转换，这是很好的参考。
java dcm转bmp源码_Dcm-files-into-BMP-or-Jpg 用c++写的一个将dicom格式的图片转化为JPG或者BMP的算法，对医学工作者有一定帮助 Special Effec...
weixin_35792236的博客
03-06 152
文件名称: Dcm-files-into-BMP-or-Jpg下载 收藏√ [5 4 3 2 1]开发工具: Visual C++文件大小: 362 KB上传时间: 2016-01-16下载次数: 0提 供 者: lubo详细说明：用c++写的一个将dicom格式的图片转化为JPG或者BMP的算法，对医学工作者有一定帮助-With c++ write a dicom format ...
dcm文件转jpg或者bmp
热门推荐
cqltbe131421的博客
03-02 1万+
医学上用到的是dcm格式的图片，怎么转其他格式呢， 这个就要用到DCMTK了 我们在安装DCMTK的安装bin目录下可以看到 这些可以用的工具 其中dcmj2png.exe就具有将dcm格式的图片转为bmp.jpg,png等格式图片的功能 我们打开cmd，进入到这个目录。或者你把这个文件拷贝到其他目录，然后用cmd进入也是 可以的。 我是拷贝到z盘的根目录下的，cd到z： 输入d
DICOM医学图像格式转换的C++实现 DIOCM转 BMP、读取DICOM图像（转载）
Hello，小高同学
04-25 1294
代码： //zoukeke 1670386358@qq.com 2020年3月17日 #include "stdafx.h" #include <opencv2/opencv.hpp> #include <windows.h> using namespace cv; struct TagValue { unsigned short tag1; unsigned short tag2; }; // 原始象素类型 enum E_SourcePixelType { ESour
Dcm转换成bmp或QImage
Summer的博客
10-24 2565
1、Dcm转换成bmpif (QFile::exists(strFilePath))//strFilePath 为转换成的bmp图像及路径 { if (!QFile::remove(strFilePath)) { DxError("Old Thumbnail Delete Failed."); return Dx
DIcOM医学图像格式转换与BMP实现详解
这证明了将DICOM图像转换为常见的Windows平台格式（如BMP）对于医学图像处理和后续操作的重要性。 4. 结束语：DICOM 3.0作为医学图像的标准存储和通信格式，对于医疗机构而言，能够确保不同设备间图像的互换性和...
DICOM医学图像到BMP格式转换的C++实现方法
"本文介绍了如何使用C++实现DICOM医学图像格式转换，主要针对DICOM标准在医学成像领域的应用和重要性进行了阐述，并探讨了如何将DICOM图像转换为常见的BMP格式。" DICOM（Digital Imaging and Communications in ...
C# 读取dcm医学图像，格式转换为bmp图像的demo
10-28
本示例是一个C#程序，用于读取DICOM图像并将其转换为更常见的BMP图像格式，这在显示和存储医学图像时非常有用。以下将详细讲解实现这个功能所需的知识点： 1. **C#编程基础**：C#是微软开发的一种面向对象的编程...
DIOCM 图像查看代码
12-26
DICOM 查询 DICOM 查询DICOM 查询DICOM 查询DICOM 查询DICOM 查询
C++，dcmtk，gdcm实现dicom图像解压缩，信息读取，支持jpeg，jpeg2000，jpeg-l
11-21
C++，dcmtk，gdcm实现dicom图像解压缩，信息读取，支持jpeg，jpeg2000，jpeg-ls等多种压缩格式，支持有损解压和无损解压
dcmtk-3.6.3
12-06
dicom开发库，c++版本，可以进行医学影像的存储，传输等。
DICOM图片转换BMP的C++代码
03-31
借助于VC，使用C++开发一个软件，可以将DICOM图片转换为BMP图片，并进行简单的图像处理
Dicom图像读取与显示
04-08
使用C++读取Dicom文件，获取Dicom文件中与图像显示相关的信息，使用OpenCV进行图像显示，代码对于Dico初学者会有很大帮助
绝对可用的dicom医学图像格式转换为BMP格式（仅含exe程序）
09-17
采用MFC和delphi混合编程形成的可执行文件，...可以支持目前大部分的dicom标准的图像的转换，但必须在exe同级目录下放置TempImage文件夹才可以。如果想修改此目录，可以给我发邮件：quanxin19841112@163.com，谢谢。
DICOM医学图像到BMP格式转换软件设计与实现分析
" dicom图像格式与bmp通用格式转换软件实现分析" 在医学成像领域，DICOM（Digital Imaging and Communications in Medicine）是一种广泛使用的标准格式，用于存储、传输和打印医学图像，如CT（计算机断层扫描）和...
dcm文件如何转化为jpg或者bmp文件
wyl_0728的博客
03-21 7190
最近看到群里很多小伙伴比较关心dcm文件怎么转jpg或者是bmp文件，其实在dcm4chee中是提供了很多的方法，但是我们真正在用的时候其实重点还是在dcm文件本身，如果是经过压缩的文件，你用传统方法就没法去结项，如果是彩色的图片，你就蒙了，那么就这些问题给大家提供几种解析方式： 一、传统转化-----使用dcm2jpg 导入jar包依赖： <dependency> ...
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值