CN108122203A - 一种几何参数的校正方法、装置、设备及*** - Google Patents
一种几何参数的校正方法、装置、设备及*** Download PDFInfo
- Publication number
- CN108122203A CN108122203A CN201611073139.4A CN201611073139A CN108122203A CN 108122203 A CN108122203 A CN 108122203A CN 201611073139 A CN201611073139 A CN 201611073139A CN 108122203 A CN108122203 A CN 108122203A
- Authority
- CN
- China
- Prior art keywords
- image
- position coordinates
- dimensional
- projection
- geometric parameter
- 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.)
- Granted
Links
- 239000011324 bead Substances 0.000 claims abstract description 122
- 238000000034 method Methods 0.000 claims abstract description 28
- 230000002285 radioactive effect Effects 0.000 claims description 30
- 238000001914 filtration Methods 0.000 claims description 17
- 238000003384 imaging method Methods 0.000 claims description 17
- 238000003708 edge detection Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000002184 metal Substances 0.000 claims description 5
- 229910052751 metal Inorganic materials 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 5
- 238000013527 convolutional neural network Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 2
- 239000010931 gold Substances 0.000 claims description 2
- 229910052737 gold Inorganic materials 0.000 claims description 2
- 238000003709 image segmentation Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 21
- 238000007408 cone-beam computed tomography Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 7
- 241000272470 Circus Species 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 230000005484 gravity Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000004308 accommodation Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本申请提供一种几何参数的校正方法、装置、设备及***，所述方法包括：获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；确定校正模体中的小球球心在二维投影图像中的物理投影位置坐标；利用物理投影位置坐标获得一组全局***几何参数；根据全局***几何参数对所有二维投影图像进行重建，得到校正模体的三维重建图像；确定三维重建图像中小球球心在获取的二维投影图像上的虚拟投影位置坐标；计算物理投影位置坐标与虚拟投影位置坐标的残差；利用残差对不同机架角度下的几何参数进行校正。本申请利用小球球心的物理投影坐标与小球球心的虚拟投影坐标之间的残差对不同机架角度的***几何参数进行补偿，修正了几何参数估计误差。
Description
技术领域
本申请涉及医疗成像技术，特别涉及一种几何参数的校正方法、装置、设备及***。
背景技术
锥形束投照计算机重组断层影像(Cone Beam Computed Tomography，CBCT)成像***一般包含一个X-ray放射源、一个探测板和一个旋转机架。其中，X-ray放射源和探测板附加到旋转机架上，随着旋转机架的旋转而旋转。CBCT成像原理是围绕三维物体旋转机架获取物体在多个角度下的投影图像，然后应用锥束投影图重建算法重建出物体三维图像。
CBCT成像***在重建图像时，需要获得成像***准确的几何***参数信息，CBCT***的几何参数是指描述***的X-ray光源、物体和探测器之间的角度和位置的几何参量，如X-ray放射源到探测板的距离，X-ray放射源到旋转轴的距离，探测板的倾斜和旋转角度等。由于X-ray放射源和探测板的重力，这些参数随着旋转机架的旋转而变化，也就是说在不同旋转机架角度下，成像***的几何***参数是不同的。不准确的CBCT几何***参数信息将导致重建图像出现几何伪影，降低重建图像质量。因此，CBCT***的几何参数校正具有提高重建图像质量，防止出现图像伪影的重要意义。
目前，通用的做法为：将一个带有若干标记点的校正模体放置在X-ray放射源和探测板之间，旋转机架得到该校正模体在不同角度的2D投影图，再根据2D投影图像中的标记点与成像***的几何参数预先设定的几何关系，计算得到成像***的几何参数。其中，具体的实现过程为：
第一种方案是：校正模体是包含24个相同大小金属小球的圆柱形模体，这些小球形成两个与圆柱模体同轴的圆环，每个圆环由均匀分布的12个小球组成，其中，要求小球的位置误差不超过25μm。之后，将校正模体放在X-ray放射源和探测板之间，使校正模体的中心与等中心对齐，校正模体的轴方向与旋转机架的旋转轴对齐。旋转机架以1度的间隔旋转，在每个机架角度分别采集一幅校正模体的2D投影图。对于每一幅投影图做如下处理：首先，确定小球投影在投影图中的位置坐标；然后，在定义的投影坐标系和探测器坐标系下，利用模体坐标系下小球位置坐标和投影图上小球位置坐标的关系，给出直接计算***几何参数的解析表达式，并计算出***几何参数；最后，通过将各个角度下的X-ray放射源的坐标拟合成圆形轨迹，重新建立一个投影***坐标系，并重新计算在此投影***坐标系下的所有角度下的几何参数。但是，这种方案中，对校正模体中小球的相对位置精度要求比较严格，小球的相对位置偏差影响最后的***几何参数的计算精度。
第二种方案：通常假设CBCT的***几何参数在不同机架角度下是一组固定的全局参数，然后利用校正模体中的标记点位置信息和其投影图像的位置信息关系，计算得到这一组全局参数。但是，这种方法的假设条件在实际中通常是不成立的，利用计算的全局参数重建三维图像会引入几何伪影，降低图像质量。
因此，如何降低校正模体精度要求并且能精确计算出成像***几何参数，是目前有待解决的技术问题。
发明内容
有鉴于此，本申请提供一种几何参数的校正方法、装置、设备及***，在降低校正模体精度要求的同时，修正了不同机架角度下的几何参数的估计误差。
具体地，本申请是通过如下技术方案实现的：
第一方面提供一种几何参数的校正方法，包括：
获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
利用所述物理投影位置坐标获得一组全局***几何参数；
根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；
利用所述残差对不同机架角度下的几何参数进行校正。
第二方面提供一种几何参数的校正装置，包括：
获取单元，用于获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
第一确定单元，用于确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
第一获得单元，用于利用所述物理投影位置坐标获得出一组全局***几何参数；
重建单元，用于根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
第二确定单元，用于确定所述三维重建图像中所述小球球心在所述获取的二维投影图像中的虚拟投影位置坐标；
第二获得单元，用于获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；
校正单元，用于利用所述残差对不同机架角度下的几何参数进行校正。
第三方面提供一种图像处理设备，包括：处理器；以及用于存储所述处理器可执行指令的存储器；其中，
所述处理器被配置为：
获取所述校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
利用所述物理投影位置坐标获得一组全局***几何参数；
根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
第四方面提供一种成像***，包括：旋转机架，附加在所述旋转机架上的放射源和探测器，其特征在于，还包括：
校正模体，用于放置在所述放射源和所述探测器之间，其中，所述校正模体包括多个金属小球，所有金属小球近似组成一条直线，相邻金属小球之间的间隔近似相等；
图像处理设备，用于获取所述校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；利用所述物理投影位置坐标获得一组全局***几何参数；根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
本申请实施例中，采集校正模体在不同机架角度下的二维投影图像，并确定在获取的二维投影图像中的小球球心物理投影的位置坐标，然后利用所述小球球心物理投影位置坐标与CBCT成像***几何参数的解析关系计算出一组全局***几何参数。然后利用这组全局几何参数对所有二维投影图像进行重建，得到校正模体的三维重建图像；然后，确定三维重建图像中小球球心的位置坐标，并计算三维重建图像中的小球球心在各个采集机架角度下投影到获取的二维投影图像的虚拟投影位置坐标。最后，计算采集的小球球心的物理投影坐标与计算的小球球心的虚拟投影坐标之间的残差，并根据该残差修正各个机架角度下的几何参数估计误差。由此可见，本申请实施例中，通过对校正模体在不同机架角度下的投影图像进行解析，并计算出一组固定不变的全局***几何参数，不需要标记校正模体中小球之间精确的相对位置，这极大的方便了算法的实现，降低了成本；同时，利用采集的小球的物理投影坐标与计算的小球虚拟投影坐标之间的残差对不同机架角度的***几何参数进行补偿，修正了***几何参数固定不变这一假设引起的几何参数估计误差。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。
附图说明
图1是本申请实施例提供的一种成像***的结构示意图；
图2是本申请实施例提供的一种校正模体的结构示意图；
图3是本申请实施例提供的一种几何参数的校正方法的流程图；
图4A是本申请实施例提供的CBCT***的表示几何***参数的示意图；
图4B是本申请实施例提供的CBCT***的表示几何***参数的另一示意图；
图5A是本申请本实施例提供的校正模体中每个小球运动的一个圆形运动轨迹的示意图；
图5B是本申请实施例提供的校正模体中每个小球的圆形运动轨迹在平板探测器上的投影为一个椭圆形轨迹的示意图；
图6是本申请实施例提供的相似三角关系的示意图；
图7是本申请实施例提供的射线源到旋转轴的距离的计算原理示意图；
图8是本申请提供的几何参数的校正装置所在设备的一种硬件结构图；
图9是本申请实施例提供的一种几何参数的校正装置的结构示意图；
图10是本申请实施例提供的一种几何参数的校正装置的另一结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
为了便于理解，在说明本申请实施例之前，先简单介绍下本申请涉及到的应用架构。
请参阅图1，是本申请实施例提供的一种成像***的结构示意图；所述成像***，比如，锥形束投照计算机重组断层影像(CBCT，Come beam computer tomography)成像***，该成像***一般包括：图像处理设备10，旋转机架11，附加在所述旋转机架11上的放射源(比如X-ray放射源)12，探测器13和校正模体14，旋转机架11可以围绕着旋转轴15旋转；其中，
校正模体14，用于放置在放射源12和探测器13之间，其中，校正模体14含有若干个金属小球，所有金属小球近似组成一条直线，相邻金属小球之间的间隔可以相等，也可以近似相等。
需要说明的是，在放置校正模体14时，应尽量使校正模体14与旋转机架11的旋转轴平行但不同轴，以保证校正模体14中的金属小球投影全部在探测器13上。
其中，本实施例提供的一种校正模体可以是一个长为约28cm直径为6mm的棍状物体，该校正模体中以嵌入21个金属小球为例，每个金属小球的直径为3mm金属小球、相邻金属小球之间的距可以为10mm，也可以近似等于10mm，当然，可以根据需要进行适应性调整。也就是说，相邻金属小球的间距要求不必特别精确。具体如图2所示，图2为本申请实施例提供的一种校正模体的结构示意图，图2中，相邻金属小球之间的间距以相等为例，具体应用中不限于此。
其中，图像处理设备10，分别与探测器和校正模体连接，可以应用到各种计算机上，用于获取所述校正模体14在不同机架角度下投影到探测器上的多幅二维投影图像；确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体14中；利用所述物理投影位置坐标获得一组全局***几何参数；根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
其中，图像处理设备10对图像处理的过程详见下述图3至图8所对应的实施例，在此不再赘述。
请参阅图3，是本申请实施例提供的一种几何参数的校正方法的流程图，所述方法包括：
步骤301：获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
其中，该步骤中的校正模体如图2所示，在此不再赘述。
该步骤中，需要先将校正模体放在放射源和探测器之间，以便于成像***能采集校正模体在不同角度下投影到探测器上的二维(即2D)投影图像。需要说明的是，在摆放校正模体时，应尽量使校正模体与旋转机架的旋转轴平行但不同轴，并且保证校正模体中的金属小球(下文统称为小球)投影全部在探测器上。
该步骤中，成像***旋转扫描一圈，就可以获得校正模体在不同机架角度下投影到探测器上的二维投影图像，其具体过程对于本领域技术人员来说，已是熟知技术，在此不再赘述。
步骤302：确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
其中，物理投影位置坐标用表示,m＝1,2,…M,n＝1,2,…,N，其中M是获取的二维投影图像的总张数，N是小球的个数。
该步骤中，一种确定方式为：
首先，对二维投影图像进行高斯滤波去除随机噪声得到滤波图像；其具体的实现过程对于本领域技术人员来说，已是熟知技术，在此不再赘述。
其次，根据所述滤波图像计算出小球投影的轮廓；
具体实现时，可以对滤波图像采用分水岭分割算法初步计算出小球投影的轮廓，并用小球周围背景的像素值填充其轮廓，得到背景图像；
再后，用滤波图像减去背景图像得到无背景图像；
最后，采用Canny边缘检测算法对无背景图像进行边缘检测，得到小球投影的边缘，并将边缘拟合成圆，圆心即为小球的物理投影位置坐标。
另一种确定方式为：先利用已训练好的卷积神经网络对所述二维投影图像进行分割，然后计算分割出图像中小球球心的物理投影位置坐标，其具体的计算过程对于本领域技术人员来说，已是熟知技术，在此不再赘述。
步骤303：利用所述物理投影位置坐标计算出一组全局***几何参数；
其中，一组全局***几何参数可以用R,D,(u0,v0),η表示，具体如图4A图4B所示，图4A和图4B分别为本实施例提供的CBCT***的表示几何***参数的示意图；
如图4A所示，R为放射源到旋转轴z轴的垂直距离；D为放射源垂直旋转轴z轴到探测器的距离；(u0,v0)为放射源垂直旋转轴z轴在探测器上的垂直投影点坐标；如图4B所示，η为探测器绕垂直投影点(u0,v0)的旋转角度。
由于运动是相互的，放射源和探测板围绕旋转轴旋转一周而校正模体静止不动，可以看作放射源和探测板静止不动而校正模体围绕旋转轴旋转一周，下文中的所有分析都是在此运动模型下进行的。此时，校正模体中的每个小球的运动为一个圆形运动轨迹，而每个小球的圆形运动轨迹在平板探测器上的投影为一个椭圆形轨迹，如图5A和5B所示。图5A为本实施例提供的校正模体中每个小球运动的一个圆形运动轨迹的示意图；图5B为本实施例提供的校正模体中每个小球的圆形运动轨迹在平板探测器上的投影为一个椭圆形轨迹的示意图。
下面分别提供计算出了一组全局***几何参数R,D,u0,v0,η的实现过程。
在图5A和图5B中，小球i的圆形轨迹与x轴和y轴的交点分别为θi1、θi2、θi3和θi4，它们在探测板投影的椭圆形轨迹分别为Ai1(ui1,vi1)、Ai2(ui2,vi2)、Ai3(ui3,vi3)和Ai4(ui4,vi4)，Ai1、Ai2、Ai3和Ai4分别为椭圆轨迹的短轴极值点和长轴极值点。
然后，再结合图6中所示的相似三角关系的示意图，可以得到如下关系式：
其中，A0为放射源垂直旋转轴z轴在探测器上的投影点。
将上述两式相加减可得：
然后有：
此外，由图6还可以很明显地得出下式：
于是可得：
由上式可知，v0与D呈线性关系，且此直线的参数可以通过计算N个小球的极值点得到。可以得到N个下面的数值对(Yi,Xi)，
之后可以用一条直线X＝v0+D·Y来拟合上述数值对(Yi,Xi)，从而可以求出几何参数v0和D。
下面来计算u0和η，首先计算投影图中每个小球形成的椭圆轨迹的中心坐标。椭圆轨迹的中心坐标可以通过上述的四个极值点的横坐标和纵坐标来得到，具体公式如下：
而后可以用u＝k+tanη·v来拟合上述所求椭圆轨迹中心点(ui0,vi0)。当求出k和tanη后，可以将前面得到的v0代入u＝k+tanη·v来得到u0，而η可由tan η的值得到。
射线源到旋转轴的距离R的计算原理示意图如图7所示。其中，下标1和2分别代表校正模体中的某两个小球，l为已知的两个小球间的距离，h为两个小球圆形轨迹的垂直距离，r1和r2分别为小球1和小球2的圆形轨迹的半径。于是有
其中，α1＝θ10-θ14，α2＝θ20-θ24，θ10和θ20为小球1和2的初始位置角度，θ14和θ24为小球1和2的一个极值点。
从前面的步骤中已知如下关系式
利用上述比例关系及已知的两个小球间距离l可以求得R，具体公式如下：
至此，计算出了一组全局***几何参数R,D,u0,v0,η。
步骤304：根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
该步骤中，根据计算出的一组全局***几何参数，利用FDK算法对所有二维投影图像进行重建，得到校正模体的3D重建图像。其中，利用FDK算法对所有二维投影图像进行重建的过程对于本领域技术人员来说已是熟知技术，在此不再赘述。
步骤305：确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
一种确定方式为：先确定三维重建图像中小球球心的三维位置坐标；具体过程为：
首先，对三维重建图像进行阈值滤波，将像素值低于预先给定阈值的三维像素点设置为0；
其次，对于三维重建图像中的每个像素(x0,y0,z0)，其中，x0中0表示迭代重建前的像素，若每个像素及其相邻的8个像素的值均不为0，则迭代计算每个像素(x0±k,y0±k,z0±k)邻域内的重心，其中，k是邻域范围的大小，且大于小球的半径，，如下式所示
i＝1,2,…,l
其中，xi，yi，zi中的i表示迭代第i次的像素。I(x,y,z)是阈值滤波后三维重建图像在(x,y,z)处的像素值，l为迭代次数。并将最终计算得到三维重建图像中的某个小球球心的三维位置坐标记为(xn,yn,zn),n＝1,2,…,N，其中N是小球的个数，本实例可以以21为例。
然后，根据全局***几何参数R,D,u0,v0,η，应用前向投影计算三维重建图像中若干小球球心在获取的二维投影图像中的虚拟投影位置坐标。
其中，若干小球球心用(xn,yn,zn)表示,n＝1,2,…,N；虚拟投影位置坐标用表示,m＝1,2,…M,n＝1,2,…,N，其中M是获取的二维投影图像的总张数，N是小球的个数。
步骤306：计算所述二维投影图像中小球球心的物理投影位置坐标与所述三维重建图像中小球球心在二维投影图像中的虚拟投影位置坐标的残差；
其中，物理投影位置坐标用表示，残差用 表示，如下式所示
残差(residuals)产生的主要原因是，在步骤303中假设***几何模型在机架的旋转过程中是静止不动的，而在实际中***几何模型随着机架旋转而变化，例如由于重力的原因探测器的位置会随机架旋转而发生变化。从另一个角度来看，步骤303中的CBCT成像***几何模型是实际成像***的简单近似模型，也就是，本实施例中的几何参数只是对***几何模型的一个近似描述。因此，残差是一个***性的误差，是可以通过建立模型来描述的。
步骤307：利用所述残差对不同机架角度下的几何参数进行校正。
该步骤中，先对所述残差进行曲线拟合，得到拟合的所述残差与机架角度的关系曲线；然后，利用拟合的所述关系曲线对不同机架角度下的几何参数进行补偿，修正不同机架角度下几何参数的误差。
下面以对放射源垂直旋转轴在探测器上的垂直投影点坐标(u0,v0)进行校正为例。
残差是一个***性的误差，是机架角度的正弦函数。残差满足如下公式
Δu(θ)＝au cos(buθ+cu)+du
Δv(θ)＝av cos(bvθ+cv)+dv
上述残差与机架角度的关系曲线可以通过对
n＝1,2,…,N进行曲线拟合得到。利用上述关系对垂直投影点坐标(u0,v0)进行校正，于是得到随机架角度变化的垂直投影点坐标：
u0(θ)＝u0+Δu(θ)
v0(θ)＝v0+Δv(θ)
所述随机架角度变化的垂直投影点坐标是对步骤303中静态***几何模型的一个修正，对机架旋转过程中***几何参数偏离全局几何参数的补偿。
本申请实施例中，采集校正模体在不同机架角度下的二维投影图像，并确定在采集的二维投影图像中小球球心的物理投影位置坐标，然后利用所述小球球心的物理投影位置坐标在一个静态的简单的几何模型下，直接获得一组全局***几何参数。然后利用这组全局几何参数对所有二维投影图像进行重建，得到校正模体的三维重建图像；然后，再确定三维重建图像中小球球心的位置坐标，并计算三维重建图像中的小球球心在各个机架角度下采集的二维投影图像上的虚拟投影位置坐标。最后，计算采集的二维图像中小球球心的物理投影位置坐标与三维重建图像中的小球球心在各个机架角度下采集的二维投影图像上的虚拟投影位置坐标的残差，并根据该残差修正各个机架角度下的几何参数。由此可见，本申请实施例中，利用校正模体在不同机架角度下的投影图像运用解析法，计算出一组固定不变的全局***几何参数，不需要标记校正模体中小球之间精确的相对位置，这极大的方便了算法的实现，降低了成本；同时，利用采集的小球的物理投影坐标与计算的小球虚拟投影坐标之间的残差对不同机架角度的***几何参数进行补偿，修正了***几何参数固定不变这一假设引起的几何参数估计误差。
与前述几何参数的校正方法的实施例相对应，本申请还提供了一种几何参数的校正装置的实施例。
本申请提供的几何参数的校正装置的实施例可以应用在设备上。装置实施例可以通过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑意义上的装置，是通过其所在设备的处理器将非易失性存储器中对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图8所示，为本申请提供的几何参数的校正装置所在设备的一种硬件结构图，除了图8所示的处理器、内存、网络接口、以及非易失性存储器之外，实施例中装置所在的设备通常根据该设备的实际功能，还可以包括其他硬件，对此不再赘述。
请参考图9，为本申请实施例提供的一种几何参数的校正装置的结构示意图，所述装置包括：获取单元91，第一确定单元92，第一获得单元93，重建单元94，第二确定单元95，第二获得单元96和校正单元97，其中，
获取单元91，用于获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
第一确定单元92，用于确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
第一获得单元93，用于利用所述物理投影位置坐标计算出一组全局***几何参数；
重建单元94，用于根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
第二确定单元95，用于确定所述三维重建图像中所述小球球心在各个机架角度下采集的二维投影图像上的虚拟投影位置坐标；
第二获得单元96，用于计算所述物理投影位置坐标与所述虚拟投影位置坐标的残差；
校正单元97，用于利用所述残差对不同机架角度下的几何参数进行校正补偿。
可选的，在另一实施例中，该实施例在上述实施例的基础上，所述第一确定单元91包括：滤波单元911，第一计算单元912，填充单元913，第二计算单元914，边缘检测单元915和第一拟合单元916，其结构示意图如图10。其中，
滤波单元911，用于对所述二维投影图像进行高斯滤波得到滤波图像；
第一计算单元912，用于根据所述滤波图像计算出小球投影的轮廓；
填充单元913，用于利用小球周围背景的像素值填充所述轮廓，得到背景图像；
第二计算单元914，用于将所述滤波图像减去背景图像得到无背景图像；
边缘检测915，用于对所述无背景图像进行边缘检测，得到小球投影的边缘；
第一拟合单元916，用于将所述边缘拟合成圆，圆心为小球的物理投影位置坐标。
可选的，在另一实施例中，该实施例在上述实施例的基础上，所述第一确定单元91包括：分割单元和位置坐标计算单元(图中未示)，其中，
分割单元，用于利用训练好的卷积神经网络对所述二维投影图像进行分割，得到所述二维投影图像的小球；
位置坐标计算单元，用于计算分割出图像中小球球心的物理投影位置坐标。
可选的，在另一实施例中，该实施例在上述实施例的基础上，所述第一获得单元93，具体用于在机架旋转过程中，利用所述物理投影位置坐标在一个静态的***几何模型下，解析计算出一组全局***几何参数，其中，所述全局***几何参数包括：放射源到旋转轴的垂直距离，放射源垂直旋转轴到探测器的距离，放射源垂直旋转轴在探测器上的垂直投影点坐标和探测器绕垂直投影轴的旋转角度。
可选的，在另一实施例中，该实施例在上述实施例的基础上，所述第二确定单元95包括：确定子单元和投影单元(图中未示)，其中，
确定子单元，用于确定三维重建图像中小球球心的三维位置坐标；
投影单元，用于根据所述三维位置坐标计算三维重建图像中小球球心在各个机架角度下采集的二维投影图像上的虚拟投影位置坐标。
可选的，在另一实施例中，该实施例在上述实施例的基础上，第二获得单元96包括：第二拟合单元和补偿单元(图中未示)，其中，
第二拟合单元，用于对所述残差和不同机架角度进行曲线拟合，得到拟合的所述残差与机架角度的关系曲线；
补偿单元，用于利用拟合的所述关系曲线对不同机架角度下的几何参数进行补偿，修正不同机架角度下几何参数的误差。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
本申请实施例还提供一种图像处理设备，包括：处理器；以及用于存储所述处理器可执行指令的存储器；其中，所述处理器被配置为：
获取所述校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
利用所述物理投影位置坐标获得一组全局***几何参数；
根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
该实施例中，处理器所执行的各个步骤的实现过程详见上述，在此不再赘述。
本申请实施例中，通过对校正模体在不同机架角度下的投影图像进行解析，并计算出一组固定不变的全局***几何参数，不需要标记校正模体中小球之间精确的相对位置，这极大的方便了算法的实现，降低了成本；同时，利用采集的小球的物理投影坐标与计算的小球虚拟投影坐标之间的残差对不同机架角度的***几何参数进行补偿，修正了***几何参数固定不变这一假设引起的几何参数估计误差。
本领域的技术人员可以清楚地了解到本发明实施例中的技术可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本发明实施例中的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
以上所述仅为本申请的较佳实施例而已，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。
Claims (15)
1.一种几何参数的校正方法，其特征在于，包括：
获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
利用所述物理投影位置坐标获得一组全局***几何参数；
根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；
利用所述残差对不同机架角度下的几何参数进行校正。
2.根据权利要求1所述的方法，其特征在于，所述确定小球球心在所述二维投影图像中的物理投影位置坐标包括：
对所述二维投影图像进行高斯滤波得到滤波图像；
根据所述滤波图像计算出小球投影的轮廓；
利用小球周围背景的像素值填充所述轮廓，得到背景图像；
将所述滤波图像减去背景图像得到无背景图像；
对所述无背景图像进行边缘检测，得到小球投影的边缘；
将所述边缘拟合成圆，圆心为小球球心的物理投影位置坐标。
3.根据权利要求1所述的方法，其特征在于，所述确定所述二维投影图像中小球球心的物理投影位置坐标包括：
利用训练好的卷积神经网络对所述二维投影图像进行图像分割；
计算分割出图像中小球球心的物理投影位置坐标。
4.根据权利要求1至3任一项所述的方法，其特征在于，利用所述物理投影位置坐标获得出一组全局***几何参数包括：
在机架旋转过程中，利用所述物理投影位置坐标在静态的***几何模型下，解析计算出一组全局***几何参数，所述全局***几何参数包括：放射源到旋转轴的垂直距离，放射源垂直旋转轴到探测器的距离，放射源垂直旋转轴在探测器上的垂直投影点坐标和探测器绕垂直投影轴的旋转角度。
5.根据权利要求4所述的方法，其特征在于，所述确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标包括：
确定三维重建图像中小球球心的三维位置坐标；
根据全局***几何参数计算所述三维重建图像中小球球心在获取的二维投影图像中的虚拟投影位置坐标。
6.根据权利要求5所述的方法，其特征在于，利用所述残差对不同机架角度下的几何参数进行校正包括：
对所述残差进行曲线拟合，得到拟合的所述残差与机架角度的关系曲线；
利用拟合的所述关系曲线对不同机架角度下的几何参数进行补偿，修正不同机架角度下几何参数的误差。
7.一种几何参数的校正装置，其特征在于，包括：
获取单元，用于获取校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
第一确定单元，用于确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
第一获得单元，用于利用所述物理投影位置坐标获得出一组全局***几何参数；
重建单元，用于根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
第二确定单元，用于确定所述三维重建图像中所述小球球心在所述获取的二维投影图像中的虚拟投影位置坐标；
第二获得单元，用于获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；
校正单元，用于利用所述残差对不同机架角度下的几何参数进行校正。
8.根据权利要求7所述的装置，其特征在于，所述第一确定单元包括：
滤波单元，用于对所述二维投影图像进行高斯滤波得到滤波图像；
第一计算单元，用于根据所述滤波图像计算出小球投影的轮廓；
填充单元，用于利用小球周围背景的像素值填充所述轮廓，得到背景图像；
第二计算单元，用于将所述滤波图像减去背景图像得到无背景图像；
边缘检测单元，用于对所述无背景图像进行边缘检测，得到小球投影的边缘；
第一拟合单元，用于将所述边缘拟合成圆，圆心为小球球心的物理投影位置坐标。
9.根据权利要求7所述的装置，其特征在于，所述第一确定单元包括：
分割单元，用于利用训练好的卷积神经网络对所述二维投影图像进行分割；
位置坐标计算单元，用于计算分割出图像中小球球心的物理投影位置坐标。
10.根据权利要求7至9任一项所述的装置，其特征在于
所述第一获得单元，具体用于在机架旋转过程中，利用所述物理投影位置坐标在静态的***几何模型下，解析计算出一组全局***几何参数，其中，所述全局***几何参数包括：放射源到旋转轴的垂直距离，放射源垂直旋转轴到探测器的距离，放射源垂直旋转轴在探测器上的垂直投影点坐标和探测器绕垂直投影轴的旋转角度。
11.根据权利要求10所述的装置，其特征在于，所述第二确定单元包括：
确定子单元，用于确定三维重建图像中小球球心的三维位置坐标；
投影单元，用于根据所述三维位置坐标计算三维重建图像中小球球心在所述获取的二维投影图像中的虚拟投影位置坐标。
12.根据权利要求11所述的装置，其特征在于，第二获得单元包括：
第二拟合单元，用于对所述残差进行曲线拟合，得到拟合的所述残差与机架角度的关系曲线；
补偿单元，用于利用拟合的所述关系曲线对不同机架角度下的几何参数进行补偿，修正不同机架角度下几何参数的误差。
13.一种图像处理设备，其特征在于，包括：处理器；以及用于存储所述处理器可执行指令的存储器；其中，
所述处理器被配置为：
获取所述校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；
确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；
利用所述物理投影位置坐标获得一组全局***几何参数；
根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；
确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；
获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
14.一种成像***，包括：旋转机架，附加在所述旋转机架上的放射源和探测器，其特征在于，还包括：
校正模体，用于放置在所述放射源和所述探测器之间，其中，所述校正模体包括多个金属小球，所有金属小球近似组成一条直线，相邻金属小球之间的间隔近似相等；
图像处理设备，用于获取所述校正模体在不同机架角度下投影到探测器上的多幅二维投影图像；确定小球球心在所述二维投影图像中的物理投影位置坐标，所述小球嵌在所述校正模体中；利用所述物理投影位置坐标获得一组全局***几何参数；根据所述全局***几何参数对所有二维投影图像进行重建，得到所述校正模体的三维重建图像；确定所述三维重建图像中所述小球球心在获取的二维投影图像上的虚拟投影位置坐标；获得所述物理投影位置坐标与所述虚拟投影位置坐标的残差；利用所述残差对不同机架角度下的几何参数进行校正。
15.根据权利要求14所述的***，其特征在于，所述校正模体与所述旋转机架的旋转轴平行且不同轴。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611073139.4A CN108122203B (zh) | 2016-11-29 | 2016-11-29 | 一种几何参数的校正方法、装置、设备及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611073139.4A CN108122203B (zh) | 2016-11-29 | 2016-11-29 | 一种几何参数的校正方法、装置、设备及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108122203A true CN108122203A (zh) | 2018-06-05 |
| CN108122203B CN108122203B (zh) | 2020-04-07 |
Family
ID=62226571
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611073139.4A Active CN108122203B (zh) | 2016-11-29 | 2016-11-29 | 一种几何参数的校正方法、装置、设备及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108122203B (zh) |
Cited By (29)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109146800A (zh) * | 2018-07-23 | 2019-01-04 | 广州华端科技有限公司 | 锥束计算机断层成像图像校正方法和*** |
| CN109549661A (zh) * | 2018-12-29 | 2019-04-02 | 北京纳米维景科技有限公司 | 一种探测器几何校正体模及校正方法 |
| CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
| CN109717889A (zh) * | 2018-12-14 | 2019-05-07 | 深圳市深图医学影像设备有限公司 | 口腔锥束ct***几何参数校正模型、方法及*** |
| CN109820531A (zh) * | 2019-01-30 | 2019-05-31 | 东软医疗***股份有限公司 | 一种锥束型ct几何参数的校正方法和装置 |
| CN109875591A (zh) * | 2019-03-14 | 2019-06-14 | 浙江大学 | 一种锥束ct***的几何标定方法 |
| CN110025327A (zh) * | 2019-05-29 | 2019-07-19 | 明峰医疗***股份有限公司 | 一种pet/ct的配准装置和方法 |
| CN110084855A (zh) * | 2019-04-19 | 2019-08-02 | 合肥中科离子医学技术装备有限公司 | 一种改进cbct几何参数标定算法 |
| CN110327064A (zh) * | 2019-06-19 | 2019-10-15 | 东软医疗***股份有限公司 | 源像距检测方法、装置及设备 |
| CN110353713A (zh) * | 2019-07-16 | 2019-10-22 | 上海联影医疗科技有限公司 | 几何校正模体、重建图像几何参数的确定方法和装置 |
| CN110766629A (zh) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | Cbct***几何校正方法、装置、计算机设备和存储介质 |
| CN110974280A (zh) * | 2019-12-23 | 2020-04-10 | 常州博恩中鼎医疗科技有限公司 | 一种适用于锥束ct的几何参数精确自动化校正方法 |
| CN111161173A (zh) * | 2019-12-18 | 2020-05-15 | 北京城市网邻信息技术有限公司 | 图像的校正信息获取、模型构建方法和装置、介质 |
| CN111248934A (zh) * | 2018-11-30 | 2020-06-09 | 辽宁开普医疗***有限公司 | 一种用于cbct***机械校正的方法及*** |
| CN111419257A (zh) * | 2020-04-07 | 2020-07-17 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描***和可读存储介质 |
| CN112233157A (zh) * | 2020-10-14 | 2021-01-15 | 首都师范大学 | 微纳米ct***转台误差的频域校正方法 |
| CN112991464A (zh) * | 2021-03-19 | 2021-06-18 | 山东大学 | 基于立体视觉的三维重建的点云误差补偿方法及*** |
| CN113286053A (zh) * | 2021-05-14 | 2021-08-20 | 广州广电运通金融电子股份有限公司 | 一种文件扫描方法、设备、介质及产品 |
| CN113643428A (zh) * | 2021-08-17 | 2021-11-12 | 北京唯迈医疗设备有限公司 | 一种适用于多自由度锥形束ct的全参数几何校准方法 |
| CN113782174A (zh) * | 2021-11-12 | 2021-12-10 | 极限人工智能有限公司 | Ct机的校验方法、装置、电子设备及可读存储介质 |
| CN113876346A (zh) * | 2021-11-16 | 2022-01-04 | 南京佗道医疗科技有限公司 | 一种倾斜图像的迭代校正方法 |
| CN113963057A (zh) * | 2021-09-07 | 2022-01-21 | 于留青 | 成像几何关系标定方法、装置、电子设备以及存储介质 |
| WO2022027545A1 (zh) * | 2020-08-07 | 2022-02-10 | 西安大医集团股份有限公司 | 图像数据的处理方法、放疗设备的等中心验证方法及*** |
| CN114460110A (zh) * | 2022-03-08 | 2022-05-10 | 中国电子科技集团公司第三十八研究所 | 一种伺服***误差补偿方法 |
| CN114712728A (zh) * | 2022-03-31 | 2022-07-08 | 西安大医集团股份有限公司 | 投影参数修正方法及计算机设备 |
| CN115409736A (zh) * | 2022-09-16 | 2022-11-29 | 深圳市宝润科技有限公司 | Dr的几何校正方法、装置、计算机设备及存储介质 |
| CN115501500A (zh) * | 2021-06-23 | 2022-12-23 | 西安大医集团股份有限公司 | 抖动偏移量确定方法、图像校准方法、设备、***及介质 |
| CN116831607A (zh) * | 2023-09-01 | 2023-10-03 | 赛诺威盛科技(北京)股份有限公司 | 用于医疗影像设备的校正装置、校正方法和医疗影像设备 |
| CN117392025A (zh) * | 2023-12-11 | 2024-01-12 | 浙江大学杭州国际科创中心 | X射线锥束ct位姿矫正方法、***及装置 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7950849B2 (en) * | 2005-11-29 | 2011-05-31 | General Electric Company | Method and device for geometry analysis and calibration of volumetric imaging systems |
| CN103729827A (zh) * | 2012-10-11 | 2014-04-16 | 南京普爱射线影像设备有限公司 | 影像增强器c形臂x射线机三维重建重叠伪影校正方法 |
| CN103969269A (zh) * | 2013-01-31 | 2014-08-06 | Ge医疗***环球技术有限公司 | 用于几何校准ct扫描仪的方法和装置 |
| CN104665862A (zh) * | 2015-02-16 | 2015-06-03 | 清华大学 | 在cbct中消除几何伪影的方法以及使用该方法的cbct*** |
| CN105769233A (zh) * | 2016-02-29 | 2016-07-20 | 江苏美伦影像***有限公司 | 一种几何校正方法 |
| CN105869130A (zh) * | 2016-04-20 | 2016-08-17 | 广州华端科技有限公司 | Ct***几何校正效果的检验方法及装置 |
| CN105931202A (zh) * | 2016-04-20 | 2016-09-07 | 广州华端科技有限公司 | 几何校正模体的校正方法和*** |
| CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct***几何参数校正模型及方法 |
| CN106127697A (zh) * | 2016-06-07 | 2016-11-16 | 哈尔滨工业大学 | 无人机机载成像高光谱几何校正方法 |
-
2016
- 2016-11-29 CN CN201611073139.4A patent/CN108122203B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7950849B2 (en) * | 2005-11-29 | 2011-05-31 | General Electric Company | Method and device for geometry analysis and calibration of volumetric imaging systems |
| CN103729827A (zh) * | 2012-10-11 | 2014-04-16 | 南京普爱射线影像设备有限公司 | 影像增强器c形臂x射线机三维重建重叠伪影校正方法 |
| CN103969269A (zh) * | 2013-01-31 | 2014-08-06 | Ge医疗***环球技术有限公司 | 用于几何校准ct扫描仪的方法和装置 |
| CN104665862A (zh) * | 2015-02-16 | 2015-06-03 | 清华大学 | 在cbct中消除几何伪影的方法以及使用该方法的cbct*** |
| CN105769233A (zh) * | 2016-02-29 | 2016-07-20 | 江苏美伦影像***有限公司 | 一种几何校正方法 |
| CN105869130A (zh) * | 2016-04-20 | 2016-08-17 | 广州华端科技有限公司 | Ct***几何校正效果的检验方法及装置 |
| CN105931202A (zh) * | 2016-04-20 | 2016-09-07 | 广州华端科技有限公司 | 几何校正模体的校正方法和*** |
| CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct***几何参数校正模型及方法 |
| CN106127697A (zh) * | 2016-06-07 | 2016-11-16 | 哈尔滨工业大学 | 无人机机载成像高光谱几何校正方法 |
Non-Patent Citations (2)
| Title |
|---|
| 池明辉: "锥束工业CT***几何参数误差校正算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 钱贝贝: "医用锥束CT射束硬化校正方法研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (40)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109146800B (zh) * | 2018-07-23 | 2019-08-13 | 广州华端科技有限公司 | 锥束计算机断层成像图像校正方法和*** |
| CN109146800A (zh) * | 2018-07-23 | 2019-01-04 | 广州华端科技有限公司 | 锥束计算机断层成像图像校正方法和*** |
| CN111248934B (zh) * | 2018-11-30 | 2023-09-29 | 辽宁开普医疗***有限公司 | 一种用于cbct***机械校正的方法及*** |
| CN111248934A (zh) * | 2018-11-30 | 2020-06-09 | 辽宁开普医疗***有限公司 | 一种用于cbct***机械校正的方法及*** |
| CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
| CN109717889A (zh) * | 2018-12-14 | 2019-05-07 | 深圳市深图医学影像设备有限公司 | 口腔锥束ct***几何参数校正模型、方法及*** |
| CN109549661A (zh) * | 2018-12-29 | 2019-04-02 | 北京纳米维景科技有限公司 | 一种探测器几何校正体模及校正方法 |
| CN109820531A (zh) * | 2019-01-30 | 2019-05-31 | 东软医疗***股份有限公司 | 一种锥束型ct几何参数的校正方法和装置 |
| CN109875591B (zh) * | 2019-03-14 | 2020-05-29 | 浙江大学 | 一种锥束ct***的几何标定方法 |
| CN109875591A (zh) * | 2019-03-14 | 2019-06-14 | 浙江大学 | 一种锥束ct***的几何标定方法 |
| CN110084855A (zh) * | 2019-04-19 | 2019-08-02 | 合肥中科离子医学技术装备有限公司 | 一种改进cbct几何参数标定算法 |
| CN110025327A (zh) * | 2019-05-29 | 2019-07-19 | 明峰医疗***股份有限公司 | 一种pet/ct的配准装置和方法 |
| CN110025327B (zh) * | 2019-05-29 | 2024-05-07 | 明峰医疗***股份有限公司 | 一种pet/ct的配准装置和方法 |
| CN110327064A (zh) * | 2019-06-19 | 2019-10-15 | 东软医疗***股份有限公司 | 源像距检测方法、装置及设备 |
| CN110353713A (zh) * | 2019-07-16 | 2019-10-22 | 上海联影医疗科技有限公司 | 几何校正模体、重建图像几何参数的确定方法和装置 |
| CN110766629A (zh) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | Cbct***几何校正方法、装置、计算机设备和存储介质 |
| CN110766629B (zh) * | 2019-10-17 | 2022-03-01 | 广州华端科技有限公司 | Cbct***几何校正方法、装置、计算机设备和存储介质 |
| CN111161173A (zh) * | 2019-12-18 | 2020-05-15 | 北京城市网邻信息技术有限公司 | 图像的校正信息获取、模型构建方法和装置、介质 |
| CN110974280A (zh) * | 2019-12-23 | 2020-04-10 | 常州博恩中鼎医疗科技有限公司 | 一种适用于锥束ct的几何参数精确自动化校正方法 |
| CN111419257A (zh) * | 2020-04-07 | 2020-07-17 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描***和可读存储介质 |
| CN111419257B (zh) * | 2020-04-07 | 2022-11-22 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描***和可读存储介质 |
| WO2022027545A1 (zh) * | 2020-08-07 | 2022-02-10 | 西安大医集团股份有限公司 | 图像数据的处理方法、放疗设备的等中心验证方法及*** |
| CN112233157A (zh) * | 2020-10-14 | 2021-01-15 | 首都师范大学 | 微纳米ct***转台误差的频域校正方法 |
| CN112233157B (zh) * | 2020-10-14 | 2022-03-08 | 俐玛精密测量技术(苏州)有限公司 | 微纳米ct***转台误差的频域校正方法 |
| CN112991464A (zh) * | 2021-03-19 | 2021-06-18 | 山东大学 | 基于立体视觉的三维重建的点云误差补偿方法及*** |
| CN113286053A (zh) * | 2021-05-14 | 2021-08-20 | 广州广电运通金融电子股份有限公司 | 一种文件扫描方法、设备、介质及产品 |
| CN115501500A (zh) * | 2021-06-23 | 2022-12-23 | 西安大医集团股份有限公司 | 抖动偏移量确定方法、图像校准方法、设备、***及介质 |
| CN113643428A (zh) * | 2021-08-17 | 2021-11-12 | 北京唯迈医疗设备有限公司 | 一种适用于多自由度锥形束ct的全参数几何校准方法 |
| CN113963057A (zh) * | 2021-09-07 | 2022-01-21 | 于留青 | 成像几何关系标定方法、装置、电子设备以及存储介质 |
| CN113782174A (zh) * | 2021-11-12 | 2021-12-10 | 极限人工智能有限公司 | Ct机的校验方法、装置、电子设备及可读存储介质 |
| CN113876346A (zh) * | 2021-11-16 | 2022-01-04 | 南京佗道医疗科技有限公司 | 一种倾斜图像的迭代校正方法 |
| CN113876346B (zh) * | 2021-11-16 | 2024-06-07 | 佗道医疗科技有限公司 | 一种倾斜图像的迭代校正方法 |
| CN114460110A (zh) * | 2022-03-08 | 2022-05-10 | 中国电子科技集团公司第三十八研究所 | 一种伺服***误差补偿方法 |
| CN114460110B (zh) * | 2022-03-08 | 2023-06-06 | 中国电子科技集团公司第三十八研究所 | 一种伺服***误差补偿方法 |
| CN114712728A (zh) * | 2022-03-31 | 2022-07-08 | 西安大医集团股份有限公司 | 投影参数修正方法及计算机设备 |
| CN115409736A (zh) * | 2022-09-16 | 2022-11-29 | 深圳市宝润科技有限公司 | Dr的几何校正方法、装置、计算机设备及存储介质 |
| CN116831607A (zh) * | 2023-09-01 | 2023-10-03 | 赛诺威盛科技(北京)股份有限公司 | 用于医疗影像设备的校正装置、校正方法和医疗影像设备 |
| CN116831607B (zh) * | 2023-09-01 | 2023-12-01 | 赛诺威盛科技(北京)股份有限公司 | 用于医疗影像设备的校正装置、校正方法和医疗影像设备 |
| CN117392025A (zh) * | 2023-12-11 | 2024-01-12 | 浙江大学杭州国际科创中心 | X射线锥束ct位姿矫正方法、***及装置 |
| CN117392025B (zh) * | 2023-12-11 | 2024-02-23 | 浙江大学杭州国际科创中心 | X射线锥束ct位姿矫正方法、***及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108122203B (zh) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108122203A (zh) | 一种几何参数的校正方法、装置、设备及*** | |
| KR100667014B1 (ko) | 단층상 및 3차원 화상작성방법 | |
| Cho et al. | Accurate technique for complete geometric calibration of cone‐beam computed tomography systems | |
| US9435901B2 (en) | Systems and methods for improving the quality of images in a PET scan | |
| US7307252B2 (en) | Detector head position correction for hybrid SPECT/CT imaging apparatus | |
| CN103384498B (zh) | 探测装置 | |
| CN104107065B (zh) | 3d图像集在不同空间之间的最佳变换 | |
| JP4415762B2 (ja) | 断層撮影装置 | |
| CN110381839A (zh) | 处置计划图像、分次内3D图像和分次内2D x射线图像的图像配准 | |
| US20050117708A1 (en) | Method and system for calibrating a source and detector instrument | |
| JP2008520257A5 (zh) | ||
| JP7154611B2 (ja) | インテリアｃｔ画像生成方法 | |
| CN104665862B (zh) | 在cbct中消除几何伪影的方法以及使用该方法的cbct*** | |
| Ferrucci et al. | Measurement of the X-ray computed tomography instrument geometry by minimization of reprojection errors—Implementation on experimental data | |
| US7459689B2 (en) | Detector head position calibration and correction for SPECT imaging apparatus using virtual CT | |
| US8755586B2 (en) | Image reconstruction including shift-variant blur compensation | |
| CN108333197B (zh) | 偏置扫描模式下工业ct***转台旋转中心标定方法 | |
| CN107192726A (zh) | 板壳物体快速高分辨三维锥束计算机层析成像方法及装置 | |
| Niebler et al. | Projection‐based improvement of 3D reconstructions from motion‐impaired dental cone beam CT data | |
| CN107233105A (zh) | 一种用于ct图像重建的修正方法及修正*** | |
| CN104507391A (zh) | 摄影控制装置、摄影***、摄影控制方法以及程序 | |
| ES2717001T3 (es) | Procedimiento de segmentación de la imagen de un objeto reconstruido mediante reconstrucción tridimensional | |
| EP3910589A1 (en) | Image registration method and apparatus, radiotherapy device and computer readable storage medium | |
| Aoi et al. | System design and development of a pinhole SPECT system for quantitative functional imaging of small animals | |
| CN113963056B (zh) | Ct图像重建方法、装置、电子设备以及存储介质 |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220325 Address after: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: Shenyang Neusoft Medical Systems Co.,Ltd. Address before: Room 117 and 119, 1000 Ziyue Road, Minhang District, Shanghai, 200241 Patentee before: Shanghai Neusoft Medical Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right |