基于盲源分离从x射线序列图像中分离图像信息的方法
【技术领域】
[0001] 本发明专利涉及X射线图像降噪，具体为基于盲源分离将图像噪声和图像信息视 为不相关分量而将图像信息分离出来的一种图像降噪算法。
【背景技术】
[0002] 自从伦琴在1895年发现X射线以来，X射线成像技术使探索物体内部结构成为可 能。典型的X射线成像系统由射线源、X射线影像探测器、CCD/CM0S、影像输出及其他光电 器件，由这些光电器件产生的热噪声、散粒噪声和量子噪声等噪声会严重降低最终输出图 像质量。
[0003] X射线图像噪声遵循泊松分布，而常用的图像降噪方法主要分为空域、时间域、频 率域几类降噪方法。
[0004] 盲源分离是一种基于统计的信号处理技术，最初的经典应用为"鸡尾酒会"。利 用盲源分离可以分析多维数据，因此，盲源分离也逐渐被应用于图像处理方面，比如图像融 合、图像增强、特征提取、伪影消除、混合图像分离、消除散乱线等。而利用盲源分离进行图 像降噪通常需要噪声的先验知识或者训练含噪声图像进行稀疏编码收缩。
【发明内容】
[0005] 技术问题：本发明专利的目的在于提供一种基于盲源分离从X射线序列图像中分 离图像信息的方法，具体采用定点独立分量分析（Fixed-Point Independent Component Analysis，FPICA)或奇异值分解（Singular Value Decomposition, SVD)，将 X 射线图像中 图像信息有效提取出来的图像降噪处理。
[0006] 技术方案：本发明是一种基于盲源分离从X射线序列图像中分离图像信息的方 法，对一个相对静止物体拍摄含有m帧图像的序列，由于图像序列中每帧图像均由相同性 质的噪声和信息组成，只有噪声的随机出现使不同图像间均存在差异，将含有m帧图像的 序列看作是图像噪声和图像信息的m种组合，利用盲源分离得到n-1帧图像噪声分量和1 帧图像信息分量，全部η个分量中标准方差值最大者即为图像信息分量，据此判断并输出 图像信息分量，m, η为整数，且m彡η彡2。
[0007] 利用盲源分离出n-1帧图像噪声分量和1帧图像信息分量，其方法是利用盲源分 离中高阶统计的定点独立分量分析或二阶统计奇异值分解进行分量分离，具体方法如下：
[0008] 对静止目标拍摄一个包含不小于2帧的X射线图像序列，序列中任一帧中的图像 信息是唯一的，噪声是随机的，因此每帧图像都不同；利用盲源分离，每帧图像可看作1个 图像信息分量和n-1个噪声分量的组合，如式（1)所示：
[0009]
[0010] 式⑴中，X = [X1, X2,…，Xm]T, Xi表示观察图像序列中任何一帧；下标i表示每帧 在图像序列中的次序，取值I, ;X表示观察图像序列矩阵;A = Iia1, a2, - ,aj，a;是分 量混合的系数向量，由系数ay组成；S = [s D s2,…，8」1'表示由图像分量s 组成的矩阵；下 标j表不分量序号，取值1，…，n ;m, η为整数，且m多η多2 ;T表不矩阵转置；
[0011] 通过定点独立分量分析（Fixed-Point Independent Component Analysis， FPICA)或奇异值分解（Singular Value Decomposition, SVD)方法均可得到矩阵W，从而得 到无限逼近分量8]的分量y ,，其关系由公式（2)所示：
[0012] Y = WX = WAS -S (2)
[0013] 式⑵中，Y = Iiy1, y2,…，yn]T，y_j表不分呙出的噪声分量或图像信息分量，下标 j表不分量的序号，Y表不由分离出的分量组成的矩阵，T表不矩阵转置；W表不分离矩阵。
[0014] 所述对一个相对静止物体拍摄含有多帧图像的序列，取不小于2帧X射线图像形 成图像序列。
[0015] 将所述的图像噪声和图像信息视为组成整幅图像的独立分量，不同帧图像视为噪 声和图像信息的不同组合。
[0016] 有益效果：X射线图像降噪目的是尽最大可能从含噪声的原图像中提取图像有用 信息，利用其进行后续图像分析等工作。
[0017] 传统的去噪方式将噪声和有用信息同等对待，认为噪声和有用信息在某个变换域 中是分布于不同区间的。然而，噪声和有用信息图像在某个区域中通常是相互叠加的。因 此，这些传统方法无法滤除设定噪声区间以外的噪声，并且可能滤除该区间的有用信息。分 离可以有效处理有用信息和噪声之间的关系，确保分离出的有用信息分量和噪声分量之间 相互独立，可以在去噪的同时，尽量保留图像有用信息，保证降噪后的图像质量。
[0018] 本发明是一种既不需要噪声先验知识也不需要对含噪声图像进行训练的图像降 噪处理方法，仅利用盲源分离的本质即可简单有效的将在图像中唯一的、固定的图像信息 从图像序列中提取出来。
【具体实施方式】
[0019] 利用X射线成像设备拍摄一个图像序列，所拍摄目标处于相对静止状态；利 用盲源分离中高阶统计的定点独立分量分析方法或者二阶统计的奇异值分解对图像序 列进行分离；计算被分离出的每个分量的标准方差（Standard Deviation, SD)或方差 (Variance)，值最大者即为所求的图像有用信息分量，而多个噪声分量的标准方差或方差 值彼此接近你，但相对图像有用信息的标准方差或方差却小。
[0020] 该方法可以利用盲源分离，将含有m帧的图像序列中每一帧作为噪声和信息的不 同组合，可将含有m帧图像的序列分离出n-1帧噪声分量和1帧图像信息分量；无需噪声先 验信息和图像训练即可将图像信息分量分离并输出，该方法不仅可简单实现X射线图像降 噪，且降噪效率且可随图像序列中帧数量增加和盲源分离的提高而提高。
[0021] 可以利用盲源分离中高阶统计的定点独立分量分析（Fixed-Point Independent Component Analysis，FPICA)或二阶统计奇异值分解（Singular Value Decomposition， SVD)进行分量分呙。
[0022] 取相对静止目标的不小于2帧X射线图像形成图像序列。
[0023] 将图像噪声和图像信息视为组成整幅图像的独立分量，不同帧图像视为噪声和图 像信息的不同组合。
[0024] 利用盲源分离使噪声和信息分开，获取标准方差（Standard Deviation, SD)或方 差（Variation)最大的分量即为图像信息分量，并将其输出。
[0025] 实例：
[0026] 对静止目标拍摄一个包含不小于2帧的X射线图像序列帧的X射线图像序列，每 一帧图像均由固定图像信息和随机图像噪声组成，每帧图像中，图像信息是唯一的，噪声随 机的，因此，序列中每帧图像各不相同，可以视为由1个有用图像信息分量和n-1个噪声分 量的组合。利用盲源分离（Blind Source S印aration，BSS)对图像序列中噪声分量和有用 信息分量分离。以下公式（1)中，X1表示图像序列中每一帧图像 ；S ]为图像分量；a u为分 量的权重系数。
[0027]
[0028] 式⑴中，X = [Xl，X2，…，xJT，X1表示观察图像序列中任何一帧；下标i表示每帧 在图像序列中的次序，取值1，;X表示观察图像序列矩阵;A = Iia1, a2, - ,aj，a;是分 量混合的系数向量，由系数ay组成；S = [s D s2,…，8」1'表示由图像分量s 组成的矩阵；下 标j表不分量序号，取值1，…，n ;m, η为整数，且m多η多2 ;T表不矩阵转置。
[0029] 通过定点独立分量分析（Fixed-Point Independent Component Analysis， FPICA)或奇异值分解（Singular Value Decomposition, SVD)方法均可得到矩阵W，从而得 到无限逼近分量8]的分量y ,，其关系由公式（2)所示：
[0030] Y = WX = WAS - S (2)
[0031] 式⑵中，Y = Iiy1, y2,…，yn]T，y_j表不分呙出的噪声分量或图像信息分量，下标 j表不分量的序号，Y表不由分离出的分量组成的矩阵，T表不矩阵转置；W表不分离矩阵。 在η个y」分量中，分别计算其标准方差（Standard Deviation, SD)或方差（Variance)，由 于图像信息分量的标准方差（方差）远大于噪声分量的标准方差（方差），各个噪声分量的 标准方差（方差）数值非常接近，因此可判断标准方差（方差）最大的为有用图像信息，可 以作为结果输出。
【主权项】
1. 一种基于盲源分离从X射线序列图像中分离图像信息的方法，其特征在于，对一个 相对静止物体拍摄含有m帧图像的序列，由于图像序列中每帧图像均由相同性质的噪声和 信息组成，只有噪声的随机出现使不同图像间均存在差异，将含有m帧图像的序列看作是 图像噪声和图像信息的m种组合，利用盲源分离得到n-1帧图像噪声分量和1帧图像信息 分量，全部η个分量中标准方差值最大者即为图像信息分量，据此判断并输出图像信息分 量，m, η为整数；且m彡η彡2。2. 根据权利要求1所述的基于盲源分离从X射线序列图像中分离图像信息的方法，其 特征在于利用盲源分离出n-1帧图像噪声分量和1帧图像信息分量，其方法是利用盲源分 离中高阶统计的定点独立分量分析或二阶统计奇异值分解进行分量分离，具体方法如下： 对静止目标拍摄一个包含不小于2帧的X射线图像序列，序列中任一帧中的图像信息 是唯一的，噪声是随机的，因此每帧图像都不同；利用盲源分离，每帧图像可看作1个图像 信息分量和n-1个噪声分量的组合，如式（1)所不：式⑴中，X = [X1, X2,…，xm]T, Xi表示观察图像序列中任何一帧；下标i表示每帧在图 像序列中的次序，取值1,…，m ;X表示观察图像序列矩阵;A = Iia1, a2,…，an]，8;是分量混 合的系数向量，由系数ay组成；S = [s D S2，…，8」1'表示由图像分量s 组成的矩阵；下标j 表不分量序号，取值1，…，n ;m, η为整数，且m多η多2 ;T表不矩阵转置； 通过定点独立分量分析FPICA或奇异值分解SVD方法均可得到矩阵W，从而得到无限逼 近分量Sj的分量y」，其关系由公式⑵所不： Y = WX = WAS - S (2) 式⑵中，Y = Iiy1, y2,…，yn]T，y_j表不分尚出的噪声分量或图像信息分量，下标j表 示分量的序号，Y表示由分离出的分量组成的矩阵，T表示矩阵转置;W表示分离矩阵。3. 根据权利要求1所述的基于盲源分离从X射线序列图像中分离图像信息的方法，其 特征在于所述对一个相对静止物体拍摄含有m帧图像的序列，取不小于2帧X射线图像形 成图像序列。4. 据权利要求1所述的基于盲源分离从X射线序列图像中分离图像信息的方法，其特 征在于将所述的图像噪声和图像信息视为组成整幅图像的独立分量，不同帧图像视为噪声 和图像信息的不同组合。
【专利摘要】本发明提出了一种基于盲源分离从X射线序列图像中分离图像信息的方法，该方法可以利用盲源分离中高阶统计的定点独立分量分析或二阶统计的奇异值分解进行分离；取相对静止目标的不小于2帧X射线图像形成图像序列；将图像噪声和图像信息视为组成整幅图像的独立分量；图像信息是稳定的，而噪声在不同帧图像中出现的随机性使不同帧图像可视为噪声和图像信息的不同组合；利用盲源分离使噪声和图像信息分开，计算标准方差最大的分量即为图像信息分量，并将其输出。该方法不需要噪声先验信息，将可将有用图像信息提取出来，提取效率随图像序列中图像数量增加而提高，也将随盲源分离的技术改进而提高，是一种实用且有潜力的X射线图像降噪方法。
【IPC分类】G06T5/00
【公开号】CN105160643
【申请号】CN201510632896
【发明人】喻春雨
【申请人】南京邮电大学
【公开日】2015年12月16日
【申请日】2015年9月29日