CN104684479B - 体动检测装置及方法 - Google Patents
体动检测装置及方法 Download PDFInfo
- Publication number
- CN104684479B CN104684479B CN201380049990.1A CN201380049990A CN104684479B CN 104684479 B CN104684479 B CN 104684479B CN 201380049990 A CN201380049990 A CN 201380049990A CN 104684479 B CN104684479 B CN 104684479B
- Authority
- CN
- China
- Prior art keywords
- contrast
- radiation image
- low frequency
- component
- frequency component
- 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
- 238000001514 detection method Methods 0.000 title claims description 44
- 238000000034 method Methods 0.000 title description 24
- 230000005855 radiation Effects 0.000 claims abstract description 131
- 238000000605 extraction Methods 0.000 claims description 19
- 238000003384 imaging method Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 5
- 210000000746 body region Anatomy 0.000 description 12
- 230000002285 radioactive effect Effects 0.000 description 12
- 241000894007 species Species 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 3
- 230000003321 amplification Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- 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/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- 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
-
- 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/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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/20172—Image enhancement details
- G06T2207/20201—Motion blur correction
-
- 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/30008—Bone
- G06T2207/30012—Spine; Backbone
-
- 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/30168—Image quality inspection
-
- 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/30196—Human being; Person
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Quality & Reliability (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Multimedia (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
无论拍摄条件或被摄体如何都能够高精度地检测放射线图像所含的体动。对比度计算部(24)算出解析点设定部(23)设定的解析点处的、进行了转换的放射线图像(G0)的高频分量的对比度及低频分量的对比度。比率计算部(25)算出高频分量的对比度与低频分量的对比度的比率。判定部(26)基于比率判定拍摄放射线图像(G0)时有无体动。显示控制部(27)将判定结果显示于显示部(40)。
Description
技术领域
本发明涉及一种使用放射线图像对拍摄时的被摄体的体动进行检测的体动检测装置及方法。另外，本发明涉及一种在放射线图像中设定用于对体动进行检测的解析点的解析点设定装置及方法。
背景技术
在对被摄体照射X射线等放射线而拍摄被摄体，并取得被摄体的放射线图像时，有时由于拍摄中被摄体的移动即体动，而取得的放射线图像会产生模糊，画质会变差。特别是在拍摄对象为颈部、胸部、腹部及腰部等放射线的照射时间较长的部位或拍摄小孩的情况下，容易发生体动，因此，放射线图像变差成为问题。
可通过进行拍摄的技师观察拍摄中的被摄体的状态并立即确认取得的放射线图像来辨别有无这种体动。但是，由于很多时候摄影室中的显示器的分辨率较低，缩小显示放射线图像，且观察环境比较亮，因此，难以观察显示的放射线图像来判定有无体动。在这种情况下，能够放大放射线图像，但放大放射线图像的操作会降低技师的工作效率。
因此，提出了基于放射线图像自动地检测体动的方法。例如，提出了如下方法：预先学习放射线图像中含的边缘的劣化作为纵向和横向两个方向的特征量，并从放射线图像提取边缘，基于所学习的特征量与所提取的边缘，对在放射线图像中是否产生了体动引起的模糊进行判定(参考专利文献1)。另外，还提出了如下方法：学习关于头部的侧面这一特定部位的特征点，并使用学习结果对取得的放射线图像的体动进行检测(参考非专利文献1)。
为了检测这种体动，需要在放射线图像上设定用于检测体动的解析点。因此，技师可以通过手动操作在显示的放射线图像上设定解析点而对体动进行检测。另外，还提出了如下方法：自动地提取具有用于体动检测的高空间频率的部位的部分像素作为解析点(参考专利文献2)。
现有技术文献
专利文献
专利文献1：美国专利第7899229号的说明书
专利文献2：日本特开2005-205086号公报
非专利文献
非专利文献1：Proc.Of SPIE Vol.6914,69140U,(2008)
发明内容
发明所要解决的课题
但是，由于专利文献1及非专利文献1所记载的方法使用了关于特定部位的学习结果，因此，无法用于拍摄未进行学习的部位的放射线图像。另外，在专利文献1所记载的方法中，由于学习边缘的劣化作为两个方向的特征点，因此，在通过与学习时不同的拍摄条件取得放射线图像的情况下、或者放射线图像中含有的被摄体等与学习的放射线图像中含有的被摄体不同的情况下，无法高精度地检测体动，因此，存在检测的稳定性低，且可应用的被摄体的范围较窄的问题。
另外，在拍摄中发生了体动的情况下，在所取得的放射线图像中，由于体动偏移而高频分量的对比度下降，因此，在专利文献2所述的方法中，在发生体动的情况下，无法高精度地提取解析点。
本发明鉴于上述情况而作出，其目的在于，无论拍摄条件或者被摄体如何都可高精度地检测放射线图像中所含的体动。
另外，本发明的另一目的在于，在放射线图像上高精度地设定用于检测体动的解析点。
用于解决课题的方案
本发明的体动检测装置，其特征在于，具备：
对比度计算单元，算出放射线图像的高频分量的对比度及低频分量的对比度；
比率计算单元，算出高频分量的对比度与低频分量的对比度的比率；及
判定单元，基于比率判定放射线图像所含的被摄体有无体动。
另外，在本发明的体动检测装置中，可以是对比度计算单元根据放射线图像生成含有高频分量的高频图像及含有低频分量的低频图像，并根据高频图像及低频图像算出高频分量的对比度及低频分量的对比度。
在该情况下，也可以是对比度计算单元在放射线图像所含的、用于算出高频分量的对比度及低频分量的对比度的一个以上的解析点，算出高频分量的对比度及低频分量的对比度。
另外，在该情况下，也可以是对比度计算单元在高频图像和低频图像的各解析点，决定以解析点为中心的解析区域所含的边缘部分的倾斜方向，并算出高频图像和低频图像各自的倾斜方向的对比度作为高频分量的对比度及低频分量的对比度。
另外，在本发明的体动检测装置中，可以是在有多个解析点的情况下，比率计算单元对应在多个解析点的各解析点决定的每个倾斜方向，算出高频分量的对比度及低频分量的对比度的比率。
在该情况下，也可以是判定单元将每个倾斜方向的比率中的值最小的比率决定为表示体动的指标值，通过将指标值与特定的阈值进行比较来判定有无体动。
另外，在本发明的体动检测装置中，也可以进一步具备解析点设定单元，对放射线图像设定解析点。
在该情况下，也可以是解析点设定单元从放射线图像提取被摄体的区域，并提取区域中的低频分量超过特定的阈值的点作为解析点。
另外，在本发明的体动检测装置中，也可以进一步具备显示控制单元，用于将有无体动的判定结果显示于显示单元。
特别是在设定解析点的情况下，也可以是显示控制单元将表示解析点的指标值的大小的信息赋予解析点，并显示放射线图像。
本发明的体动检测方法，其特征在于，
算出放射线图像的高频分量的对比度及低频分量的对比度，
算出高频分量的对比度与低频分量的对比度的比率，
基于比率判定放射线图像所含的被摄体有无体动。
另外，也可以将本发明的体动检测方法作为用于使计算机执行的程序而进行提供。
本发明的解析点设定装置，其特征在于，具备：
被摄体区域提取单元，从放射线图像提取作为对象的被摄体的区域；
解析点提取单元，提取上述被摄体的区域中的空间频率的低频分量超过特定的阈值的点作为用于进行检测体动用的处理的解析点。
另外，在本发明的解析点设定装置中，也可以基于上述被摄体的区域中的像素值的标准差来决定阈值。
本发明的体动检测装置，其特征在于，具备：
本发明的解析点设定装置；
对比度计算单元，在上述解析点算出所述放射线图像的高频分量的对比度及低频分量的对比度；
比率计算单元，算出所述高频分量的对比度与所述低频分量的对比度的比率；
判定单元，基于该比率判定上述放射线图像所含的被摄体有无体动。
本发明的解析点设定方法，其特征在于，
从放射线图像提取作为对象的被摄体的区域，
提取上述被摄体的区域中的空间频率的低频分量超过特定的阈值的点作为用于进行检测体动用的处理的解析点。
本发明的另一体动检测方法，其特征在于，
通过本发明的解析点设定方法在上述放射线图像上设定解析点，
在上述解析点算出上述放射线图像的高频分量的对比度及低频分量的对比度，
算出上述高频分量的对比度与上述低频分量的对比度的比率，
基于该比率判定上述放射线图像所含的被摄体有无体动。
此外，本发明的解析点设定方法和另一体动检测方法也可以作为用于由计算机执行的程序提供。
发明效果
根据本发明的体动检测装置及方法，算出放射线图像的高频分量的对比度及低频分量的对比度，算出高频分量的对比度与低频分量的对比度的比率，并基于比率判定放射线图像所含的被摄体有无体动。在此，在拍摄中产生了体动的情况下，在所取得的放射线图像中，高频分量的对比度降低，但低频分量的对比度几乎不降低。另外，根据拍摄条件或被摄体的种类不同，放射线图像整体的对比度会变动，但其变动表现为低频分量的对比度的差。因此，能够使用高频分量的对比度与低频分量的对比度的比率作为表示有无体动的指标。由于本发明不使用学习结果，而是使用不依赖于拍摄条件或被摄体的种类的高频分量的对比度及低频分量的对比度，因此，能够不受拍摄条件或被摄体的种类的影响地稳定地对有无体动进行判定。
另外，根据本发明的解析点设定装置及方法，从放射线图像提取作为对象的被摄体的区域，并提取被摄体的区域中的空间频率的低频分量超过特定的阈值的点作为用于进行检测体动的处理的解析点。在此，在拍摄中产生了体动的情况下，在所取得的放射线图像中，高频分量的对比度降低，但低频分量的对比度几乎不降低。因此，根据本发明，能够不受体动影响地在放射线图像中高精度地设定用于进行检测体动的处理的解析点，其结果是，可以高精度地检测体动。
附图说明
图1是表示应用了本发明的实施方式的解析点设定装置及体动检测装置的放射线图像拍摄系统的结构的示意图。
图2是表示本实施方式中进行的处理的流程图。
图3是用于对多分辨率转换进行说明的图。
图4是表示解析点设定部的结构的概略框图。
图5是表示放射线图像的例子的图。
图6是表示设定的解析点的例子的图。
图7是表示以解析点为中心的8个方向的图。
图8是用于对解析区域及水平方向的线进行说明的图。
图9是表示判定为存在体动的情况下的判定结果的显示画面的例子的图。
图10是表示判定为存在体动的情况下的判定结果的显示画面的例子的图。
图11是表示判定为存在体动的情况下的判定结果的显示画面的例子的图。
图12是表示在横轴上绘制对比度比的值，在纵轴上绘制各对比度比的值中判定为存在体动的图像数量及判定为没有体动的图像数量的直方图的图。
具体实施方式
以下，参照附图详细地对本发明的实施方式进行说明。图1是表示应用了本发明的实施方式的解析点设定装置及体动检测装置的放射线图像拍摄系统的结构的示意图。本实施方式的放射线图像拍摄系统包含：拍摄装置10，通过拍摄被摄体而取得被摄体的放射线图像G0；信号处理部20，具备基于放射线图像G0而检测体动的、本实施方式的解析点设定装置及体动检测装置；输入部30，用于对信号处理部20提供各种指示；显示部40，显示通过拍摄而取得的放射线图像等。
拍摄装置10具备：射线球管12，对被摄体H照射放射线R；拍摄控制部13，对射线球管12的驱动进行控制；及载置被摄体H的拍摄台14。在拍摄台14上设有输出被摄体H的放射线检测信号的放射线检测器11。放射线检测器11的矩阵状配置的每个像素输出与照射放射线的能级对应的放射线检测信号，并对该检测信号进行A/D转换处理，作为表示被摄体H的放射线图像G0的数字图像信号而输出。
此外，作为放射线检测器11，可以使用：如日本特开平7-72253号公报所记载的那样，层叠有接受放射线的照射而发出可见光的闪烁器和对该可见光进行检测的固态光检测元件的放射线检测器；或如日本特开2010-206067号公报所记载的那样，其具有放射线光导层，上述放射线光导层接收放射线的照射而输出与其能量对应的电信号。
信号处理部20具备：输入表示放射线图像G0的数字图像信号的预处理部21、多分辨率转换部22、解析点设定部23、对比度计算部24、比率计算部25、判定部26及显示控制部27。
输入部30例如由键盘31及鼠标32等构成，用于将放射线技师等使用者的各种指示输入信号处理部20。
显示部40例如由液晶显示器或CRT显示器等构成，并根据需要显示有无体动的判定结果、所拍摄的被摄体的放射线图像等。
上述的信号处理部20、输入部30及显示部40例如可以由普通的个人计算机等的计算机系统构成。
接着，对放射线图像的拍摄进行说明。在拍摄放射线图像时，放射线检测器11载置在拍摄装置10的拍摄台14的上方，并在其上方放置被摄体H。在该状态下，通过操作拍摄控制部13而驱动射线球管12，放射线R透过被摄体H而照射到放射线检测器11上。当拍摄结束时，从放射线检测器11取得表示放射线图像G0的数字图像信号。该放射线图像G0可在显示部40中重放显示。
放射线图像G0被输入信号处理部20。以下，对信号处理部20进行的处理进行说明。图2是表示在信号处理部20中进行的处理的流程图。对于输入了信号处理部20的放射线图像G0，首先在预处理部21中接受：对放射线的照射不均及放射线检测器11的检测不均等引起的信号值的变动进行校正的处理；为了提高放射线图像G0的画质而对浓度、对比度、频率分量等进行校正的图像处理；及根据需要而适当进行的其它处理(预处理，步骤ST1)。
接着，多分辨率转换部22对预处理后的放射线图像G0进行多分辨率转换，并提取出放射线图像G0的高频分量及低频分量(步骤ST2)。图3是用于对多分辨率转换进行说明的图。首先，多分辨率转换部22通过σ＝1的高斯滤光片对放射线图像G0进行滤波处理，将放射线图像G0缩小至1/2而生成高斯分量G1。高斯分量G1为将放射线图像G0缩小至1/2而成的分量。接着，多分辨率转换部22进行例如三次B样条插值等插值运算，而将高斯分量G1放大至与放射线图像G0相同的尺寸，并从放射线图像G0减去放大后的高斯分量G1，从而生成最高频段的拉普拉斯分量L0。在本实施方式中，将该拉普拉斯分量L0作为高频分量使用。此外，在本实施方式中，为了方便，将最高频段称为第零频段。
接着，多分辨率转换部22通过σ＝1的高斯滤光片对高斯分量G1进行滤波处理，将高斯分量G1进一步缩小至1/2而生成高斯分量G2，将高斯分量G2放大至与高斯分量G1相同的尺寸，并从高斯分量G1减去放大后的高斯分量G2，从而生成第一频段的拉普拉斯分量L1。进一步，反复进行上述处理直至生成所需的频段的拉普拉斯分量，由此，生成多个频段的拉普拉斯分量Lj(j＝0～n)。
在本实施方式中，反复进行上述处理直至获得第三频段的拉普拉斯分量L3，并以形成与第0频段的拉普拉斯分量L0相同的尺寸的方式将拉普拉斯分量L3放大8倍，从而生成拉普拉斯分量L3′。而且，将该拉普拉斯分量L3′作为低频分量使用。
在此，高斯分量的各像素的信号值表示该像素的浓度，拉普拉斯分量的各像素的信号值表示该像素中的该频段的频率分量的大小。
此外，也可以通过使用小波变换等其它多分辨率转换的方法，生成频段不同的多个频段图像而将最高频段的频段图像作为高频分量使用，将与上述相同的第三频段的频段图像放大8倍的图像作为低频分量使用。另外，低频分量的频段不限于第三频段，只要是比最高频段低的频段，就能够将任意频段的频率分量作为低频分量使用。另外，也可以使用放射线图像G0本身作为高频分量。
接着，解析点设定部23在实施了预处理的放射线图像G0中设定后述的用于算出对比度的解析点(步骤ST3)。具体来说，使用低频分量L3′设定解析点。以下，对解析点的设定进行说明。图4是表示解析点设定部23的结构的概略框图。如图4所示，解析点设定部23具备：被摄体区域提取部51及解析点提取部52。被摄体区域提取部51为了排除放射线图像G0中含有的皮肤线、照射野框的影响，将检测体动的对象限定在放射线图像G0中含有的、包含骨部的被摄体区域内，而从放射线图像G0提取被摄体区域。在此，为了判定在拍摄中被摄体是否产生了体动，需要提取放射线图像G0中的高对比度的区域。因此，在本实施方式中，被摄体区域提取部51提取包含骨部的区域作为被摄体区域。
具体来说，采用使用关于作为对象的被摄体的学习结果的方法、关于被摄体的规则库的方法、图切割方法等各种分割技术，从放射线图像G0提取被摄体区域。例如，如图5所示，在从包含腰椎的放射线图像G0中提取腰椎部分的情况下，使用对腰椎的形状进行了学习的学习结果，从放射线图像G0提取包含腰椎的、预定大小的区域作为被摄体区域A0。
解析点提取部52从被摄体区域A0中提取作为判定有无体动的处理的基准的解析点。解析点优选是从容易判定有无体动的部分提取。具体来说，优选在骨部的边缘提取解析点。因此，解析点提取部52从多分辨率转换部22生成的低频分量L3′中的被摄体区域A0提取像素值比较高的像素作为解析点。具体来说，解析点提取部52将低频分量L3′的被摄体区域A0内的各像素的像素值与阈值Th1进行比较，提取像素值超过阈值Th1的像素位置作为解析点。另外，阈值Th1例如通过Th1＝σ×α算出。在此，σ为低频分量L3′的被摄体区域A0内的像素值的标准差，α为系数。系数α根据低频分量L3′的各像素的像素值适当设定即可。图6表示这样设定的解析点的例子。
接着，对比度计算部24算出解析点的高频分量的对比度及低频分量的对比度(步骤ST4)。在此，在对被摄体照射放射线的微小的时间(数十毫秒～数百毫秒)内，体动表现为向一个方向的偏移。在发生体动的情况下，放射线图像G0在与体动的方向相同的方向上变差，而在与体动的方向垂直的方向上不变差。因此，对比度计算部24首先算出各解析点的边缘的倾斜方向。具体来说，如图7所示，对比度计算部24在低频分量L3′上设定以解析点P0为中心的8个方向。例如，在将从纸面的下方朝向上方的方向设为基准(0度)的情况下，设定-67.5度、-45度、-22.5度、0度、22.5度、45度、67.5度及90度8个方向。另外，将8个方向分别设为θ1～θ8。另外，设定的方向不限于8个方向，也可以设定2个方向、4个方向或16个方向等任意数量。
而且，在每一方向上分别设定以解析点P0为中心的特定尺寸的解析区域。解析区域例如是9×9像素的矩形区域。而且，在解析区域内假想地设定有多条线，并算出各条线的像素值的最大值与最小值的差值作为各条线的对比度。在此，线的方向与设定了该解析区域的方向一致。例如，在对图7所示的方向θ1设定解析区域的情况下，设定于该解析区域的线的方向为θ1。
图8是用于对解析区域及水平方向的线进行说明的图。如图8所示，在对解析点P0的某方向θ1设定的解析区域A2内设定有多条线。另外，使线的数量与解析区域A2的垂直方向的像素数量相同即可。即，若解析区域A2的尺寸为9×9像素，则线的数量为9即可。但是，也可以适当稀疏化而减少线的数量。
而且，对比度计算部24算出全部线的对比度的平均值，作为在低频分量L3′中设定该解析区域的方向的对比度。而且，对比度计算部24对于在全部8个方向θ1～θ8上设定的各解析区域，算出对比度的平均值作为关于各解析区域的对比度CLθi(i＝1～8)。此外，对比度计算部24将对算出了所算出的8个对比度CLθi中值最高的解析区域设定的方向决定为该解析点的边缘的倾斜方向。
同样，对于高频分量L0也与低频分量L3′相同地，对比度计算部24算出全部的解析点中关于全部8个方向θ1～θ8的8个对比度CHθi(i＝1～8)。而且，将对算出所算出的8个对比度CHθi中值最高的解析区域设定的方向决定为该解析点的边缘的倾斜方向。
另外，因为在各解析点中，无论是低频分量L3′还是高频分量L0，对比度变化的方向均相同，因此，各解析点的边缘的倾斜方向一致。因此，算出在各解析点中相同边缘倾斜方向的、各一个高频分量的对比度CHθi及低频分量的对比度CLθi。
在此，在拍摄中产生了体动的情况下，在所取得的放射线图像G0中高频分量的对比度降低，但低频分量的对比度几乎不降低。另外，虽然根据拍摄条件或被摄体的种类不同，而放射线图像整体的对比度会变动，但该变动表现为低频分量的对比度的差。因此，高频分量的对比度与低频分量的对比度的比率可以作为表示有无体动的指标使用。
因此，比率计算部25算出各解析点的高频分量的对比度CHθi及低频分量的对比度CLθi的比率，更具体来说，算出高频分量的对比度CHθi相对于低频分量的对比度CLθi的比率CRθi(＝CHθi/CLθi)(步骤ST5)。此外，算出每一个方向θi的比率CRθi的平均值CRmθi，并算出最低的平均值CRmθi，作为用于判定有无体动的指标值(步骤ST6)。
另外，用于判定有无体动的指标值的计算不限于上述方法，也可以算出各解析点的全部8个方向的高频分量的对比度CHθi及低频分量的对比度CLθi的比率CRθi，并对每个方向θi算出全部解析点的高频分量的对比度CHθi及低频分量的对比度CLθi的比率CRθi的平均值CRmθi，算出该平均值中最低的平均值CRmL，作为用于判定有无体动的指标值。
如上所述，在拍摄中产生了体动的情况下，在所取得的放射线图像G0中，高频分量的对比度降低，但低频分量的对比度几乎不降低。因此，在没有发生体动的情况下，指标值为接近1的值，在产生了体动的情况下，指标值比1小，并且体动越大则值越小。
因此，判定部26将指标值与阈值Th2进行比较，对指标值是否在阈值Th2以上进行判定(步骤ST7)。在指标值为阈值Th2以上的情况下，判定部26判定为没有体动(步骤ST8)，在低于阈值Th2的情况下，判定部26判定为存在体动(步骤ST9)。另外，阈值Th2为比1小的预定的值。
显示控制部27将判定部26的判定结果显示于显示部40(步骤ST10)，并结束处理。图9是表示判定为存在体动的情况下的判定结果的显示画面的图。如图9所示，在判定结果的显示画面中，对解析点标记了表示放射线图像G0、体动的方向及大小的箭头。此外，标记有箭头的解析点的数量可以是将对各解析点算出的比率CRθi按升序排列时的前预定数量(图9中为4个)，也可以是指标值低于阈值Th2的全部解析点。另外，箭头也可以根据比率CRθi的大小而变更颜色，如图9所示，也可以根据比率CRθi的大小变更线型。此外，在本实施方式中，由于比率CRθi越小则体动越大，因此，箭头用于表示比率CRθi越小则体动越大。另外，如图10所示，也可以对应比率CRθi的大小变更箭头的长度。在此，在图10中，比率越小则箭头越长。通过这样显示判定结果，可容易地确认放射线图像G0中产生了体动的部分及体动的大小。
另外，即使如上所述标记箭头的情况下，有时也想要通过目视在放射线图像G0中确认体动。但是，由于放射线图像的摄影室中的显示部40的分辨率较小，且观察环境也比较亮，不容易通过目测确认体动。因此，可以设为如图11所示那样在放射线图像G0上设定预定尺寸的关注区域W，并将关注区域W的图像放大而与放射线图像G0并排显示。另外，关注区域W的位置可以以包含标记有箭头的解析点的方式预先设定，也可以由使用者使用输入部30移动到放射线图像G0上的所期望的位置。
此外，也可以将指标值显示于显示部40。在该情况下，由于仅显示数值的话难以判定体动的大小，因此，也可以使用到目前为止判定有无体动的结果，使人容易了解目前解析中的放射线图像G0中有多大程度的体动。例如，如图12所示，在横轴上绘制指标值即对比度的比率(对比度比)的值，在纵轴上绘制各对比度比的值中判定为存在体动的图像数量及判定为没有体动的图像数量的直方图，也可以在该直方图中通过箭头指示出目前解析中的放射线图像G0的指标值。另外，在图12中，直方图中的阈值Th2左侧的区域表示存在体动，阈值右侧的区域表示没有体动。
另一方面，在判定为没有体动的情况下，显示“没有体动”的消息，或者发出声音即可。另外，在判定为存在体动的情况下，也可以显示“发现体动。请重新拍摄”的消息，或者发出声音。
这样，在本实施方式中，基于放射线图像G0的高频分量的对比度及低频分量的对比度的比率，对放射线图像G0所含有的被摄体有无体动进行判定。在此，在拍摄中产生了体动的情况下，在所取得的放射线图像G0中，高频分量的对比度降低，但低频分量的对比度几乎不降低。另外，根据拍摄条件或被摄体的种类不同，而放射线图像整体的对比度会变动，但该变动表现为低频分量的对比度的差。因此，高频分量的对比度与低频分量的对比度的比率可以作为表示有无体动的指标使用。在本实施方式中，由于不是如上述专利文献1等方法那样使用学习结果，而使用不依赖于拍摄条件或被摄体的种类的高频分量的对比度及低频分量的对比度来判定有无体动，因此，能够不受拍摄条件或被摄体的种类的影响地稳定地对有无体动进行判定。
另外，在本实施方式中，从放射线图像G0提取被摄体区域A0，并提取被摄体区域A0中低频分量L3′超过阈值Th1的像素，作为用于体动检测的解析点。在此，在拍摄中产生了体动的情况下，在所取得的放射线图像G0中，高频分量的对比度降低，但低频分量的对比度几乎不降低。因此，根据本实施方式，能够不受体动影响地在放射线图像中高精度地设定用于进行检测体动的处理的解析点，其结果是，能够高精度地检测体动。
另外，在上述实施方式中，在解析点设定部23中设定用于判定有无体动的解析点，也可以设为将放射线图像G0显示于显示部40，而使用者通过手动操作设定解析点。另外，也可以将解析点设定成放射线图像G0上的固定位置。
另外，在上述实施方式中，在多分辨率转换部22中对换放射线图像G0进行多分辨率转而取得放射线图像G0的高频分量及低频分量，但也可以对包含解析点的预定范围的区域的放射线图像G0进行傅立叶变换而求出关于解析点的功率谱，并算出功率谱中预定的高频段的值及低频段的值，作为解析点高频分量的对比度及低频分量的对比度。在该情况下，高频分量的对比度与低频分量的对比度的比率实质上是功率谱的高频分量与低频分量的比率。
另外，在上述实施方式中，在使用放射线检测器11拍摄被摄体的放射线图像的拍摄装置10中，使用所取得的放射线图像对体动进行判定，当然，在将被摄体的放射线图像信息存储记录于如日本特开平8-266529号公报、日本特开平9-24039号公报等所示的作为放射线检测体的存储型荧光体片材，并通过从存储型荧光体片材以光电形式读取而取得放射线图像，并使用该放射线图像判定体动时也适用。
Claims (15)
1.一种体动检测装置，其特征在于，具备：
对比度计算单元，算出放射线图像的高频分量的对比度及低频分量的对比度；
比率计算单元，算出所述高频分量的对比度与所述低频分量的对比度的比率；及
判定单元，基于该比率判定所述放射线图像所含的被摄体有无体动。
2.如权利要求1所述的体动检测装置，其特征在于，
所述对比度计算单元根据所述放射线图像生成含有所述高频分量的高频图像及含有所述低频分量的低频图像，并根据该高频图像及该低频图像算出所述高频分量的对比度及所述低频分量的对比度。
3.如权利要求2所述的体动检测装置，其特征在于，
所述对比度计算单元在所述放射线图像所含的、用于算出所述高频分量的对比度及所述低频分量的对比度的一个以上的解析点，算出所述高频分量的对比度及所述低频分量的对比度。
4.如权利要求3所述的体动检测装置，其特征在于，
所述对比度计算单元在所述高频图像和所述低频图像各自的所述解析点，决定以该解析点为中心的解析区域所含的边缘部分的倾斜方向，并算出所述高频图像和所述低频图像各自的该倾斜方向的对比度作为所述高频分量的对比度及所述低频分量的对比度。
5.如权利要求4所述的体动检测装置，其特征在于，
在有多个所述解析点的情况下，所述比率计算单元对应在该多个解析点的各所述解析点决定的每个所述倾斜方向，算出所述高频分量的对比度及所述低频分量的对比度的比率。
6.如权利要求5所述的体动检测装置，其特征在于，
所述判定单元将每个所述倾斜方向的所述比率中的值最小的比率决定为表示所述体动的指标值，通过将该指标值与特定的阈值进行比较来判定有无所述体动。
7.如权利要求3或4所述的体动检测装置，其特征在于，
还具备解析点设定单元，对所述放射线图像设定所述解析点。
8.如权利要求7所述的体动检测装置，其特征在于，
所述解析点设定单元从所述放射线图像提取所述被摄体的区域，并提取该区域中的所述低频分量超过特定的阈值的点作为解析点。
9.如权利要求1至5中任一项所述的体动检测装置，其特征在于，
还具备显示控制单元，将有无所述体动的判定结果显示于显示单元。
10.如权利要求6所述的体动检测装置，其特征在于，
还具备显示控制单元，将有无所述体动的判定结果显示于显示单元。
11.如权利要求10所述的体动检测装置，其特征在于，
所述显示控制单元将表示所述解析点的所述指标值的大小的信息赋予该解析点，并显示所述放射线图像。
12.一种体动检测方法，其特征在于，
算出放射线图像的高频分量的对比度及低频分量的对比度，
算出所述高频分量的对比度与所述低频分量的对比度的比率，
基于该比率判定所述放射线图像所含的被摄体有无体动。
13.一种体动检测装置，其特征在于，具备：
被摄体区域提取单元，从放射线图像提取作为对象的被摄体的区域；
解析点提取单元，提取所述被摄体的区域中的空间频率的低频分量超过特定的阈值的点作为用于进行检测体动用的处理的解析点；
对比度计算单元，在所述解析点算出所述放射线图像的高频分量的对比度及低频分量的对比度；
比率计算单元，算出所述高频分量的对比度与所述低频分量的对比度的比率；
判定单元，基于该比率判定所述放射线图像所含的被摄体有无体动。
14.如权利要求13所述的体动检测装置，其特征在于，
基于所述被摄体的区域中的像素值的标准差来决定所述阈值。
15.一种体动检测方法，其特征在于，
从放射线图像提取作为对象的被摄体的区域，
提取所述被摄体的区域中的空间频率的低频分量超过特定的阈值的点作为用于进行检测体动用的处理的解析点，
在所述解析点算出所述放射线图像的高频分量的对比度及低频分量的对比度，
算出所述高频分量的对比度与所述低频分量的对比度的比率，
基于该比率判定所述放射线图像所含的被摄体有无体动。
Applications Claiming Priority (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012214137A JP5839710B2 (ja) | 2012-09-27 | 2012-09-27 | 解析点設定装置および方法、並びに体動検出装置および方法 |
| JP2012214136 | 2012-09-27 | ||
| JP2012-214137 | 2012-09-27 | ||
| JP2012-214136 | 2012-09-27 | ||
| PCT/JP2013/005557 WO2014050045A1 (ja) | 2012-09-27 | 2013-09-20 | 体動検出装置および方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104684479A CN104684479A (zh) | 2015-06-03 |
| CN104684479B true CN104684479B (zh) | 2017-11-14 |
Family
ID=50741843
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201380049990.1A Active CN104684479B (zh) | 2012-09-27 | 2013-09-20 | 体动检测装置及方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US9619893B2 (zh) |
| EP (1) | EP2901932B1 (zh) |
| JP (1) | JP5839710B2 (zh) |
| CN (1) | CN104684479B (zh) |
| WO (1) | WO2014050045A1 (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6789620B2 (ja) | 2015-10-08 | 2020-11-25 | キヤノン株式会社 | 画像処理装置及びその制御方法、コンピュータプログラム |
| WO2017115432A1 (ja) * | 2015-12-28 | 2017-07-06 | 株式会社島津製作所 | 放射線装置 |
| JP6671267B2 (ja) * | 2016-09-08 | 2020-03-25 | 富士フイルム株式会社 | 画像処理装置、方法およびプログラム |
| TWI656507B (zh) * | 2017-08-21 | 2019-04-11 | 瑞昱半導體股份有限公司 | 電子裝置 |
| CN109427044B (zh) * | 2017-08-25 | 2022-02-25 | 瑞昱半导体股份有限公司 | 电子装置 |
| JP7098341B2 (ja) * | 2018-01-30 | 2022-07-11 | キヤノン株式会社 | 制御装置、放射線撮影システム、制御方法及びプログラム |
| US11980493B2 (en) * | 2020-07-31 | 2024-05-14 | Ricoh Company, Ltd. | Image processing apparatus, imaging system, image processing method, and storage medium |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5617143A (en) * | 1993-03-29 | 1997-04-01 | Canon Kabushiki Kaisha | Movement detection device and encoding apparatus using the same |
| CN101129266A (zh) * | 2006-08-22 | 2008-02-27 | 柯尼卡美能达医疗印刷器材株式会社 | X线图像处理系统 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06284399A (ja) * | 1993-03-29 | 1994-10-07 | Canon Inc | 動き検出装置 |
| JPH0772253A (ja) | 1993-09-01 | 1995-03-17 | Fuji Photo Film Co Ltd | 放射線画像信号読出方法およびそれに用いられる放射線検出器 |
| JPH08266529A (ja) | 1995-03-31 | 1996-10-15 | Fuji Photo Film Co Ltd | 骨塩定量分析方法 |
| JPH0924039A (ja) | 1995-07-10 | 1997-01-28 | Fuji Photo Film Co Ltd | 骨塩定量分析方法および装置 |
| JP3564247B2 (ja) * | 1996-11-26 | 2004-09-08 | 松下電器産業株式会社 | 画像動き補正装置 |
| JP2001299733A (ja) * | 2000-04-27 | 2001-10-30 | Konica Corp | Ｐｃｉ放射線画像処理装置、ｐｃｉ放射線画像検出処理装置、ｐｃｉ放射線画像出力装置及びｐｃｉ画像診断支援装置 |
| JP2005205086A (ja) | 2004-01-26 | 2005-08-04 | Hitachi Medical Corp | Ｘ線診断装置 |
| JP2006034521A (ja) * | 2004-07-26 | 2006-02-09 | Konica Minolta Medical & Graphic Inc | 医用画像処理装置及び医用画像処理方法 |
| US7899229B2 (en) | 2007-08-06 | 2011-03-01 | Hui Luo | Method for detecting anatomical motion blur in diagnostic images |
| JP2010020606A (ja) | 2008-07-11 | 2010-01-28 | Good Medix:Kk | 医療情報仲介サーバ、そのプログラム及び医療情報仲介システム。 |
| JP2010206067A (ja) | 2009-03-05 | 2010-09-16 | Fujifilm Corp | 放射線検出器 |
| JP2010246862A (ja) * | 2009-03-24 | 2010-11-04 | Japan Health Science Foundation | 医用画像作成装置及びプログラム |
-
2012
- 2012-09-27 JP JP2012214137A patent/JP5839710B2/ja active Active
-
2013
- 2013-09-20 WO PCT/JP2013/005557 patent/WO2014050045A1/ja active Application Filing
- 2013-09-20 CN CN201380049990.1A patent/CN104684479B/zh active Active
- 2013-09-20 EP EP13841697.9A patent/EP2901932B1/en active Active
-
2015
- 2015-03-24 US US14/666,468 patent/US9619893B2/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5617143A (en) * | 1993-03-29 | 1997-04-01 | Canon Kabushiki Kaisha | Movement detection device and encoding apparatus using the same |
| CN101129266A (zh) * | 2006-08-22 | 2008-02-27 | 柯尼卡美能达医疗印刷器材株式会社 | X线图像处理系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2901932A1 (en) | 2015-08-05 |
| JP2014064857A (ja) | 2014-04-17 |
| EP2901932A4 (en) | 2016-12-21 |
| CN104684479A (zh) | 2015-06-03 |
| US9619893B2 (en) | 2017-04-11 |
| US20150193948A1 (en) | 2015-07-09 |
| EP2901932B1 (en) | 2018-09-19 |
| WO2014050045A1 (ja) | 2014-04-03 |
| JP5839710B2 (ja) | 2016-01-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104684479B (zh) | 体动检测装置及方法 | |
| US9230322B2 (en) | Image processor, treatment system, and image processing method | |
| US9947101B2 (en) | Radiographic image analysis device and method, and recording medium having program recorded therein | |
| JP5828649B2 (ja) | 画像処理装置、画像処理方法、及びコンピュータプログラム | |
| US20110305381A1 (en) | Radiological image area extracting apparatus, radiological image area extraction program, and radiographic apparatus | |
| US8831323B2 (en) | Method and apparatus for measuring activity of a tracer | |
| US9595116B2 (en) | Body motion detection device and method | |
| JP2018130336A (ja) | 放射線撮影装置、放射線撮影システム、放射線撮影方法、及びプログラム | |
| WO2011092982A1 (ja) | 動態画像処理システム及びプログラム | |
| US20160210740A1 (en) | Method and system for spine position detection | |
| WO2018110623A1 (ja) | 足部判定システム及び方法 | |
| US9867586B2 (en) | Stereo X-ray tube based suppression of outside body high contrast objects | |
| US20200240934A1 (en) | Tomography apparatus and controlling method for the same | |
| JP6184925B2 (ja) | 放射線画像解析装置、方法、及びプログラム | |
| JP6373952B2 (ja) | 放射線画像解析装置および方法並びにプログラム | |
| JP4416823B2 (ja) | 画像処理装置、画像処理方法、及びコンピュータプログラム | |
| US20160300369A1 (en) | Iterative reconstruction with system optics modeling using filters | |
| Jabri et al. | Quantitative assessment of image quality enhancement due to unsharp-mask processing in x-ray fluoroscopy | |
| EP3154026B1 (en) | Image processing apparatus, control method thereof, and computer program | |
| CN117015799A (zh) | 检测x射线图像中的异常 | |
| JP2010172560A (ja) | 放射線画像撮影装置及び画像処理装置 | |
| EP2285282B1 (en) | Spatial orientation method of an immobilized subject | |
| JP2005109908A (ja) | 画像処理装置および画像処理プログラム | |
| Soares et al. | Assessing the impact of correlated noise in digital mammography: a virtual clinical trial | |
| JP6120913B2 (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |