CN102376097A - 非对称检测器ct迭代重建方法 - Google Patents
非对称检测器ct迭代重建方法 Download PDFInfo
- Publication number
- CN102376097A CN102376097A CN2010102620547A CN201010262054A CN102376097A CN 102376097 A CN102376097 A CN 102376097A CN 2010102620547 A CN2010102620547 A CN 2010102620547A CN 201010262054 A CN201010262054 A CN 201010262054A CN 102376097 A CN102376097 A CN 102376097A
- Authority
- CN
- China
- Prior art keywords
- projection
- data
- iteration
- positive
- negative
- 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.)
- Pending
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种非对称检测器CT迭代重建方法,其特征在于包括以下步骤:从CT扫描机获得的生数据经过空气校正、非线性校正、串扰校正、散焦校正以及定标校正获得投影数据;判断上述投影数据为断层扫描数据还是螺旋扫描数据;如为断层扫描数据,则对断层扫描数据进行分组;计算投影权重,得到权值矩阵;根据上述权值矩阵进行OSEM迭代正、反投影;判断是否满足迭代终止条件,如满足,则终止迭代正、反投影,得到初始CT图像;对上述初始CT图像进行后处理,得到重建的CT图像。本发明方法直接使用不完整投影数据进行迭代反投影重建,可以减轻伪影,抑制由插值造成的数据不一致性带来的噪声水平差异,主要适用于四排以下CT检测器。
Description
技术领域
本发明涉及医学图像处理技术,具体的说是一种非对称检测器CT迭代重建方法。
背景技术
在CT扫描中,为了增大扫描视野或减少辐射剂量,通常在检测器中心通道的一侧增加或者减少部分检测器模块,从而形成非对称检测器阵列,如图1所示,减少了部分检测器;相对于对称检测器结构,图中的非对称结构相当于有部分投影数据缺失,导致重建视野缩小。如图2所示,扇束几何重组为平行束几何,其中FOV_In为[-γm+α,γm-α]所能够覆盖的视野范围,FOV_Out为完整扫描视野与FOV_In之间的视野部分。该图仅为示意图,一个通道和它的对面通道不一定能完全对应,要补偿缺失部分的数据,目前的处理方法都是通过插值,利用其对面通道数据(Complementarydata)反投影,但这种方法在数据不一致性较大时会引入伪影,需要寻找适当的平滑过渡函数来抑制。
在统计计算中,最大期望(EM,Expectation-Maximization)算法是在概率模型中寻找参数最大似然估计的算法,其中概率模型依赖于无法观测的隐藏变量。最大期望算法经过两个步骤交替进行计算,第一步是计算期望(E),也就是将隐藏变量像能够观测到的一样包含在内从而计算最大似然的期望值;另外一步是最大化(M),也就是最大化在E步上找到的最大似然的期望值从而计算参数的最大似然估计。M步上找到的参数然后用于另外一个E步计算,这个过程不断交替进行。
传统的MLEM(Maximum Likelihood Expectation Maximization)方法计算式为:
其中:
μ(k)(j):第j个像素在第k次迭代中的衰减系数
P(i):第i个检测器单元的测量投影值
w(i,j):第j个像素对第i个投影的贡献权重
MLEM方法旨在寻找与测量的投影数据具有最大似然性(ML)的估计解,其迭代过程是由EM算法来实现的。由于是以统计规律为基础,MLEM重建法具有很好的抗噪声能力,尤其是在处理统计性差的数据时,更能显示出它相对于解析法的优越性,但是这种方法仍然存在迭代法的运算量大,运算时间长等缺点。
发明内容
针对现有技术中在CT扫描中心通道两侧检测器不对称的情况下(去掉一侧的部分通道),无法获得完整的一圈反投影数据的不足之处,本发明要解决的技术问题是提供一种可以减轻伪影和噪声水平差异的非对称检测器CT迭代重建方法。
为解决上述技术问题,本发明采用的技术方案是:
本发明一种非对称检测器CT迭代重建方法包括以下步骤:
从CT扫描机获得的生数据经过空气校正、非线性校正、串扰校正、散焦校正以及定标校正获得投影数据;
判断上述投影数据为断层扫描数据还是螺旋扫描数据;
如为断层扫描数据,则对断层扫描数据进行分组;
计算投影权重,得到权值矩阵;
根据上述权值矩阵进行OSEM迭代正、反投影;
判断是否满足迭代终止条件,如满足,则终止迭代正、反投影,得到初始CT图像;
对上述初始CT图像进行后处理,得到重建的CT图像。
如果不满足迭代终止条件,则返回根据上述权值矩阵进行迭代正、反投影步骤。
如为螺旋扫描数据,则对螺旋扫描数据经z方向滤波后重组到建像平面后再进行分组;
接续计算投影权重,得到权值矩阵步骤。
所述分组是对一圈投影数据中0~90度范围内的投影数据进行间隔方式的子集分组,对子集中的每一个角度做正、反投影的同时进行与该角度相隔90度、180度以及270度的其他三个角度的正、反投影。
所述OSEM迭代正、反投影为:在每一次迭代过程中将投影数据分成L个子集,每一个子集对重建图像各像素点值校正后,重建图像被更新一次,所有的子集运算一遍。
所述OSEM迭代正、反投影为直接使用不完整投影数据进行迭代正、反投影重建。
本发明具有以下有益效果及优点:
1.本发明方法直接使用不完整投影数据进行迭代反投影重建,可以减轻伪影,抑制由插值造成的数据不一致性带来的噪声水平差异,主要适用于四排以下CT检测器。
附图说明
图1为现有技术中减少部分检测器的非对称检测器阵列;
图2为现有技术中非对称平行束几何示意图;
图3为本发明方法流程图;
图4为本发明方法涉及的加权函数示意图;
图5为本发明方法涉及的平行束几何示意图
图6为本发明方法中FBP缺失数据一圈反投影;
图7为本发明方法中FBP缺失数据插值成半圈反投影;
图8为本发明方法中OSEM缺失数据一圈反投影。
具体实施方式
如图3所示,本发明非对称检测器CT迭代重建方法包括以下步骤:
从CT扫描机获得的生数据(Raw Data)经过空气校正、非线性校正、串扰校正、散焦校正以及定标校正获得投影数据;
判断上述投影数据为断层扫描数据还是螺旋扫描数据;
如为断层扫描数据,则对断层扫描数据进行分组;
计算投影权重,得到权值矩阵;
根据上述权值矩阵进行OSEM迭代正、反投影;
判断是否满足迭代终止条件,如满足,则终止迭代正、反投影,得到初始CT图像;
对上述初始CT图像进行后处理,得到重建的CT图像。
如果不满足迭代终止条件,则返回根据上述权值矩阵进行迭代正、反投影步骤;
如为螺旋扫描数据,则对螺旋扫描数据经z方向(即三维坐标系中的z轴方向)滤波后重组到建像平面后再进行分组;
接续计算投影权重,得到权值矩阵步骤。
1.投影数据的重组及分组
对于断层扫描数据,不需要进行重组,只需进行分组;而对于螺旋扫描数据,通过传统的z方向滤波的方法将不同z位置的投影数据重组到建像平面。
本发明采用正交的子集排序方法,将0~90度范围内的投影数据进行间隔方式的子集分组(将投影划分为N个子集,每隔M/N个投影选取一个投影,其中M为0~90度范围内的投影数),每次进行四个相互正交角度(即与该角度相隔90度、180度以及270度)的迭代正、反投影。
2.OSEM迭代正、反投影
为了加快收敛速度,减少运算时间,提高图像质量,人们提出了很多快速算法,其中OSEM(Ordered Subsets Expectation Maximization有序子集期望最大化)是很有应用前景的一种快速迭代重建算法,它是在最大似然期望法MLEM的基础上发展起来的。
所述OSEM迭代正、反投影为直接使用不完整投影数据进行迭代正、反投影重建,替代通过插值补充数据后进行滤波反投影重建的方法,可以抑制由数据不一致性带来的噪声水平差异,减轻伪影。
OSEM方法属于“块”迭代法中的一种,在每一次迭代过程中将投影数据分成L个子集,每一个子集对重建图像各像素点值校正以后,重建图像便被更新一次,所有的子集运算一遍,称为一次迭代过程。在MLEM方法一次迭代过程中,使用所有的投影数据对重建图像每一个像素点的值进行校正,重建图像被更新一次,而在OSEM方法一次迭代过程中重建图像被更新L次,所以OSEM方法具有加快收敛的作用。
k=0,初始化
重复下面步骤直到迭代到指定迭代次数为止。
w(i,j)为投影权重,μl(j)为像素衰减值;
对每一个像素计算:
C)
式(3)对i的求和是在第l个子集内进行的。
3.计算投影权重
本发明采用传统的函数加权方式:
射线γi进入物体并在经过的所有点处采样该物体(如图4所示),由于每一点并非都与图像网格点重合,所以就要以每一采样点为中心的加权函数h(x,y)来加权计算采样点值,落在加权函数h范围中的网格点xj都通过加权函数适当加权而结合进来。
w(i,j)=∫h(X(si)-X(xj),Y(si)-Y(xj))dsi (4)
采用贝塞尔-凯塞函数作为加权函数,其方程为:
式中:Im为第一类修正贝塞尔函数,m为贝塞尔函数的阶数,r为函数半径,d为距函数中心点距离(从像素xj到射线Pi的最短距离),α为控制函数形状的参数。
权值的闭合形式为:
非对称反投影几何示意图如图5所示。被扫描物体的阴影部分在投影P1下未被照射到,而在P2下被照射到,如果采取一圈滤波反投影的方式重建,这部分CT值由于滤波采样的不均衡偏低。图6是一个典型的水模非对称反投影的例子,去掉中心通道左侧部分通道,同时水模偏心放置,此时模体中的一部分不能在全部投影下都被完全照射到,因此造成CT值明显偏低。
要解决此问题,可以将数据重组到完整半圈,需进行插值,用对面通道数据补偿缺失数据,同时要在截断处需做平滑处理,否则由于数据的不一致性会引入明显的伪影,如图7所示:
如果用OSEM迭代重建方法进行缺失数据的一圈反投影,由于它是使用EM算法求参数的最大似然估计,而不是投影采样贡献累加,投影不均衡而造成部分CT值偏低的现象消失,同时也不需要进行插值,可以一定程度上抑制数据不一致带来的伪影及噪声差异,非对称OSEM迭代反投影图像结果如图8所示。
图像后处理包括图像定标、去环、去金属伪影、骨硬化校正等。
Claims (6)
1.一种非对称检测器CT迭代重建方法,其特征在于包括以下步骤:
从CT扫描机获得的生数据经过空气校正、非线性校正、串扰校正、散焦校正以及定标校正获得投影数据;
判断上述投影数据为断层扫描数据还是螺旋扫描数据;
如为断层扫描数据,则对断层扫描数据进行分组;
计算投影权重,得到权值矩阵;
根据上述权值矩阵进行OSEM迭代正、反投影;
判断是否满足迭代终止条件,如满足,则终止迭代正、反投影,得到初始CT图像;
对上述初始CT图像进行后处理,得到重建的CT图像。
2.按权利要求1所述的非对称检测器CT迭代重建方法,其特征在于:
如果不满足迭代终止条件,则返回根据上述权值矩阵进行迭代正、反投影步骤。
3.按权利要求1所述的非对称检测器CT迭代重建方法,其特征在于:
如为螺旋扫描数据,则对螺旋扫描数据经z方向滤波后重组到建像平面后再进行分组;
接续计算投影权重,得到权值矩阵步骤。
4.按权利要求1或3所述的非对称检测器CT迭代重建方法,其特征在于:所述分组是对一圈投影数据中0~90度范围内的投影数据进行间隔方式的子集分组,对子集中的每一个角度做正、反投影的同时进行与该角度相隔90度、180度以及270度的其他三个角度的正、反投影。
5.按权利要求1所述的非对称检测器CT迭代重建方法,其特征在于:所述OSEM迭代正、反投影为:在每一次迭代过程中将投影数据分成L个子集,每一个子集对重建图像各像素点值校正后,重建图像被更新一次,所有的子集运算一遍。
6.按权利要求1所述的非对称检测器CT迭代重建方法,其特征在于:所述OSEM迭代正、反投影为直接使用不完整投影数据进行迭代正、反投影重建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102620547A CN102376097A (zh) | 2010-08-25 | 2010-08-25 | 非对称检测器ct迭代重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102620547A CN102376097A (zh) | 2010-08-25 | 2010-08-25 | 非对称检测器ct迭代重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102376097A true CN102376097A (zh) | 2012-03-14 |
Family
ID=45794646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102620547A Pending CN102376097A (zh) | 2010-08-25 | 2010-08-25 | 非对称检测器ct迭代重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102376097A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103705263A (zh) * | 2012-10-01 | 2014-04-09 | 西门子公司 | 不对称探测的断层造影成像 |
WO2014172938A1 (zh) * | 2013-04-27 | 2014-10-30 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
CN104821002A (zh) * | 2014-02-05 | 2015-08-05 | 西门子公司 | 计算机断层成像中图像数据的迭代重建 |
CN105054959A (zh) * | 2015-08-12 | 2015-11-18 | 沈阳东软医疗***有限公司 | 空气校正方法、装置及设备 |
CN105615912A (zh) * | 2015-12-17 | 2016-06-01 | 沈阳东软医疗***有限公司 | 一种ct扫描方法和*** |
CN105701778A (zh) * | 2016-01-11 | 2016-06-22 | 赛诺威盛科技(北京)有限公司 | 一种从ct图像中去除金属伪影的方法 |
CN106572833A (zh) * | 2014-08-04 | 2017-04-19 | 株式会社日立制作所 | X射线ct装置、数据处理装置以及投影数据生成方法 |
CN109171790A (zh) * | 2018-09-26 | 2019-01-11 | 沈阳东软医疗***有限公司 | 一种ct扫描数据的处理方法、装置及ct机 |
CN110337673A (zh) * | 2017-01-05 | 2019-10-15 | 通用电气公司 | 基于深度学习的对在断层摄影重建中使用的数据的估计 |
CN110559001A (zh) * | 2019-09-03 | 2019-12-13 | 明峰医疗***股份有限公司 | 一种ct扫描仪散焦辐射强度分布测量方法 |
CN110998651A (zh) * | 2017-08-24 | 2020-04-10 | 爱克发有限公司 | 一种用于生成对象的增强断层扫描图像的方法 |
CN113168721A (zh) * | 2018-11-28 | 2021-07-23 | 皇家飞利浦有限公司 | 用于重建对象的图像的*** |
CN115115738A (zh) * | 2022-08-29 | 2022-09-27 | 威海市博华医疗设备有限公司 | 一种对肺癌影像成像的修正方法及装置 |
CN115294232A (zh) * | 2022-09-30 | 2022-11-04 | 浙江太美医疗科技股份有限公司 | 重建算法的识别方法、装置、电子设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1416781A (zh) * | 2002-12-30 | 2003-05-14 | 北京质子科技开发有限公司 | 正电子断层扫描中热源高分辨快速图像迭代重建方法 |
US20040174960A1 (en) * | 2003-03-05 | 2004-09-09 | Jiang Hsieh | Progressive updating approach for volumetric CT image reconstruction |
US20070075248A1 (en) * | 2005-09-30 | 2007-04-05 | Cardiovascular Imaging Technologies L.L.C. | Computer program, method, and system for hybrid ct attenuation correction |
CN101034479A (zh) * | 2006-03-10 | 2007-09-12 | Ge医疗***环球技术有限公司 | 图像重建方法和x射线ct设备 |
-
2010
- 2010-08-25 CN CN2010102620547A patent/CN102376097A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1416781A (zh) * | 2002-12-30 | 2003-05-14 | 北京质子科技开发有限公司 | 正电子断层扫描中热源高分辨快速图像迭代重建方法 |
US20040174960A1 (en) * | 2003-03-05 | 2004-09-09 | Jiang Hsieh | Progressive updating approach for volumetric CT image reconstruction |
US20070075248A1 (en) * | 2005-09-30 | 2007-04-05 | Cardiovascular Imaging Technologies L.L.C. | Computer program, method, and system for hybrid ct attenuation correction |
CN101034479A (zh) * | 2006-03-10 | 2007-09-12 | Ge医疗***环球技术有限公司 | 图像重建方法和x射线ct设备 |
Non-Patent Citations (3)
Title |
---|
H.MALCOLM HUDSON ET AL.: "Accelerated Image Reconstruction Using Ordered Subsets of Projection Data", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》, vol. 13, no. 4, 31 December 1994 (1994-12-31), pages 601 - 609 * |
曾理等: "锥束CT的图像分块OSEM重建算法", 《计算机工程与应用》, vol. 43, no. 32, 11 November 2007 (2007-11-11) * |
陈英茂等: "图像重建-有序子集最大期望值法", 《中华核医学杂志》, vol. 22, no. 6, 25 December 2002 (2002-12-25), pages 379 - 381 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103705263B (zh) * | 2012-10-01 | 2015-12-02 | 西门子公司 | 不对称探测的断层造影成像 |
CN103705263A (zh) * | 2012-10-01 | 2014-04-09 | 西门子公司 | 不对称探测的断层造影成像 |
US9439617B2 (en) | 2012-10-01 | 2016-09-13 | Siemens Aktiengesellschaft | Tomographic imaging with asymmetric detection |
WO2014172938A1 (zh) * | 2013-04-27 | 2014-10-30 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
CN104821002A (zh) * | 2014-02-05 | 2015-08-05 | 西门子公司 | 计算机断层成像中图像数据的迭代重建 |
CN106572833B (zh) * | 2014-08-04 | 2019-09-24 | 株式会社日立制作所 | X射线ct装置、数据处理装置以及投影数据生成方法 |
CN106572833A (zh) * | 2014-08-04 | 2017-04-19 | 株式会社日立制作所 | X射线ct装置、数据处理装置以及投影数据生成方法 |
CN105054959A (zh) * | 2015-08-12 | 2015-11-18 | 沈阳东软医疗***有限公司 | 空气校正方法、装置及设备 |
CN105054959B (zh) * | 2015-08-12 | 2018-06-15 | 沈阳东软医疗***有限公司 | 空气校正方法、装置及设备 |
US10238360B2 (en) | 2015-08-12 | 2019-03-26 | Shenyang Neusoft Medical Systems Co., Ltd. | Air calibration |
CN105615912A (zh) * | 2015-12-17 | 2016-06-01 | 沈阳东软医疗***有限公司 | 一种ct扫描方法和*** |
CN105615912B (zh) * | 2015-12-17 | 2020-06-12 | 东软医疗***股份有限公司 | 一种ct扫描方法和*** |
CN105701778A (zh) * | 2016-01-11 | 2016-06-22 | 赛诺威盛科技(北京)有限公司 | 一种从ct图像中去除金属伪影的方法 |
CN110337673A (zh) * | 2017-01-05 | 2019-10-15 | 通用电气公司 | 基于深度学习的对在断层摄影重建中使用的数据的估计 |
CN110337673B (zh) * | 2017-01-05 | 2023-10-20 | 通用电气公司 | 基于深度学习的在断层摄影重建中估计数据的方法和*** |
CN110998651A (zh) * | 2017-08-24 | 2020-04-10 | 爱克发有限公司 | 一种用于生成对象的增强断层扫描图像的方法 |
CN109171790A (zh) * | 2018-09-26 | 2019-01-11 | 沈阳东软医疗***有限公司 | 一种ct扫描数据的处理方法、装置及ct机 |
CN109171790B (zh) * | 2018-09-26 | 2022-03-15 | 东软医疗***股份有限公司 | 一种ct扫描数据的处理方法、装置及ct机 |
CN113168721A (zh) * | 2018-11-28 | 2021-07-23 | 皇家飞利浦有限公司 | 用于重建对象的图像的*** |
US11954769B2 (en) | 2018-11-28 | 2024-04-09 | Koninklijke Philips N.V. | System for reconstructing an image of an object |
CN110559001A (zh) * | 2019-09-03 | 2019-12-13 | 明峰医疗***股份有限公司 | 一种ct扫描仪散焦辐射强度分布测量方法 |
CN115115738A (zh) * | 2022-08-29 | 2022-09-27 | 威海市博华医疗设备有限公司 | 一种对肺癌影像成像的修正方法及装置 |
CN115115738B (zh) * | 2022-08-29 | 2022-11-15 | 威海市博华医疗设备有限公司 | 一种对肺癌影像成像的修正方法及装置 |
CN115294232A (zh) * | 2022-09-30 | 2022-11-04 | 浙江太美医疗科技股份有限公司 | 重建算法的识别方法、装置、电子设备和存储介质 |
CN115294232B (zh) * | 2022-09-30 | 2023-01-06 | 浙江太美医疗科技股份有限公司 | 重建算法的识别方法、装置、电子设备和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102376097A (zh) | 非对称检测器ct迭代重建方法 | |
CN103732147B (zh) | X射线计算机断层摄影装置以及图像重建方法 | |
US10896486B2 (en) | Denoising method and system for preserving clinically significant structures in reconstructed images using adaptively weighted anisotropic diffusion filter | |
CN103279929B (zh) | 一种基于余弦积分的ct图像金属轨迹预测和伪影去除方法 | |
CN101404088B (zh) | Ct图像重建的方法及*** | |
CN103366389A (zh) | Ct图像重建方法 | |
CN102103757B (zh) | 锥束图像重建方法及装置 | |
Levkovilz et al. | The design and implementation of COSEN, an iterative algorithm for fully 3-D listmode data | |
CN102973291B (zh) | C型臂半精确滤波反投影断层成像方法 | |
CN102831627A (zh) | 一种基于gpu多核并行处理的pet图像重建方法 | |
CN104751429B (zh) | 一种基于字典学习的低剂量能谱ct图像处理方法 | |
JP2013085965A (ja) | 円軌道コーンビームコンピュータ断層撮影(ct)においてアーチファクトを大幅に軽減するための方法およびシステム | |
CN102456227A (zh) | Ct图像重建方法及装置 | |
US20120275673A1 (en) | Truncation correction imaging enhancement method and system | |
CN102521853A (zh) | 一种反投影权重锥束ct重建算法 | |
CN103479379B (zh) | 一种倾斜螺旋扫描的图像重建方法及装置 | |
CN103679706A (zh) | 一种基于图像各向异性边缘检测的ct稀疏角度重建方法 | |
CN102331433A (zh) | 大尺寸工业长管道管壁的外部螺旋锥束ct扫描成像方法 | |
CN111696166A (zh) | 基于fdk型预处理矩阵的圆周锥束ct快速迭代重建方法 | |
CN102938154B (zh) | 一种基于粒子滤波的动态pet图像重建方法 | |
CN102488528B (zh) | 一种层析成像几何参数的校准方法 | |
CN106803237A (zh) | 一种基于极限学***均图像去噪方法 | |
CN109712213B (zh) | Pet图像重建方法、***、可读存储介质和设备 | |
US8121246B2 (en) | Radiographic apparatus and arithmetic processing program | |
CN116188615A (zh) | 一种基于正弦域和图像域的稀疏角度ct重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120314 |