CN110786875B - 放射线图像处理装置及放射线图像处理方法 - Google Patents
放射线图像处理装置及放射线图像处理方法 Download PDFInfo
- Publication number
- CN110786875B CN110786875B CN201910712355.6A CN201910712355A CN110786875B CN 110786875 B CN110786875 B CN 110786875B CN 201910712355 A CN201910712355 A CN 201910712355A CN 110786875 B CN110786875 B CN 110786875B
- Authority
- CN
- China
- Prior art keywords
- group
- metal
- marks
- detector
- image processing
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 230000005855 radiation Effects 0.000 claims abstract description 58
- 239000003550 marker Substances 0.000 claims description 27
- 230000005484 gravity Effects 0.000 claims description 25
- 241000562569 Riodinidae Species 0.000 abstract description 84
- 229910052751 metal Inorganic materials 0.000 abstract description 63
- 239000002184 metal Substances 0.000 abstract description 63
- 238000003384 imaging method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 238000000034 method Methods 0.000 description 12
- 238000002372 labelling Methods 0.000 description 5
- 230000010365 information processing Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 229910052774 Proactinium Inorganic materials 0.000 description 2
- 239000011133 lead Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 1
- 229910052721 tungsten Inorganic materials 0.000 description 1
- 239000010937 tungsten Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T7/00—Details of radiation-measuring instruments
- G01T7/005—Details of radiation-measuring instruments calibration techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- 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/44—Constructional features of apparatus for radiation diagnosis
-
- 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
-
- 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
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
- G06T7/74—Determining position or orientation of objects or cameras using feature-based methods involving reference images or patches
-
- 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/30204—Marker
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Optics & Photonics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
提供一种放射线图像处理装置及放射线图像处理方法。本发明的课题在于,对拍摄时的管球的位置进行高精度推定。本发明的放射线图像处理装置(6)具备:组处理部(62b),根据拍摄图像(P)中的多个金属标识(51a~51d、52a~52d)的图像各自的面积将多个金属标识分类为离检测器(4)相对远的第1组和离检测器(4)相对近的第2组;标识分类部(62c),针对每个经分类的组,根据拍摄图像(P)的图像平面上的多个金属标识(51a~51d、52a~52d)的相对位置对多个金属标识进行分类;及对处理部(62d),将相对位置匹配的第1组的金属标识(51a~51d)和第2组的金属标识(52a~52d)配对选择。
Description
本申请主张基于2018年8月2日申请的日本专利申请第2018-145986号的优先权。其申请的全部内容通过参考援用于本说明书中。
技术领域
本发明涉及一种推定作为放射线源的管球的位置的放射线图像处理装置及放射线图像处理方法。
背景技术
进行断层合成摄影的放射线图像摄影装置通过变更作为放射线源的管球的位置的同时进行拍摄来获取多个图像的信号。该放射线图像摄影装置与被摄体一并拍摄具有成对的金属标识的体模,并根据显示于图像的1对金属标识的位置来推定拍摄时的管球的位置(例如,参考专利文献1)。
专利文献1:日本特开2006-181252号公报
如专利文献1中记载的以往的放射线图像摄影装置那样根据1对金属标识的位置来推定管球位置的结构中,有可能因载置被测体的顶板倾斜等而产生金属标识未被拍摄到拍摄图像中的状态。其结果,存在无法对拍摄时的管球的位置进行高精度推定这一问题。
发明内容
本发明是鉴于上述情况而完成的,其目的在于,提供一种能够对拍摄时的管球的位置进行高精度推定的放射线图像处理装置及放射线图像处理方法。
为了实现上述目的,本发明中的放射线图像处理装置具备:检测器,以与被保持为能够移动的作为放射线源的管球对置的方式配置,且将通过从该管球照射的放射线拍摄的拍摄图像转换成图像数据;及体模,设置于所述管球与所述检测器之间,且在离所述检测器相对远的位置具有作为第1组的多个标识,在离所述检测器相对近的位置具有作为第2组的多个标识,所述放射线图像处理装置还具备:判定机构,根据所述拍摄图像中的多个标识的图像各自的面积,判定该拍摄图像中的多个标识相当于所述第1组的标识和所述第2组的标识中的哪一个;及对处理机构,根据通过所述判定机构判定的所述第1组和所述第2组各自中的多个标识的相对位置关系,将所述第1组的标识和所述第2组的标识配对选择。
本发明的放射线图像处理装置中,通过具有如上所述的判定机构及对处理机构,能够在相对于检测器的远近方向上划分拍摄图像中的多个标识的基础上配对选择。由此,通过根据从拍摄图像中的多个标识中配对选择的标识的位置来推定拍摄时的管球的位置,例如即使在体模倾斜的情况下,也能够选择适当的对,从而能够准确地推定拍摄时的管球的位置。
上述放射线图像处理装置中,优选还具备根据通过对处理机构配对选择的标识的位置坐标推定管球的位置的管球位置推定处理机构。由此,能够对拍摄时的管球的位置进行高精度推定。并且,优选以如下区域的中心点为重心,并根据离重心的位置对多个标识进行分类,该区域被在拍摄图像的图像平面上的一方向上位于最外部的标识和在与该一方向正交的方向上位于最外部的标识包围。由此,拍摄图像中,即使多个标识的图像的一部分缺失,也能够进行对选择。并且,优选多个标识包含能够由第1组的标识和第2组的标识构成至少2个对的数量的标识。由此,能够利用最小限度的标识来进行对选择。
本发明中的放射线图像处理方法为放射线图像处理装置的放射线图像处理方法,所述放射线图像处理装置具备:检测器,以与被保持为能够移动的作为放射线源的管球对置的方式配置,且将通过从该管球照射的放射线拍摄的拍摄图像转换成图像数据;及体模,设置于所述管球与所述检测器之间,且在离所述检测器相对远的位置具有作为第1组的多个标识,在离所述检测器相对近的位置具有作为第2组的多个标识,所述放射线图像处理方法包括:
判定工序,根据所述拍摄图像中的多个标识的图像各自的面积,判定该拍摄图像中的多个标识相当于所述第1组的标识和所述第2组的标识中的哪一个;及对处理工序,根据通过所述判定工序判定的所述第1组和所述第2组各自中的多个标识的相对位置关系,将所述第1组的标识和所述第2组的标识配对选择。
本发明的放射线图像处理方法中,通过在相对于检测器的远近方向上划分拍摄图像中的多个标识的基础上配对选择,由此根据从拍摄图像中的多个标识中配对选择的标识的位置来推定拍摄时的管球的位置,例如即使在体模倾斜的情况下,也能够选择适当的对,从而能够准确地推定拍摄时的管球的位置。
上述放射线图像处理方法中,优选还具备根据通过对处理工序配对选择的标识的位置坐标推定管球的位置的管球位置推定处理工序。由此,能够对拍摄时的管球的位置进行高精度推定。并且,优选以如下区域的中心点为重心,并根据离重心的位置对多个标识进行分类,该区域被在拍摄图像的图像平面上的一方向上位于最外部的标识和在与该一方向正交的方向上位于最外部的标识包围。由此,拍摄图像中,即使多个标识的图像的一部分缺失,也能够进行对选择。并且,优选多个标识包含能够由第1组的标识和第2组的标识构成至少2个对的数量的标识。由此,能够利用最小限度的标识来进行对选择。
发明效果
根据本发明,能够对拍摄时的管球的位置进行高精度推定。
附图说明
图1是表示具备本发明的实施方式所涉及的放射线图像处理装置的放射线图像摄影装置的整体结构的示意图。
图2是表示放射线图像摄影装置所具备的体模的示意图。
图3是表示放射线图像处理装置的功能结构的框图。
图4的(A)~(C)是对放射线图像处理装置所具备的组处理机构进行的组处理进行说明的示意图。
图5的(A)及(B)是对放射线图像处理装置所具备的标识分类机构进行的标识分类进行说明的示意图。
图6的(A)~(C)是对获取放射线图像处理装置中的管球的位置的处理进行说明的示意图。
图7是对放射线图像处理装置进行的处理的流程进行说明的流程图。
符号说明
1-放射线图像摄影装置,2-管球,3-位置变更机构,4-检测器,5-体模,50-金属标识,51a、51b、51c、51d-金属标识(标识),52a、52b、52c、52d-金属标识(标识),6-放射线图像处理装置,60-信号获取部,61-图像生成部,62-标识信息处理部,62a-标签处理部,62b-组处理部(组处理机构),62c-标识分类部(标识分类机构),62d-对处理部(对处理机构),63-管球位置推定处理部,7-摄影控制部,T-被摄体,P-拍摄图像,G1、G2-重心,S1-制作二值化图像的工序,S2-加标签的工序,S3-计算各区域的面积的工序,S4-根据面积判定标识的远近的工序,S5-决定标识对的工序,S6-推定管球坐标的工序。
具体实施方式
以下,根据附图对将本发明具体化的实施方式进行说明。
参考图1~图7,对具备本发明的实施方式所涉及的放射线图像处理装置6的放射线图像摄影装置1进行说明。
首先,参考图1~图3,对具备本发明的实施方式所涉及的放射线图像处理装置6的放射线图像摄影装置1的结构进行说明。图1是表示具备放射线图像处理装置6的放射线图像摄影装置1的整体结构的示意图。图2是表示放射线图像摄影装置1所具备的体模5的示意图。图3是表示放射线图像处理装置6的功能结构的框图。
图1所示的放射线图像摄影装置1为进行医用断层合成摄影的装置,其通过变更作为放射线源的管球2的位置的同时拍摄被摄体T来获取多个图像数据。具体而言,放射线图像摄影装置1具备管球2、位置变更机构3、检测器4、体模5、放射线图像处理装置6、摄影控制部7等。
管球2根据来自摄影控制部7的信号施加高电压,由此产生放射线(X射线)并朝向检测器4照射该放射线。该管球2被位置变更机构3保持为能够移动。位置变更机构3根据来自摄影控制部7的信号变更管球2的位置。
检测器4为以与管球2对置的方式配置,且将通过从管球2照射的放射线拍摄的拍摄图像转换成图像数据的FPD(Flat Panel Detector,平板探测器)。即,检测器4将放射线转换成电信号,并且读取所转换的电信号作为图像的信号来向放射线图像处理装置6输出该图像的信号。另外,检测器4具备多个转换元件(省略图示)和配置于这些多个转换元件上的像素电极(省略图示)。并且,多个转换元件及像素电极以规定的周期(像素间距)配置。
体模5设置于管球2与检测器4之间,且与被摄体T一并进行拍摄以推定管球2的位置。如图2所示,体模5由树脂等构成,其在内部具有由铝、金、铅、钨等金属构成的多个金属标识50、51a、51b、51c、51d、52a、52b、52c、52d。金属标识51a及金属标识52a在相对于检测器4的远近方向上划分配置而成对。金属标识51b及金属标识52b在相对于检测器4的远近方向上划分配置而成对。金属标识51c及金属标识52c在相对于检测器4的远近方向上划分配置而成对。金属标识51d及金属标识52d在相对于检测器4的远近方向上划分配置而成对。
在此,成对的金属标识在远近方向上彼此至少远离70mm而配置。并且,成对的金属标识在从远近方向观察时(俯视体模5时)配置于不重合的位置。
返回到图1进行说明。放射线图像处理装置6为对通过检测器4获得的图像的信号进行处理的装置。该放射线图像处理装置6具有CPU(Central Processi ng Unit,中央处理单元)、GPU(Graphics Processing Unit,图形处理单元)或构成为用于图像处理的FPGA(Field-Programmable Gate Array,现场可编程门阵列)等处理器(省略图示)。
如图3所示,放射线图像处理装置6通过执行放射线图像处理程序,实现信号获取部60、图像生成部61、标识信息处理部62、管球位置推定处理部63等各种功能。
信号获取部60获取通过检测器4(参考图1)获得的、拍摄被摄体T和体模5而得的图像的信号,并向图像生成部61输出所获取的图像的信号。
图像生成部61根据由信号获取部60输出的图像的信号,例如生成根据预先决定的像素值的阈值进行二值化的拍摄图像P(参考图4的(A))。图4的(A)中示出了在拍摄图像P中描绘出9个金属标识的状况,且省略了被摄体T的描绘。
标识信息处理部62通过图像识别处理通过图像生成部61生成的拍摄图像P(参考图4的(A))中的金属标识51a~51d、52a~52d(参考图4的(A))的信息。具体而言,标识信息处理部62具备标签处理部62a、组处理部62b、标识分类部62c、对处理部62d等。
标签处理部62a分别对拍摄图像P(参考图4的(A))中的多个金属标识51a~51d、52a~52d(参考图4的(A))进行加标签处理,从而彼此区分各金属标识的图像和其他图像。
组处理部62b根据通过标签处理部62a进行加标签的拍摄图像P中的多个金属标识51a~51d、52a~52d(参考图4的(A))各自的面积,将拍摄图像P中的金属标识51a~51d、52a~52d(参考图4的(A))分类为离检测器4(参考图1)相对远的第1组和离检测器4(参考图1)相对近的第2组。
具体而言,组处理部62b在计算拍摄图像P中的多个金属标识51a~51d、52a~52d(参考图4的(A))各自的面积之后,计算所计算出的面积的最大值和最小值的平均值。并且,组处理部62b在拍摄图像P中将具有大于平均值的面积的金属标识51a~51d(参考图4的(B))判定为离检测器4(参考图1)相对远(在体模5内位于上部。)而分类为第1组,并且将具有小于平均值的面积的金属标识52a~52d(参考图4的(C))判定为离检测器4(参考图1)相对近(在体模5内位于下部。)而分类为第2组。
标识分类部62c针对每个经分类的组根据拍摄图像P的图像平面上的多个金属标识51a~51d(参考图5的(A))、52a~52d(参考图5的(B))的相对位置来对该多个金属标识51a~51d(参考图5的(A))、52a~52d(参考图5的(B))进行分类。
具体而言,标识分类部62c以如下区域(被金属标识51a~51d包围的区域)的中心点为第1组的金属标识51a~51d(参考图5的(A))的坐标的重心G1(参考图5的(A)),并根据离该重心G1(参考图5的(A))的位置对第1组的金属标识51a~51d(参考图5的(A))进行分类,该区域被拍摄图像P(参考图5的(A))的图像平面上的一方向(x方向)上位于最外部的金属标识和在与该一方向正交的方向(y方向)上位于最外部的金属标识包围。
如图5的(A)所示,金属标识51a被分类为位于坐标的重心G1的左上方。金属标识51b被分类为位于坐标的重心G1的右上方。金属标识51c被分类为位于坐标的重心G1的右下方。金属标识51d被分类为位于坐标的重心G1的左下方。
返回到图3进行说明。标识分类部62c以如下区域(被金属标识52a~52d包围的区域)的中心点为第2组的金属标识52a~52d(参考图5的(B))的坐标的重心G2(参考图5的(B)),并根据离该重心G2(参考图5的(B))的位置对第2组的金属标识52a~52d(参考图5的(B))进行分类,该区域被在拍摄图像P(参考图5的(B))的图像平面上的一方向(x方向)上位于最外部的金属标识和在与该一方向正交的方向(y方向)上位于最外部的金属标识包围。
如图5的(B)所示,金属标识52a被分类为位于坐标的重心G2的左上方。金属标识52b被分类为位于坐标的重心G2的右上方。金属标识52c被分类为位于坐标的重心G2的右下方。金属标识52d被分类为位于坐标的重心G2的左下方。
返回到图3进行说明。对处理部62d将相对位置匹配的第1组的金属标识51a~51d(参考图5的(A))和第2组的金属标识52a~52d(参考图5的(B))配对选择。
具体而言,对处理部62d将位于坐标的重心G1(参考图5的(A))的左上方的金属标识51a(参考图5的(A))和位于坐标的重心G2(参考图5的(B))的左上方的金属标识52a(参考图5的(B))配对选择。并且,对处理部62d将位于坐标的重心G1(参考图5的(A))的右上方的金属标识51b(参考图5的(A))和位于坐标的重心G2(参考图5的(B))的右上方的金属标识52b(参考图5的(B))配对选择。
并且,对处理部62d将位于坐标的重心G1(参考图5的(A))的右下方的金属标识51c(参考图5的(A))和位于坐标的重心G2(参考图5的(B))的右下方的金属标识52c(参考图5的(B))配对选择。而且,对处理部62d将位于坐标的重心G1(参考图5的(A))的左下方的金属标识51d(参考图5的(A))和位于坐标的重心G2(参考图5的(B))的左下方的金属标识52d(参考图5的(B))配对选择。
管球位置推定处理部63根据通过对处理部62d配对选择的金属标识51a~51d(参考图5的(A))、52a~52d(参考图5的(B))的位置坐标来推定管球2(参考图1)的位置。金属标识的位置坐标例如以FPD(检测器4)的中心为坐标的基准。
接着,参考图6以成对的金属标识51a、52a为例对获取管球2的位置信息的处理进行说明。图6的(A)表示管球2及金属标识51a、52a的位置关系的示意图。图6的(B)是拍摄图像P的示意图。图6的(C)是以矢量图示出管球2、金属标识51a、52a及拍摄图像P中的金属标识51a、52a的位置关系例。另外,为了便于说明,图6的(A)中未图示体模5。
如图6的(A)所示,在将管球2配置于规定的位置并从倾斜方向照射放射线的情况下,透过金属标识51a及金属标识52a的放射线在检测器4上到达不同的点。因此,如图6的(B)所示,拍摄图像P中的金属标识51a、52a描绘于不同的位置。若将管球2的位置设为S,将金属标识51a的位置设为M1,将金属标识52a的位置设为M2,将拍摄图像P中的金属标识51a的位置设为I1,将拍摄图像P中的金属标识52a的位置设为I2,则能够获得如图6的(C)所示的矢量图。
如图6的(C)所示,管球2、金属标识51a、52a及拍摄图像P中的金属标识51a、52a成为外分点的关系。即,拍摄图像P中的金属标识51a为对由管球2和金属标识51a构成的线段SM1以t1:(1-t1)的比率进行外分的点。并且,拍摄图像P中的金属标识52a为对由管球2和金属标识52a构成的线段SM2以t2:(1-t2)的比率进行外分的点。根据上述关系,获得以下式(1)及式(2)。
[数1]
在此,将管球2的位置S的位置坐标定义为(x,y,Sd)。并且,将金属标识51a的位置M1的位置坐标定义为(Pa,Pb,Pd+Ps)。并且,将金属标识52a的位置M2的位置坐标定义为(Pa,Pb,Pd)。并且,将拍摄图像P中的金属标识51a的位置I1的位置坐标定义为(a1,b1,0)。并且,将拍摄图像P中的金属标识52a的位置I2的位置坐标定义为(a2,b2,0)。
另外,x为管球2的X方向上的坐标。并且,y为管球2的Y方向上的坐标。并且,Pa为金属标识51a、52a的X方向上的坐标。并且,Pb为金属标识51a、52a的Y方向上的坐标。并且,Sd为检测器4至管球2为止的Z方向上的距离(SID:Source Image receptor Distance)。并且,Pd为检测器4至金属标识52a为止的Z方向上的距离。并且,Ps为金属标识51a、52a彼此之间的Z方向上的距离。
根据管球2的位置坐标、金属标识51a、52a的位置坐标、拍摄图像P中的金属标识51a、52a的位置坐标、上述式(1)及式(2)获得以下式(3)~(8)。
[数2]
x*t1+Pa*(1-t1)=a1…(3)
x*t2+Pa*(1-t2)=a2…(4)
y*t1+Pa*(1-t1)=b1…(5)
y*t2+Pa*(1-t2)=b2…(6)
Sd*t1+(Pd+Ps)*(1-t1)=0…(7)
Sd*t2+Pd*(1-t2)=0…(8)
在上述式(3)~式(8)中,Sd、Pd、Ps分别为已知的值,因此未知值和式的数量变得相等,从而能够获取管球2的位置信息。具体而言,根据式(7)及式(8)获得以下式(9)及式(10)。
[数3]
在此,若将上述式(9)的解设为t1=α,将上述式(10)的解设为t2=β,则根据上述式(1)及式(2)获得以下式(11)。并且,根据上述式(3)及式(4)获得以下式(12)。
[数4]
上述式(11)的a1及a2通过获取拍摄图像P中的金属标识51a、52a的X坐标的值而得。并且,上述式(12)的b1及b2通过获取拍摄图像P中的金属标识51a、52a的Y坐标的值而得。因此,放射线图像处理装置6能够根据拍摄图像P中的金属标识51a、52a的坐标值和上述式(11)及式(12)来获取拍摄拍摄图像P时的管球2的位置信息。
(管球的位置的推定)
接着,参考图7对放射线图像处理装置6进行的处理进行说明。图7是对放射线图像处理装置6进行的管球位置推定处理的流程进行说明的流程图。
如图7所示,放射线图像处理装置6依次执行制作二值化图像的工序S1、加标签的工序S2、计算各区域的面积的工序S3、根据面积判定标识的远近的工序S4、决定标识对的工序S5及推定管球坐标的工序S6。
制作二值化图像的工序S1为图像生成部61根据通过检测器4检测出的图像的信号生成经二值化的图像P的工序。
加标签的工序S2为标签处理部62a分别对拍摄图像P中的金属标识51a~51d、52a~52d进行加标签来进行彼此区分的工序。在此,例如,将经二值化的拍摄图像P的中央部分的矩形区域作为关心区域(关注区域),并针对该关心区域的加标签处理,由此提取具有相同的像素值的像素(相邻像素具有相同的像素值的像素)相连的区域来分配相同的标签,从而区分包括金属标识的体模的图像区域和被摄体T的图像区域。也可以根据构成关心区域的像素的像素值的平均值提取特征量来识别被摄体T和金属标识。
计算各区域的面积的工序S3为由组处理部62b计算通过标签处理部62a进行加标签的拍摄图像P中的多个金属标识51a~51d、52a~52d各自的面积的工序。该计算各区域的面积的工序S3中,组处理部62b还一并计算所计算出的面积的最大值和最小值的平均值。
根据面积判定标识的远近的工序S4为由组处理部62b将上述计算出的平均值作为阈值,并将具有大于平均值的面积的拍摄图像P中的金属标识51a~51d判定为离检测器4相对远(在图2的体模5内位于上部。)而分类为第1组,并且将具有小于平均值的面积的拍摄图像P中的金属标识52a~52d判定为离检测器4相对近(在图2的体模5内位于下部。)而分类为第2组的工序。
决定标识对的工序S5为在由标识分类部62c针对每个通过组处理部62b进行分类的组根据多个金属标识51a~51d、52a~52d的xy坐标平面上的相对位置对多个金属标识51a~51d、52a~52d进行分类的基础上,由对处理部62d将相对位置匹配的第1组的金属标识51a~51d和第2组的金属标识52a~52d配对选择的工序。
具体而言,如参考图5的(A)和图5的(B)进行说明,例如选择金属标识51a和金属标识52a的对、金属标识51b和金属标识52b的对、金属标识51c和金属标识52c的对、金属标识51d和金属标识52d的对中的任一对作为相对位置匹配的金属标识。
另外,体模5中,作为在相对于检测器4的远近方向上划分配置而成对的多个金属标识,在此,例示了构成4个对的数量,但并不限定于此。例如,即使考虑因体模5倾斜等而一部分金属标识未被拍摄到拍摄图像中的情况,管球位置的推定中只要具备构成至少2个对的数量的金属标识即可。并且,标识并不限定于金属制,只要X射线的吸收量多,则可以是任何原材料。
推定管球坐标的工序S6为由管球位置推定处理部63通过参考图6的(A)、图6的(B)及图6的(C)进行说明的方法,根据通过对处理部62d配对选择的金属标识51a~51d、52a~52d的位置坐标来推定管球2的位置的工序。
(实施方式的效果)
本发明的实施方式中,能够获得如下效果。
如上所述,本实施方式中,放射线图像处理装置6具备:检测器4,以与保持为能够移动的作为放射线源的管球2对置的方式配置,且将通过从管球2照射的放射线拍摄的拍摄图像P转换成图像数据;及体模5,设置于管球2与检测器4之间,且具有相对于检测器4在远近方向上划分配置的多个金属标识51a~51d、52a~52d,所述放射线图像处理装置还具备:组处理部62b,根据拍摄图像P中的多个金属标识51a~51d、52a~52d的图像各自的面积,将多个金属标识51a~51d、52a~52d分类为离检测器4相对远的第1组和离检测器4相对近的第2组;标识分类部62c,针对每个经分类的组,根据拍摄图像P的图像平面上的多个金属标识51a~51d、52a~52d的相对位置对多个金属标识51a~51d、52a~52d进行分类;及对处理部62d,将相对位置匹配的第1组的金属标识51a~51d和第2组的金属标识52a~52d配对选择。
由此,能够在相对于检测器4的远近方向上划分拍摄图像P中的多个金属标识51a~51d、52a~52d的基础上配对选择。由此,能够根据拍摄图像P中的多个对的金属标识51a~51d、52a~52d的位置来推定拍摄时的管球2的位置,因此例如即使在体模5倾斜的情况下,也能够准确地推定拍摄时的管球2的位置。即,能够高精度地进行拍摄时的管球2的位置的推定。
并且,本实施方式中,放射线图像处理装置6具备根据通过对处理部62d配对选择的金属标识51a~51d、52a~52d的位置坐标来推定管球2的位置的管球位置推定处理部63。
并且,本实施方式中,标识分类部62c以如下区域的中心点为重心G1、G2,并根据离重心G1、G2的位置对多个金属标识51a~51d、52a~52d进行分类,该区域被在拍摄图像P的图像平面上的一方向上位于最外部的金属标识和在与该一方向正交的方向上位于最外部的金属标识包围。
并且,本实施方式中,多个金属标识51a~51d、52a~52d包含能够由第1组的金属标识51a~51d和第2组的金属标识52a~52d构成至少2个对的数量的金属标识。
并且,本实施方式中的放射线图像处理方法中,放射线图像处理装置6具备:检测器4,以与保持为能够移动的作为放射线源的管球2对置的方式配置,且将通过从管球2照射的放射线拍摄的拍摄图像P转换成图像数据;及体模5,设置于管球2与检测器4之间,且具有相对于检测器4在远近方向上划分配置的多个金属标识51a~51d、52a~52d,放射线图像处理方法包括:组处理工序,根据拍摄图像P中的多个金属标识51a~51d、52a~52d的图像各自的面积,将多个金属标识51a~51d、52a~52d分类为离检测器4相对远的第1组和离检测器4相对近的第2组;标识分类工序,针对每个经分类的组,根据拍摄图像P的图像平面上的多个金属标识51a~51d、52a~52d的相对位置对多个金属标识51a~51d、52a~52d进行分类;及对处理工序,将相对位置匹配的第1组的金属标识51a~51d和第2组的金属标识52a~52d配对选择。
并且,本实施方式中,包括根据通过对处理工序配对选择的金属标识51a~51d、52a~52d的位置坐标来推定管球2的位置的管球位置推定处理工序。
并且,本实施方式中,标识分类工序中以如下区域的中心点为重心G1、G2,并根据离重心G1、G2的位置对多个金属标识51a~51d、52a~52d进行分类,该区域被在拍摄图像P的图像平面上的一方向上位于最外部的金属标识和在与该一方向正交的方向上位于最外部的金属标识包围。
并且,本实施方式中,多个金属标识51a~51d、52a~52d包含能够由第1组的金属标识51a~51d和第2组的金属标识52a~52d构成至少2个对的数量的金属标识51a~51d、52a~52d构成。
(变形例)
另外,应认为上述实施方式在所有方面是例示,而并不限制本发明。本发明的范围由技术方案来表示,而不是上述实施方式的说明,而且还包括与技术方案等同的含义及范围内的所有变更(变形例)。
例如,上述实施方式中,体模5具有金属标识50,但由于金属标识50确认放射线图像处理装置6中的处理是否准确,或者用于图像的重构等,因此体模5也可以不具有金属标识50。
Claims (6)
1.一种放射线图像处理装置,其具备:检测器,以与被保持为能够移动的作为放射线源的管球对置的方式配置,且将通过从该管球照射的放射线拍摄的拍摄图像转换成图像数据;及体模,设置于所述管球与所述检测器之间,且在离所述检测器相对远的位置具有作为第1组的多个标识,在离所述检测器相对近的位置具有作为第2组的多个标识,所述放射线图像处理装置的特征在于,还具备:
判定机构,根据所述拍摄图像中的多个标识的图像各自的面积,判定该拍摄图像中的多个标识相当于所述第1组的标识和所述第2组的标识中的哪一个;
对处理机构,根据通过所述判定机构判定的所述第1组和所述第2组各自中的多个标识的相对位置关系,将所述第1组的标识和所述第2组的标识配对选择;及
管球位置推定处理机构,根据通过所述对处理机构配对选择的标识的位置坐标来推定所述管球的位置。
2.根据权利要求1所述的放射线图像处理装置,其特征在于,
以如下区域的中心点为重心,并根据离该重心的位置对所述多个标识进行分类,该区域被在所述拍摄图像的图像平面上的一方向上位于最外部的标识和在与该一方向正交的方向上位于最外部的标识包围。
3.根据权利要求1或2所述的放射线图像处理装置,其特征在于,
所述多个标识包含能够由所述第1组的标识和所述第2组的标识构成至少2个对的数量的标识。
4.一种放射线图像处理装置的放射线图像处理方法,所述放射线图像处理装置具备:检测器,以与被保持为能够移动的作为放射线源的管球对置的方式配置,且将通过从该管球照射的放射线拍摄的拍摄图像转换成图像数据;及体模,设置于所述管球与所述检测器之间,且在离所述检测器相对远的位置具有作为第1组的多个标识,在离所述检测器相对近的位置具有作为第2组的多个标识,所述放射线图像处理方法的特征在于,包括:
判定工序,根据所述拍摄图像中的多个标识的图像各自的面积,判定该拍摄图像中的多个标识相当于所述第1组的标识和所述第2组的标识中的哪一个;
对处理工序,根据通过所述判定工序判定的所述第1组和所述第2组各自中的多个标识的相对位置关系,将所述第1组的标识和所述第2组的标识配对选择;及
管球位置推定处理工序,根据通过所述对处理工序配对选择的标识的位置坐标来推定所述管球的位置。
5.根据权利要求4所述的放射线图像处理方法,其特征在于,
以如下区域的中心点为重心,并根据离该重心的位置对所述多个标识进行分类,该区域被在所述拍摄图像的图像平面上的一方向上位于最外部的标识和在与该一方向正交的方向上位于最外部的标识包围。
6.根据权利要求4或5所述的放射线图像处理方法,其特征在于,
所述多个标识包含能够由所述第1组的标识和所述第2组的标识构成至少2个对的数量的标识。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018-145986 | 2018-08-02 | ||
JP2018145986A JP6988732B2 (ja) | 2018-08-02 | 2018-08-02 | 放射線画像処理装置および放射線画像処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110786875A CN110786875A (zh) | 2020-02-14 |
CN110786875B true CN110786875B (zh) | 2024-03-12 |
Family
ID=69427011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910712355.6A Active CN110786875B (zh) | 2018-08-02 | 2019-08-02 | 放射线图像处理装置及放射线图像处理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11353606B2 (zh) |
JP (1) | JP6988732B2 (zh) |
CN (1) | CN110786875B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005021345A (ja) * | 2003-07-01 | 2005-01-27 | Toshiba Corp | X線立体再構成処理装置、x線撮影装置、x線立体再構成処理方法及びx線立体撮影補助具 |
CN1573321A (zh) * | 2003-06-10 | 2005-02-02 | 株式会社岛津制作所 | X光照相设备 |
JP2006043047A (ja) * | 2004-08-03 | 2006-02-16 | Toshiba Corp | X線立体撮影装置 |
EP1760457A2 (de) * | 2005-07-13 | 2007-03-07 | Carl Zeiss Industrielle Messtechnik GmbH | Verfahren und eine Anordnung zum Kalibrieren einer Messanordnung |
CN102253061A (zh) * | 2011-04-19 | 2011-11-23 | 东南大学 | 立式锥束计算机断层成像校准系统及应用该系统的方法 |
CN103635139A (zh) * | 2011-07-06 | 2014-03-12 | 富士胶片株式会社 | X射线成像设备和用于其的校准方法 |
CN106667515A (zh) * | 2012-11-05 | 2017-05-17 | 上海联影医疗科技有限公司 | Ct扫描仪散焦校正方法 |
Family Cites Families (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2487044C (fr) * | 2002-05-23 | 2010-05-11 | Centre National De La Recherche Scientifique (C.N.R.S.) | Dispositif de stereoradiographie et procedure d'utilisation |
JP4415762B2 (ja) * | 2003-06-10 | 2010-02-17 | 株式会社島津製作所 | 断層撮影装置 |
US6888924B2 (en) * | 2003-06-25 | 2005-05-03 | General Electric Company | Method, apparatus, and medium for calibration of tomosynthesis system geometry using fiducial markers with non-determined position |
JP4711612B2 (ja) * | 2003-09-11 | 2011-06-29 | 株式会社日立メディコ | 定位脳手術支援システム |
US7016456B2 (en) * | 2003-10-31 | 2006-03-21 | General Electric Company | Method and apparatus for calibrating volumetric computed tomography systems |
US7177386B2 (en) * | 2004-03-15 | 2007-02-13 | Varian Medical Systems Technologies, Inc. | Breathing synchronized computed tomography image acquisition |
JP4479503B2 (ja) | 2004-12-28 | 2010-06-09 | 株式会社島津製作所 | 断層撮影装置 |
US7950849B2 (en) * | 2005-11-29 | 2011-05-31 | General Electric Company | Method and device for geometry analysis and calibration of volumetric imaging systems |
US8121368B2 (en) * | 2008-02-19 | 2012-02-21 | The Board Of Trustees Of The Leland Stanford Junior University | 3D real-time tracking of human anatomy using combined kV and MV imaging |
EP2245986B1 (de) * | 2008-08-22 | 2013-10-16 | BrainLAB AG | Pyramidenförmige Röntgenmarkervorrichtung |
JP5493360B2 (ja) * | 2009-01-08 | 2014-05-14 | オムロン株式会社 | X線検査方法、x線検査装置およびx線検査プログラム |
CA2655001C (en) * | 2009-02-20 | 2015-11-24 | Queen's University At Kingston | Marker localization using intensity-based registration of imaging modalities |
JP2011019707A (ja) * | 2009-07-15 | 2011-02-03 | Fujifilm Corp | X線撮影装置、x線撮影装置の制御方法、及びプログラム |
US9082177B2 (en) * | 2009-11-25 | 2015-07-14 | Dental Imaging Technologies Corporation | Method for tracking X-ray markers in serial CT projection images |
US8363919B2 (en) * | 2009-11-25 | 2013-01-29 | Imaging Sciences International Llc | Marker identification and processing in x-ray images |
WO2012007036A1 (en) * | 2010-07-14 | 2012-01-19 | Brainlab Ag | Method and system for determining an imaging direction and calibration of an imaging apparatus |
US8948842B2 (en) * | 2011-01-21 | 2015-02-03 | Headwater Partners Ii Llc | Radiation treatment with multiple imaging elements |
DE102011007794B4 (de) * | 2011-04-20 | 2019-05-23 | Siemens Healthcare Gmbh | Verfahren zur geometrisch korrekten Zuordnung von 3D-Bilddaten eines Patienten |
DE102012200893A1 (de) * | 2012-01-23 | 2013-07-25 | Siemens Aktiengesellschaft | Optische Justagevorrichtung |
US20130229495A1 (en) * | 2012-03-01 | 2013-09-05 | Ali-Reza Bani-Hashemi | Method for calibrating an imaging system |
US9183354B2 (en) * | 2012-08-15 | 2015-11-10 | Musc Foundation For Research Development | Systems and methods for image guided surgery |
CN104783824B (zh) * | 2014-01-20 | 2020-06-26 | 上海联影医疗科技有限公司 | X射线成像系统的校正方法 |
US10420626B2 (en) * | 2015-06-30 | 2019-09-24 | Canon U.S.A., Inc. | Fiducial markers, systems, and methods of registration |
US10004564B1 (en) * | 2016-01-06 | 2018-06-26 | Paul Beck | Accurate radiographic calibration using multiple images |
JP6732489B2 (ja) * | 2016-03-18 | 2020-07-29 | 株式会社東芝 | X線撮影機器のアライメント調整支援装置、方法及びプログラム |
CN107468350B (zh) * | 2016-06-08 | 2020-12-08 | 北京天智航医疗科技股份有限公司 | 一种三维图像专用标定器、手术定位系统及定位方法 |
US10258302B2 (en) * | 2017-04-13 | 2019-04-16 | Apn Health, Llc | Rapid 3D cardiac parameter mapping |
EP3421086B1 (en) * | 2017-06-28 | 2020-01-15 | OptiNav Sp. z o.o. | Determination of geometrical information about a medical treatment arrangement comprising a rotatable treatment radiation source unit |
EP3843038B1 (en) * | 2019-12-23 | 2023-09-20 | HTC Corporation | Image processing method and system |
-
2018
- 2018-08-02 JP JP2018145986A patent/JP6988732B2/ja active Active
-
2019
- 2019-08-02 CN CN201910712355.6A patent/CN110786875B/zh active Active
- 2019-08-02 US US16/530,132 patent/US11353606B2/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1573321A (zh) * | 2003-06-10 | 2005-02-02 | 株式会社岛津制作所 | X光照相设备 |
JP2005021345A (ja) * | 2003-07-01 | 2005-01-27 | Toshiba Corp | X線立体再構成処理装置、x線撮影装置、x線立体再構成処理方法及びx線立体撮影補助具 |
JP2006043047A (ja) * | 2004-08-03 | 2006-02-16 | Toshiba Corp | X線立体撮影装置 |
EP1760457A2 (de) * | 2005-07-13 | 2007-03-07 | Carl Zeiss Industrielle Messtechnik GmbH | Verfahren und eine Anordnung zum Kalibrieren einer Messanordnung |
CN102253061A (zh) * | 2011-04-19 | 2011-11-23 | 东南大学 | 立式锥束计算机断层成像校准系统及应用该系统的方法 |
CN103635139A (zh) * | 2011-07-06 | 2014-03-12 | 富士胶片株式会社 | X射线成像设备和用于其的校准方法 |
CN106667515A (zh) * | 2012-11-05 | 2017-05-17 | 上海联影医疗科技有限公司 | Ct扫描仪散焦校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110786875A (zh) | 2020-02-14 |
JP6988732B2 (ja) | 2022-01-05 |
US20200064502A1 (en) | 2020-02-27 |
JP2020018692A (ja) | 2020-02-06 |
US11353606B2 (en) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3992851A1 (en) | Image classification method, apparatus and device, storage medium, and medical electronic device | |
CN107909010B (zh) | 一种道路障碍物检测方法与装置 | |
EP2803015B1 (en) | Two stage categorization of objects in images | |
CN106254860B (zh) | 照相机的校准装置、照相机系统、以及照相机的校准方法 | |
JP5549345B2 (ja) | 画像収集装置に用いる天空検出装置及び方法 | |
US8699786B2 (en) | 3D model generating apparatus, method and CRM by line pattern imaging | |
CN111103317A (zh) | 扫描晶片的系统及方法 | |
JP6736414B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
KR20140142368A (ko) | 결함 판정 장치, 방사선 촬상 시스템, 및 결함 판정 방법 | |
CN110546644A (zh) | 识别装置、识别方法以及识别程序 | |
US20150371381A1 (en) | Image region mapping device, 3d model generating apparatus, image region mapping method, and image region mapping program | |
CN110765903A (zh) | 行人重识别方法、装置及存储介质 | |
Führ et al. | Camera self-calibration based on nonlinear optimization and applications in surveillance systems | |
CN102013007B (zh) | 用于检测面部的装置和方法 | |
CN110786875B (zh) | 放射线图像处理装置及放射线图像处理方法 | |
CN105551042B (zh) | 一种扫描床标记点位置的确定方法和装置 | |
CN104966283A (zh) | 图像分层配准方法 | |
CN114266718A (zh) | 医用图像处理装置及方法、学习模型的学习方法 | |
KR101991452B1 (ko) | 유두위치 검출 방법, 유방 이미지 표시방법 및 유두위치 검출 장치 | |
CN105758337B (zh) | 一种获取透镜平面与图像传感器平面之间夹角的方法 | |
EP3933388A1 (en) | Feature point recognition system and recognition method | |
CN107784623B (zh) | X射线成像设备的图像处理方法及装置 | |
CN115423804A (zh) | 影像标定方法及装置、影像处理方法 | |
CN114240984A (zh) | 一种圆形标志点边缘提取方法及其应用 | |
JP5908816B2 (ja) | 骨塩量計測装置および方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |