CN107481208A - Dicom医学影像动态线性调窗方法 - Google Patents
Dicom医学影像动态线性调窗方法 Download PDFInfo
- Publication number
- CN107481208A CN107481208A CN201710820652.3A CN201710820652A CN107481208A CN 107481208 A CN107481208 A CN 107481208A CN 201710820652 A CN201710820652 A CN 201710820652A CN 107481208 A CN107481208 A CN 107481208A
- Authority
- CN
- China
- Prior art keywords
- mrow
- window
- image
- mfrac
- mtd
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 230000006870 function Effects 0.000 claims abstract description 22
- 238000013507 mapping Methods 0.000 claims abstract description 14
- 230000000694 effects Effects 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 2
- 210000001015 abdomen Anatomy 0.000 description 4
- 239000003814 drug Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供一种DICOM医学影像动态线性调窗方法,包括以下步骤:步骤一,DICOM图像数据读取,根据DICOM文件的标准格式读入图像的数据元素,数据元素包括标签信息、类型、长度、数据信息,并保存于内存之中;步骤二,设置DICOM图像窗口的窗宽和窗位,标记DICOM文件中的窗宽、窗位,通过改进的灰度映射窗函数G(V)用于动态调窗。本发明的优点在于采用改进的灰度映射窗函数G(V),映射灰度图像,实现图像显示更加细致、清晰,并保证动态调窗图像的显示的实时性。
Description
技术领域
本发明涉及医学图像处理,特别涉及一种基于应用的DICOM医学影像动态线性调窗方法。
背景技术
DICOM(Digital Imaging and Communications in Medicine)即医学数字成像和通信,是医学图像和相关信息的国际标准(ISO 12052)。它定义了质量能满足临床需要的可用于数据交换的医学图像格式。
DICOM被广泛应用于放射医疗,心血管成像以及放射诊疗诊断设备(X射线,CT,核磁共振,超声等),并且在眼科和牙科等其它医学领域得到越来越深入广泛的应用。在数以万计的在用医学成像设备中,DICOM是部署最为广泛的医疗信息标准之一。当前大约有百亿级符合DICOM标准的医学图像用于临床使用。
自从1985年DICOM标准第一版发布以来,DICOM给放射学实践带来了革命性的改变,X光胶片被全数字化的工作流程所代替。就像Internet成为信息传播应用的全新平台,DICOM使“改变临床医学面貌”的高级医学图像应用成为可能。比如在急诊科中,心脏负荷测试,乳腺癌的检查,DICOM为医生和病人服务,是医学成像有效工作的标准。
在医学图像生成后,在其他平台上显示医学图像并获得与图像相关信息是十分重要的。由于医学图像中所存储的数据是相对值,并且其灰度级一般为4096灰度级或者更高,在实际中,这远远高出显示器所能显示的最大灰度级范围(一般为256灰度级)。为了逼真的显示医学图像,需要将其灰度级映射到显示器设备能够显示的灰度级,并允许医护人员实时地根据需要。动态调节所显示的灰度范围,这就是对医学图像进行的调窗预处理技术。
目前针对DICOM格式医学图像显示的研究主要是线性调窗和动态线性调窗的快速实现,以及将DICOM格式医学图像转化成位图后的图像增强。然而,使用当前常用灰度映射函数的时候,由于DICOM标准中并没有对窗位的最大,最小值和窗宽的最大值作限制,当采用线性窗函数调节窗宽窗位时候,经常会出现像素值超出实际的像素值范围的现象,有必要将原有算法进行改进,改善图像的显示质量和效果,以及显示的实时性。
发明内容
本发明的目的在于提供一种DICOM医学影像动态线性调窗方法,以实现图像显示更加细致、清晰,增强图像的优点,并保证动态调窗图像的显示的实时性。
为实现上述目的,本发明采用的技术方案是:
一种DICOM医学影像动态线性调窗方法,包括以下步骤,
步骤二,设置DICOM图像窗口的窗宽和窗位,标记DICOM文件中的窗宽、窗位,通过灰度映射窗函数G(V)用于动态调窗。
步骤三,基于步骤二的灰度映射窗函数G(V)计算得出的结果,以及步骤一中的标签信息填充位图结构体,并将构造的位图显示出来;
所述步骤二中,先在整个图像内找到灰度最大值max和最小值min,设置窗宽W=max-min,窗位C=(max-min)/2,这样可以使图像不同灰度阶层的对比度达到最大,达到了一种直方图均衡化的效果,所述灰度最大值max和最小值min不是唯一确定的,而是根据调整后的截距和斜率来制约动态调整范围;
所述灰度映射窗函数G(V)方程如下:
其中,V表示图像数据,G(V)为图像的灰度值,W为窗宽(显示图像的范围),C为窗位(显示区域的中心位置),gmax和gmin分别是图像灰度值的最大值和最小值,S(n)和I(n)分别为窗函数的斜率和截距的绝对值。当S(n)和I(n)都取1的时候,上式(1.1)即变为基本灰度映射窗函数,这样做的好处是将输入值与调整后的窗函数的斜率与截距相关联,这样有效改善了原有算法的显示效果,当我们调整对话框来改变窗宽和窗位时候,图像的显示的实时性取得了很好的改善。可以看出,该方法不将gmax和gmin分别设为图像中最大最小灰度值,而是根据调整后的截距和斜率来制约动态调整范围。
与相关技术相比,本发明提供的具有以下有益效果:
本发明DICOM医学影像动态线性调窗方法使图像显示更加细致,增强图像的目的,当我们调整对话框来改变窗宽和窗位时候,图像的显示的实时性取得了很好的改善。
附图说明
图1本发明DICOM文件的组成;
图2为本发明实施例所述的采用基本窗函数线性调窗的腹部成像图;
图3为本发明实施例所述的采用改进窗函数线性调窗的腹部成像图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本发明的DICOM医学影像动态线性调窗方法,包括以下步骤:
步骤一,DICOM图像数据读取。
DICOM文件的组成和格式如图1,根据DICOM文件的标准格式读入图像的数据元素,数据元素包括标签信息、类型、长度、数据信息,并保存于内存之中。
步骤二,设置图像窗口的窗宽和窗位。
在研究DICOM解析数据发现,很多DICOM文件中未标记窗宽窗位,于是实验设定,在整个图像内找到灰度最大值max和最小值min,设置窗W=max-min,窗位C=(max-min)/2,这样可以使图像不同灰度阶层的对比度达到最大,达到了一种直方图均衡化的效果。但同时,我们不将gmax和gmin分别设为图像中最大最小灰度值,而是根据调整后的截距和斜率来制约动态调整范围。
步骤三,采用灰度映射窗函数G(V),映射灰度图像。
灰度映射窗函数G(V)方程如下:
V表示图像数据,G(V)为图像的灰度值,W为窗宽(显示图像的范围),C为窗位(显示区域的中心位置),gmax和gmin分别是图像灰度值的最大值和最小值,S(n)和I(n)分别为窗函数的斜率和截距的绝对值,当S(n)和I(n)都取1的时候,上式(1.1)即变为基本灰度映射窗函数,这样做的好处是将输入值与调整后的窗函数的斜率与截距相关联,这样有效改善了原有算法的显示效果,当我们调整对话框来改变窗宽和窗位时候,图像的显示的实时性取得了很好的改善。
对比试验效果
选取一个腹部的CT图像进行试验,gmin的值为3500,gmax的值为3900,S(n)=255/400,I(n)=255*35/4,通过本发明改进灰度映射函数得到的图像(图3),采用基本灰度映射函数(见式1.2)得到的图象(图2)与显示效果比较发现,本发明调窗预处理技术图像不仅对腹部能够更清晰地显示,同时能够对体表边缘进行显示,展现出更多的图像细节,图像效果有了很大的改善。
使用本发明线性调窗可以获得更好的图像显示效果,也为后期的图像处理打下基础,保证了动态调窗的实时性,为医学图像的研究及临床应用提供了方便。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (1)
1.一种DICOM医学影像动态线性调窗方法,其特征在于,包括以下步骤:
步骤一,DICOM图像数据读取,根据DICOM文件的标准格式读入图像的数据元素,数据元素包括标签信息、类型、长度、数据信息,并保存于内存之中;
步骤二,设置DICOM图像窗口的窗宽和窗位,标记DICOM文件中的窗宽、窗位,通过灰度映射窗函数G(V)用于动态调窗。
步骤三,基于步骤二的灰度映射窗函数G(V)计算得出的结果,以及步骤一中的标签信息填充位图结构体,并将构造的位图显示出来;
所述步骤二中,先在整个图像内找到灰度最大值max和最小值min,设置窗宽W=max-min,窗位C=(max-min)/2,这样可以使图像不同灰度阶层的对比度达到最大,达到了一种直方图均衡化的效果,所述灰度最大值max和最小值min不是唯一确定的,而是根据调整后的截距和斜率来制约动态调整范围;
所述灰度映射窗函数G(V)方程如下:
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>V</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>V</mi>
<mo>&le;</mo>
<mi>int</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mo>-</mo>
<mi>W</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mi>I</mi>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&lsqb;</mo>
<mi>V</mi>
<mo>-</mo>
<mi>int</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mo>-</mo>
<mi>W</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mi>I</mi>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
<mo>&rsqb;</mo>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mi>g</mi>
<mi>max</mi>
</msub>
<mi>W</mi>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>e</mi>
<mi>l</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>V</mi>
<mo>&GreaterEqual;</mo>
<mi>int</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>C</mi>
<mo>+</mo>
<mi>W</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mi>I</mi>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,V表示图像数据,G(V)为图像的灰度值,W为窗宽(显示图像的范围),C为窗位(显示区域的中心位置),gmax和gmin分别是图像灰度值的最大值和最小值,S(n)和I(n)分别为窗函数的斜率和截距的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710820652.3A CN107481208A (zh) | 2017-09-13 | 2017-09-13 | Dicom医学影像动态线性调窗方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710820652.3A CN107481208A (zh) | 2017-09-13 | 2017-09-13 | Dicom医学影像动态线性调窗方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107481208A true CN107481208A (zh) | 2017-12-15 |
Family
ID=60585082
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710820652.3A Withdrawn CN107481208A (zh) | 2017-09-13 | 2017-09-13 | Dicom医学影像动态线性调窗方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107481208A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
CN110223261A (zh) * | 2019-05-20 | 2019-09-10 | 上海联影医疗科技有限公司 | 医学图像处理方法和系统、存储介质及计算机设备 |
CN111803104A (zh) * | 2020-07-20 | 2020-10-23 | 上海市第六人民医院 | 一种医学影像显示方法、介质及电子设备 |
CN112053328A (zh) * | 2020-08-18 | 2020-12-08 | 深圳蓝韵医学影像有限公司 | 数字x射线尘肺图像的处理方法及系统 |
CN112562829A (zh) * | 2021-02-19 | 2021-03-26 | 南京景三医疗科技有限公司 | 一种自适应生成dicom影像默认窗宽窗位的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102485179A (zh) * | 2010-12-02 | 2012-06-06 | 深圳迈瑞生物医疗电子股份有限公司 | 影像窗宽/窗位调节方法、装置及医学影像设备 |
CN107145730A (zh) * | 2017-04-27 | 2017-09-08 | 广州慧扬健康科技有限公司 | Dicom医学影像文件读取适配系统 |
-
2017
- 2017-09-13 CN CN201710820652.3A patent/CN107481208A/zh not_active Withdrawn
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102485179A (zh) * | 2010-12-02 | 2012-06-06 | 深圳迈瑞生物医疗电子股份有限公司 | 影像窗宽/窗位调节方法、装置及医学影像设备 |
CN107145730A (zh) * | 2017-04-27 | 2017-09-08 | 广州慧扬健康科技有限公司 | Dicom医学影像文件读取适配系统 |
Non-Patent Citations (1)
Title |
---|
顾巍: "虚拟肝脏手术规划系统的关键技术研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
CN110223261A (zh) * | 2019-05-20 | 2019-09-10 | 上海联影医疗科技有限公司 | 医学图像处理方法和系统、存储介质及计算机设备 |
CN110223261B (zh) * | 2019-05-20 | 2022-03-01 | 上海联影医疗科技股份有限公司 | 医学图像处理方法和系统、存储介质及计算机设备 |
CN111803104A (zh) * | 2020-07-20 | 2020-10-23 | 上海市第六人民医院 | 一种医学影像显示方法、介质及电子设备 |
CN112053328A (zh) * | 2020-08-18 | 2020-12-08 | 深圳蓝韵医学影像有限公司 | 数字x射线尘肺图像的处理方法及系统 |
CN112053328B (zh) * | 2020-08-18 | 2024-06-28 | 深圳蓝影医学科技股份有限公司 | 数字x射线尘肺图像的处理方法及系统 |
CN112562829A (zh) * | 2021-02-19 | 2021-03-26 | 南京景三医疗科技有限公司 | 一种自适应生成dicom影像默认窗宽窗位的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107481208A (zh) | Dicom医学影像动态线性调窗方法 | |
Seeram et al. | Digital radiography | |
US8577110B2 (en) | Device, method and computer readable recording medium containing program for separating image components | |
Seibert | Tradeoffs between image quality and dose | |
Malkus et al. | A method to extract image noise level from patient images in CT | |
US20230306657A1 (en) | Noise suppression using deep convolutional networks | |
US10182783B2 (en) | Visualization of exposure index values in digital radiography | |
Carter et al. | Digital Radiography and PACS E-Book: Digital Radiography and PACS E-Book | |
US20050069186A1 (en) | Medical image processing apparatus | |
Senthil et al. | Enhancement Sushisen algorithms in images analysis technologies to increase computerized tomography images | |
CN111915523A (zh) | Dr图像亮度自适应调整方法及系统 | |
Seeram et al. | Digital image processing concepts | |
JP2002044413A (ja) | 放射線画像処理方法および放射線画像処理装置 | |
Carroll | Digital Radiography in Practice | |
KR20140134903A (ko) | 의료 영상 화질 개선 방법 및 그 장치 | |
Kim et al. | Refinement of image quality in panoramic radiography using a generative adversarial network | |
Collins et al. | Comparison of measurements from photographed lateral cephalograms and scanned cephalograms | |
Svenson et al. | Can adaptive post-processing of storage phosphor plate panoramic radiographs provide better image quality? A comparison of anatomical image quality of panoramic radiographs before and after adaptive processing | |
JPWO2007108346A1 (ja) | 放射線画像処理方法、放射線画像処理装置及び放射線画像撮影システム | |
Poggi | The Digital Image and Its Processing: Quality Controls in Mammography: An Overview | |
KR20020072638A (ko) | 흉부 단순 엑스선 영상에서 늑골의 경계를 자동으로추출하는 방법 | |
US20050047635A1 (en) | Method, apparatus and program for image processing | |
CN101536034A (zh) | 用于放射图像的色调等级变换 | |
NL7909276A (nl) | Stralingsbeeldverwerking. | |
Allen | Guidelines for Image Evaluation |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20171215 |
|
WW01 | Invention patent application withdrawn after publication |