CN101801273B - 放射线图像处理装置以及放射线图像处理方法 - Google Patents
放射线图像处理装置以及放射线图像处理方法 Download PDFInfo
- Publication number
- CN101801273B CN101801273B CN2007801005937A CN200780100593A CN101801273B CN 101801273 B CN101801273 B CN 101801273B CN 2007801005937 A CN2007801005937 A CN 2007801005937A CN 200780100593 A CN200780100593 A CN 200780100593A CN 101801273 B CN101801273 B CN 101801273B
- Authority
- CN
- China
- Prior art keywords
- grid
- mentioned
- image processing
- images
- image
- 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.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 102
- 238000003672 processing method Methods 0.000 title claims description 6
- 238000012545 processing Methods 0.000 claims abstract description 65
- 239000000203 mixture Substances 0.000 claims description 52
- 238000000034 method Methods 0.000 claims description 20
- 238000001914 filtration Methods 0.000 claims description 11
- 230000000737 periodic effect Effects 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 4
- 238000000926 separation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 208000031481 Pathologic Constriction Diseases 0.000 description 4
- 210000001215 vagina Anatomy 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 238000002835 absorbance Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000218691 Cupressaceae Species 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000000586 desensitisation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/06—Diaphragms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4035—Arrangements for generating radiation specially adapted for radiation diagnosis the source being combined with a filter or grating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5252—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明的放射线图像处理装置,具备分离单元、调整单元和除去单元。由于采用通过被检测体进行摄影所得到的放射线图像来生成补正图像,因此不必预先获取用于格栅补正的数据而可以由一次摄影来获取各图像。而且,将通过被检测体进行摄影所得到的放射线图像分离成格栅图像和非格栅图像,并基于非格栅图像调整在实际空间上的格栅图像的强度,并生成补正图像,因此以一次摄影可以实现高精度的补正处理。
Description
技术领域
本发明涉及对通过用于除去散射放射线的格栅(grid)而摄影的放射线图像进行处理的放射线图像处理装置以及放射线图像处理程序。
背景技术
放射线图像处理装置将透过被检测体的放射线能量蓄积到被称为“平板型放射线检测器(FPD:Flat Panel Detector)”的图像传感器上,并转换成电信号而进行可视化。此时,FPD蓄积透过了被检测体的直线的放射线能量,并且还蓄积透过被检测体时散射的能量低的放射线。所以,因散射放射线而导致所摄影的放射线图像的对比度降低等,使图像质量降低。通常,为了除去这些散射放射线,通过交替地配置使用放射线吸收率高的铅等的放射线吸收层和使用放射线吸收率低的铝或纸等的放射线透过层而构成的格栅来进行摄影。
但是,通过格栅摄影被检测体的情况下,由放射线吸收层的影响而花纹模样的像(格栅像)会重叠到作为摄影图像的放射线图像上。以往,为了减轻该格栅像,有一种摇动系统中一边摇动格栅一边进行摄影。但是,实现该摇动系统时,存在装置的大型化,成本提高等问题。
因此,为了除去具有由格栅引起的偏差的格栅像,有一种方法中基于构成图像传感器的多个像素的灵敏度偏差,根据不通过被检测体而通过格栅预先摄影的图像数据,来创建用于补正被检测体的放射线图像的补正用图像,并采用该补正图像来除去由格栅引起的图像成分(例如参照专利文献1)。除此以外,为了对应于多种类型的格栅或各种各样的图像传感器,提出了基于格栅相关的格栅频率与图像传感器的抽样频率的关系,进行滤波处理的方法(例如参照专利文献2-4)等。
专利文献1:日本特开2005-324069号公报
专利文献2:日本特许第2507659号公报(第1、3-5页,图1-3、7)
专利文献3:日本特开2004-344670号公报(第7页,图7,表1)
专利文献4:日本特开2003-260053号公报(第6页,表1)
但是,在上述的专利文献1的情况下,必须如上所述那样不通过被检测体而通过格栅对成为用于除去格栅像的基准的补正用图像进行预先摄影。实际上,每次使用时预先获取不通过被检测体的补正用图像是烦杂的,在实用上是不现实的。而且,为了得到精度高的补正结果,必须要使获取补正用图像时、以及实际摄影被检测体时的摄影环境或条件等完全一致,但在实际的使用状况下要使摄影环境或条件完全一致是在很多情况下实用上有困难的。
而且,在上述的专利文献1的情况下,假定为格栅纹在整个图像中稳定,通过空间滤波从作为补正对象的对象图像中提取格栅纹成分。实际上,由于格栅内的放射线吸收层的倾斜随图像位置而不同,有时图像中央部和端部格栅纹的频带是不同的,并不一定是稳定的。从而,由于采用实际上频带不稳定的补正用图像,在假定频带在整个图像中稳定的前提下进行补正,因此其补正结果的精度必然很低。而且,通过预先生成补正用图像,可以降低灵敏度的偏差,但是由于没有考虑创建补正图像时的补正本身的精度的偏差,因此并未必能实现高精度的补正。
发明内容
本发明,鉴于这样的事情,其目的是提供一种放射线图像处理装置以及放射线图像处理程序,不用预先获取用于格栅补正的数据而以一次摄影来可以实现精度高的补正处理。
本发明,为了达成这样的目的,采取以下的构成。
即,本发明的放射线图像处理装置,对通过用于除去散射放射线的格栅进行摄影所得到的放射线图像进行处理,上述放射线图像处理装置的特征在于,具备:分离单元,其将上述的放射线图像分离成包含上述格栅的格栅像成分的格栅图像和包含其它成分的非格栅图像;调整单元,其基于上述非格栅图像,对在实际空间上的上述格栅图像的强度进行调整,并生成调整图像;和除去单元,其从上述非格栅图像中减去上述调整图像,并生成除去格栅影响后所得到的补正图像。
根据本发明的放射线图像处理装置,分离单元将放射线图像分离成包含格栅的格栅像成分的格栅图像和包含其它成分的非格栅图像,而调整单元基于非格栅图像、对在实际空间上的格栅图像的强度进行调整,并生成调整图像。然后,除去单元从非格栅图像减去调整图像,并生成除去格栅影响后所得到的补正图像。这样,由于采用通过被检测体进行摄影所得到的放射线图像来生成补正图像,因此可以不通过被检测体,并且不必通过格栅对成为除去格栅像的基准的补正用图像进行预先摄影,而以一次摄影就可以获取各图像,没有必要使摄影环境或条件等的完全一致,进而不受限制。而且,将通过被检测体进行摄影所得到的放射线图像分离成格栅图像和非格栅图像,并基于非格栅图像对在实际空间上的格栅图像的强度进行调整,并生成补正图像,因此利用基于调整的格栅图像的各像素中的强度来考虑了各像素的偏差,可以实现高精度的补正处理。其结果,不必预先获取用于格栅补正的数据而以一次摄影即可实现高精度的补正处理。
上述的分离单元的一例,是在空间频域上分离格栅图像和非格栅图像。由于在空间频域上分离,因此可以更可靠地分离格栅图像成分和其它成分。当然,分离单元的其它一例,也可以在实际空间上分离格栅图像和非格栅图像。
而且,分离单元,提取放射线图像的周期性的图形(pattern),由此分离图像成分。特别是,格栅像的成分,以具有某种特定范围的频带的周期性的图形为多,因此,采用使相当于格栅像的特定的频带通过的带通滤波处理,来分离图像成分是理想的。而且,利用格栅像重叠的放射线图像在实际空间上成为周期性的图形这一点,分离单元也可以对于与放射线图像在实际空间上的强度有关的分布图,提取出与三角函数类似的成分,由此分离图像成分。
而且,当基于非格栅图像对在实际空间上的格栅图像的强度进行调整之际,按照使从非格栅图像减去调整图像后所得到的图像中不包含与格栅图像相同的成分的方式,调整单元调整格栅图像的强度,并生成调整图像。具体地,当非格栅图像的格栅位置的局部的强度变化不为单调时,按照使格栅图像的格栅位置的强度等于非格栅图像的格栅位置的局部的强度变化的方式,调整单元调整格栅图像的强度,并生成调整图像。或者,按照使格栅图像的格栅位置的强度等于非格栅图像的格栅位置的强度与对非格栅图像平滑化后所得到的图像的格栅位置的强度之差的方式,调整单元调整格栅图像的强度,并生成调整图像。通过这样的调整,在从非格栅图像减去调整图像后所得到的图像(即补正图像)中不包含与格栅图像相同的成分,可以实现更高精度的补正。
而且,本发明的放射线图像处理方法,用于执行通过用于除去散射放射线的格栅进行摄影所得到的放射线图像的处理,上述放射线图像处理方法的特征在于,包括下述步骤:分离步骤,将上述的放射线图像分离成包含上述格栅的格栅像成分的格栅图像和包含其它成分的非格栅图像;调整步骤,基于上述非格栅图像,对在实际空间上的上述格栅图像的强度进行调整,并生成调整图像;和除去步骤,从上述非格栅图像中减去上述调整图像,并生成除去格栅影响后所得到的补正图像。
根据本发明的放射线图像处理程序,由于采用通过被检测体进行摄影所得到的放射线图像来生成补正图像,不必预先获取用于格栅补正的数据,而以一次摄影就可以获取各图像。而且,将通过被检测体进行摄影所得到的放射线图像分离成格栅图像和非格栅图像,并基于非格栅图像对在实际空间上的格栅图像的强度进行调整,并生成补正图像,因此以一次摄影就可以实现高精度的补正处理。
发明效果
根据本发明的放射线图像处理装置以及放射线图像处理程序,由于采用通过被检测体进行摄影所得到的放射线图像来生成补正图像,不必预先获取用于格栅补正的数据,而以一次摄影就可以获取各图像。而且,将通过被检测体进行摄影所得到的放射线图像分离成格栅图像和非格栅图像,并基于非格栅图像对在实际空间上的格栅图像的强度进行调整,并生成补正图像,因此以一次摄影就可以实现高精度的补正处理。
附图说明
图1是实施例的放射线图像处理装置的方框图。
图2是平板型放射线检测器(FPD)的检测面的示意图。
图3(a)~(c)是分离的说明中所采用的各图像的概要图。
图4是表示一系列调整的流程的流程图。
图5(a)~(c)是调整的说明中的各分布图(profile)的概要图。
图中:
4-图像处理部,31-格栅,41-分离部,42-调整部,43-除去部。
具体实施方式
以下,参照附图对本发明的实施例进行说明。图1是实施例的放射线图像处理装置的方框图,图2是平板型放射线检测器(FPD)的检测面的示意图,图3是分离的说明中所采用的各图像的概要图,图4是表示一系列调整流程的流程图,图5是调整的说明中所采用的各分布图的概要图。
本实施例的放射线图像处理装置,如图1所示,具备:装载被检测体M的顶板1;朝向被检测体M并照射放射线(例如X射线)的放射线源2(例如X射线管);检测从放射线源2照射并透过被检测体M的放射线的平板型放射线检测器(以下略记为“FPD”)3;基于由FPD3检测出的放射线进行图像处理的图像处理部4;对由图像处理部4进行了各种图像处理所得的放射线图像进行显示的显示部5。显示部5由监视器和电视机等的显示单元来构成。而且,在FPD3的检测面侧配置格栅31。
图像处理部4,由中央运算处理器(CPU)等构成。另外,将用于进行各种图像处理的程序等写入到以ROM(只读存储器)等为代表的存储介质并进行存储,并从该存储介质读出程序等之后由图像处理部4的CPU执行,从而进行对应该程序的图像处理。特别是,图像处理4的后述的分离部41、调整部42、除去部43,执行与放射线图像的分离、格栅图像的强度调整、格栅影响的除去相关的程序,由此分别进行对应该程序的放射线图像的分离、格栅图像的强度调整、格栅影响的除去。与放射线图像的分离、格栅图像的强度调整、格栅影响的除去相关的程序,相当于本发明的放射线图像处理程序。
图像处理部4具备:分离部41,其将放射线图像分离成包含格栅31的格栅像成分的格栅图像和包含其它成分的非格栅图像;调整部42,其基于非格栅图像,对在实际空间上的格栅图像的强度进行调整,并生成调整图像;除去部43,其从非格栅图像减去调整图像,并生成除去了格栅影响的补正图像。分离部41,相当于本发明的分离单元,调整部42,相当于本发明的调整单元,除去部43,相当于本发明的除去单元。
FPD3,如图2所示,在其检测面上将可感应放射线的多个检测元件d排列成二维矩阵状来构成。检测元件d,将透过被检测体M的放射线变换成电信号后姑且蓄积,并读出该蓄积的电信号,由此检测放射线。将由各个检测元件分别检测出的电信号,变换成对应电信号的像素值,并对分别与检测元件d的位置对应的像素分配其像素值,由此输出放射线图像,并将放射线图像送到图像处理部4的分离部41(参照图1)。
对实际的摄影以及各图像有关的数据的流向进行说明。按下省略图示的摄影按钮,使放射线源2发生放射线,将该放射线朝向被检测体M照射,并与其联动而开始摄影。发生的放射线透过被检测体M,通过格栅31入射到FPD3,FPD3检测放射线并输出放射线图像来进行摄影。如上所述,在FPD3的检测面侧配设格栅31,因此当放射线透过被检测体M之际发生的散射放射线被格栅31除去。所以,在从FPD3输出的放射线图像上重叠了格栅31的影响所产生的格栅像。
该重叠格栅像的放射线图像被送入图像处理部4的分离部41,分离部41把放射线图像分离成包含格栅像成分的格栅图像,和包含其它成分的非格栅图像。在本实施例中用频率分析来分离成分。即,就格栅像的成分而言,由于具有某种特定范围的频带的周期的图形很多,由通过相当格栅像的特定频带的带通滤波处理来分离成分。
当格栅31是标称格栅(按标准设计的格栅)的情况下,只要采用作为已知的标称格栅的格栅31的格栅密度Dg、以及同样已知的FPD3的图像传感器的抽样间距Ps及抽样频率fs,以fg=Ps·Dg·fs的公式就可以求出与格栅相关的格栅频率fg。另一方面,奈奎斯特频率fn是抽样频率fs的二分之一,即用fn=fs/2来表示。由抽样定理,如果fg>fn,则格栅频率fg为格栅像的频率,如果fg≤fn,则2fn-fg(即fs-fg)为格栅像的频率。然后,利用带通滤波处理只通过格栅像的频率,并将通过的格栅像的频率成分作为格栅图像,将剩余的频率成分作为非格栅图像,来分离放射线图像。
当格栅31不是标称格栅(即在整个图像中不稳定)、并且如上所述那样在图像中央部和端部中格栅纹的频带不同的情况下,则不采用上述的抽样定理,而在空间频域上分离格栅图像和非格栅图像。具体地,对放射线图像进行傅立叶变换求出功率谱,将功率谱表示极大值的范围的频带,作为格栅像的频率,来分离格栅图像和非格栅图像。
更具体地说明。将放射线图像的纵横作为x方向、y方向,将各像素值作为f(x、y),则在空间频域上的傅立叶变换后的F(μ、v),采用f(x、y)表示为如下的公式(1)。
F(μ、v)=∑∑f(x、y)·e-j2π(μx+vy)/N....(1)
其中,上述的公式(1)中的j为复数,∑是y=0~N-1为止的总和,x=0~N-1为止的总和。
由上述公式(1)对放射线图像f(x、y)进行傅立叶变换的空间频域上的图像F(μ、v)为如图3(a)所示。即,图3(a)是将放射线图像f(x、y)的功率谱进行可视化的图像F(μ、v)。将图3(a)的图像分割为预先设定的规定大小的块,在各个块中以块内的频谱强度的平均置换的频谱数据表示在图3(b)中。对图3(b)的频谱数据进行扫描,并提取表示极大值位置的块的结果表示在图3(c)中。在图3(c)中,将图中的涂黑所示的位置的块作为极大值来提取。
就(傅立叶变换后的)空间频域上的图像F(μ、v)而言,如果没有格栅的影响,则强度从图3的图像中心向端部单调地减少。如果不是单调减少而在有的地方出现极大值,则该地方(块)应该被认为是格栅像的频带。所以,从图3的图像中心进行扫描。与后述的格栅图像的强度调整中的“单调”一样,在本说明书中,设强度按照没有取极值的方式变化的情况为“单调”、强度按照取极值的方式变化的情况为不“单调”。
通过带通滤波处理将傅立叶变换的成分分离成相当于被提取的块位置的频带的成分(图3(c)的涂黑的块)、以及相当于除此以外的位置的频带的成分(图3(c)中的白部分),并分别对它们进行傅立叶反变换而进行图像化,由此生成格栅图像和非格栅图像。傅立叶反变换后的f(x、y),采用F(μ、v)表示为下述的公式(2)。
f(x、y)=∑∑F(μ、v)·ej2π(μx+vy)/N....(2)
其中,上述的公式(2)中的∑是v=0~N-1为止的总和,μ=0~N-1为止的总和。
这样,当格栅31为标称格栅的情况下,由抽样定理求出格栅像的频率,并根据使该格栅像的频率通过的带通滤波处理来分离图像成分,当格栅31不是标称格栅的情况下,将相当于在傅立叶变换的空间频域上功率谱为最大值的块位置的频率作为格栅像的频率,并根据使该格栅像的频率通过的带通滤波处理来分离图像成分,由此在空间频域上分离格栅图像和非格栅图像。
另外,对于分离的手法,并不限定于上述的内容。例如,也可以由分离为比预先设定的规定的值更大频带、和除此以外的频带,由此分离格栅图像和非格栅图像。也并不限定于空间频域上分离如上所述的格栅图像和非格栅图像,也可以利用重叠了格栅像的放射线图像在实际空间上成为周期性的图形,对于与放射线图像在实际空间上的强度相关的分布图,提取出与三角函数类似的成分(例如正弦波、余弦波、以及他们的积/和,或者与其它函数的积/和),由此分离图像成分。也就是,也可以从放射线图像在实际空间上的分布图中的、采用与三角函数类似的成分而与格栅像正交(格栅像的断面)的分布图中提取类似的成分(极值),将被提取的成分作为格栅图像。而且,对于频率分析并不限定于傅立叶变换,也可以采用小波变换或伽柏滤波等的其它频率分析。由分离部41进行的放射线图像的分离,相当于本发明的分离步骤。
将在分离部41生成的格栅图像和非格栅图像送入调整部42,并且仅将非格栅图像送入除去部43。调整部42,为了使从非格栅图像中减去后述的调整图像而得到的图像不包含与格栅图像相同的成分,调整部42调整格栅图像的强度,并生成调整图像。在本实施例中,调整部42以图4的流程进行调整。
(步骤S1)极值?
向与格栅像正交的方向(格栅像的断面方向)扫描格栅图像,调查是否为关注像素的强度表示极值的像素。具体地,对于与格栅像正交(格栅像的断面)的强度有关的分布图,在格栅图像与非格栅图像之间进行比较。图5(a)是格栅图像的分布图,图5(b)是与图5(a)相同位置的不单调时的非格栅图像的分布图,图5(c)是与图5(a)同样位置的单调时的非格栅图像的波形。
(步骤S2)下一像素
在图5(a)所示的分布图中,如果不是关注像素的强度表示极值的像素,就沿着与格栅像正交的方向推进到下一个像素。然后,返回到步骤S1,直到找到表示极值的像素为止,重复进行步骤S1、S2,沿着与格栅像正交的方向进行扫描。
(步骤S3)单调?
在图5(a)所示的分布图中,如果是关注像素的强度表示极值的像素,则如图5(a)所示,在强度表示极大或者极小的位置(即格栅图像的格栅位置)中,调查非格栅图像是否在局部单调变化。也就是,调查非格栅图像的格栅位置的局部强度变化是否单调。对于非格栅图像是否局部单调变化,是通过对格栅图像的强度的极大与极小所反转的位置(图5(a)中的A、C)、与其内部的极值的位置(图5(a)中的B)所分别对应的点的强度进行比较来决定的。如果将格栅图像的A、B、C点所对应的非格栅图像,如图5(b)、图5(c)所示那样设为A’、B’、C’,则将A’的强度以及C’的强度的平均值与B’的强度进行比较,当产生超过规定的范围的差异的情况下,判断为不单调。
如图5(b)所示,当格栅图像的B’成为极值的情况下,由于A’的强度以及C’的强度的平均值与B’的强度之差,超过规定的范围,因此判断为不单调。该B’的强度,可以被认为受格栅图像的B的强度影响而成为极值。也就是,由分离部41,从包含格栅像成分的格栅图像中分离出非格栅图像,但是被格栅图像的B的强度所影响的成分很小而且包含在非格栅图像中,因此可以认为非格栅图像的格栅位置的局部的强度变化是不单调的。如果不单调,则进入下一步骤S4。
相反,如图5(c)所示,当非格栅图像的B’没有成为极值的情况下,强度小的A’的强度以及强度大的C’的强度的平均值,由于A’以及C’的强度相抵消,以强度为中等的B’的强度没有变化。所以,A’的强度以及C’的强度的平均值与B’的强度之差,就收敛在规定的范围内,并判断为单调。如果是单调,则跳过下一步骤S4而进入再下一个步骤S5。
(步骤S4)调整
当在步骤S3中判断为非格栅图像的格栅位置的局部的强度变化不为单调的情况下,按照使格栅图像的格栅位置(即B)的强度等于非格栅图像的格栅位置(即B’)的局部的强度变化(A’的强度以及C’的强度的平均值与B’的强度之差)的方式,将格栅图像的强度调整为常数倍。
(步骤S5)所有的像素?
对一个关注像素进行步骤S1~S4(或者如果为单调则进行步骤S1~S3),并调查是否对所有的像素进行了扫描。当对所有的像素进行了扫描的情况下,认为对所有的像素进行了调整部42的格栅图像的强度调整进而结束一系列的调整。当没有对所有的像素进行扫描的情况下,返回到步骤S1,反复进行步骤S1~S4(或者如果为单调则进行步骤S1~S3),并沿着与格栅像正交的方向进行扫描,直到对所有的像素进行为止。
这样,进行步骤S1~S5来完成一系列的调整。另外,对于调整的手法,并不限定上述的内容。虽然对非格栅图像的格栅位置的局部的强度变化是否为单调进行了调查,但是除此以外,在调整部42中,创建对非格栅图像进行平滑化所得的图像,并且按照使格栅图像的格栅位置(B)的强度等于非格栅图像的格栅位置(B’)的强度与上述的平滑化了的图像的格栅位置的强度之差的方式,调整部42也可以调整格栅图像的强度,并生成调整图像。由调整部42进行的格栅图像的强度调整,相当于本发明的调整步骤。
由调整部42调整后所生成的调整图像被送入除去部43,并且如上所述那样将在分离部41中生成的非格栅图像送入除去部43。除去部43,从非格栅图像中减去调整图像,并生成除去了格栅影响的补正图像。由除去部43进行的格栅影响的除去,相当于本发明的除去步骤。
由除去部43除去后所生成的补正图像被送入显示部5,显示部5显示输出补正图像。而且,将由除去部43除去后所生成的补正图像,写入到RAM(随机存储器)等代表性的存储介质并进行存储,而根据适当必要也可以读出,或者也可以印刷输出到打印机等代表性的印刷单元上。
根据本实施例的放射线图像处理装置,分离部41将放射线图像分离成包含格栅31的格栅像的成分的格栅图像和包含其它成分的非格栅图像,调整部42基于非格栅图像,调整实际空间上的格栅图像的强度,并生成调整图像。然后,除去部43,从非格栅图像减去调整图像,并生成除去了格栅影响的补正图像。这样,由于是采用通过被检测体M而摄影的放射线图像来生成的补正图像,因此不通过被检测体M且也没有必要通过格栅31预先摄影成为用于除去格栅像的基准的补正图像,以一次摄影就可以获取各图像,没有必要使摄影环境或条件等完全一致,进而不受限制。而且,由于将通过被检测体M摄影的放射线图像分离为格栅图像和非格栅图像,并基于非格栅图像,调整实际空间上的格栅图像的强度,并生成补正图像,因此,利用基于调整的格栅图像的各像素中的强度来考虑了各像素的偏差(例如,FPD3的图像传感器的灵敏度的偏差或补正的偏差等),可以实现高精度的补正处理。其结果,不必预先获取用于格栅补正的数据而以一次摄影实现精度很高的补正处理。
在本实施例中,当格栅31不是标称格栅的情况下,分离部41,在空间频域上分离格栅图像和非格栅图像。由于是在空间频域上进行分离,因此可以更可靠地分离格栅像的成分和其它成分。当然,如对于放射线图像的实际空间上的强度相关的分布图,提取出与三角函数类似的成分的变形实施例所述那样,也可以在实际空间上分离格栅图像和非格栅图像。
分离部41,通过提取放射线图像的周期性的图形来分离图像成分。如上所述格栅像的成分,在很多情况下为具有某种特定范围的频带的周期性的图形,因此最好由通过相当于格栅像的特定频带的带通滤波处理来分离图像成分。
当基于非格栅图像、在调整实际空间上的格栅图像的强度之际,为了使从非格栅图像减去了调整图像的图像中不包含与格栅图像同样的成分,调整部42调整格栅图像的强度,并生成调整图像。具体地,当非格栅图像的格栅位置的局部的强度变化不为单调时,按照使格栅图像的格栅位置的强度等于非格栅图像的格栅位置的局部的强度变化的方式,调整部42调整格栅图像的强度,并生成图像。通过这样的调整,从非格栅图像减去了调整图像而所得到的图像(即补正图像)中不包含与格栅图像同样的成分,可以实现精度更高的补正。
本发明,并不限定上述的实施方式,可以进行下述的变形实施。
(1)在上述的实施例中,放射线采用X射线为例进行说明,当然也适用于X射线以外的放射线(例如γ射线等)。
(2)在上述的实施例中,放射线图像处理装置,用于医用等,是在图1所示的顶板1上放置被检测体进行摄影的构造,但是并不限定于此。例如,在用于工业用等的非破坏检查装置上,将被检测体(此时检查对象物为被检测体)搬运到传送带上进行摄影的构造也可以,用于医用的X射线CT装置等那样的构造也可以。而且,将进行摄影的装置作为外部装置以另外装置构成,简单地具备分离单元、调整单元、除去单元的装置也可以。
(3)格栅像容易显现的情况,是在采用由与各个像素对应配置的检测元件构成的平板型放射线检测器(FPD)那样的图像传感器的情况。所以,本发明适用于FPD,可以实现高精度的补正,特别有用,但是图像传感器并不限定于FPD。只要是通常采用的图像传感器,本发明都可以适用。
Claims (9)
1.一种放射线图像处理装置,对通过用于除去散射放射线的格栅进行摄影所得到的放射线图像进行处理,上述放射线图像处理装置的特征在于,具备:
分离单元,其将上述的放射线图像分离成包含上述格栅的格栅像成分的格栅图像和包含其它成分的非格栅图像;
调整单元,其基于上述非格栅图像,对在实际空间的上述格栅图像的强度进行调整,并生成调整图像;和
除去单元,其从上述非格栅图像中减去上述调整图像,并生成除去格栅影响后所得到的补正图像。
2.如权利要求1所述的放射线图像处理装置,其特征在于,
上述分离单元,在空间频域上分离上述格栅图像和非格栅图像。
3.如权利要求1或2所述的放射线图像处理装置,其特征在于,
上述分离单元,提取上述放射线图像中的周期性的图形,由此分离图像成分。
4.如权利要求3所述的放射线图像处理装置,其特征在于,
上述分离单元,采用使特定的频带通过的带通滤波处理,来分离图像成分。
5.如权利要求3所述的放射线图像处理装置,其特征在于,
上述分离单元,对于与上述放射线图像在实际空间上的强度有关的分布图,提取出与三角函数类似的成分,由此分离图像成分。
6.如权利要求1或2所述的放射线图像处理装置,其特征在于,
按照使从上述非格栅图像减去上述调整图像后所得到的图像中不包含与上述格栅图像相同的成分的方式,上述调整单元调整格栅图像的强度,并生成调整图像。
7.如权利要求6所述的放射线图像处理装置,其特征在于,
当上述非格栅图像的格栅位置的局部的强度变化不为单调时,按照使上述格栅图像的格栅位置的强度等于非格栅图像的格栅位置的局部的强度变化的方式,上述调整单元调整格栅图像的强度,并生成调整图像。
8.如权利要求6所述的放射线图像处理装置,其特征在于,
按照使上述格栅图像的格栅位置的强度等于上述非格栅图像的格栅位置的强度与对非格栅图像平滑化后所得到的图像的格栅位置的强度之差的方式,上述调整单元调整格栅图像的强度,并生成调整图像。
9.一种放射线图像处理方法,用于执行通过用于除去散射放射线的格栅进行摄影所得到的放射线图像的处理,上述放射线图像处理方法的特征在于,包括下述步骤:
分离步骤,将上述的放射线图像分离成包含上述格栅的格栅像成分的格栅图像和包含其它成分的非格栅图像;
调整步骤,基于上述非格栅图像,在实际空间上对上述格栅图像的强度进行调整,并生成调整图像;和
除去步骤,从上述非格栅图像中减去上述调整图像,并生成除去格栅影响后所得到的补正图像。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2007/069262 WO2009044452A1 (ja) | 2007-10-02 | 2007-10-02 | 放射線画像処理装置および放射線画像処理プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101801273A CN101801273A (zh) | 2010-08-11 |
CN101801273B true CN101801273B (zh) | 2013-05-15 |
Family
ID=40525885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007801005937A Active CN101801273B (zh) | 2007-10-02 | 2007-10-02 | 放射线图像处理装置以及放射线图像处理方法 |
Country Status (6)
Country | Link |
---|---|
US (1) | US9655570B2 (zh) |
EP (1) | EP2196148B1 (zh) |
JP (1) | JP5212376B2 (zh) |
KR (1) | KR101095270B1 (zh) |
CN (1) | CN101801273B (zh) |
WO (1) | WO2009044452A1 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5361575B2 (ja) * | 2009-07-01 | 2013-12-04 | 株式会社東芝 | X線画像撮影装置及び画像処理装置 |
US8831325B2 (en) | 2009-12-29 | 2014-09-09 | Shimadzu Corporation | Radiographic image processing apparatus and radiographic image processing program |
JP5407937B2 (ja) * | 2010-03-02 | 2014-02-05 | 株式会社島津製作所 | X線撮影装置 |
EP2700359B1 (en) * | 2011-04-22 | 2015-07-15 | Shimadzu Corporation | X-ray diagnostic apparatus and x-ray diagnostic program |
JP6128463B2 (ja) * | 2013-11-06 | 2017-05-17 | 富士フイルム株式会社 | 放射線画像処理装置および方法並びにプログラム |
JP2018038500A (ja) * | 2016-09-06 | 2018-03-15 | キヤノン株式会社 | 放射線撮影装置、放射線撮影システム、放射線撮影方法、及びプログラム |
JP6815818B2 (ja) * | 2016-10-17 | 2021-01-20 | キヤノン株式会社 | 放射線撮影システム及び放射線撮影方法 |
EP3427664A1 (en) | 2017-07-13 | 2019-01-16 | Koninklijke Philips N.V. | A device for scatter correction in an x-ray image and a method for scatter correction in an xray image |
JP6897585B2 (ja) * | 2018-01-24 | 2021-06-30 | コニカミノルタ株式会社 | 放射線画像処理装置、散乱線補正方法及びプログラム |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1853570A (zh) * | 2005-04-18 | 2006-11-01 | 通用电气公司 | Ct成像方法和系统 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2507659B2 (ja) | 1990-04-04 | 1996-06-12 | 株式会社東芝 | 放射線画像読取装置 |
US5661818A (en) * | 1995-01-27 | 1997-08-26 | Eastman Kodak Company | Method and system for detecting grids in a digital image |
US6333990B1 (en) * | 1998-06-02 | 2001-12-25 | General Electric Company | Fourier spectrum method to remove grid line artifacts without changing the diagnostic quality in X-ray images |
US6269176B1 (en) * | 1998-12-21 | 2001-07-31 | Eastman Kodak Company | Method for x-ray antiscatter grid detection and suppression in digital radiography |
US20010033638A1 (en) * | 2000-02-04 | 2001-10-25 | Hitoshi Inoue | Image acquisition method and apparatus |
JP4746761B2 (ja) | 2001-05-01 | 2011-08-10 | キヤノン株式会社 | 放射線画像処理装置、放射線画像処理方法、記憶媒体、及びプログラム |
JP3825989B2 (ja) * | 2001-05-01 | 2006-09-27 | キヤノン株式会社 | 放射線画像処理装置、画像処理システム、放射線画像処理方法、コンピュータ読出可能な記憶媒体、及びコンピュータプログラム |
US7142705B2 (en) | 2001-05-01 | 2006-11-28 | Canon Kabushiki Kaisha | Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program |
JP3825988B2 (ja) | 2001-05-01 | 2006-09-27 | キヤノン株式会社 | 放射線画像処理装置、画像処理システム、放射線画像処理方法、記録媒体、及びプログラム |
US7050618B2 (en) * | 2002-02-08 | 2006-05-23 | Eastman Kodak Company | Method for antiscatter stationary grid artifacts detection and attenuation in digital radiographic images |
JP3631215B2 (ja) | 2002-03-12 | 2005-03-23 | キヤノン株式会社 | 放射線画像処理装置、放射線画像処理システム、放射線撮影システム、放射線撮影装置、放射線画像処理方法、コンピュータ可読記憶媒体、及びプログラム |
JP4596748B2 (ja) * | 2003-05-07 | 2010-12-15 | キヤノン株式会社 | 放射線画像撮影装置及び放射線画像撮影装置における再構成方法 |
JP4363095B2 (ja) * | 2003-07-01 | 2009-11-11 | コニカミノルタエムジー株式会社 | 医用画像処理装置及び医用画像処理システム |
JP3903027B2 (ja) * | 2003-08-07 | 2007-04-11 | キヤノン株式会社 | 放射線画像処理方法及び装置並びにグリッドの選別方法及び装置 |
JP2004344670A (ja) | 2004-06-18 | 2004-12-09 | Canon Inc | X線画像処理装置及び処理方法 |
US7479969B2 (en) * | 2005-02-24 | 2009-01-20 | Agfa Healthcare, N.V. | Grid detection method |
US7796792B2 (en) * | 2005-06-29 | 2010-09-14 | Agfa Healthcare, N.V. | Method of identifying disturbing frequencies originating from the presence of an anti-scatter grid during acquisition of a radiation image |
JP4393436B2 (ja) | 2005-08-10 | 2010-01-06 | キヤノン株式会社 | 放射線画像処理装置、画像処理システム、放射線画像処理方法、記憶媒体及びプログラム |
-
2007
- 2007-10-02 CN CN2007801005937A patent/CN101801273B/zh active Active
- 2007-10-02 JP JP2009535910A patent/JP5212376B2/ja active Active
- 2007-10-02 KR KR1020107004757A patent/KR101095270B1/ko active IP Right Grant
- 2007-10-02 WO PCT/JP2007/069262 patent/WO2009044452A1/ja active Application Filing
- 2007-10-02 EP EP07829002.0A patent/EP2196148B1/en not_active Not-in-force
- 2007-10-02 US US12/680,773 patent/US9655570B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1853570A (zh) * | 2005-04-18 | 2006-11-01 | 通用电气公司 | Ct成像方法和系统 |
JP2006297083A (ja) * | 2005-04-18 | 2006-11-02 | General Electric Co <Ge> | 分散型x線源及び補間に基づく再構成を用いたct撮像の方法及びシステム |
US7227923B2 (en) * | 2005-04-18 | 2007-06-05 | General Electric Company | Method and system for CT imaging using a distributed X-ray source and interpolation based reconstruction |
Also Published As
Publication number | Publication date |
---|---|
EP2196148A4 (en) | 2013-06-26 |
US20100215243A1 (en) | 2010-08-26 |
WO2009044452A1 (ja) | 2009-04-09 |
CN101801273A (zh) | 2010-08-11 |
KR101095270B1 (ko) | 2011-12-20 |
KR20100040324A (ko) | 2010-04-19 |
US9655570B2 (en) | 2017-05-23 |
JPWO2009044452A1 (ja) | 2011-01-27 |
JP5212376B2 (ja) | 2013-06-19 |
EP2196148A1 (en) | 2010-06-16 |
EP2196148B1 (en) | 2015-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101801273B (zh) | 放射线图像处理装置以及放射线图像处理方法 | |
TWI327375B (en) | Method and apparatus for setting black level in an imager using both optically black and tied pixels | |
JP2017143436A5 (zh) | ||
ATE500688T1 (de) | Verfahren und vorrichtung zur verringerung der effekte von dunkelstrom und defekten pixeln in einer abbildungseinrichtung | |
JP5698577B2 (ja) | 画像処理装置、画像処理方法並びに画像処理プログラム | |
CN103312996B (zh) | 图像处理设备、图像处理方法和放射线系统 | |
CN103824258B (zh) | 图像处理装置以及图像处理方法 | |
US6542189B1 (en) | Frequency compensated multiple sampling method for a video digital camera | |
CN105021282B (zh) | 快照式偏振高光谱相机及成像方法 | |
JP2003233818A (ja) | 周期的パターン判定方法および装置並びにプログラム | |
JP5818475B2 (ja) | 画像処理装置、画像処理方法並びに画像処理プログラム | |
KR101716111B1 (ko) | 이물질 검출 시스템 및 방법 | |
JP2014176565A (ja) | 画像処理装置、放射線撮影装置、画像処理方法、コンピュータプログラム及び記憶媒体 | |
JP6159811B2 (ja) | 解像度および焦点改善 | |
JP5635389B2 (ja) | 画像処理装置、画像処理プログラム、及びx線画像診断装置 | |
JP2018186502A5 (zh) | ||
JP4968057B2 (ja) | 放射線検出器のライン状異常画像素子の検出方法と放射線撮像装置 | |
JP7208723B2 (ja) | 画像処理装置およびその制御方法、放射線撮影システム、プログラム | |
JP4712436B2 (ja) | デジタル信号における周期的変動の抑制 | |
US6995387B2 (en) | Suppression of periodic variations in a digital signal | |
JP6570594B2 (ja) | 画像処理装置、放射線撮影装置、画像処理方法、コンピュータプログラム及び記憶媒体 | |
JP2005167773A (ja) | 画像処理方法及び装置 | |
Landon | A pipeline for digital restoration of deteriorating photographic negatives | |
TW202238460A (zh) | 放射線圖像處理方法、機械學習方法、學習完成模型、機械學習之前處理方法、放射線圖像處理模組、放射線圖像處理程式及放射線圖像處理系統 | |
JP5456404B2 (ja) | 欠陥画素検出方法および画像処理装置 |
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 |