CN112562829A - 一种自适应生成dicom影像默认窗宽窗位的方法 - Google Patents
一种自适应生成dicom影像默认窗宽窗位的方法 Download PDFInfo
- Publication number
- CN112562829A CN112562829A CN202110188159.0A CN202110188159A CN112562829A CN 112562829 A CN112562829 A CN 112562829A CN 202110188159 A CN202110188159 A CN 202110188159A CN 112562829 A CN112562829 A CN 112562829A
- Authority
- CN
- China
- Prior art keywords
- image
- optimal
- gray
- histogram
- level
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F9/00—Arrangements for program control, e.g. control units
- G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
- G06F9/44—Arrangements for executing specific programs
- G06F9/451—Execution arrangements for user interfaces
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Computer Interaction (AREA)
- Radiology & Medical Imaging (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种自适应生成DICOM影像默认窗宽窗位的方法。该方法包括获取医学图像数据;计算图像灰度的最大值和最小值;根据图像灰度的最大值和最小值构建图像的灰度直方图;计算灰度直方图中的所有波峰区域;计算组织区域对应最佳直方图波峰区域;根据最佳直方图波峰区域来计算图像最佳显示灰度值区间,并根据图像最佳显示灰度值区间计算最佳窗宽和最佳窗位;将计算出的最佳窗宽和最佳窗位写入DICOM图像文件中,作为默认窗宽窗位。本发明不需要复杂的计算,即可高效地得到DICOM影像最佳窗宽和最佳窗位,并写入到DICOM文件中,作为默认窗宽窗位,使得图像显示明显更加清晰且符合阅图习惯。
Description
技术领域
本发明涉及医学图像的生成技术领域,具体涉及一种自适应生成DICOM影像默认窗宽窗位的方法。
背景技术
DICOM 文件是指按照DICOM (Digital Imaging Communications in Medicine)标准而存储的医学文件。DICOM 文件主要组成部分就是数据集合。这不仅包括医学图像,还包括许多和医学图像有关的信息,如病人姓名、图像大小等。窗宽窗位就是其中一个非常重要的信息。
窗宽是CT图像上显示的CT值范围,在此CT值范围内的组织和病变均以不同的模拟灰度显示。而CT值高于此范围的组织和病变,无论高出程度有多少,均以白影显示,不再有灰度差异;反之,低于此范围的组织结构,不论低的程度有多少,均以黑影显示,也无灰度差别。增大窗宽,则图像所示CT值范围加大,显示具有不同密度的组织结构增多,但各结构之间的灰度差别减少。减小窗宽,则显示的组织结构减少,然而各结构之间的灰度差别增加。
窗位是窗的中心位置,同样的窗宽,由于窗位不同,其所包括CT值范围的CT值也有差异。
DICOM图像浏览软件(包括医院的PACS及其他后处理软件)一般都是通过读取DICOM的默认窗宽窗位或者基于默认窗宽窗位来优化窗宽窗位参数来显示图像,因此DICOM图像的默认窗宽窗位很大程度上决定了图像是否能默认清晰显示,也就是决定了用户是否需要进行手动调节窗宽窗位。因此,给生成的DICOM图像设置合适的窗宽窗位尤其重要。。
发明内容
本发明的目的是针对现有技术存在的不足,提供一种自适应生成DICOM影像默认窗宽窗位的方法。
为实现上述目的,本发明提供了一种自适应生成DICOM影像默认窗宽窗位的方法,包括:
步骤1:获取医学图像数据;
步骤4:计算所述灰度直方图中的所有波峰区域;
步骤5:计算组织区域对应最佳直方图波峰区域;
步骤6:根据步骤5中求得的最佳直方图波峰区域来计算图像最佳显示灰度值区间,并根据图像最佳显示灰度值区间计算最佳窗宽和最佳窗位;
步骤7:将步骤6中计算出的最佳窗宽和最佳窗位写入DICOM图像文件中,作为默认窗宽窗位。
进一步的,所述医学图像数据包括医疗设备采集获取的原始图像数据和对原始图像数据处理后产生的图像数据。
进一步的,所述步骤4中计算灰度直方图中的所有波峰区域具体包括:
计算数组M1中两个级数之差,其中,,其中,和分别
为数组M1中第i个和i+1个级数,i为大于1的自然数,如时,其中,为设定阈值,则
级数即判断为相邻级数,并从数组M1中去除,以此得到新的数组M2,数组M2中灰度级数
即组成灰度直方图中所有的波峰区域。
进一步的,所述步骤5中计算组织区域对应最佳直方图波峰区域具体包括:
进一步的,所述最佳窗宽和最佳窗位的计算方式如下:
进一步的,所述b的计算方式如下:
有益效果:本发明通过对图像进行直方图分析,估计组织所属灰度区间,不需要复杂的计算,即可高效地得到DICOM影像最佳窗宽和最佳窗位,并写入到DICOM文件中,作为默认窗宽窗位,使得图像显示明显更加清晰且符合阅图习惯。
附图说明
图1是本发明实施例的自适应生成DICOM影像默认窗宽窗位的方法的流程示意图;
图2为窗宽窗位自适应前冠脉CT图像显示效果;
图3为图2中冠脉CT图像对应的直方图及窗宽窗位的示意图;
图4为图2中冠脉CT图像经过本发明窗宽窗位自适应后的显示效果图;
图5为图4中冠脉CT图像对应的直方图及窗宽窗位的示意图;
图6为经过本发明设置的DICOM默认窗宽窗位查看结果。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,本实施例在以本发明技术方案为前提下进行实施,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围。
如图1所示,本发明实施例提供了一种自适应生成DICOM影像默认窗宽窗位的方法,包括:
步骤1:获取医学图像数据。具体来说,上述医学图像的原始图像数据,可以由医疗图像采集设备(如CT等)采集获得的原始图像数据,也可以是在医学图像处理过程中,对原始图像数据处理后产生的图像数据,或其他图像来源。
步骤3:根据图像灰度的最大值和最小值构建图像的灰度直方图。直
方图的组数可为某一默认值,一般选择128或256。而直方图的组距由图像灰度的
最大值和最小值及组数共同决定。构建图像的灰度直方图,需要遍历整
个图像区域,为减少计算量,加快计算速度,可通过等距采样遍历图像。通过实验,通过一定
程度的等距采样绘制直方图,并不会明显影响最终计算得到的最佳窗宽窗位数值,比如采
样距离为3。组距的计算方式如下:
步骤4:计算灰度直方图中的所有波峰区域。可采用阈值分割方法提取连续的直方图灰度区间,具体如下:
根据整个灰度直方图统计的图像像素点总数来确定一个阈值,将频数小于阈值的灰度级予以剔除。作为优选实施例,阈值=直方图统计的图像像素点总数*0.0015。
频数小于阈值的灰度级为不具有统计意义的信息,排除在窗宽之外,频数在阈值以上
的灰度级为具有统计意义的信息。
计算数组M1中两个级数之差,其中,,其中,和分别
为数组M1中第i个和i+1个级数,i为大于1的自然数,如时,其中,为设定阈值,则
级数即判断为相邻级数,并从数组M1中去除,以此得到新的数组M2,数组M2中灰度级数
即组成灰度直方图中所有的波峰区域。直方图的组数为128时,可采用5,直方图的
组数为256时,可采用10。
步骤5:计算组织区域对应最佳直方图波峰区域。计算组织区域对应最佳直方图波峰区域具体包括:
步骤6:根据步骤5中求得的最佳直方图波峰区域来计算图像最佳显示灰度值区间,并根据图像最佳显示灰度值区间计算最佳窗宽和最佳窗位。具体的,最佳窗宽和最佳窗位的计算方式如下:
步骤7:将步骤6中计算出的最佳窗宽和最佳窗位写入DICOM图像文件中,作为默认窗宽窗位。如图6所示,窗宽对应在DICOM文件中的Tag ID是(0028,1051),窗位对应在DICOM文件中的Tag ID是(0028,1050),通过将最佳窗宽(WW)和最佳窗位(WL)写入对应的Tag中,随着DICOM文件一起保存,即生成了以自适应生成DICOM影像为默认窗宽窗位的DICOM图像。新生成的DICOM图像对应的窗宽窗位tag可通过第三方DICOM图像浏览器查看。
如图2所示,在自适应窗宽窗位之前的冠脉CTA图像显示比较模糊,其灰度直方图及窗宽窗位如图3所示;如图4所示,在经过本发明自适应窗宽窗位算法的处理后,冠脉CTA图像显示明显更加清晰且符合阅图习惯,其灰度直方图及窗宽窗位如图5所示。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,其它未具体描述的部分,属于现有技术或公知常识。在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
2.根据权利要求1所述的自适应生成DICOM影像默认窗宽窗位的方法,其特征在于,所述医学图像数据包括医疗设备采集获取的原始图像数据和对原始图像数据处理后产生的图像数据。
6.根据权利要求5所述的自适应生成DICOM影像默认窗宽窗位的方法,其特征在于,所述步骤5中计算组织区域对应最佳直方图波峰区域具体包括:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188159.0A CN112562829B (zh) | 2021-02-19 | 2021-02-19 | 一种自适应生成dicom影像默认窗宽窗位的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110188159.0A CN112562829B (zh) | 2021-02-19 | 2021-02-19 | 一种自适应生成dicom影像默认窗宽窗位的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112562829A true CN112562829A (zh) | 2021-03-26 |
CN112562829B CN112562829B (zh) | 2021-05-14 |
Family
ID=75035949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110188159.0A Active CN112562829B (zh) | 2021-02-19 | 2021-02-19 | 一种自适应生成dicom影像默认窗宽窗位的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112562829B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114663653A (zh) * | 2022-05-16 | 2022-06-24 | 南京佗道医疗科技有限公司 | 一种医学影像感兴趣区域的窗位窗宽计算方法 |
CN115568873A (zh) * | 2022-12-07 | 2023-01-06 | 佗道医疗科技有限公司 | 一种ct影像自动调整方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105862A (zh) * | 2007-08-02 | 2008-01-16 | 宁波大学 | 一种医学图像窗口参数的自适应调整方法 |
CN104376568A (zh) * | 2014-11-28 | 2015-02-25 | 成都影泰科技有限公司 | 一种基于格式的dicom医学图像处理方法 |
CN107481208A (zh) * | 2017-09-13 | 2017-12-15 | 许人灿 | Dicom医学影像动态线性调窗方法 |
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
CN110223261A (zh) * | 2019-05-20 | 2019-09-10 | 上海联影医疗科技有限公司 | 医学图像处理方法和***、存储介质及计算机设备 |
-
2021
- 2021-02-19 CN CN202110188159.0A patent/CN112562829B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105862A (zh) * | 2007-08-02 | 2008-01-16 | 宁波大学 | 一种医学图像窗口参数的自适应调整方法 |
CN104376568A (zh) * | 2014-11-28 | 2015-02-25 | 成都影泰科技有限公司 | 一种基于格式的dicom医学图像处理方法 |
CN107481208A (zh) * | 2017-09-13 | 2017-12-15 | 许人灿 | Dicom医学影像动态线性调窗方法 |
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
CN110223261A (zh) * | 2019-05-20 | 2019-09-10 | 上海联影医疗科技有限公司 | 医学图像处理方法和***、存储介质及计算机设备 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114663653A (zh) * | 2022-05-16 | 2022-06-24 | 南京佗道医疗科技有限公司 | 一种医学影像感兴趣区域的窗位窗宽计算方法 |
CN114663653B (zh) * | 2022-05-16 | 2022-08-12 | 南京佗道医疗科技有限公司 | 一种医学影像感兴趣区域的窗位窗宽计算方法 |
CN115568873A (zh) * | 2022-12-07 | 2023-01-06 | 佗道医疗科技有限公司 | 一种ct影像自动调整方法 |
CN115568873B (zh) * | 2022-12-07 | 2023-08-22 | 佗道医疗科技有限公司 | 一种ct影像自动调整方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112562829B (zh) | 2021-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0780791B1 (en) | Automatic tone scale adjustment using image activity measures | |
US8861886B2 (en) | Enhanced visualization for medical images | |
CN112562829B (zh) | 一种自适应生成dicom影像默认窗宽窗位的方法 | |
US20030161516A1 (en) | Method for enhancing the contrast of an image | |
CN114677391B (zh) | 一种脊椎图像分割方法 | |
JPH03206572A (ja) | 階調変換自動化装置 | |
CN113989407B (zh) | Ct影像中肢体部位识别模型训练方法及*** | |
Liu et al. | Automatic lung segmentation based on image decomposition and wavelet transform | |
CN116152505A (zh) | 基于x光数据的骨骼目标识别分割方法 | |
US20030161548A1 (en) | Noise reduction method | |
US20030161518A1 (en) | Gradation processing method | |
CN114972067A (zh) | 一种x光小牙片图像增强方法 | |
CN114693682B (zh) | 一种基于图像处理的脊椎特征识别方法 | |
US20030161519A1 (en) | Multiscale gradation processing method | |
EP1345170B1 (en) | Noise reduction method. | |
CN101897592B (zh) | X射线影像中金属异物的提取方法和装置 | |
US20030161497A1 (en) | Method for enhancing the contrast of an image | |
CN112052854B (zh) | 一种实现自适应对比度增强的医疗图像可逆信息隐藏方法 | |
CN116993764B (zh) | 一种胃部ct智能分割提取方法 | |
CN116402816B (zh) | 一种体检ct影像数据的管理方法及*** | |
CN116128774B (zh) | 一种胃部螺旋ct数据增强处理方法 | |
US20030161517A1 (en) | Method of normalising a digital signal representation of an image | |
EP1347413A1 (en) | Method for enhancing the contrast of an image. | |
CN117252825A (zh) | 基于口腔全景片图像的龋齿识别方法和装置 | |
Kalhor et al. | Assessment of histogram-based medical image contrast enhancement techniques; an implementation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |