CN102104784A - 一种大数据量的像素集窗宽窗位调节方法 - Google Patents
一种大数据量的像素集窗宽窗位调节方法 Download PDFInfo
- Publication number
- CN102104784A CN102104784A CN 201010157645 CN201010157645A CN102104784A CN 102104784 A CN102104784 A CN 102104784A CN 201010157645 CN201010157645 CN 201010157645 CN 201010157645 A CN201010157645 A CN 201010157645A CN 102104784 A CN102104784 A CN 102104784A
- Authority
- CN
- China
- Prior art keywords
- window
- image
- window width
- interest
- width
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及一种大数据量的像素集窗宽窗位调节方法,针对DICOM格式数字医学影像通过设置其DIB颜色对照表,使得灰阶影像在普通RGB显示器上以加窗方式显示;在实时调节影像的窗宽窗位的过程中,通过取模运算过滤响应频率;影像兴趣区截取缩放;建立调窗变换速查表,查调窗变换速查表遍历兴趣区像素集获得显示灰阶值等方法来降低像素集调窗过程中的计算量,使得即使对于数据量非常大的影像文件,也能做到实时,连续,平滑地调节窗宽窗位。
Description
技术领域
本发明涉及一种计算机处理数字影像窗宽窗位调节的方法,尤其涉及大数据量的像素集调节窗宽窗位过程中的速度优化方法。
背景技术
计算机处理影像的技术帮助人们方便地调节影像达到所需的效果,对于医学放射检查DICOM格式的影像,由于其反映的是人体组织在放射射线照射下的密度成像,所以其像素是只包含灰阶信息,不包含RGB彩色信息的影像数据。
而DICOM格式的影像数据不同于普通RGB图像的图像位存储方式,其通常使用8位,16位,24位位深来表示影像灰阶层次,灰阶层次通常达到4096灰度级甚至更高,而普通RGB显示器只能显示256级灰阶,所以必须通过加窗方法来查看不同灰阶的影像效果。
一张放射X光数字影像的像素通常达到千万像素级别,而传统的调节窗宽窗位算法利用如下像素转换方法:
上式各符号含义如下:
M(pix[I,j]):座际为i,j的像素在当前窗宽窗位条件下在显示器上显示的灰阶值;
Mmax: 当前显示器能够显示的最大灰阶值;
W: 当前窗宽值;
Pix[i,j]: 座标为i,j像素的像素值;
C: 当前窗位值;
使用该传统算法则在调节窗宽窗位过程中针对每个像素平均需要进行5次乘法(除法),4次加法(减法),2次比较运算,而当用户每秒移动鼠标调节窗宽窗位的过程中,每秒钟需要将影像文件所有的像素进行20~30轮以上的运算,这个计算量是非常巨大的,对于大数据量的影像像素集难于做到实时连续平滑的处理响应。
发明内容
本发明的目的是建立一种快速,平滑,优化的大数据量像素集窗宽窗位调节方法。
一种大数据量的像素集窗宽窗位调整方法,包括如下处理步骤:
第一步:针对DICOM影像文件,打开文件,装载像素数据到内存。
第二步:构建DIB颜色对照表,设置颜色对照表中Red=Green=Blue,使得影像以灰阶方式显示。
第三步:在用户执行窗宽窗位调节的计算过程中包含如下子处理步骤。
1.通过取模运算过滤减少调窗响应频率,通过针对一个配置参数的取模(Mod)运算降低对于用户连续操作的响应频率来降低运算负担,由于适应人眼的帧率是15~24帧/秒,而调窗过程中所涉及的是不发生图像位移,缩放的固定像素位置灰阶变换,是人眼对于变换帧率相对不敏感情况,所以取模运算配置参数最低可以设置到10帧/秒,人眼主观感受仍然是可以接受的流畅感觉,如图1所示。
2.对当前影像兴趣区进行截取缩放运算,使得需要进行调窗计算处理的像素集仅包括当前用户界面所正显示中的像素集,并且按照当前用户界面所正显示的宽,高值缩放影像,重构影像对象,通过这一步骤的处理,在窗宽窗位调节过程中忽略掉当前用户界面兴趣区之外的像素,如图2所示。
3.建立调窗变换速查表:通过计算两次调窗响应过程中的鼠标位置座标差,得到当前新的窗宽窗位值C,W;并通过计算当前窗宽变换斜率RSetp,递增斜率步长的方式来建立长度为最大像素值的调窗变换速查表。如图3所示。
4.遍历当前兴趣区像素集,通过查找获取下标为Pix[i,j]像素值的调窗变换速查表Table的对应显示灰阶值来获取得到当前窗宽窗位值下整个兴趣区像素集的灰阶显示数据。如图4所示。
附图说明
本发明包括如下附图:
图1是通过取模运算过滤调窗响应频率原理图;
图2是当前影像兴趣区截取缩放减少参与计算像素量原理图;
图3建立调窗变换速查表原理图;
图4查找调窗变换速查表获得像素灰阶显示值原理图。
具体实施方法
下面结合附图对本发明做进一步详细说明。
本发明是一种快速,平滑,优化的大数据量像素集窗宽窗位调节方法;针对DICOM格式灰阶影像文件,调节窗宽窗位运算进行速度优化,执行如下处理步骤;
第一步:打开DICOM格式灰阶影像文件:读取其元数据字段和影像数据字段,构建Hashtable数据结构对象,保存每个数据元素标签(TAG)与值类型(Value Representation),值的长度(Value Field)和数值(Value)的对应关系;并且读取像素存储区域到内存像素集数组对象;
第二步:构建DIB颜色对应表:普通RGB显示器通过256级色阶的RGB三种颜色组合来表示彩色,但是对于放射照射所产生的组织密度灰阶影像,则只显示其灰阶效果,通过遍历DIB对象的颜色表,将颜色对应表的R,G,B分量设置为相等的数值,这样,在将只包含灰阶影像信息的DICOM影像显示在普通的RGB显示器上显示时就会表现出256级灰阶的灰度图显示效果。
第三步:响应用户的鼠标调窗操作其中包含如下子方法步骤:
1.通过取模运算过滤减少调窗响应频率:由于人眼正常运动图像帧率适应频率是15~24帧/秒,而在执行窗宽窗位的调节过程中,用户眼前所呈现的是一幅位置固定,但每个像素点的灰阶值在连续改变的动态影像,没有像素移动或者缩放操作,动态图像效果表现为是人眼相对不敏感的固定像素位置灰阶值变化过程,所以在配置较低的机器运行环境下,可以通过取模参数配置将每秒响应帧率降低到10帧/秒。
在用户鼠标移动调节窗宽,窗位的过程中,系统通过取模预算(MOD)一个配置参数,间隔地过滤掉参数阈值之外的鼠标移动响应事件,使得在连续的鼠标移动事件中,系统可以忽略掉一部分调窗计算量,避免由于系统计算负担过大而导致的系统停滞不响应现象,给用户的调窗操作体验是连续平滑的,如图1所示。
2.对当前影像兴趣区进行截取缩放运算:一个数字X光放射检查所产生的影像通常有一千万,甚至几千万像素,如果整体遍历全部像素则是一个非常庞大的计算数据量。而用户在某一时刻正在用户界面上所查看的图像的像素集则只与用户当前阅片窗口的界面大小有关,通常是百万,数百万像素级别,所以减少连续调窗过程中的计算数据量的重要方法之一就是忽略当前用户兴趣区之外的影像数据,只处理当前用户兴趣区范围之内的数据,并且通过影像截取缩放运算构建一个与当前用户兴趣区DIB座标一致,大小相同的DICOM影像对象,并将这个影像对象作为调窗运算的数据对象,从而避免针对整个DICOM影像文件进行整体运算的计算负担,如图2所示。
3.建立调窗变换速查表:首先通过获取鼠标移动连续调窗事件响应过程中的两次响应事件之间的鼠标位置座标差值,与当前影像窗宽窗位值相加(减),得到新的当前窗宽窗位值C,W;
通过获取当前显示器的最大灰阶级数Mmax除以当前灰阶影像的最大灰阶值PMax得到当前窗宽变换斜率RSetp;
构建一个长度为PMax的调窗变换速查表,表中数值从下标为零开始递增,递增步长为RSetp;
而对于像素值Pix[i,j]其值小于最小窗宽,或大于最大窗宽的情况则分别赋值0或者Mmax,如图3所示。
4.查调窗变换速查表遍历兴趣区像素集设置显示灰阶值:遍历兴趣区像素集的所有像素,以每一个像素的实际灰阶值Pix[i,j]作为查找调窗变换速查表Table的下标索引,并将该下标值所对应的调窗变换表数值Table[Pix[i,j]]更新目标显示DIB中的显示器灰阶值用于显示调窗后的影像效果,如图4所示。
Claims (5)
1.一种大数据量的像素集窗宽窗位调节方法,针对影像像素数据进行实时窗宽窗位调节,通过减少调窗过程中的计算量负担,使调窗操作平滑连续,其特征在于其包含如下实现步骤:
(1)取模运算过滤事件响应频率;
(2)影像兴趣区截取缩放;
(3)建立调窗变换速查表;
(4)查调窗变换速查表获取显示灰阶值等方法来降低像素集调窗过程中的计算负担来实现连续平滑的调窗效果。
2.根据权利要求1所述一种大数据量的影像像素集的窗宽窗位调节方法,其特征在于所述的取模运算过滤事件响应频率处理,通过取模运算(MOD)配置参数,间隔地过滤掉参数阈值之外的鼠标移动响应事件,使得在连续的鼠标移动事件中,系统可以忽略掉一部分调窗计算量,避免由于系统计算负担过大而导致的系统停滞不响应现象。
3.根据权利要求1所述一种大数据量的影像像素集的窗宽窗位调节方法,其特征在于所述进行影像兴趣区截取缩放处理步骤忽略当前用户兴趣区之外的影像数据,只处理当前用户兴趣区范围之内的数据,并且通过影像截取缩放运算构建一个与当前用户兴趣区DIB座标-致,大小相同的DICOM影像对象,并将这个影像对象作为调窗运算的数据对象,从而避免针对整个DICOM影像文件进行整体运算的计算负担。
4.根据权利要求1所述一种大数据量的影像像素集的窗宽窗位调节方法,其特征在于所述的建立调窗变换速查表步骤,通过获取当前显示器的最大灰阶级数Mmax除以当前灰阶影像的最大灰阶值PMax得到当前窗宽变换斜率RSetp;并构建一个长度为PMax的调窗变换速查表,表中数值从下标为零开始递增,递增步长为RSetp;而对于像素值P[i,j]其值小于最小窗宽(C-W/2),或大于最大窗宽(C+W/2)的情况则分别赋值0或者Mmax。
5.根据权利要求1所述一种大数据量的影像像素集的窗宽窗位调节方法,其特征在于所述的查调窗变换速查表遍历兴趣区像素集设置显示灰阶值:遍历兴趣区像素集的所有像素,以每一个像素的实际灰阶值作为查找调窗变换速查表的下标索引,并用该下标值所对应的调窗变换表数值更新目标显示DIB中的显示器灰阶值,用于显示调窗后的影像效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010157645 CN102104784A (zh) | 2010-04-28 | 2010-04-28 | 一种大数据量的像素集窗宽窗位调节方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010157645 CN102104784A (zh) | 2010-04-28 | 2010-04-28 | 一种大数据量的像素集窗宽窗位调节方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102104784A true CN102104784A (zh) | 2011-06-22 |
Family
ID=44157252
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010157645 Pending CN102104784A (zh) | 2010-04-28 | 2010-04-28 | 一种大数据量的像素集窗宽窗位调节方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102104784A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102760294A (zh) * | 2012-08-01 | 2012-10-31 | 赵建峰 | 一种高精度影像高度替代联动显示方法 |
CN103914852A (zh) * | 2014-03-14 | 2014-07-09 | 兰州交通大学 | 基于cuda的dicom医学影像动态非线性调窗方法 |
CN104333738A (zh) * | 2014-10-29 | 2015-02-04 | 李建 | 一种基于dicom的医学高清动态及静态图像存贮系统 |
CN104851083A (zh) * | 2015-05-29 | 2015-08-19 | 上海京颐科技股份有限公司 | 医学影像的调窗方法及装置 |
CN104905808A (zh) * | 2014-03-12 | 2015-09-16 | 上海联影医疗科技有限公司 | 定位片图像的实时显示方法和定位片图像的实时显示系统 |
CN107145730A (zh) * | 2017-04-27 | 2017-09-08 | 广州慧扬健康科技有限公司 | Dicom医学影像文件读取适配系统 |
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
CN110928469A (zh) * | 2018-09-19 | 2020-03-27 | 上海西门子医疗器械有限公司 | 设置最佳显示窗宽和窗位的方法 |
US10885700B2 (en) * | 2018-06-27 | 2021-01-05 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for fusing image data |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1220859A (zh) * | 1997-09-01 | 1999-06-30 | 通用电器横河医疗系统株式会社 | 图象显示方法、图象显示装置和基于图象的诊断装置 |
JP2000067231A (ja) * | 1998-08-18 | 2000-03-03 | Hitachi Medical Corp | 画像表示装置 |
CN1547161A (zh) * | 2003-12-08 | 2004-11-17 | 西安理工大学 | 彩色多窗ct图像的自动生成方法 |
CN1722838A (zh) * | 2004-07-15 | 2006-01-18 | 三星电子株式会社 | 使用基础层的可伸缩性视频编码方法和设备 |
CN101105862A (zh) * | 2007-08-02 | 2008-01-16 | 宁波大学 | 一种医学图像窗口参数的自适应调整方法 |
-
2010
- 2010-04-28 CN CN 201010157645 patent/CN102104784A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1220859A (zh) * | 1997-09-01 | 1999-06-30 | 通用电器横河医疗系统株式会社 | 图象显示方法、图象显示装置和基于图象的诊断装置 |
JP2000067231A (ja) * | 1998-08-18 | 2000-03-03 | Hitachi Medical Corp | 画像表示装置 |
CN1547161A (zh) * | 2003-12-08 | 2004-11-17 | 西安理工大学 | 彩色多窗ct图像的自动生成方法 |
CN1722838A (zh) * | 2004-07-15 | 2006-01-18 | 三星电子株式会社 | 使用基础层的可伸缩性视频编码方法和设备 |
CN101105862A (zh) * | 2007-08-02 | 2008-01-16 | 宁波大学 | 一种医学图像窗口参数的自适应调整方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102760294A (zh) * | 2012-08-01 | 2012-10-31 | 赵建峰 | 一种高精度影像高度替代联动显示方法 |
CN104905808A (zh) * | 2014-03-12 | 2015-09-16 | 上海联影医疗科技有限公司 | 定位片图像的实时显示方法和定位片图像的实时显示系统 |
CN103914852A (zh) * | 2014-03-14 | 2014-07-09 | 兰州交通大学 | 基于cuda的dicom医学影像动态非线性调窗方法 |
CN104333738A (zh) * | 2014-10-29 | 2015-02-04 | 李建 | 一种基于dicom的医学高清动态及静态图像存贮系统 |
CN104851083A (zh) * | 2015-05-29 | 2015-08-19 | 上海京颐科技股份有限公司 | 医学影像的调窗方法及装置 |
CN107145730A (zh) * | 2017-04-27 | 2017-09-08 | 广州慧扬健康科技有限公司 | Dicom医学影像文件读取适配系统 |
CN109035203A (zh) * | 2018-06-25 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | 医学图像处理方法、装置、设备及存储介质 |
US10885700B2 (en) * | 2018-06-27 | 2021-01-05 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for fusing image data |
US11200727B2 (en) | 2018-06-27 | 2021-12-14 | Shanghai United Imaging Healthcare Co., Ltd. | Method and system for fusing image data |
CN110928469A (zh) * | 2018-09-19 | 2020-03-27 | 上海西门子医疗器械有限公司 | 设置最佳显示窗宽和窗位的方法 |
CN110928469B (zh) * | 2018-09-19 | 2023-08-22 | 上海西门子医疗器械有限公司 | 设置最佳显示窗宽和窗位的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102104784A (zh) | 一种大数据量的像素集窗宽窗位调节方法 | |
CN104408752B (zh) | 一种基于混合色调映射算法的高动态范围图像压缩方法 | |
CN102741879B (zh) | 由单眼图像产生深度图的方法及其系统 | |
DE112018007721T5 (de) | Aufnehmen und modifizieren von 3D-Gesichtern unter Verwendung neuronaler Bild-und Zeitverfolgungsnetze | |
US20210341998A1 (en) | Gaze-point determining method, contrast adjusting method, and contrast adjusting apparatus, virtual reality device and storage medium | |
CN107886503A (zh) | 一种消化道解剖位置识别方法及装置 | |
CN106875361A (zh) | 一种基于深度卷积神经网络去除泊松噪声的方法 | |
CN106897981A (zh) | 一种基于引导滤波的低照度图像增强方法 | |
CN109993698A (zh) | 一种基于生成对抗网络的单幅图像超分辨纹理增强方法 | |
CN103024300A (zh) | 一种高动态范围图像显示方法及装置 | |
CN105469407A (zh) | 一种基于改进的引导滤波器的人脸图像图层分解方法 | |
CN101930595A (zh) | 图像处理方法和图像处理设备 | |
CN112561846A (zh) | 训练图像融合模型的方法、装置和电子设备 | |
CN107610675A (zh) | 一种基于动态电平的图像处理方法及装置 | |
CN110634147A (zh) | 基于双边引导上采样的图像抠图方法 | |
CN110349085A (zh) | 一种基于生成对抗网络的单幅图像超分辨特征增强方法 | |
CN109727316B (zh) | 虚拟现实图像的处理方法及其系统 | |
Liu et al. | Learning noise-decoupled affine models for extreme low-light image enhancement | |
CN103106671B (zh) | 一种基于视觉注意机制的图像感兴趣区域检测方法 | |
CN103413339A (zh) | 十亿像素高动态范围图像重建与显示的方法 | |
Qian et al. | An adaptive enhancement algorithm based on visual saliency for low illumination images | |
Everingham et al. | Head-mounted mobility aid for low vision using scene classification techniques | |
CN106023057B (zh) | 一种用于皮下静脉显影仪的控制处理系统及成像方法 | |
CN113298744B (zh) | 一种端到端的红外与可见光图像融合方法 | |
Wang et al. | Learning pixel-adaptive weights for portrait photo retouching |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
DD01 | Delivery of document by public notice |
Addressee: Liang Wei Document name: the First Notification of an Office Action |
|
DD01 | Delivery of document by public notice |
Addressee: Liang Wei Document name: Notification that Application Deemed to be Withdrawn |
|
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110622 |