CN108492341B - 一种基于像素顶点的平行束投影方法 - Google Patents
一种基于像素顶点的平行束投影方法 Download PDFInfo
- Publication number
- CN108492341B CN108492341B CN201810109560.9A CN201810109560A CN108492341B CN 108492341 B CN108492341 B CN 108492341B CN 201810109560 A CN201810109560 A CN 201810109560A CN 108492341 B CN108492341 B CN 108492341B
- Authority
- CN
- China
- Prior art keywords
- detector
- projection
- pixel
- calculating
- center
- 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
Images
Classifications
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
Abstract
本发明属于图像投影与图像重建领域,公开了一种基于像素顶点的平行束投影方法,以图像中心为原点建立二维坐标系,计算每个像素中心的坐标,以探测器阵列中心为原点建立一维坐标系,计算每个探测器的坐标;在确定的投影角度下,计算每个像素中心投影到探测器上的坐标;计算投影的上底和下底的长度;计算像素四个顶点投影到探测器上的坐标;根据坐标判断该像素点投影到哪些探测器上,计算对应的面积即为权值。本发明具有投影精度高的优点,可用于断层成像迭代重建。本发明将像素的四个顶点投影到探测器上,在一般角度下得到形状为梯形的投影,计算对应探测器的权值。和基于距离驱动投影模型相比,本发明的投影结果误差更小。
Description
技术领域
本发明属于图像投影与图像重建技术领域,尤其涉及一种基于像素顶点的平行束投影方法。
背景技术
目前,业内常用的现有技术是这样的:光学投影断层成像(Optical TomographyProjection,OPT)技术主要用于1-10mm生物样本的成像,具有分辨率高、设备成本低、操作方便,既能进行结构成像也能进行分子成像等优点,很好的适用于生命科学领域小尺度生物样本的成像需求。OPT利用光在弱散射介质中沿直线传播的特点,在光束穿透生物样本得到各角度投影图后,利用投影图进行三维重建得到样本的三维图像。重建算法分为解析法和迭代法两大类,迭代法重建中需要反复使用投影和反投影操作。常见的投影模型有以下三种:1.基于像素驱动模型:沿投影角度将像素中心投影到探测器上,执行线性外插将像素点对投影数据的贡献按比例分配到相应的探测器上;2.基于射线驱动模型:射线沿投影角度经过探测器中心,图像的每行或每列与射线相交,通过插值运算可得到这些数值,然后将所有的行或列的值累加起来,构成该射线对应的探测器像素接收到的投影值;3.基于距离驱动模型:距离驱动中,不再使用插值系数作为权值,而是将图像像素和探测器像素一同投影在横轴上,然后计算它们的重合区域,以此作为该图像像素和探测器像素之间的加权权值。在上述几种投影模型中,基于像素驱动的模型在反投影操作中广泛使用,如果应用于正向投影则需要外插运算,从而会引入高频伪影,导致重建结果质量很差,因此很少在正向投影中使用。基于射线的模型在正向投影中广泛使用,但该模型如果应用于反向投影,则也会因为外插运算会引入高频伪影,同样地,会导致重建结果质量很差。基于距离投影模型将像素的两条边界中点沿射线方向投影到探测器,得到形状为矩形的投影,而两条射线两侧的小部分面积会被忽略掉,与真实投影结果相比仍有误差,因此投影的精度仍然不够高。
综上所述,现有技术存在的问题是:
(1)基于像素驱动的模型在反投影操作中广泛使用,如果应用于正向投影则会引入高频伪影,很少在正向投影中使用。
(2)基于射线的模型在正向投影中广泛使用,但模型如果应用于反向投影,引入高频伪影。
(3)基于距离投影模型投影的精度不够高。
解决上述技术问题的难度和意义:
建立和实际情况相同的投影模型是解决上述技术问题的关键,如果能够解决上述问题,对迭代重建结果的质量会有较大提升。
发明内容
针对现有技术存在的问题,本发明提供了一种基于像素顶点平行束投影方法。
本发明是这样实现的,一种基于像素顶点平行束投影方法,所述基于像素顶点平行束投影方法包括:以图像中心为原点建立二维坐标系,计算每个像素中心的坐标,以探测器阵列中心为原点建立一维坐标系,计算每个探测器的坐标;在确定的投影角度下,计算每个像素中心投影到探测器上的坐标;计算投影的上底和下底的长度;计算像素四个顶点投影到探测器上的坐标;根据坐标判断该像素点投影到哪些探测器上,计算对应的面积即为权值。
进一步,所述基于像素顶点平行束投影方法包括以下步骤:
(1)以图像中心为原点建立二维坐标系I,计算每个像素中心的坐标(xi,yi),i=1,2,…N,N为图像中像素的个数,以探测器阵列中心为原点建立一维坐标系D,计算每个探测器边界的坐标ti,i=1,2,…M,M为探测器的个数;
(2)在确定的投影角度θ下,计算每个像素中心投影到探测器上的坐标Pi:
Pi=xi·cosθ-yi·sinθ+r;
其中r是坐标系I和D的原点在横向上的距离:r=OI-OD;
(3)计算投影的上底和下底的长度w1和w2:
其中X是像素的边长,θ是投影角度;
(4)计算像素的四个顶点投影到探测器上的坐标:
(5)根据像素四个顶点投影之后的坐标判断该像素投影到哪些探测器上,将投影按面积分配到探测器上,每个探测器分配到的面积为像素对该探测器的权值贡献,投影和探测器的相对位置有十种情况,因此权值的计算也对应如下情况:w=dh;
(6)重复(2)至(5),完成所有投影角度的投影。
进一步,所述(3)中计算投影上底和下底的长度时,将像素的两条对角线投影到探测器上,长度较长的一条为下底,长度较短的一条为上底;投影角度为0°,90°,180°或270°,两条对角线投影长度相等,此时像素的投影为矩形;投影角度为45°,135°,225°或315°时,其中一条对角线垂直于探测器,投影长度为0,像素的投影为等腰三角形。
进一步,所述(5)中判断像素投影到哪些探测器上及对应的权值计算方法具体如下:某一探测器的左边界坐标为Dbe,右边界坐标为Dend;
1)如果Dend≤Pube满足,则像素没有投影到该探测器上;
6)如果Ptbe≤Dbe,Dend≤Ptend两式同时满足,则分配到该探测器的投影形状为矩形,使用w=dh计算权值;
12)如果Dbe≥Puend满足,则像素没有投影到该探测器上。
本发明的另一目的在于提供一种应用所述基于像素顶点平行束投影方法的光学投影断层成像系统。
综上所述,本发明的优点及积极效果为:本发明将像素的四个顶点投影到探测器上,在一般角度下得到形状为梯形的投影,计算对应探测器的权值。而基于距离驱动的投影模型是将像素的两条边界中点投影到探测器,得到形状为矩形的投影。和基于距离驱动投影模型相比,本发明的投影结果误差更小。
附图说明
图1是本发明实施例提供的基于像素顶点平行束投影方法流程图。
图2是本发明实施例提供的投影模型、像素点投影的示意图;
图中:(a)为本发明投影模型的示意图;(b)为一个像素点投影的详细示意图。
图3是本发明实施例提供的探测器和投影相对位置示意图;
图中:(a)至(b)为探测器和投影相对位置的十种情况示意图;
图4是本发明实施例提供的待投影的原始图像示意图。
图5是本发明实施例提供的原始图像的前向投影结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
光学投影断层成像(Optical Tomography Projection,OPT)技术主要用于1-10mm生物样本的成像,具有分辨率高、设备成本低、操作方便,既能进行结构成像也能进行分子成像等优点,很好的适用于生命科学领域小尺度生物样本的成像需求。将图像像素的四个顶点投影到探测器上,得到形状为梯形的投影(投影角度为0°,90°,180°,270°时为矩形,投影角度为45°,135°,225°,315°时为三角形),将投影按面积分配到探测器上,每个探测器分配到的面积即为像素对该探测器的权值贡献。
如图1所示,本发明实施例提供的基于像素顶点平行束投影方法包括以下步骤:
S101:以图像中心为原点建立二维坐标系,计算每个像素中心的坐标,以探测器阵列中心为原点建立一维坐标系,计算每个探测器的坐标;
S102:在确定的投影角度下,计算每个像素中心投影到探测器上的坐标;
S103:计算投影的上底和下底的长度;
S104:计算像素四个顶点投影到探测器上的坐标;
S105:根据坐标判断该像素点投影到哪些探测器上,计算对应的面积即为权值。
下面结合附图对本发明的应用原理作进一步的描述。
在此实施例中,以Shepp–Logan模型图像作为原始图像,图像大小为256×256,每隔1°进行一次投影,共进行360次,得到的投影图像大小为256×360,图像边长X设置为1,坐标系I和D的原点在横向上的距离r设置为0,探测器尺寸设置为1,探测器个数设置为256。
本发明基于像素顶点的平行束投影模型,具体步骤包括如下:
步骤1:以图像中心为原点建立二维坐标系I,计算每个像素中心的坐标(xi,yi)(i=1,2,…,65536,65536为图像中像素的个数),以探测器阵列中心为原点建立一维坐标系D,计算每个探测器边界的坐标ti(i=1,2,…,256,256为探测器的个数)。
步骤2:在确定的投影角度θ下,计算每个像素中心投影到探测器上的坐标Pi:
Pi=xi·cosθ-yi·sinθ+r;
其中r是坐标系I和D的原点在横向上的距离r=OI-OD。
步骤3:计算投影的上底和下底的长度w1和w2:
其中X是像素的边长,θ是投影角度。
步骤(3)中计算投影上底和下底的长度时,将像素的两条对角线投影到探测器上,长度较长的一条为下底,长度较短的一条为上底。如果投影角度为0°,90°,180°或270°,两条对角线投影长度相等,此时像素的投影为矩形;如果投影角度为45°,135°,225°或315°时,其中一条对角线垂直于探测器,投影长度为0,此时像素的投影为等腰三角形。
步骤4:计算像素四个顶点投影到探测器上的坐标:
步骤5:根据像素四个顶点投影之后的坐标判断该像素投影到哪些探测器上,将投影按面积分配到探测器上,每个探测器分配到的面积即为像素对该探测器的权值贡献。投影和探测器的相对位置有十种情况,因此权值的计算也对应如下十种情况:
(5e)w=dh;
假设某一探测器的左边界坐标为Dbe,右边界坐标为Dend,
1)如果Dend≤Pube满足,则像素没有投影到该探测器上;
2)如果Dbe≤Pube,Pube<Dend≤Ptbe两式同时满足,则分配到该探测器的投影形状为三角形,使用步骤(5a)中的公式计算权值;
3)如果Dbe≤Pube,Ptbe<Dend≤Ptend两式同时满足,则分配到该探测器的投影形状为直角梯形,使用步骤(5b)中的公式计算权值;
4)如果Dbe≤Pube,Ptend<Dend<Puend两式同时满足,则分配到该探测器的投影形状为四边形,使用步骤(5c)中的公式计算权值;
5)如果Pube<Dbe<Ptbe,Ptend<Dend<Puend两式同时满足,则分配到该探测器的投影形状为六边形,使用步骤(5d)中的公式计算权值;
6)如果Ptbe≤Dbe,Dend≤Ptend两式同时满足,则分配到该探测器的投影形状为矩形,使用步骤(5e)中的公式计算权值;
7)如果Ptbe≤Dbe<Ptend,Ptend<Dend<Puend两式同时满足,则分配到该探测器的投影形状为五边形,使用步骤(5f)中的公式计算权值;
8)如果Pube<Dbe<Ptbe,Ptbe<Dend≤Puend两式同时满足,则分配到该探测器的投影形状为五边形,使用步骤(5g)中的公式计算权值;
9)如果Pube<Dbe<Ptbe,Dend≥Puend两式同时满足,则分配到该探测器的投影形状为四边形,使用步骤(5h)中的公式计算权值;
10)如果Ptbe≤Dbe<Ptend,Dend≥Puend两式同时满足,则分配到该探测器的投影形状为直角梯形,使用步骤(5i)中的公式计算权值;
11)如果Ptend≤Dbe<Puend,Dend≥Puend两式同时满足,则分配到该探测器的投影形状为三角形,使用步骤(5j)中的公式计算权值;
12)如果Dbe≥Puend满足,则像素没有投影到该探测器上。
步骤6:重复步骤(2)至(5),完成360个投影角度的投影。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于像素顶点的平行束投影方法,其特征在于,所基于像素顶点的平行束投影方法包括:以图像中心为原点建立二维坐标系,计算每个像素中心的坐标,以探测器阵列中心为原点建立一维坐标系,计算每个探测器的坐标;在确定的投影角度下,计算每个像素中心投影到探测器上的坐标;计算投影的上底和下底的长度;计算像素四个顶点投影到探测器上的坐标;根据坐标判断像素顶点投影到哪些探测器上,计算对应的面积即为权值;
所述基于像素顶点平行束投影方法包括以下步骤:
(1)以图像中心为原点建立二维坐标系I,计算每个像素中心的坐标(xi,yi),i=1,2,…N,N为图像中像素的个数,以探测器阵列中心为原点建立一维坐标系D,计算每个探测器边界的坐标ti,i=1,2,…M,M为探测器的个数;
(2)在确定的投影角度θ下,计算每个像素中心投影到探测器上的坐标Pi:
Pi=xi·cosθ-yi·sinθ+r;
其中r是坐标系I和D的原点在横向上的距离:r=OI-OD;
(3)计算投影的上底和下底的长度w1和w2:
其中X是像素的边长,θ是投影角度;
(4)计算像素的四个顶点投影到探测器上的坐标:
(5)根据像素四个顶点投影之后的坐标判断该像素投影到哪些探测器上,将投影按面积分配到探测器上,每个探测器分配到的面积为像素对该探测器的权值贡献,投影和探测器的相对位置有十种情况,因此权值的计算也对应如下情况:w=dh;
(6)重复(2)至(5),完成所有投影角度的投影。
2.如权利要求1所述的基于像素顶点的平行束投影方法,其特征在于,所述(3)中计算投影上底和下底的长度时,将像素的两条对角线投影到探测器上,长度较长的一条为下底,长度较短的一条为上底;投影角度为0°,90°,180°或270°,两条对角线投影长度相等,此时像素的投影为矩形;投影角度为45°,135°,225°或315°时,其中一条对角线垂直于探测器,投影长度为0,像素的投影为等腰三角形。
3.如权利要求1所述的基于像素顶点的平行束投影方法,其特征在于,所述(5)中判断像素投影到哪些探测器上及对应的权值计算方法具体如下:某一探测器的左边界坐标为Dbe,右边界坐标为Dend;
1)如果Dend≤Pube满足,则像素没有投影到该探测器上;
6)如果Ptbe≤Dbe,Dend≤Ptend两式同时满足,则分配到该探测器的投影形状为矩形,使用w=dh计算权值;
12)如果Dbe≥Puend满足,则像素没有投影到该探测器上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810109560.9A CN108492341B (zh) | 2018-02-05 | 2018-02-05 | 一种基于像素顶点的平行束投影方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810109560.9A CN108492341B (zh) | 2018-02-05 | 2018-02-05 | 一种基于像素顶点的平行束投影方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108492341A CN108492341A (zh) | 2018-09-04 |
CN108492341B true CN108492341B (zh) | 2022-02-25 |
Family
ID=63344545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810109560.9A Active CN108492341B (zh) | 2018-02-05 | 2018-02-05 | 一种基于像素顶点的平行束投影方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108492341B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146071B (zh) * | 2019-05-16 | 2021-04-23 | 成都精位科技有限公司 | 定位维度提升方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7693318B1 (en) * | 2004-01-12 | 2010-04-06 | Pme Ip Australia Pty Ltd | Method and apparatus for reconstruction of 3D image volumes from projection images |
US7778392B1 (en) * | 2004-11-02 | 2010-08-17 | Pme Ip Australia Pty Ltd | Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs) |
CN102279970A (zh) * | 2010-06-13 | 2011-12-14 | 清华大学 | 基于gpu的螺旋锥束ct重建方法 |
CN103310471A (zh) * | 2012-03-09 | 2013-09-18 | 株式会社日立医疗器械 | Ct 图像生成装置及方法、ct 图像生成系统 |
CN103976753A (zh) * | 2013-02-08 | 2014-08-13 | 株式会社日立医疗器械 | Ct图像生成装置和ct图像生成方法 |
CN104599312A (zh) * | 2013-10-30 | 2015-05-06 | 北京莫高丝路文化发展有限公司 | 球幕多通道投影仿真方法 |
CN105118030A (zh) * | 2015-08-11 | 2015-12-02 | 上海联影医疗科技有限公司 | 医学图像金属伪影的校正方法及装置 |
CN105832358A (zh) * | 2016-03-16 | 2016-08-10 | 西安电子科技大学 | 一种基于系统校准的旋转双平板pet系统的成像方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104899827B (zh) * | 2015-05-26 | 2018-04-10 | 大连理工大学 | 基于固定分辨率条件下的离散Radon投影和Mojette投影转换方法 |
-
2018
- 2018-02-05 CN CN201810109560.9A patent/CN108492341B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7693318B1 (en) * | 2004-01-12 | 2010-04-06 | Pme Ip Australia Pty Ltd | Method and apparatus for reconstruction of 3D image volumes from projection images |
US7778392B1 (en) * | 2004-11-02 | 2010-08-17 | Pme Ip Australia Pty Ltd | Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs) |
CN102279970A (zh) * | 2010-06-13 | 2011-12-14 | 清华大学 | 基于gpu的螺旋锥束ct重建方法 |
CN103310471A (zh) * | 2012-03-09 | 2013-09-18 | 株式会社日立医疗器械 | Ct 图像生成装置及方法、ct 图像生成系统 |
CN103976753A (zh) * | 2013-02-08 | 2014-08-13 | 株式会社日立医疗器械 | Ct图像生成装置和ct图像生成方法 |
CN104599312A (zh) * | 2013-10-30 | 2015-05-06 | 北京莫高丝路文化发展有限公司 | 球幕多通道投影仿真方法 |
CN105118030A (zh) * | 2015-08-11 | 2015-12-02 | 上海联影医疗科技有限公司 | 医学图像金属伪影的校正方法及装置 |
CN105832358A (zh) * | 2016-03-16 | 2016-08-10 | 西安电子科技大学 | 一种基于系统校准的旋转双平板pet系统的成像方法 |
Non-Patent Citations (5)
Title |
---|
3D forward and back-projection for X-ray CT using separable footprints;Y Long 等;《IEEE Transactions on Medical Imaging》;20101130;第29卷(第11期);第1839-1850页 * |
FDK算法中一种新的插值方法;翟静 等;《中国体视学与图像分析》;20070930;第12卷(第3期);第221-225页 * |
Locality Preserving Projection Based on Endmember Extraction for Hyperspectral Image Dimensionality Reduction and Target Detection;Wang, Yiting 等;《APPLIED SPECTROSCOPY》;20160930;第70卷(第9期);第1573-1581页 * |
Novel registration for microcomputed tomography and bioluminescence imaging based on iterated optimal projection;Ma, Xibo 等;《JOURNAL OF BIOMEDICAL OPTICS》;20130228;第18卷(第2期);第026013:1-9页 * |
Pixel purity vertex component analysis;Rob Heylen 等;《2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)》;20171204;第213-216页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108492341A (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107993263B (zh) | 环视系统自动标定方法、汽车、标定装置及存储介质 | |
CN105913417B (zh) | 基于透视投影直线的几何约束位姿方法 | |
CN110163918A (zh) | 一种基于射影几何的线结构光标定方法 | |
CN101577002A (zh) | 应用于目标检测的鱼眼镜头成像系统标定方法 | |
CN104778716B (zh) | 基于单幅图像的卡车车厢体积测量方法 | |
CN104807405B (zh) | 一种基于光线角度标定的三维坐标测量方法 | |
CN114018932A (zh) | 基于矩形标定物的路面病害指标测量方法 | |
CN114359405A (zh) | 一种离轴沙姆3d线激光相机的标定方法 | |
CN114998448A (zh) | 一种多约束双目鱼眼相机标定与空间点定位的方法 | |
CN114283203A (zh) | 一种多相机系统的标定方法及系统 | |
CN108492341B (zh) | 一种基于像素顶点的平行束投影方法 | |
CN109506629B (zh) | 一种水下核燃料组件检测装置旋转中心标定的方法 | |
CN113096058B (zh) | 空间目标多源数据参数化仿真与MixCenterNet融合检测方法 | |
Pan et al. | Method of using non-specific target to complete structured light plane calibration | |
CN113658144A (zh) | 路面病害几何信息的确定方法、装置、设备及介质 | |
CN117611438A (zh) | 一种基于单目图像的2d车道线到3d车道线的重构方法 | |
CN107240133A (zh) | 一种立体视觉映射模型建立方法 | |
CN115239801B (zh) | 一种对象定位方法及装置 | |
TWI509564B (zh) | 三維成像的投影方法 | |
CN113804105B (zh) | 一种加速器射野几何与透视成像几何耦合校准装置及方法 | |
CN106156471A (zh) | 一种多层次流场可视化方法 | |
CN115409898A (zh) | 一种基于特殊环形标定板的高精度相机标定方法及设备 | |
CN114998496A (zh) | 一种基于场景航拍图及稀疏点云的正射影像快速生成方法 | |
CN111783024A (zh) | 一种中性原子图像的局部三维磁层离子通量分布反演方法 | |
WO2020143352A1 (zh) | 安检设备及其图像检测方法 |
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 |