CN101877125A - 一种基于小波域统计信号的图像融合处理方法 - Google Patents
一种基于小波域统计信号的图像融合处理方法 Download PDFInfo
- Publication number
- CN101877125A CN101877125A CN2009102430677A CN200910243067A CN101877125A CN 101877125 A CN101877125 A CN 101877125A CN 2009102430677 A CN2009102430677 A CN 2009102430677A CN 200910243067 A CN200910243067 A CN 200910243067A CN 101877125 A CN101877125 A CN 101877125A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- coefficient
- initial pictures
- wavelet coefficient
- fusion
- 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
Images
Landscapes
- Image Analysis (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种基于小波域统计信号处理的图像融合方法,属于图像处理的技术领域。首先对一组初始图像进行小波分解;对低频尺度系数采用加权平均的方法得到融合后的低频尺度系数;对于高频小波系数,在三个方向上建立小波域隐马尔可夫树模型,使用期望最大算法估计融合后的高频小波系数;对融合后低频尺度系数和高频小波系数进行小波逆变换,得到融合图像。本发明在处理高频小波系数时,设各初始图像的高频小波系数为“模糊”的真实图像的高频小波系数与非高斯噪声之和,并使用了连续的“模糊因子”,这简化了图像处理方法的运算,同时能够达到较好的融合效果。
Description
技术领域
本发明涉及一种基于小波域统计信号的图像融合处理方法,属于图像处理的技术领域。
背景技术
图像融合处理方法利用从同一场景获得的多幅初始图像之间信息的冗余性和互补性,最终得到一幅综合之前多幅初始图像优点的内容丰富的融合图像。其中,数码成像时对多聚焦图像的融合是融合处理方法的主要应用之一。
在图像融合处理方法中,基于小波变换的融合处理方法是其中很重要的一类,并已经得到广泛应用。由于人眼视觉***是在不同尺度上处理接收到的视觉信号,而且对细节信息非常敏感,基于小波变换的图像融合处理方法就是针对上述特性,将图像分解到多个尺度上,并对低频能量和高频细节分别进行处理,从而达到比其他融合处理方法更好的融合效果。其中,基于小波域统计信号的融合处理方法,将待融合的一组初始图像建模为被“模糊”的真实图像与非高斯噪声之和,再利用估计理论估计出真实图像(即融合图像)。
尽管小波变换通常被认为是“去相关”变换,但Crouse证实了小波系数之间存在着一定的相关性,它表现为尺度内的聚集性和尺度间的持续性。小波域隐马尔可夫树模型能够充分反映小波系数在尺度内和尺度间的相关性,利用这个模型的树状结构,使用期望最大的迭代算法,可以更准确估计出融合图像。
Blum提出初始图像的高频小波系数形成模型的“模糊因子”为离散值-1,0,1,这使得期望最大算法中迭代时更新“模糊因子”的步骤略为复杂。同时,在某些场景下,“模糊因子”的离散取值显得较为单一。
发明目的
本发明的目的是提出一种基于小波域统计信号的图像融合处理方法,根据待融合的初始图像的成像特点,使用连续的“模糊因子”,并且在期望最大算法更新“模糊因子”时使用显式更新方法,以此减少图像处理的运算量。
本发明提出的基于小波域统计信号的图像融合处理方法,包括以下步骤:
(1-1)对待融合的一组初始图像中的每个图像分别进行N层小波分解,分别得到每个初始图像的低频尺度系数和每个初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数;
(1-2)对初始图像的上述低频尺度系数进行加权融合,对显著性高的图像尺度系数赋予大的权值,得到融合的低频尺度系数;
(1-3)设各初始图像的高频小波系数为“模糊”的真实图像的高频小波系数与非高斯噪声之和,对上述初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数分别建立隐马尔可夫树模型,使用期望最大的迭代算法对真实图像的高频小波系数进行估计,得到融合的高频小波系数;
(1-4)将上述融合的低频尺度系数和融合的高频小波系数进行小波逆变换,得到融合图像。
上述方法中,步骤(1-2)对初始图像的低频尺度系数进行加权融合的方法,包括以下各步骤:
(2-1)设zi为第i幅初始图像的低频子带位于坐标(u,v)处的尺度系数,则初始图像的尺度系数的显著性为: 其中,
(2-2)根据上述尺度系数的显著性,得到坐标(u,v)处各幅初始图像的加权值为:
(2-3)根据上述各幅初始图像的加权值,得到坐标(u,v)处的融合低频尺度系数为:
上述方法中,步骤(1-3)得到融合的高频小波系数的方法,包括以下各步骤:
(3-1)设各初始图像在X方向上的高频小波系数为 其中,X=HL、LH或HH,zi,k(j)表示第i幅初始图像在X方向上的第k棵小波树上第j个节点处的小波系数,k=1,…,K,K为初始图像在X方向上的小波树总数,j=1,…,P,P为每棵小波树上的节点总数,一般P=(4N-1)/3,wk(j)表示真实图像在X方向上的第k棵小波树上第j个节点处的小波系数,即待估计的融合后小波系数,βi,k(j)为“模糊因子”,-1≤βi,k(j)≤1,εi,k(j)为加性非高斯噪声,由零均值和M态高斯混合模型拟合得到;
(3-2)使用期望最大算法推导出迭代公式,迭代估计上述参数,得到真实图像在X方向上的小波系数wk(j),即融合后的小波系数。
本发明提出的基于小波域统计信号的图像融合处理方法,具有以下优点:
(1)本发明方法在参数估计时,使用了小波域隐马尔可夫树模型,考虑了小波系数尺度内和尺度间的相关性,因此估计的参数更准确,得到的融合图像效果更好。
(2)本发明方法中的“模糊因子”采用连续值,更符合待融合的初始图像成像的实际情况,这使得对初始图像的高频小波系数建模更加准确,因此能够得到较好的融合图像。
(3)本发明方法中的期望最大算法在更新“模糊因子”时,有显式表达,使得计算简化,减少了图像处理的运算量。
附图说明
图1是本发明的图像融合处理方法的原理框图。
图2是本发明方法中建立的高频小波系数的小波树树状模型的平面图。其中,S1是一棵小波树的根节点,ρ(i)是节点i的父节点。
图3是单棵小波树的树状结构图。其中,S1是这棵小波树的根节点,ρ(i)是节点i的父节点。
图4是本发明方法中对高频小波系数进行融合的方法流程图。
图5和图6分别是两幅待融合的初始图像。其中图5在右边大表面清晰聚焦,左边小表面失焦。图6相反。
图7是利用本发明方法对图5和图6的图像进行融合后的融合图像。
具体实施方式:
本发明提出的基于小波域统计信号的图像融合处理方法,其原理框图如图1所示,包括以下步骤:
(1-1)对待融合的一组初始图像中的每个图像分别进行N层小波分解,分别得到每个初始图像的低频尺度系数和每个初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数,一般N=3或4,图1所示,N=3;
(1-2)对初始图像的上述低频尺度系数进行加权融合,对显著性高的图像尺度系数赋予大的权值,得到融合的低频尺度系数;
(1-3)设各初始图像的高频小波系数为“模糊”的真实图像的高频小波系数与非高斯噪声之和,对上述初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数分别建立隐马尔可夫树模型,使用期望最大的迭代算法对真实图像的高频小波系数进行估计,得到融合的高频小波系数;
(1-4)将上述融合的低频尺度系数和融合的高频小波系数进行小波逆变换,得到融合图像。
本发明方法中,对所有低频尺度系数进行加权融合的方法,包括以下各步骤:
(2-1)设zi为第i幅初始图像的低频子带位于坐标(u,v)处的尺度系数,则初始图像的尺度系数的显著性为: 其中,
(2-2)根据上述尺度系数的显著性,得到坐标(u,v)处各幅初始图像的加权值为:
(2-3)根据上述各幅初始图像的加权值,得到坐标(u,v)处的融合低频尺度系数为:
本发明方法中,得到融合的高频小波系数的方法,包括以下各步骤:
(3-1)设各初始图像在X方向上的高频小波系数为zi,k(j)=βi,k(j)wk(j)+εi, k(j),其中,X=HL、LH或HH,HL为水平方向,LH为垂直方向,HH为对角方向,zi,k(j)表示第i幅初始图像在X方向上的第k棵小波树上第j个节点处的小波系数,k=1,…,K,K为初始图像在X方向上的小波树总数,j=1,…,P,P为每棵小波树上的节点总数,一般P=(4N-1)/3,wk(j)表示真实图像在X方向上的第k棵小波树上第j个节点处的小波系数,即待估计的融合后小波系数,βi,k(j)为“模糊因子”,由于初始图像可以认为由真实图像部分模糊形成,变换到小波域则表现为对真实图像的小波系数的正负加权,所以,这里βi,k(j)可以在[-1,1]之间连续取值,εi,k(j)为加性非高斯噪声,由零均值和M态高斯混合模型拟合得到;
(3-2)使用期望最大算法推导出迭代公式,迭代估计上述参数,得到真实图像在X方向上的小波系数wk(j),即融合后的小波系数。
上述得到融合的高频小波系数的方法中,对水平方向、垂直方向和对角方向上的小波系数分别建立模型,如图2所示。其中,S1是一棵小波树的根节点,一般选取最细尺度(J=1)上的小波系数作为小波树的根节点;节点ρ(i)与节点i互为父、子节点。在这个树状模型中,尺度J(1≤J<N)上的每一个节点都在尺度J+1上有四个子节点。设尺度J上的节点在子带中的坐标为(u0,v0),则它在尺度J+1上的子节点坐标分别为(2u0-1,2v0-1),(2u0,2v0-1),(2u0-1,2v0)和(2u0,2v0)。
上述期望最大的迭代算法在估计真实图像在X方向上的小波系数时,还要涉及到一系列的参数,包括:状态概率 表示噪声εi,k(j)为状态m的概率,同时由于εi,k(j)wk(j)是确定值,所以也表示zi,k(j)为状态m的概率,这里认为在X方向上对任意的i,k,是一样的,因此可简写为Ps(j)(m),其中m=1,…,M;状态转移概率 表示父节点ρ(j)处状态为n时其子节点j处状态为m的概率,其中n=1,…,M;σm 2(j)表示εi,k(j)为状态m时的高斯方差。上述所有的未知参数构成了待估计的模型参数集: 在期望最大算法中,上述模型参数集被迭代更新,从而得到真实图像在X方向上的小波系数wk(j),即融合后的小波系数。
上述期望最大算法中,对模型参数进行估计的具体步骤包括:
(1)初始化模型参数集Φ;
(2)利用模型参数集Φ中的参数,使用前向后向算法计算P(Si,k(j)=m|z,Φ)和P(Si,k(j)=m,Si,k(ρ(j))=n,|z,Φ),其中,z表示所有初始图像在X方向上的所有小波系数;
(3)根据(2)步得到的概率值,逐步更新模型参数集,并记为Φ′;
(4)若||Φ′-Φ||充分小,则循环停止;否则,Φ=Φ′,回到(2)。
上述初始化模型参数集,包括以下步骤:
(1)状态概率Ps(j)(m)=1/M,其中m=1,…,M;
(2)状态转移概率 其中m=1,…,M,n=1,…,M;
(3)将初始图像的小波系数的加权平均值作为真实图像小波系数wk(j)的初始值,即 其中,权值λi,k(j)由初始图像小波系数的显著性来决定。小波系数的显著性为 。其中 (u,v)表示zi,k(j)在子带中的坐标。这时,权值
(4)“模糊因子”βi,k(j)=1,其中i=1,…,q,k=1,…,K,j=1,…,P,即每个初始图像的小波系数都可以完全反映真实图像的小波系数;
(5)高斯方差σm 2(j)的选取。令 这里取γ=10。同时,各状态方差要满足
上述使用前向后向算法计算P(Si,k(j)=m|z,Φ)和P(si,k(j)=m,Si,k(ρ(j))=n,|z,Φ)时,假设了小波树之间彼此独立,所以需要对每个初始图像在X方向上的每棵小波树分别计算概率值。对于单棵小波树的计算,描述时可省略zi,k(j),βi,k(j)和wk(j)的下标,同时设小波分解的层数为N,J=1为最细尺度,J=N为最粗尺度,ρ(j)是节点j的父节点,c(j)是节点j的子节点的集合。前向后向算法包括以下步骤:
(1)设Tj是以j节点为根的子树。若Ti是Tj的子树,则Tj\i是Tj除去Ti。由此定义:
Bj(m)≡f(Tj|S(j)=m,Φ)
Bj,ρ(j)(m)≡f(Tj|S(ρ(j))=m,Φ)
Bρ(j)\j(m)≡f(Tρ(j)\j|S(ρ(j))=m,Φ)
Aj(m)≡p(S(j)=m,T1\j|Φ)
对于M态混合高斯模型,设z(j)为状态m的概率密度函数为:
(2)前向算法中,初始化最细尺度(J=1)上节点j处的似然函数:
由此向粗尺度循环计算,对尺度J上的节点有:
J=J+1
直至J=N时停止;
(3)后向算法中,初始化最粗尺度(J=N)根节点处:
A1(m)=PS(1)(m)
由此向细尺度循环计算,对尺度J上的节点有:
J=J-1
直至J=1时停止;
(4)由上述求出的A和B值,最终得到条件概率:
上述更新模型参数集的过程,包括以下步骤:
(1)更新状态概率PS(j)(m)和状态转移概率aj,ρ(j) m,n
(2)更新高斯方差σm 2(j)
(3)更新“模糊因子”βi,k(j)
设 则
(4)更新wk(j)值
以图5和图6两幅初始图像为例,使用本发明方法进行融合,得到图7的融合结果。在两幅初始图像中,两个表面分别聚焦,而得到的融合图像在两个表面都清晰聚焦,得到了比较好的融合结果。
Claims (3)
1.一种基于小波域统计信号的图像融合处理方法,其特征在于该方法包括以下步骤:
(1-1)对待融合的一组初始图像中的每个图像分别进行N层小波分解,分别得到每个初始图像的低频尺度系数和每个初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数;
(1-2)对初始图像的上述低频尺度系数进行加权融合,对显著性高的图像尺度系数赋予大的权值,得到融合的低频尺度系数;
(1-3)设各初始图像的高频小波系数为“模糊”的真实图像的高频小波系数与非高斯噪声之和,对上述初始图像在水平方向、垂直方向和对角方向上的N层高频小波系数分别建立隐马尔可夫树模型,使用期望最大的迭代算法对真实图像的高频小波系数进行估计,得到融合的高频小波系数;
(1-4)将上述融合的低频尺度系数和融合的高频小波系数进行小波逆变换,得到融合图像。
2.如权利要求1所述的方法,其提特征在于所述的步骤(1-2)对初始图像的低频尺度系数进行加权融合的方法,包括以下各步骤:
(2-1)设zi为第i幅初始图像的低频子带位于坐标(u,v)处的尺度系数,则初始图像的尺度系数的显著性为: 其中,
(2-2)根据上述尺度系数的显著性,得到坐标(u,v)处各幅初始图像的加权值为:
(2-3)根据上述各幅初始图像的加权值,得到坐标(u,v)处的融合低频尺度系数为:
3.如权利要求1所述的方法,其提特征在于所述的步骤(1-3)得到融合的高频小波系数的方法,包括以下各步骤:
(3-1)设各初始图像在X方向上的高频小波系数为zi,k(j)=βi,k(j)wk(j)+εi,k(j),其中,X=HL、LH或HH,zi,k(j)表示第i幅初始图像在X方向上的第k棵小波树上第j个节点处的小波系数,k=1,…,K,K为初始图像在X方向上的小波树总数,j=1,…,P,P为每棵小波树上的节点总数,一般P=(4N-1)/3,wk(j)表示真实图像在X方向上的第k棵小波树上第j个节点处的小波系数,即待估计的融合后小波系数,βi,k(j)为“模糊因子”,-1≤βi,k(j)≤1,εi,k(j)为加性非高斯噪声,由零均值和M态高斯混合模型拟合得到;
(3-2)使用期望最大算法推导出迭代公式,迭代估计上述参数,得到真实图像在X方向上的小波系数wk(j),即融合后的小波系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910243067 CN101877125B (zh) | 2009-12-25 | 2009-12-25 | 一种基于小波域统计信号的图像融合处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910243067 CN101877125B (zh) | 2009-12-25 | 2009-12-25 | 一种基于小波域统计信号的图像融合处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101877125A true CN101877125A (zh) | 2010-11-03 |
CN101877125B CN101877125B (zh) | 2013-02-13 |
Family
ID=43019670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910243067 Expired - Fee Related CN101877125B (zh) | 2009-12-25 | 2009-12-25 | 一种基于小波域统计信号的图像融合处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101877125B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306381A (zh) * | 2011-06-02 | 2012-01-04 | 西安电子科技大学 | 基于Beamlet与小波变换的图像融合方法 |
CN102567977A (zh) * | 2011-12-31 | 2012-07-11 | 南京理工大学 | 一种基于小波的红外偏振图像自适应融合方法 |
CN102908122A (zh) * | 2012-09-17 | 2013-02-06 | 广州市伟迈机电科技有限公司 | 数字化仪成像及图像拼接处理方法 |
CN103198456A (zh) * | 2013-03-21 | 2013-07-10 | 西安电子科技大学 | 基于方向波域隐马尔可夫树模型的遥感图像融合方法 |
CN103679670A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 一种基于改进模型的pcnn多源图像融合方法 |
CN109447930A (zh) * | 2018-10-26 | 2019-03-08 | 太原科技大学 | 小波域光场全聚焦图像生成算法 |
CN113191210A (zh) * | 2021-04-09 | 2021-07-30 | 杭州海康威视数字技术股份有限公司 | 一种图像处理方法、装置及设备 |
WO2022088976A1 (zh) * | 2020-10-29 | 2022-05-05 | 北京达佳互联信息技术有限公司 | 图像处理方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1254770C (zh) * | 2004-05-13 | 2006-05-03 | 上海交通大学 | 基于期望值最大与离散小波框架的图像融合方法 |
CN100395777C (zh) * | 2005-12-20 | 2008-06-18 | 上海电力学院 | 基于最优小波滤波器组的多传感器图像融合方法 |
-
2009
- 2009-12-25 CN CN 200910243067 patent/CN101877125B/zh not_active Expired - Fee Related
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102306381A (zh) * | 2011-06-02 | 2012-01-04 | 西安电子科技大学 | 基于Beamlet与小波变换的图像融合方法 |
CN102567977A (zh) * | 2011-12-31 | 2012-07-11 | 南京理工大学 | 一种基于小波的红外偏振图像自适应融合方法 |
CN102567977B (zh) * | 2011-12-31 | 2014-06-25 | 南京理工大学 | 一种基于小波的红外偏振图像自适应融合方法 |
CN102908122A (zh) * | 2012-09-17 | 2013-02-06 | 广州市伟迈机电科技有限公司 | 数字化仪成像及图像拼接处理方法 |
CN103679670B (zh) * | 2012-09-25 | 2016-08-31 | 中国航天科工集团第二研究院二〇七所 | 一种基于改进模型的pcnn多源图像融合方法 |
CN103679670A (zh) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | 一种基于改进模型的pcnn多源图像融合方法 |
CN103198456B (zh) * | 2013-03-21 | 2015-06-03 | 西安电子科技大学 | 基于方向波域隐马尔可夫树模型的遥感图像融合方法 |
CN103198456A (zh) * | 2013-03-21 | 2013-07-10 | 西安电子科技大学 | 基于方向波域隐马尔可夫树模型的遥感图像融合方法 |
CN109447930A (zh) * | 2018-10-26 | 2019-03-08 | 太原科技大学 | 小波域光场全聚焦图像生成算法 |
CN109447930B (zh) * | 2018-10-26 | 2021-08-20 | 太原科技大学 | 小波域光场全聚焦图像生成算法 |
WO2022088976A1 (zh) * | 2020-10-29 | 2022-05-05 | 北京达佳互联信息技术有限公司 | 图像处理方法及装置 |
CN113191210A (zh) * | 2021-04-09 | 2021-07-30 | 杭州海康威视数字技术股份有限公司 | 一种图像处理方法、装置及设备 |
CN113191210B (zh) * | 2021-04-09 | 2023-08-29 | 杭州海康威视数字技术股份有限公司 | 一种图像处理方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN101877125B (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101877125B (zh) | 一种基于小波域统计信号的图像融合处理方法 | |
CN110378844B (zh) | 基于循环多尺度生成对抗网络的图像盲去运动模糊方法 | |
CN104217406B (zh) | 一种基于剪切波系数处理的sar图像降噪方法 | |
CN107133923B (zh) | 一种基于自适应梯度稀疏模型的模糊图像非盲去模糊方法 | |
CN105913393A (zh) | 一种自适应小波阈值图像去噪算法及装置 | |
CN105279740A (zh) | 一种基于稀疏正则化的图像去噪方法 | |
Petrovic et al. | Enforcing integrability for surface reconstruction algorithms using belief propagation in graphical models | |
CN102567973A (zh) | 基于改进的形状自适应窗口的图像去噪方法 | |
CN104008537A (zh) | 结合cs-ct-chmm的噪声图像融合新方法 | |
WO2023197601A1 (zh) | 一种基于梯度场的点云修复方法 | |
CN112884671A (zh) | 一种基于无监督生成对抗网络的模糊图像复原方法 | |
CN112232134A (zh) | 一种基于沙漏网络结合注意力机制的人体姿态估计方法 | |
CN111754447A (zh) | 基于多状态上下文隐马尔科夫模型的红外和可见光图像融合方法 | |
CN116563146A (zh) | 基于可学习曲率图的图像增强方法与*** | |
CN112712019A (zh) | 一种基于图卷积网络的三维人体姿态估计方法 | |
CN102810202B (zh) | 基于分数阶差分加权的图像多步残差反馈迭代滤波方法 | |
CN109544567B (zh) | 基于两阶段区域合并策略的极化合成孔径雷达图像分割方法 | |
CN114998107A (zh) | 一种图像盲超分辨网络模型、方法、设备及存储介质 | |
CN113780389B (zh) | 基于一致性约束的深度学习半监督密集匹配方法及*** | |
CN107341449A (zh) | 一种基于云块特征变化的静止气象卫星降水估算方法 | |
CN106934398A (zh) | 基于超像素聚类和稀疏表示的图像去噪方法 | |
CN102136134B (zh) | 基于mrf先验的sar图像去斑方法 | |
CN102314687B (zh) | 一种红外序列图像中的小目标检测方法 | |
CN116167947A (zh) | 一种基于噪声水平估计的图像降噪方法 | |
CN115913995A (zh) | 一种基于卡尔曼滤波修正的云服务动态QoS预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130213 Termination date: 20161225 |