CN1259890C - 图像处理装置和方法 - Google Patents
图像处理装置和方法 Download PDFInfo
- Publication number
- CN1259890C CN1259890C CNB2003101018998A CN200310101899A CN1259890C CN 1259890 C CN1259890 C CN 1259890C CN B2003101018998 A CNB2003101018998 A CN B2003101018998A CN 200310101899 A CN200310101899 A CN 200310101899A CN 1259890 C CN1259890 C CN 1259890C
- Authority
- CN
- China
- Prior art keywords
- image
- mentioned
- difference image
- difference
- shade
- 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.)
- Expired - Fee Related
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 title claims description 65
- 230000008859 change Effects 0.000 claims abstract description 24
- 230000009466 transformation Effects 0.000 claims description 45
- 238000002059 diagnostic imaging Methods 0.000 claims description 22
- 239000004744 fabric Substances 0.000 claims description 17
- 239000008280 blood Substances 0.000 claims description 4
- 210000004369 blood Anatomy 0.000 claims description 4
- 238000012856 packing Methods 0.000 claims 1
- 230000008034 disappearance Effects 0.000 abstract description 19
- 238000006243 chemical reaction Methods 0.000 abstract description 8
- 230000002542 deteriorative effect Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 33
- 101100021996 Arabidopsis thaliana CYP97C1 gene Proteins 0.000 description 25
- 101100510695 Arabidopsis thaliana LUT2 gene Proteins 0.000 description 25
- 238000010586 diagram Methods 0.000 description 12
- 230000009471 action Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000002372 labelling Methods 0.000 description 6
- 239000000654 additive Substances 0.000 description 4
- 230000000996 additive effect Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000004195 computer-aided diagnosis Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000008676 import Effects 0.000 description 3
- 210000004072 lung Anatomy 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 0 **(C(CCC1=CC2=C1)=CC13)=C1C3=CC2=C Chemical compound **(C(CCC1=CC2=C1)=CC13)=C1C3=CC2=C 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- XILPGHLPTFYDNG-UHFFFAOYSA-N C(C1)CC1C1CCCC1 Chemical compound C(C1)CC1C1CCCC1 XILPGHLPTFYDNG-UHFFFAOYSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/24—Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
-
- 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/10116—X-ray image
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提供一种图像处理方法,其中,基于由用户选择的显示模式,在生成的差分图像中选择了仅显示阴影增大(消失)时,将差分图像进行灰阶变换,使得只有阴影增大(消失)的区域的亮度与其他区域的亮度不同,通过显示灰阶变换后的差分图像,使阴影变化的状态按状态单独作为图像来显示,使得能够容易地观察阴影变化。另外,抽取出过去图像和当前图像中的负图像/正图像信息,进而输入差分图像上的阴影和诊断上的判断基准的对应关系。阴影定义规定,对于图像诊断上病变部的加深或减轻,其在差分图像上对应于高浓度(高亮度值)或者低浓度(低亮度值)之一。
Description
技术领域
本发明涉及图像处理方法,特别是适合用于显示医疗图像中的被检测体的时效变化的图像处理方法。
背景技术
近年,为响应关于医用图像的数字化的国际标准DICOM(digitalinformation and communication in medicine),日本放射线装置工业会(JIRA)制定了MIPS(medical image processing system)标准,进行医用图像的数字化。
在医用图像领域中,通过计算机分析医用X射线图像、CT扫描图像、MRI图像等的CAD(Computer Aided Diagnosis)的使用正在发展,其中特别关注取得拍摄相同部位的图像的时效差别的时效差分方式的CAD。
在该时效差分方式的CAD中,例如对不同时间点拍摄的1组胸部纯X射线图像进行图像分析。并且,在各个图像中寻求解剖学上的相同位置,使当前或过去的图像之一发生变形,进行每个像素的差分处理(例如参照专利文献1)。
由这样的差分处理而得到的差分图像的亮度值对应于当前和过去的图像间的图像信号的变化。即,如果过去图像与当前图像之间没有变化,则其差分值为0,如果产生某些变化,则呈现与该变化对应的亮度值的变化。
图22是说明不同时间的相同部位中的图像信号的变化例的图,用1维图来表示过去图像的信号和当前图像的信号,及由此得到的差分图像的信号。
在图22A的过去图像的信号中,如果在大致平坦轮廓的部分出现了用图22B的区域A所示那样的表示阴影的信号,则从过去图像的信号中减去当前图像的信号,取出其差分,生成差分图像的信号后,得到图22C所示那样的轮廓。
另外,图22D是以产生上述变化时的胸部X射线图像为对象的差分图像的例子,如该图所示,描绘出与新产生的阴影对应的亮度值低的区域A。这样,通过显示差分图像,就能容易地观察2幅图像间产生的变化了。
如上述那样,使用差分图像后,具有更容易地观察在不同时间点产生的图像的变化的效果。
然而,在图像间的变化中,并非总只发生阴影的增大或消失的一种,有时两者都发生。在这样的情况下,存在如下问题,即读影者(医师)必须对差分图像显示的各个阴影的变化十分注意。
而且还存在这样的问题,即差分图像所显示的阴影的个数较多时,容易发生漏看阴影变化的情况。
还存在这样的问题,即在观察中,不仅是各个阴影的变化,还必须把握整体进程,但是,如果阴影个数较多,就很难判断这样的进程。
在放射线图像中,大像素值在显示上是与高亮度对应还是与低亮度对应,有时是因图像生成装置不同而不同的。例如,像近年正进行实用化的FPD(Flat Panel Detector)那样,在将透过被检测体的X射线直接或间接地变换成电信号、生成数字图像的装置中,作为正图像而生成图像数据。在该正图像中,具有高X射线透过率的诸如肺野等部位的像素值变大,大像素值在显示上与高亮度对应。因此,在例如胸部X射线图像中,对应于肺野的像素值变大了。
另一方面,因为通常放射线图像都是用负图像进行观察的,所以用胶片扫描器等将X射线胶片数字化后的图像数据与肺野对应的像素值变小了。
因此,使负图像和正图像混合进行差分处理后,差分图像中的阴影出现方式不统一,是极难进行判断的。
发明内容
本发明是鉴于上述问题点而做出的,其目的在于,在显示差分图像时,即使是发生多个阴影变化的情况,也能够容易且准确地观察每个阴影变化,同时能够容易且准确地进行时效观察中的判断。
进而,本发明的目的在于提供一种生成即使是如上述那样在图像的正/负不统一的图像间生成也具有一贯性的差分图像的方法和装置,特别是提供一种可以输出在差分图像上的阴影的显示状态上具有一贯性的差分图像的方法和装置。
本发明提供一种图像处理方法,其特征在于,包括:输入步骤,输入在不同时间点拍摄的至少2幅医用图像;差分图像生成步骤,在由上述输入步骤输入的2幅医用图像间进行位置重合,差分处理对应的坐标点的图像信号并生成差分图像;差分图像存储步骤,将由上述差分图像生成步骤生成的差分图像存储到存储介质中;基准差分值确定步骤,作为灰阶变换处理的基准值确定由上述差分图像存储步骤所存储的差分图像的基准差分值;状态选择步骤,使用户选择成为由上述差分图像生成步骤所生成的差分图像的读影(interpret)的指标的阴影变化的状态;差分图像处理步骤,根据由上述基准差分值确定步骤确定的基准差分值和由上述状态选择步骤选择的阴影变化的状态,来对上述差分图像进行灰阶变换处理;以及输出步骤,输出由上述差分图像处理步骤处理过的差分图像。
本发明的其他目的、特征以及优点,通过以附图为参照的下面的说明将会弄明白。此外,在附图中对相同或相似的结构附加相同的参照标号。
附图说明
图1是表示本发明的第1实施方式、表示图像处理装置的基本结构的一例的图。
图2是表示本发明的第1实施方式、说明图像处理装置的各部分的详细动作(图像处理方法)的一例的流程图。
图3是表示本发明的第1实施方式、表示差分图像的直方图的一例的图。
图4是表示本发明的第1实施方式、表示显示单元显示的监控器画面的例子的图。
图5是表示本发明的第1实施方式、表示对应各显示模式的灰阶变换特性的一例的图。
图6是表示本发明的第1实施方式、表示显示单元显示的实际图像的例子的图。
图7是表示本发明的第2实施方式、表示图像处理装置的基本结构的一例的图。
图8是表示本发明的第2实施方式、说明图像处理装置的各部分的详细动作(图像处理方法)的一例的流程图。
图9是表示本发明的第2实施方式、说明分析单元的更详细的动作的流程图。
图10是表示本发明的第2实施方式、表示显示单元显示的监控画面的例子的图。
图11是表示本发明的第3实施方式、表示对应各显示模式的灰阶变换特性的一例的图。
图12是说明作为第4实施方式的医用图像处理装置的结构的框图。
图13是表示作为第4实施方式的医用图像处理装置的内部结构的框图。
图14是说明作为第4实施方式的医用图像处理装置的功能结构的框图。
图15是说明查找表的特性的图。
图16是说明查找表的组合的图。
图17是说明差分图像中的阴影定义的图。
图18是说明位置重合单元的动作的流程图。
图19是图像信号的轮廓变化的说明图。
图20是说明作为第5实施方式的医用图像处理装置的功能结构的框图。
图21是说明第5实施方式的变形例的功能结构的框图。
图22是表示以往的技术、说明不同时间点的相同部位中的图像信号的变化例的图。
具体实施方式
下面参照附图详细说明本发明的实施方式。
<第1实施方式>
下面参照附图说明本发明的第1实施方式。
图1是图示适用本实施方式的图像处理装置的基本结构的一例的图。
首先说明整体的动作。由图像输入单元1输入的至少2幅图像在差分处理单元2中进行处理,生成差分图像。上述差分图像暂时存储在存储器单元3中,在适当灰阶变换单元4中实施必要的处理,在显示单元5中输出显示。此处显示单元5是例如高精细的CRT监控器、液晶监控器等的显示装置,据此进行医师的诊断。
另外,在以后的说明中,作为上述差分处理的基础的1对图像也同样存储在存储器单元3中,在适当显示单元5中输出显示,与上述差分图像一同提供给诊断用。
下面,基于图2的流程图,说明本实施方式中的图像处理装置的各部分的详细动作(图像处理方法)的一例。
[步骤S100]
首先,在步骤S 100中,由图像输入单元1输入1对图像。在图1中,图像输入单元1是输入由无图示的摄影装置拍摄的多幅图像的装置,从例如硬盘那样的外部存储装置输入提供给差分处理的1对图像。
在此,作为差分处理的对象的1对图像是在不同时间点由预定的摄影装置进行拍摄、预先存储好的,但也不限定于此。例如,也可以是1幅图像已经拍摄好、存储在外部存储装置中了,而当前时间点的另1幅图像拍摄后直接经图像输入单元1输入。
而且,图像输入单元1也可以是输入保存在经网络连接的外部的存储装置中的图像的接口装置。
[步骤S200]、[步骤S300]
在上述步骤S 100中输入的1对图像在差分处理单元2中进行位置重合及差分处理。据此,生成差分图像,存储在存储器单元3中。关于上述1对图像的位置重合及差分处理的细节,可以使用公知的技术,例如日本特开平7-37074号公报所揭示的方法。
具体来说,例如,在数字化的1对图像中,其中一幅有非线性变形,在该变形了的图像与另一幅不变形的图像间进行相减处理,进行两图像的位置重合,生成上述差分图像。
另外,在本实施方式中,使表示产生的阴影变化的值如图12所示那样。具体来说,就是使负值与阴影增大的情况相对应,使正值与阴影消失的情况相对应。但是,这样的关系在本发明中并非必须是这样,也可以是相反的关系。
[步骤S400]
基准值确定单元6分析存储在存储器单元3中的差分图像,确定该差分图像中的基准值,存储在无图示的内部存储器中。在此,上述基准值,在差分图像中是对应于2幅图像间不发生变化的部分的像素的亮度的值,在差分图像中通常为0。
但是,由于存在位置重合误差、2幅图像的摄影条件不同,所以在1对图像间均匀地存在亮度差时,上述基准值也可能不为0。因此,在本实施方式中,基于差分图像的直方图来确定上述基准值。
图3是表示差分图像的直方图的一例的图。
在时效差分处理中,由于进行2幅图像间的位置重合,所以与无变化的部分对应的亮度值的频率最高。因此,基准值确定单元6寻求直方图30中频率最高的亮度值Cs,将该亮度值Cs作为上述基准值Cs进行存储。
[步骤S500]
显示模式选择单元7,在显示上述差分图像时,可由来自用户的输入选择以什么样的阴影变化作为对象进行显示。作为2幅图像间的阴影变化,考虑有
(1)无变化
(2)阴影增大
(3)阴影消失
三种状态,用户通过显示模式选择单元7,在下面的3种显示模式中选择其一,
(1)显示全部的阴影变化
(2)仅显示阴影增大
(3)仅显示阴影消失
在上述显示模式的选择中,使用如下用户接口,即用户在可对话性地输入的画面上,可以唯一地选择上述3种的显示模式。
图4是表示由显示选择单元7进行选择时,显示单元5所显示的监控器画面的例子的图。
在该例中,上述差分图像是显示在相当于显示单元5的监控器中并且与显示该差分图像的画面相同的画面,上述显示模式成为可选择的状态。但是,也可以用与差分图像不同的监控器及与图4所示那样的画面分开设置的开关等来选择上述显示模式。
在图4的例子中,画面上显示3个按钮41至43,只能选择其中之一。据此,只选择显示全部阴影或显示增大或消失的阴影,该选择结果被输出到灰阶变换单元4。另外,在图4的画面内显示的箭头是鼠标指针。
[步骤S600]
灰阶变换单元4基于从显示模式选择单元7输入的显示模式,变换存储在存储器单元3中的差分图像,输出到显示单元5中。下面说明对应于各显示模式的变换方法。
图5是图示灰阶变换单元4中的对应各显示模式的输入亮度和输出亮度的关系(灰阶变换特性)的一例的图。
如果显示模式是显示全部阴影变化的全阴影变化显示模式,则灰阶变换单元4基于图5(a)的关系变换差分图像并输出。如果显示模式是仅显示阴影增大的阴影增大显示模式,则灰阶变换单元4基于图5(b)的关系变换差分图像并输出。如果显示模式是仅显示阴影消失的阴影消失显示模式,则基于图5(c)的关系变换差分图像并输出。
在该图中,Smin、Smax表示输入、即差分图像的亮度的最小值和最大值。Dmin、Dmax表示输出、即显示的输出图像的亮度的最小值和最大值。Cd代表输出图像中的亮度值的中间值,它与上述的差分图像中的基准值Cs相对应。
因此,通过上述方法进行灰阶变换后,如果是上述全阴影变化显示模式,则只将差分图像在输出图像的亮度范围内做线性映射。另外,如果是上述阴影增大显示模式,则将比基准值Cs值大的亮度值全部映射为中间值Cd。进而,如果是上述阴影消失显示模式,则将比基准值Cs值小的亮度值全部映射为中间值Cd。
在此,因为中间值Cd在原始的差分图像中对应于未发生变化的部分的亮度值,所以在上述阴影增大显示模式中,只有增大的阴影变化被观察到,在上述阴影消失显示模式中,只有消失的阴影变化被观察到。
[步骤S700]
显示单元5,如上述那样,输入并显示在灰阶变换单元4中进行灰阶变换后的差分图像。作为显示的形式,既可以例如图4所示那样,将过去图像或当前图像的一者与经上述灰阶变换后的差分图像(上述输出图像)一并显示,也可以一并将上述过去图像、上述当前图像及经上述灰阶变换后的差分图像同时显示。但是,显示的形式不限于此,可以采取各种各样的形式。
图6是表示显示单元5所显示的实际图像的一例的图,图6A对应于上述全阴影变化显示模式,图6B对应于上述阴影增大显示模式,图6C对应于上述阴影消失显示模式。通过如上述那样进行灰阶变换,可以如图6所示,单独显示阴影的变化进行观察。
如以上那样,在本实施方式中,如果仅显示阴影增大(消失),则将阴影消失(增大)的区域的亮度映射为与基准值Cs对应的中间值Cd,处理差分图像进行显示,使得只有阴影增大(消失)的区域与其他区域相区别,因此,即使是在差分图像中阴影增大的区域和消失的区域混合在一起、包含多种阴影变化的情况,也可以分别将阴影变化的状态按状态作为图像进行显示。据此,能够容易地观察阴影变化,还能容易地进行病灶等的时效观察中的判断。因此,读影者不会漏看阴影的增大或消失,能容易地进行正确的读影(interpret),同时,能够容易且准确地把握差分图像整体的状态变化。
另外,在本实施方式中,如果仅显示阴影增大(阴影消失),则将阴影消失(增大)的区域的亮度映射为与基准值Cs对应的中间值Cd,变换成无阴影变化的状态的亮度进行输出,但如果只是将阴影增大(消失)的区域与其他区域相区别,就未必一定要将上述阴影消失(增大)的区域变换成无阴影变化的状态的亮度。例如,可以使上述阴影消失(增大)的区域的对比度相对降低。
<第2实施方式>
下面说明本发明的第2实施方式。在本实施方式中,实施在第1实施方式中说明的灰阶变换后,通过对发生阴影变化的部分进行分析处理,来显示阴影发生了怎样程度的增大或消失。这样,因为本实施方式只有施加了上述灰阶变换后的处理与上述第1实施方式不同,所以对与上述第1实施方式相同的部分,用与图1至图6相同的符号来表示,省略其详细说明。
图7是表示适用本实施方式的图像处理装置的基本结构的一例的图。
在该图中,本实施方式的图像处理装置除追加了分析单元8外,与图1所示的第1实施方式的图像处理装置的结构相同。
由灰阶变换单元4生成的输出图像输入到显示单元5及显示单元8中。
下面,基于图8所示的流程图,说明本实施方式的图像处理装置的各部分的详细动作(图像处理方法)的一例。但是,在图8中,关于步骤S800以外的各步骤,因为与上述第1实施方式(图2)相同,故省略说明。
[步骤S800]
在灰阶变换单元4中进行灰阶变换后的差分图像(输出图像)由分析单元8进行2值化,对上述阴影增大显示模式和上述阴影消失模式分别进行分析。具体来说,即对各显示模式分别计算表示阴影变化的区域的个数、面积等的特征量。
基于下面的图9所示的流程图,说明分析单元8的更详细的动作。
[步骤S801]
分析单元8从显示模式选择单元7输入当前的显示模式。如果显示模式是上述全阴影变化显示模式,则就此结束处理。
此时,来自分析单元8的输出变为无效,通过分配某特别的值,可以由后续的显示单元5中识别来自分析单元8的输出无效。另外,显示单元5识别显示模式,如果识别为显示模式是上述全阴影变化显示模式,也可以使来自分析单元8的输入变为无效(忽视)。
[步骤S802]、[步骤S803]
而且,分析单元8,如果当前的显示模式是上述阴影增大显示模式,则反转已输入的灰阶变换后的差分图像的灰阶,进入后续的处理。据此,在后述的2值化处理中,可以将发生阴影变化的部分统一分配为白像素。
[步骤S804]
分析单元8设定用于对输入的灰阶变换后的差分图像进行2值化的阈值。作为上述阈值的设定方法,可以使用公知技术,例如p平铺显示法。上述p平铺显示法中的面积比率的设定是根据对阴影变化的敏感度而预先设定的。
[步骤S805]
分析单元8基于在上述步骤S804中设定的阈值,将灰阶变换后的差分图像2值化,生成2值图像。据此,得到2值图像,在该2值图像中,对发生了增大或消失的阴影变化的部分分配1(白),对其他部分分配0(黑)。
[步骤S806]
对在上述步骤S805中生成的2值图像进行标记处理,对孤立了的阴影变化进行区域分割,计数分割后的区域的个数。
[步骤S807]
通过在上述步骤S806中进行的标记处理,求出分割后的各区域的面积,计算其合计。
[步骤S808]
输出在上述步骤S806、S807中求得的区域的个数和各区域的面积的合计。
图10是表示通过本实施方式中的以上处理、在显示单元5中显示的监控画面的一例的图。
在该例中,被选择的显示模式的发生阴影变化的区域的个数和面积被显示在处于差分图像显示区域(显示灰阶变换后的差分图像的区域)1001下方的显示区域1002内。据此,在各阴影变化的状态中,就可以定量地把握发生了怎样程度的变化了。
如以上那样,在本实施方式中,因为计算并显示发生阴影变化的区域的个数和面积,所以,除了获得上述第1实施方式的效果外,还能使读影者定量地把握阴影变化的状态,更容易且准确地把握差分图像整体的状态变化。
另外,在以上说明中,是合计并显示通过标记处理所分割了的各区域的面积。也可以是对每个区域分别计算面积,并将其在显示区域1002内一并显示。
<第3实施方式>
下面,说明本发明的第3实施方式。另外,本实施方式只有在进行上述灰阶变换时的处理与上述第1及第2实施方式不同,所以,对与上述第1实施方式相同的部分,用与图1至图10相同的符号来表示,省略其详细说明。
在上述第1实施方式中,用灰阶变换单元4进行了灰阶变换(参照图5),使得对输入亮度值的输出亮度值的特性分区地成为线性,但灰阶变换的方法不必限于此。
图11是图示本实施方式的灰阶变换特性的一例的图。
图11A是显示模式为上述阴影增大显示模式(仅显示阴影的增大的显示模式)时的特性。
如特性1101a所示,对于基准值Cs以下的输入,将灰阶改变为非线形,对于超过基准值Cs的输入,改变为中间值Cd。
通过这样处理,对于阴影变化值小的区域可以抑制输出,因此可以抑制由于位置重合不齐而产生的噪音。据此,可以更容易地观察阴影变化。
另外,图11B是显示模式为上述阴影消失显示模式(仅显示阴影的小时的显示模式)时的特性。
如特性1102a所示,对于基准值Cs以上的输入,将灰阶改变为非线形,对于比基准值Cs小的输入,改变为中间值Cd。
通过这样处理,对于与先前同样的阴影变化值小的区域,可以抑制输出,因此能够更容易地观察阴影变化。
进而,也可以进行灰阶变换,使之如特性1101b、1102b(图11的虚线)所示那样。这样处理后,例如显示模式是上述阴影增大显示模式时,消失的阴影的区域也能够以低对比度来确认其存在,因此,能容易地同时把握阴影的增大和消失的关系。
另外,在上述的第1至第3实施方式中,如图4所示,是从3种显示模式中选择其一,显示1幅灰阶变换后的差分图像(输出图像)的,但显示单元5的显示区域较大时,也可以选择结合了上述全阴影变化显示模式和另外2种显示模式的显示形式。
另外,也可以设置可由用户改变在上述第3实施方式中说明的灰阶变换的特性的用户接口。即,例如,根据选择的特性对差分图像进行灰阶变换,使得可由用户从图5所示的特性和图11所示的特性中选择所希望的特性。
第4实施方式
在图12中,作为本发明的第4实施方式的医用图像处理装置110经由局域网(LAN)100与医用图像生成装置130、140、150及医用数据服务器160连接。医用图像生成装置130、140是例如CT扫描装置、MRI装置,医用图像生成装置150是例如X射线摄影装置。
用医用图像生成装置130、140、150生成的医用图像被直接传送给医用图像处理装置,或者暂时保存在医用数据服务器160中后,传送给医用图像处理装置110。直接传送给医用图像处理装置110的医用图像被保存在医用图像处理装置110的存储介质中。
医用图像处理装置110具有医用图像用高精细监控器120,可高精细地显示医用图像。另外,也可以使用具有更便宜的PC用监控器170的医用图像处理装置180。
在图13中,医用图像处理装置110在总线200上连接CPU210、RAM220、ROM230、通信接口240、输入装置260,医用图像用高精细监控器120和打印机250等的输出装置经由适当的接口连接在总线200上。输入装置包括键盘、指示装置等。
CPU210用于医用图像处理装置110整体及输出装置的控制,其控制程序存储在ROM230中。通信接口240控制由LAN100进行的通信,在适当医用图像生成装置130、140、150及医疗数据服务器160之间收发医用图像及其他数据。
图14是说明医用图像处理装置的功能结构的框图,图15是表示图14的医用图像处理装置中使用的查找表的特性的框图。图16是表示图14的医用图像处理装置中使用的查找表的组合的表。
在图14中,经由局域网100,将不同时期拍摄的人体的相同部位的一组图像的当前图像及过去图像的文件从医用数据服务器160输入到图像输入单元1中。在此,各图像文件由放射线图像数据和保持附带信息的标记信息构成,该附带信息中包含表示图像数据是正图像还是负图像的信息。
表示图像数据是正图像还是负图像的信息,在例如作为医用图像文件标准的DICOM格式的文件中,称为“Photometric Interpretation”的属性值在“MONOCHROME I”时定义为表示负图像,而在“MONOCHROME II”时定义为表示正图像。
图像输入单元1在向差分处理单元2输出图像数据,还抽取表示已输出给差分处理单元2的过去图像及当前图像是负图像还是正图像的信息,将判定信号S0输入到控制单元10中。
另一方面,接受诸如来自计算机的外部存储装置或用户的直接指示的无图示的外部输入,从阴影定义输入单元11输入差分图像上的阴影与诊断上的判断基准的对应关系。阴影定义规定为,对于图像诊断上病变部的加深或减轻,其在差分图像上与高浓度(高亮度值)或低浓度(低亮度值)之一相对应。
图18是表示阴影定义的例子的图,在该图中,类型A定义为,在差分图像上,将阴影的增大作为低浓度(黑色)的区域、将阴影的减少作为高浓度(白色)的区域来表现。另一方面,类型B定义为,在差分图像上,将阴影的增大作为高浓度(白色)的区域、将阴影的减少作为低浓度(黑色)的区域来表现。在此,低浓度或高浓度在以差分图像上未发生变化的部分的亮度值作为基准时,表示为比该基准值低的值或高的值。
差分处理单元2,具有从图像输入单元1分别输入过去及当前的图像的查找表LUT1、LUT2,在这些查找表LUT1、LUT2中进行必要的变换处理。查找表LUT1、LUT2的输出被输入到电平修正单元203中,转换一幅图像的浓度,使得过去或当前的图像的直方图中心与另一幅图像的直方图中心相一致。这样,如果使用电平修正后的图像执行位置重合,就能提高位置重合的精度。
电平修正单元203的输出被输入到位置重合单元20中,进行坐标变换,使得对应像素的坐标相一致。
在位置重合单元20中,为了对当前或过去的图像进行位置重合,使得解剖学上的对应点相一致,执行图18的流程图所示的以下步骤。
步骤S201:检测当前或过去的图像中的被检测体的位置,进入到步骤S202中。
步骤S202:在当前或过去的图像的大致对应的位置上设定多个关心区域ROI,进入到步骤S203中。
步骤S203:对每个ROI进行当前或过去的图像的匹配处理。在匹配处理中,计算出例如互相关系数,将匹配度最高的位置作为转换向量来求。
步骤S204:将在步骤S203中求得的转换向量进行2维内插,计算出用于使两图像的对应坐标相一致的坐标变换的参数。
步骤S205:基于在步骤S204中求得的参数,对当前或过去的图像进行坐标变换。
步骤S206:输出未进行坐标变换的图像和经过坐标变换后的图像。
下面,将与位置重合单元20的过去图像、当前图像对应的输出分别输入到查找表LUT3、查找表LUT4中,使之进行必要的变换。
查找表LUT3、查找表LUT4的输出被输入到相加运算单元202中,执行用于差分的相加运算。
在查找表LUT1至查找表LUT4上连接控制单元9。在控制单元10上连接图像输入单元1及阴影定义输入单元11。从图像输入单元1向控制单元10输入表示过去及当前图像是负图像还是正图像的判定信号S0。从阴影定义输入单元11向控制单元10输入表示判定方法是类型A还是类型B的信息。
查找表LUT1、查找表LUT2将两图像统一为负图像或者正图像,查找表LUT3、查找表LUT4将其中一幅图像进行符号反转,使得通过相加运算单元202中的相加运算来执行差分。
相加运算单元202生成差分图像,输出给图像输出单元5。
查找表LUT1至查找表LUT4可自由选择如图15A-图15C所示那样的3种特性。图15A将输入原样输出,是无变换的查找表(TypeI),图15B是将输入进行灰阶反转(正图像/负图像变换)的查找表(TypeII),图15C是将输入进行符号反转的查找表(TypeIII)。
在图16的表A中,在将判定方法设定为类型A时,并且过去图像是负图像,当前图像是正图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别设置为TypeI、TypeII、TypeI、TypeIII。
通过图19说明过去及当前图像的图像数据是如何变化的。如图19A所示,在该例中,正图像的当前图像中存在负图像的过去图像中所不存在的病灶阴影。判定方法为类型A时,由于查找表LUT2被设置为TypeII,所以当前图像进行灰阶反转,进行平均值的修正后成为图19B所示那样。
进而经位置重合进行相加运算时,在LUT4中对当前图像进行符号反转,因此当前图像的数据轮廓成为如图19C所示那样,通过相加运算计算差分。
即,从负图像的过去图像中减去由正图像灰阶反转为负图像的当前图像,因而,在差分图像上,阴影增大时作为低浓度(黑色)的区域显现,阴影减少时作为高浓度(白色)的区域显现。
如果判定方法被设定为类型A,过去图像和当前图像是负图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeI、TypeI、TypeI、TypeIII。因此,在查找表LUT1、查找表LUT2中,当前图像和过去图像都不进行灰阶反转,在查找表LUT3、查找表LUT4中,对当前图像进行图像信号的符号反转,所以,在相加单元202,中实质性地从负图像的过去图像中减去负图像的当前图像,其结果,当阴影增大时在差分图像上作为低浓度(黑色)的区域显现,阴影减少时在差分图像上作为高浓度(白色)的区域显现。
如果判定方法被设定为类型A,过去图像是正图像、当前图像是负图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeII、TypeI、TypeI、TypeIII。于是,过去图像从正图像灰阶反转为负图像后,减去当前图像。因此,在差分图像上,在阴影增大时作为低浓度(黑色)的区域显现,在阴影减少时作为高浓度(白色)的区域显现。
如果判定方法被设定为类型A,过去图像和当前图像都是正图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeI、TypeI、TypeIII、TypeI。于是,从正图像的当前图像中减去正图像的过去图像,所以,在差分图像上,阴影增大时作为低浓度(黑色)的区域显现,阴影减少时作为高浓度(白色)的区域显现。
如果过去图像和当前图像都是负图像,通过从过去图像中减去当前图像来形成差分图像。如果过去图像和当前图像都是正图像,通过从当前图像中减去过去图像来形成差分图像。所以,不论哪种情况,在差分图像上,都是在阴影增大时作为低浓度(黑色)的区域显现,而在阴影减少时作为高浓度(白色)的区域显现。
在图16的表B中,如果判定方法被设定为类型B,过去图像和当前图像都是负图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeI、TypeI、TypeIII、TypeI。因此,从负图像的当前图像中减去负图像的过去图像,在差分图像上,当阴影增大时作为高浓度(白色)的区域显现,当阴影减少时作为低浓度(黑色)的区域显现。
如果判定方法被设定为类型B,过去图像是正图像、当前图像是负图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeI、TypeII、TypeI、TypeIII。因而,从正图像的过去图像中减去由负图像灰阶反转为正图像的当前图像,所以,在差分图像上,当阴影增大时作为高浓度(白色)的区域显现,当阴影减少时作为低浓度(黑色)的区域显现。
如果判定方法被设定为类型B,过去图像是负图像、当前图像是正图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeII、TypeI、TypeI、TypeIII。据此,从由负图像灰阶反转为正图像的过去图像中减去正图像的当前图像,在差分图像上,当阴影增大时作为高浓度(白色)的区域显现,当阴影减少时作为低浓度(黑色)的区域显现。
如果判定方法被设定为类型B,过去图像和当前图像都是正图像时,查找表LUT1、查找表LUT2、查找表LUT3、查找表LUT4分别被设置为TypeI、TypeI、TypeI、TypeIII。因此,从正图像的过去图像中减去正图像的当前图像,在差分图像上,当阴影增大时作为高浓度(白色)的区域显现,当阴影减少时作为低浓度(黑色)的区域显现。
过去图像和当前图像都是负图像时,通过从当前图像中减去过去图像,形成差分图像,过去图像和当前图像都是正图像时,通过从过去图像中减去当前图像,在差分图像上,当阴影增大时作为高浓度(白色)的区域显现,当阴影减少时作为低浓度(黑色)的区域显现。
如以上说明,通过本实施方式,不论当前和过去图像的属性如何组合,阴影的增加和减少都能以预先确定的1种形式来表现。
<第5实施方式>
在可显示DICOM格式的图像的显示装置中,对于PhotometricInterpretation的值为“MONOCHROME I”的图像数据,显示灰阶反转,而对于Photometric Interpretation的值为“MONOCHROME II”的图像数据,不进行灰阶反转,就原样进行显示。
在第4实施方式所示的医用图像处理装置中,在过去图像和当前图像具有相同的图像特性时,通过改变进行图像信号的符号反转的图像来改变相减运算的顺序。但在本实施方式中,固定相减运算的顺序,通过对差分图像改变提供的Photometric Interpretation的值而取得同样的结果。
图21是说明作为本发明的第5实施方式的医用图像处理装置的功能结构的框图。因为医用图像处理装置的结构与图12及图13所示相同,故省略其说明。
在图18中,关于与第4实施方式相同的部分使用相同的符号,省略说明。差分处理单元2所包含的查找表LUT3、查找表LUT4设定不变的特性,而不通过控制单元10改变特性。标记生成单元9生成从差分处理单元2输出的差分图像数据所附带的附带信息。在该附带信息中包含表示从差分处理单元2输出了的差分图像数据是正图像还是负图像的标记信息。
在标记生成单元9中,根据从图像输入单元1输入的过去图像数据的附带信息和当前图像的附带信息以及与被输入到阴影定义输入单元11的阴影的定义的组合,生成从差分处理单元2输出的差分图像数据所附带的附带信息。因此,从差分处理单元2输出的差分图像,作为附加了由标记生成单元9生成的附带信息的图像文件,而被输出给图像输出单元5。
在第5实施方式中,说明将查找表LUT3固定为TypeI、将查找表LUT4固定为TypeIII时的例子。
如果判定方法被设定为类型A,过去图像和当前图像都是负图像时,通过控制单元10将查找表LUT1和查找表LUT2都设定为TypeI。因此,在查找表LUT1和查找表LUT2中,当前图像和过去图像都不进行灰阶反转。通过被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,因此在相加运算单元202中,实质性地从负图像的过去图像中减去负图像的当前图像。控制单元3控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时在差分图像上作为低浓度(黑色)的区域显现,阴影减少时在差分图像上作为高浓度(白色)的区域显现。
如果判定方法被设定为类型A,过去图像是负图像、当前图像是正图像时,查找表LUT1、查找表LUT2分别被设置为TypeI、TypeII。因为查找表LUT2被设置为TypeII,所以将当前图像进行灰阶反转。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以,在相加运算单元202中,实质性地从负图像的过去图像中减去负图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时在差分图像上作为低浓度(黑色)的区域显现,阴影减少时在差分图像上作为高浓度(白色)的区域显现。
如果判定方法被设定为类型A,过去图像是正图像、当前图像是负图像时,查找表LUT1、查找表LUT2分别被设置为TypeII、TypeI。因为查找表LUT1被设置为TypeII,所以将过去图像进行灰阶反转。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以,在相加运算单元202中,实质性地从负图像的过去图像中减去负图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时在差分图像上作为低浓度(黑色)的区域显现,阴影减少时在差分图像上作为高浓度(白色)的区域显现。
另一方面,如果判定方法被设定为类型B,过去图像是正图像、当前图像是负图像时,查找表LUT1、查找表LUT2分别被设置为TypeI、TypeII。因为查找表LUT2被设置为TypeII,所以将当前图像进行灰阶反转。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以,在相加运算单元202中实质性地从正图像的过去图像中减去正图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时作为高浓度(白色)的区域显现,阴影减少时作为低浓度(黑色)的区域显现。
如果判定方法被设定为类型B,过去图像是负图像、当前图像是正图像时,查找表LUT1、查找表LUT2分别被设置为TypeII、TypeI。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以,从由负图像灰阶反转为正图像的过去图像中减去正图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时作为高浓度(白色)的区域显现,阴影减少时作为低浓度(黑色)的区域显现。
同样地,如果判定方法被设定为类型B,过去图像和当前图像都是正图像时,查找表LUT1、查找表LUT2都被设置为TypeI。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以从正图像的过去图像中减去正图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME II”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时作为高浓度(白色)的区域显现,阴影减少时作为低浓度(黑色)的区域显现。
如果判定方法被设定为类型A,过去图像和当前图像都是正图像时,查找表LUT1、查找表LUT2都被设置为TypeI。利用被固定为TypeI的查找表LUT3和被固定为TypeIII的查找表LUT4,对当前图像进行图像信号的符号反转,所以,在相加运算单元202中,实质性地从正图像的过去图像中减去正图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME I”给予Photometric Interpretation的值。
通过将“MONOCHROME I”给予Photometric Interpretation的值,在可显示DICOM格式的图像的显示装置中显示经灰阶反转的差分图像。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时在差分图像上作为低浓度(黑色)的区域显现,阴影减少时在差分图像上作为高浓度(白色)的区域显现。
如果判定方法被设定为类型B,过去图像和当前图像都是负图像时,因为查找表LUT1、查找表LUT2都被设置为TypeI,所以从负图像的过去图像中减去负图像的当前图像。控制单元10控制标记生成单元9,使得对于该差分图像,将“MONOCHROME I”给予Photometric Interpretation的值。
因此,如果在可显示DICOM格式的图像的显示装置中显示该差分图像,则阴影增大时在差分图像上作为高浓度(白色)的区域显现,阴影减少时在差分图像上作为低浓度(黑色)的区域显现。
在本实施方式中,说明了用固定的查找表LUT3和LUT4对位置重合单元201的输出进行符号变换,用相加运算单元202对变换结果进行相加运算的方法。但如图21、20所示那样,设置一个对位置重合单元201的输出,进行从过去图像中减去当前图像的处理的相减运算单元205,也是一样的。
这样,在本实施例中,固定相减运算的顺序,基于从阴影定义输入单元11输入的阴影定义和过去图像及当前图像的正图像/负图像信息,改变给予差分图像的DICOM格式的Photometric Interpretation的值,由此,即使是在正图像和负图像混合的图像间的差分图像中,也可以统一阴影的含义。
<其他实施方式>
另外,本发明可适用于以上实施方式的图像处理装置之外的装置,例如可适用于由多个设备构成的***、或者由一个设备构成的装置。
而且,通过任意的装置、***来实现上述实施方式的功能的方法都构成本发明,使普通计算机执行这样的方法的软件的程序也构成本发明。
并且,当然,将存储了上述软件的程序代码的存储介质提供给***或者装置,由该***或者装置的计算机(或CPU、MPU)读出并执行存储在存储介质中的程序代码也能达成本发明。
此时,从存储介质中读出的程序代码自身实现上述实施方式的功能,存储了该程序代码的存储介质构成本发明。
作为用于提供程序代码的存储介质,可以使用例如软盘、硬盘、光盘、光磁盘、CD-ROM、CD-R、磁带、非易失性的存储卡、ROM等。
另外,不仅由计算机执行读出的程序代码来实现上述实施方式的功能。当然,基于该程序代码的指示,由在计算机上运行的OS等进行实际处理的一部分或全部,通过该处理来实现上述实施方式的功能的情况也包含在内。
而且,当然,从存储介质中读出的程序代码在写入插在计算机中的功能扩展板或连接在计算机上的功能扩展单元所具有的存储器中后,基于该程序代码的指示,由该功能扩展板或功能扩展单元所具有的CPU等进行实际处理的一部分或全部,通过该处理来实现上述实施方式的功能的情况也包含在内。
通过本发明,对于从图像的正图像/负图像不统一的图像间生成的差分图像,可以生成具有一贯性的差分图像,可以输出差分图像上的阴影的显示形式具有一贯性的差分图像。
本发明也适用于一种程序或存储了该程序的存储介质。
当然,本发明也可以适用于包括多个装置(例如,射线发生装置、射线照相装置、图像处理装置、接口装置等)的***和具有上述装置的功能的单个装置。当本发明适用于包括多个装置的***时,上述装置经由诸如电学装置、光学装置、机械装置和/或其他装置和另一个连接。
而且,本发明还适用于包括网络(LAN和/或WAN等)的诊断援助***。
如上所述,本发明因此可以实现上述的目的。
本发明并不限于上述实施形式,在不脱离本发明的精神和范围内可以进行各种变更和修正。因此,为了公开本发明的范围,附加以下的权利要求。
Claims (9)
1.一种图像处理方法,其特征在于,包括:
输入步骤,输入在不同时间点拍摄的至少2幅医用图像;
差分图像生成步骤,在由上述输入步骤输入的2幅医用图像间进行位置重合,差分处理对应的坐标点的图像信号并生成差分图像;
差分图像存储步骤,将由上述差分图像生成步骤生成的差分图像存储到存储介质中;
基准差分值确定步骤,作为灰阶变换处理的基准值确定由上述差分图像存储步骤所存储的差分图像的基准差分值;
状态选择步骤,使用户选择用作由上述差分图像生成步骤所生成的差分图像的读影的指标的阴影变化的状态;
差分图像处理步骤,根据由上述基准差分值确定步骤确定的基准差分值和由上述状态选择步骤选择的阴影变化的状态,来对上述差分图像进行灰阶变换处理;以及
输出步骤,输出由上述差分图像处理步骤处理过的差分图像。
2.根据权利要求1所述的图像处理方法,其特征在于,
在上述差分基准值确定步骤中,根据上述差分图像的差分值的直方图,来确定上述基准差分值。
3.根据权利要求1所述的图像处理方法,其特征在于,
在上述差分图像处理步骤中,将处于未由上述状态选择步骤选择的状态的区域的差分值变换成上述基准差分值。
4.根据权利要求1所述的图像处理方法,其特征在于,
在上述差分图像处理步骤中,对上述差分图像进行灰阶变换,使得处于未由上述状态选择步骤选择的状态的区域成为相对较低的对比度。
5.根据权利要求1所述的图像处理方法,其特征在于,
还包括:
2值图像生成步骤,对由上述差分图像处理步骤处理过的差分图像实施2值化并生成2值图像;以及
计数步骤,在由上述2值图像生成步骤生成的2值图像中,计算具有预定面积的区域的数量,
其中,在上述输出步骤中,将上述计数步骤的处理结果与上述差分图像一同输出。
6.根据权利要求1所述的图像处理方法,其特征在于,
使用户通过上述状态选择步骤选择的阴影变化的状态是全阴影显示状态、阴影增大状态、阴影消失状态中的至少一种状态。
7.根据权利要求1所述的图像处理方法,其特征在于,
上述差分图像生成步骤包括:
控制步骤,基于第1图像和第2图像的图像属性决定灰阶处理方法,使得用预定的显示方法显示阴影变化区域;
灰阶处理步骤,按照上述控制步骤的决定,对上述第1图像和/或上述第2图像进行预定的灰阶处理;以及
运算步骤,从经上述灰阶处理步骤处理过的上述第1图像和上述第2图像计算出差分图像。
8.一种图像处理装置,其特征在于,包括:
输入装置,输入在不同时间点拍摄的至少2幅医用图像;
差分图像生成装置,在由上述输入装置输入的2幅医用图像间进行位置重合,差分处理对应的坐标点的图像信号并生成差分图像;
差分图像存储装置,将由上述差分图像生成装置生成的差分图像存储到存储介质中;
基准差分值确定装置,作为灰阶变换处理的基准值确定由上述差分图像存储装置所存储的差分图像的基准差分值;
状态选择装置,使用户选择用作由上述差分图像生成装置所生成的差分图像的读影的指标的阴影变化的状态;
差分图像处理装置,根据由上述基准差分值确定装置确定的基准差分值和由上述状态选择装置选择的阴影变化的状态,来对上述差分图像进行灰阶变换处理;以及
输出装置,输出由上述差分图像处理装置处理过的差分图像。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002309841A JP3793137B2 (ja) | 2002-10-24 | 2002-10-24 | 画像処理装置、画像処理方法、プログラム、および記憶媒体 |
JP309841/2002 | 2002-10-24 | ||
JP134022/2003 | 2003-05-13 | ||
JP2003134022A JP4447850B2 (ja) | 2003-05-13 | 2003-05-13 | 画像処理装置、画像処理方法、及びコンピュータプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1496714A CN1496714A (zh) | 2004-05-19 |
CN1259890C true CN1259890C (zh) | 2006-06-21 |
Family
ID=32109502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101018998A Expired - Fee Related CN1259890C (zh) | 2002-10-24 | 2003-10-23 | 图像处理装置和方法 |
Country Status (2)
Country | Link |
---|---|
US (2) | US7403646B2 (zh) |
CN (1) | CN1259890C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678477A (zh) * | 2012-09-25 | 2014-03-26 | 波音公司 | 全景可视化文档区别 |
Families Citing this family (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8265728B2 (en) * | 2003-11-26 | 2012-09-11 | University Of Chicago | Automated method and system for the evaluation of disease and registration accuracy in the subtraction of temporally sequential medical images |
JP2005253685A (ja) * | 2004-03-11 | 2005-09-22 | Konica Minolta Medical & Graphic Inc | 画像診断支援装置及び画像診断支援プログラム |
JP4537143B2 (ja) * | 2004-07-30 | 2010-09-01 | キヤノン株式会社 | 画像処理装置及びその方法、撮像装置、プログラム |
US9563980B2 (en) * | 2005-11-18 | 2017-02-07 | Autodesk, Inc. | Grip manipulatable shadows in 3D models |
US7769216B2 (en) * | 2005-12-29 | 2010-08-03 | Hologic, Inc. | Facilitating comparison of medical images |
US7630531B2 (en) * | 2006-01-31 | 2009-12-08 | Mevis Medical Solutions, Inc. | Enhanced navigational tools for comparing medical images |
US8425418B2 (en) * | 2006-05-18 | 2013-04-23 | Eigen, Llc | Method of ultrasonic imaging and biopsy of the prostate |
JP2008000536A (ja) * | 2006-06-26 | 2008-01-10 | Fujifilm Corp | 画像表示装置 |
US8064664B2 (en) * | 2006-10-18 | 2011-11-22 | Eigen, Inc. | Alignment method for registering medical images |
US7804989B2 (en) * | 2006-10-30 | 2010-09-28 | Eigen, Inc. | Object recognition system for medical imaging |
US20080161687A1 (en) * | 2006-12-29 | 2008-07-03 | Suri Jasjit S | Repeat biopsy system |
US8175350B2 (en) | 2007-01-15 | 2012-05-08 | Eigen, Inc. | Method for tissue culture extraction |
US7856130B2 (en) * | 2007-03-28 | 2010-12-21 | Eigen, Inc. | Object recognition system for medical imaging |
JP2009082227A (ja) * | 2007-09-27 | 2009-04-23 | Fujifilm Corp | 医用レポートシステム、医用レポート作成装置、及び医用レポート作成方法 |
JP2009086765A (ja) * | 2007-09-27 | 2009-04-23 | Fujifilm Corp | 医用レポートシステム、医用レポート作成装置、及び医用レポート作成方法 |
JP4528322B2 (ja) * | 2007-09-28 | 2010-08-18 | 富士フイルム株式会社 | 画像表示装置、画像表示方法、および画像表示プログラム |
US8571277B2 (en) | 2007-10-18 | 2013-10-29 | Eigen, Llc | Image interpolation for medical imaging |
US7942829B2 (en) * | 2007-11-06 | 2011-05-17 | Eigen, Inc. | Biopsy planning and display apparatus |
JP2010057902A (ja) * | 2008-08-06 | 2010-03-18 | Toshiba Corp | レポート作成支援装置、レポート作成支援システム、及び医用画像参照装置 |
JP5318646B2 (ja) * | 2009-04-17 | 2013-10-16 | ルネサスエレクトロニクス株式会社 | 差分画像生成装置、差分画像生成方法、及びプログラム |
FR2946171B1 (fr) * | 2009-05-29 | 2011-07-15 | Groupe Des Ecoles De Telecommunications Get Ecole Nationale Superieure Des Telecommunications Enst | Procede de quantification de l'evolution de pathologies impliquant des changements de volumes de corps, notamment de tumeurs |
JP5576631B2 (ja) * | 2009-09-09 | 2014-08-20 | キヤノン株式会社 | 放射線撮影装置、放射線撮影方法、及びプログラム |
JP5610761B2 (ja) * | 2009-12-16 | 2014-10-22 | キヤノン株式会社 | X線画像処理装置、x線画像処理システム、x線画像処理方法、及びコンピュータプログラム |
US8675933B2 (en) | 2010-04-30 | 2014-03-18 | Vucomp, Inc. | Breast segmentation in radiographic images |
US9076197B2 (en) | 2010-04-30 | 2015-07-07 | Vucomp, Inc. | Probability density function estimation |
WO2012006318A1 (en) | 2010-07-07 | 2012-01-12 | Vucomp, Inc. | Marking system for computer-aided detection of breast abnormalities |
WO2012046842A1 (ja) * | 2010-10-07 | 2012-04-12 | 株式会社東芝 | 読影依頼システム、読影依頼仲介装置、制御プログラム及び読影結果評価方法 |
US8792748B2 (en) * | 2010-10-12 | 2014-07-29 | International Business Machines Corporation | Deconvolution of digital images |
JP6016403B2 (ja) * | 2012-03-27 | 2016-10-26 | キヤノン株式会社 | 画像処理装置、画像処理方法 |
US10268662B2 (en) * | 2012-09-10 | 2019-04-23 | The Boeing Company | Panoptic visualization of a document according to the structure thereof |
EP2720192B1 (en) * | 2012-10-12 | 2015-04-08 | General Electric Company | Method, system and computer readable medium for liver diagnosis |
US20150301323A1 (en) * | 2012-11-15 | 2015-10-22 | Shimadzu Corporation | System for setting analysis target region |
JP6318739B2 (ja) * | 2014-03-17 | 2018-05-09 | コニカミノルタ株式会社 | 画像処理装置、およびプログラム |
WO2016174926A1 (ja) * | 2015-04-30 | 2016-11-03 | 富士フイルム株式会社 | 画像処理装置及び画像処理方法及びプログラム |
US10716544B2 (en) | 2015-10-08 | 2020-07-21 | Zmk Medical Technologies Inc. | System for 3D multi-parametric ultrasound imaging |
KR101798083B1 (ko) | 2016-02-04 | 2017-12-12 | 삼성전자주식회사 | 단층 영상 처리 장치, 방법, 및 그 방법에 관련된 기록매체 |
AT16426U1 (de) * | 2016-03-02 | 2019-08-15 | Mathias Zirm Dr | Ferndiagnose-Unterstützungsverfahren |
JP6656207B2 (ja) * | 2017-05-31 | 2020-03-04 | キヤノン株式会社 | 情報処理装置、その制御方法、及びプログラム |
JP7119974B2 (ja) * | 2018-12-14 | 2022-08-17 | コニカミノルタ株式会社 | 医用画像表示装置及び医用画像表示システム |
US11635802B2 (en) * | 2020-01-13 | 2023-04-25 | Sony Interactive Entertainment Inc. | Combined light intensity based CMOS and event detection sensor for high speed predictive tracking and latency compensation in virtual and augmented reality HMD systems |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5289374A (en) * | 1992-02-28 | 1994-02-22 | Arch Development Corporation | Method and system for analysis of false positives produced by an automated scheme for the detection of lung nodules in digital chest radiographs |
US5359513A (en) | 1992-11-25 | 1994-10-25 | Arch Development Corporation | Method and system for detection of interval change in temporally sequential chest images |
JP3002937B2 (ja) * | 1993-09-13 | 2000-01-24 | 富士写真フイルム株式会社 | エネルギーサブトラクション画像処理方法 |
JPH07299053A (ja) * | 1994-04-29 | 1995-11-14 | Arch Dev Corp | コンピュータ診断支援方法 |
US5982953A (en) * | 1994-09-02 | 1999-11-09 | Konica Corporation | Image displaying apparatus of a processed image from temporally sequential images |
US5999652A (en) * | 1995-05-25 | 1999-12-07 | Lockheed Martin Corporation | Plume or combustion detection by time sequence differentiation of images over a selected time interval |
DE19613342A1 (de) | 1996-04-03 | 1997-10-09 | Philips Patentverwaltung | Automatisches Bildauswertungsverfahren |
US5987345A (en) * | 1996-11-29 | 1999-11-16 | Arch Development Corporation | Method and system for displaying medical images |
US5982927A (en) * | 1996-12-19 | 1999-11-09 | Cognex Corporation | Methods and apparatuses for in-line solder paste inspection |
US6061476A (en) * | 1997-11-24 | 2000-05-09 | Cognex Corporation | Method and apparatus using image subtraction and dynamic thresholding |
DE69833128T2 (de) * | 1997-12-10 | 2006-08-24 | Koninklijke Philips Electronics N.V. | Bildung eines zusammengesetzten bildes aus aufeinanderfolgenden röntgenbildern |
US6118906A (en) * | 1998-02-03 | 2000-09-12 | Eastman Kodak Company | Sharpening system adjusted for measured noise of photofinishing images |
FR2777374B1 (fr) * | 1998-04-10 | 2000-05-12 | Commissariat Energie Atomique | Procede de recalage de deux images differentes d'un meme objet |
US6240201B1 (en) * | 1998-07-24 | 2001-05-29 | Arch Development Corporation | Computerized detection of lung nodules using energy-subtracted soft-tissue and standard chest images |
JP2000048036A (ja) | 1998-07-28 | 2000-02-18 | Canon Inc | 画像処理装置およびその方法 |
JP2001299733A (ja) * | 2000-04-27 | 2001-10-30 | Konica Corp | Pci放射線画像処理装置、pci放射線画像検出処理装置、pci放射線画像出力装置及びpci画像診断支援装置 |
JP3792529B2 (ja) | 2001-04-17 | 2006-07-05 | シロキ工業株式会社 | ワイヤ式ウィンドレギュレータ |
-
2003
- 2003-10-03 US US10/677,377 patent/US7403646B2/en not_active Expired - Fee Related
- 2003-10-23 CN CNB2003101018998A patent/CN1259890C/zh not_active Expired - Fee Related
-
2007
- 2007-01-16 US US11/623,565 patent/US7447346B2/en not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678477A (zh) * | 2012-09-25 | 2014-03-26 | 波音公司 | 全景可视化文档区别 |
Also Published As
Publication number | Publication date |
---|---|
US7403646B2 (en) | 2008-07-22 |
CN1496714A (zh) | 2004-05-19 |
US20070110307A1 (en) | 2007-05-17 |
US20040081342A1 (en) | 2004-04-29 |
US7447346B2 (en) | 2008-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1259890C (zh) | 图像处理装置和方法 | |
CN1264111C (zh) | 医疗用图像处理***和方法 | |
CN101032159A (zh) | 图像处理设备及其方法,图像处理程序 | |
CN1284073C (zh) | 信息显示***及其信息处理装置、指示装置和标记显示法 | |
CN100341045C (zh) | 图像显示装置和图像处理方法 | |
CN100350431C (zh) | 增强以批模式处理的肖像图像的方法和*** | |
CN1191549C (zh) | 用于医学图象的迭代滤波器结构 | |
CN1284063C (zh) | 信息处理方法和信息处理装置 | |
CN101053521A (zh) | 医用图像显示装置 | |
CN1228965C (zh) | 图像处理***、投影仪及黑白扩展处理法 | |
CN1909671A (zh) | 产生用户偏爱数据的方法和装置 | |
CN1639725A (zh) | 图像像素编码方法、图像处理方法和针对一个或多个图像像素再现的对象进行定性识别的图像处理方法 | |
CN1630886A (zh) | 三维图象的比较程序、比较方法及比较装置 | |
CN1738441A (zh) | 使用屏幕的一部分的边的梯形失真校正 | |
CN1671176A (zh) | 修正图像变形的图像处理装置、修正摄影图像变形的摄影装置 | |
CN1294530C (zh) | 彩色图象处理方法、彩色图象处理器、彩色显示器、以及用于实施彩色图象处理方法的计算机程序 | |
CN1892800A (zh) | 图像显示装置和图像显示方法 | |
CN1691740A (zh) | 放大显示装置和放大图像控制装置 | |
CN1734444A (zh) | 图像处理装置及其控制方法 | |
CN1601369A (zh) | 图像处理***、投影机以及图像处理方法 | |
CN1444907A (zh) | 图像处理装置和超声波诊断装置 | |
CN1592915A (zh) | 用户可定义的图像基准点 | |
CN1735135A (zh) | 分辨率转换设备和方法 | |
CN1400807A (zh) | 图象处理方法和设备,图象处理***以及存储介质 | |
CN1167036C (zh) | 图象处理方法和记录了图象处理程序的记录媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20060621 Termination date: 20191023 |