技术领域
本说明书涉及对非均匀介质进行非侵入性超声表征的方法和系统，尤其适用于医学成像或非破坏性测试，并且更普遍而言，适用于可以使用超声成像的所有领域。
背景技术
在声成像领域，试图用超声波主动探测未知介质来对其进行表征。这尤其是医学成像中回波描记术的原理。
声成像系统的分辨率可以被定义为辨别对象的微小细节的能力。原则上，声成像系统受衍射限制，并且其理论分辨率由λ/2(其中λ是声在介质中的波长)或检测器的有限角孔径给出。然而，实际上，分辨率通常由于声速在传播介质中的变化而降低。
具体而言，通常在声成像中，介质被认为是均匀的，其具有恒定的声速c
在过去的几十年中，已经开发出各种方法来校正像差(这主要发生在医学成像领域)，以便在超声图像的分辨率由于介质的非均匀性而降低的情况下获得具有最佳分辨率的超声图像。然而，在复数像差的情况下，波前畸变会根据图像的区域而变化。符合等晕域条件的这些区域通常无法在理论上确定，因此构成了应用现有技术的基本限制。
如图1A至图1C所示，传统的回波描记方法使用压电换能器11的阵列10，这些压电换能器11可以独立地发射和接收超声脉冲。每个换能器的位置由向量u标识。当将这样的阵列面对要研究的介质20放置时，可以以各种方式对介质20进行声穿透。
对要研究的介质进行声穿透的第一种方法是从阵列中的一个换能器发射超声脉冲，该换能器的位置通过u
已知还有其他方法可以对要研究的介质进行声穿透，这些技术称为“波束成形”，指的是通过应用延迟线对波进行相干控制。在这些技术中，已知有聚焦发射，这是常规的回波描记术广泛使用的一项技术。如图1B(左图)所示，仅需对换能器11施加适当的一组延迟以校正行进时间，从而使所有脉冲一起到达目标焦点r
最近开发的另一种波束成形技术包括用一系列平面波对介质进行声穿透。图1C示出了这种“平面波”回波描记法的原理，该原理例如在G.Montaldo等人的文章《Coherentplane-wave compounding for very high frame rate ultrasonography and transientelastography》(IEEE Trans.Ultrason.,Ferroelect.Freq.Control 56 489-506,2009)中进行了描述。在发射时向每个信号施加延迟(图1C，左图)，以形成相对于换能器阵列10以角度θ
受1950年代天文学中使用的光学像差校正技术的启发，为了校正超声成像中的相位畸变，引入了自适应聚焦技术(例如，参见D.J.Philips,S.W.Smith,O.T.Von Ramm andF.L.Thurstone,"A phase compensation technique for B-Modeechoencephalography".Ultrasound in Medicine,1,345-404,1975)。
如图2所示，该校正技术适用于所有常规的超声成像方法(图1A至图1C)。该技术假设存在薄的畸变层22，其在换能器11上引起简单的时间偏移。因此，可以通过在波束成形中通常使用的延迟中添加其他延迟定律来补偿这些时间偏移。如图2所示，该校正技术包括在不考虑像差的情况下介质20的激发(25)、像差定律的估计(26–27)以及在考虑畸变器的情况下新的发射/接收(28–29)。通过使用发射器和检测器之间的反馈回路，可以实时迭代系统，以找到用于校正像差的最佳时移。已经提出了几种方法，用于基于超声散斑(也就是说，由欠分辨散射体的随机分布进行的超声反射所产生的信号)来计算要施加到阵列中的每个换能器的延迟。一种方法是基于相邻换能器接收的信号之间的相关性(例如，参见M.O'Donnell and S.W.Flax，"Phase aberration measurements in medical ultrasound:Human studies",Ultrasonic Imaging 10,1-11,1988)。另一种方法是基于超声散斑的对比度优化((Nock,G.E.Trahey and S.W.Smith"Phase aberration correction inmedical ultrasound using speckle brightness as a quality factor",J.Acoust.Soc.Am,Vol.85,5,pp 1819-1833,1989)。然而，这些校正技术是基于薄的畸变层22的假设，该薄的畸变层22位于尽可能靠近换能器阵列10的位置。因此，对于体积性畸变器，这些方法仅在单个等晕域上有效。
等晕域的概念如图3A和图3B所示。在换能器阵列10与焦平面之间存在畸变器22。在图3A的示例中，向换能器阵列10的换能器11施加延迟，使得可以校正在畸变器22的输出处的声波前部(圆柱形波)。这导致在焦平面的给定点r
除了校正场限于单个等晕域外，由于波场的每个频谱分量可能会经历不同的相位畸变，上述像差校正技术通常不适用于超宽带信号(例如用于超声成像的信号)。为了对超声成像中的像差进行超宽带校正，已经开发出更复杂的方法，例如从时间反转概念衍生的逆滤波或匹配滤波方法。这些方法包括通过合适的时空滤波器对发射和接收的信号进行卷积，该时空滤波器使得可以优化波在介质中的聚焦并校正其像差。然而，要了解此过滤器，要求存在强大的散射体(C.Prada and M.Fink,"Eigenmodes of the time reversaloperator",Wave Motion,20,151-163,1994)或位于要成像的介质内的尤其具有侵入性的检测器阵列(例如参见M.Tanter,J.-L.Thomas and M.Fink."Time reversal and theinverse filter",J.Acoust.Soc.Am.108,223-234,2000)。
为了迭代地校正超声散斑中的像差，一种最新的方法直接受益于反射矩阵及其相关矩阵。因此，在J-L.Robert等人的文章("Green's function estimation in speckleusing the decomposition of the time reversal operator:Application toaberration correction in medical imaging",J.Acoust.Soc.Am.,123,2(2008))中，测量在输入处的聚焦发射基与输出处的换能器基之间的反射矩阵R
在本说明书的上下文中，提出了一种原创方法，该方法既不基于假星的产生，也不基于迭代方法。所提出的方法还使得能够同时校正散斑图像的各个等晕域中的像差，还可以包括更多会引起入射波的镜面反射的确定性对象(界面、大散射体等)。
发明内容
根据第一方面，本说明书涉及一种用于对非均匀介质进行非侵入性超声表征的方法，该方法包括：
-通过至少一个第一换能器阵列在所述非均匀介质的区域中产生一系列入射超声波的步骤；
-记录在输入处的发射基和输出处的换能器基之间定义的实验性反射矩阵的步骤；
-基于所述实验性反射矩阵确定至少一个第一宽带反射矩阵的步骤，所述第一宽带反射矩阵在输入处和输出处以聚焦基在第一谱带中定义；
-确定在所述聚焦基和观察基之间在所述第一谱带中定义的第一畸变矩阵的步骤，所述第一畸变矩阵是由在所述聚焦基和像差校正基之间定义的所述第一宽带反射矩阵与针对模型介质在与上述相同的基之间定义的参考反射矩阵的相位共轭矩阵的逐项矩阵乘积直接获得或在基变换之后获得；
-基于所述第一畸变矩阵确定所述非均匀介质中的所述声穿透区的物理参数的至少一个映射的步骤。
如将在说明书的其余部分更详细地描述的，所述非均匀介质的物理参数可以包括例如与所述声穿透区中包含的至少一个等晕域相关联的像差定律、表征所述介质的局部反射的参数、表征在所述介质中的声速的参数、多重散射率。
在本说明书的含义内的非均匀介质包括任何包含散射体的介质，该散射体相对于周围介质表现出声阻抗的中断，所述周围介质也可能表现出在空间上不均匀的声速(例如多层介质)。在本说明书中描述的方法和系统尤其使得可以对这种非均匀介质进行非侵入性深度超声表征。
超声波波前通过非均匀介质传播的像差对应于来自该介质的波前与在理想情况下预期的波前之间的差异。这些像差可以例如与声速在所述介质中的空间波动相关联。例如，在回波描记术的情况下，所述声波在到达目标器官之前穿过各种类型的组织。由于软组织具有不同的机械性能，所述声速远非均匀。这导致声波波前的像差，该像差导致所得到的超声图像产生畸变，并因此导致其分辨率和对比度下降。这些像差可能会损及医疗检查。
最佳像差校正基是使所述声穿透区中包含的等晕域的尺寸最大化的基。例如，在薄的畸变器位于所述换能器阵列附近的情况下，像差校正基可以是换能器基，或者在平移不变的多层介质的情况下，像差校正基可以是平面波基。在空间不均匀的体积性畸变器的一般情况下，无论所考虑的像差校正基是什么，都存在多个等晕域。因此研究所述畸变矩阵使得可以确定与这些区域中的每一个相关联的像差定律。
所述聚焦基可以包含在相对于所述换能器阵列位于所述介质的给定深度处的单个焦平面中，或者可以通过同时考虑位于不同深度的焦点，以体积来定义。
记录了所述实验性反射矩阵的发射基可以等同于所述换能器基、所述平面波基(“傅立叶基”)或任何聚焦基。
可以在各个基之间定义反射矩阵或所述畸变矩阵，通过在时域中应用延迟线或通过在频域中与传播矩阵进行矩阵乘积来执行基变换。在说明书的其余部分，术语“聚焦”反射矩阵或“聚焦”畸变矩阵分别是指在输入处和输出处以聚焦基定义的反射矩阵或畸变矩阵。
记录的所述实验性反射矩阵是“实数”矩阵(也就是说，是由时域中的实系数组成的矩阵)，每个换能器记录的电信号是实数。然而，所述宽带聚焦反射矩阵是“复数”矩阵(也就是说，由复数分析信号形成的系数组成的矩阵)。可以通过仅保留所记录的电信号频谱的正频率来获得形成所述宽带聚焦反射矩阵的复数分析信号。这些复数分析信号的幅角与所述介质反射的波的相位相吻合。
模型介质或参考介质是例如声速在其中恒定的均匀介质。在理论上可以针对该参考介质建立作为传播矩阵的参考反射矩阵(或更简单而言，“参考矩阵”)(也就是说，针对该参考介质用每个焦平面上的虚拟平面镜获得的超声波的反射矩阵)。根据传播介质的理论知识程度，所述参考介质可以采用更复杂的形式(例如多层介质等)。
在聚焦基和像差校正基之间定义的第一宽带反射矩阵与针对模型介质定义的参考反射矩阵的相位共轭矩阵的逐项矩阵乘积或“哈达玛积”，等于从所述第一宽带反射矩阵中的每个元素的相位减去所述参考反射矩阵中的相应元素的相位。具体而言，一个矩阵或第一矩阵的相位共轭(在本说明书中更简称为“共轭”)的元素是复数，该复数的模与所述第一矩阵的相同，但相位相反。因此，从所述第一宽带反射矩阵的每个元素的相位减去预期的弹道分量(由所述模型介质定义)，这使得可以针对每个输入聚焦发射隔离由所述介质反射的场的畸变分量。申请人已经表明，分析所述畸变矩阵使得尤其可以辨别在所述声穿透区中包含的等晕域，并且可以在观察平面中确定与每个等晕域相关联的像差定律。
根据一个或多个示例性实施例，所述观察基是所述像差校正基。可以试图确定所述聚焦基和所述像差校正基之间的第一畸变矩阵，例如以量化和/或校正像差。
在这种情况下，所述第一畸变矩阵的确定可以包括：
-在输入处或输出处对所述第一宽带聚焦反射矩阵进行基变换，以获得在所述聚焦基和所述像差校正基之间确定的所述第一宽带反射矩阵；并且
-对所述第一宽带反射矩阵与在相同的基之间定义的参考反射矩阵的相位共轭矩阵进行逐项矩阵乘积。
根据一个或多个示例性实施例，所述观察基是聚焦基，并且参考“聚焦的”第一畸变矩阵。可以试图确定“聚焦”畸变矩阵(也就是说，在输入处和输出处以聚焦基定义的矩阵)，例如以确定扩散函数(或“点扩散函数”，PSF)的映射，扩散函数被定义为在输入处和输出处的聚焦基之间的空间脉冲响应。基于所述PSF，例如可以对在所研究介质的声穿透区中的聚焦质量(即像差水平)、声速或多重散射率进行局部量化。
可以通过从在所述聚焦基和所述像差校正基之间定义的所述第一畸变矩阵的基变换来获得所述聚焦的第一畸变矩阵。
或者，可以基于所述第一宽带聚焦反射矩阵，借助于所述第一宽带聚焦反射矩阵的每个行和/或列与以相同的基定义的参考反射矩阵的相同行和/或列之间的空间相关性来直接构建所述聚焦的第一畸变矩阵。基于所述聚焦的第一畸变矩阵，仅通过基变换就可以确定在所述聚焦基和所述像差校正基之间定义的所述第一畸变矩阵。
根据一个或多个示例性实施例，所述第一宽带聚焦反射矩阵的确定包括：
-在输入处和输出处对以频域表示的所述实验性反射矩阵进行基变换，以形成聚焦反射矩阵；
-在所述聚焦反射矩阵的所述第一谱带上进行相干加和，以形成所述第一宽带聚焦反射矩阵。
或者，可以通过将时移应用于所述实验性反射矩阵来获得所述第一宽带聚焦反射矩阵。然而，这种方法不能分析不同谱带上的实验数据，这在像差的频率分散的情况下可能是不利的。
根据一个或多个示例性实施例，所述非均匀介质的物理参数的至少一个映射的确定包括：
-确定所述第一畸变矩阵在所述聚焦基中的不变量，以在所述聚焦基中识别至少一个第一等晕域；
-对于识别出的每个第一等晕域，确定至少一个第一像差定律在所述校正基中的映射。
为了基于所述第一畸变矩阵来确定所述等晕域和相关的像差定律，可以采用几种方法。根据一个或多个示例性实施例，对这些不变量的确定包括所述第一畸变矩阵的奇异值分解、所述第一畸变矩阵归一化后(也就是说，其每个元素的模已经被归一化，但是其相位被保留)的奇异值分解或者所述第一畸变矩阵的相关矩阵归一化后(也就是说，所述第一畸变矩阵的相关矩阵的每个元素的模已经被归一化)的奇异值分解。
特别假定存在非均匀介质进行的镜面反射，申请人已经表明，可以通过所述第一畸变矩阵的奇异向量的线性组合来获得非均匀介质关于整个声穿透区的局部反射率的映射。
申请人还已经表明，所述畸变矩阵的奇异值分解使得可以过滤信号子空间(以其行和/或其列之间具有实质相关性为特征的矩阵)中的噪声子空间(其行和列之间不具有相关性的随机矩阵)，所述噪声子空间既包含实验噪声，又包含由聚焦基上游发生的多重散射事件引起的反射场的非相干贡献。
根据一个或多个示例性实施例，超声表征方法还包括在所述观察基中确定由所述一个或多个第一像差定律校正的第一反射矩阵。尤其是在假定样品发生随机散射式反射时，这使得可以确定所述介质在像差校正后的局部反射率的映射或“图像”。
根据一个或多个示例性实施例，基于所述声穿透区的所述校正后的第一反射矩阵，可以确定像差校正后的畸变矩阵，该校正后的畸变矩阵在所述像差校正基中与在所述像差校正基中确定的所述校正后的反射矩阵和所述参考反射矩阵的相位共轭矩阵的逐项矩阵乘积相对应。所述像差校正后的畸变矩阵使得可以通过确定其在聚焦基中的不变量，在识别出的所述第一等晕域中完善校正。还使得可以在聚焦基中识别至少一个第二等晕域，并且对于识别出的每个第二等晕域，确定第二像差定律在所述像差校正基中的映射。因此，取决于声穿透区中包含的等晕域数量，可以根据需要多次迭代该方法，以获得介质在像差校正后的局部反射率的映射或“图像”。当假定随机散射或中间反射(也就是说，镜面反射和随机散射性反射的混合反射)时，该迭代过程尤其适用。
根据一个或多个示例性实施例，对所述非均匀介质的物理参数的至少一个映射的确定包括针对识别出的至少一个第一等晕域，在所述聚焦基中确定第一扩散函数(PSF)。扩散函数对应于在像差校正平面中测量的像差定律的空间傅里叶变换。其在每个等晕域内在空间上都是不变的。根据一个示例，所述扩散函数可以从所述聚焦畸变矩阵(也就是说，在输入处和输出处以聚焦基所定义的)的奇异值分解获得。
根据一个或多个示例性实施例，所述非均匀介质的物理参数的至少一个映射的确定包括确定在所述声穿透区中的声速。这涉及例如研究所述畸变矩阵的第一奇异值作为速度模型的函数的变化，所述速度模型用于对参考矩阵进行建模。最佳速度模型是使畸变矩阵的第一奇异值最大化的模型。一种可选方案是研究聚焦畸变矩阵。具体而言，这使得可以研究扩散函数作为用于计算参考矩阵的速度模型的函数的变化。最佳速度模型是使波长归一化后的扩散函数的宽度最小化的模型。
根据一个或多个示例性实施例，所述非均匀介质的物理参数的至少一个映射的确定包括确定在所述声穿透区中的多重散射率。例如，这涉及研究像差已经被预先校正的聚焦畸变矩阵。超出焦点后获得的非相干强度(多重散射)的平均值与在焦点处的强度(所述畸变矩阵的中心行)的比值可以给出局部多重散射率。
根据一个或多个示例性实施例，所述表征方法还包括识别和/或消除反射场的镜面分量和/或在换能器阵列与非均匀介质的各个界面或层之间产生的多重反射。为此，可以在输入处和输出处将畸变矩阵投影到傅立叶基上。在此基础上，所述场的镜面反射分量和多重反射分量看起来是精确的反射角和入射角对。因此，它们可以很容易被过滤，仅保留所述反射场的随机散射分量(散斑)。然后，对随机散射分量的这种区分允许直接获得像差定律，该像差定律将被应用在输入处和输出处以校正所述反射矩阵并获得所述介质的最佳图像，这在镜面分量占主导时是不可能的。另外，在傅立叶平面中过滤所述畸变矩阵使得可以消除多重反射(探头与人体的表面层之间的多次回波)，该多重反射污染例如在浅深度处的超声图像。
根据一个或多个示例性实施例，所述表征方法还包括：
-基于所述实验性反射矩阵确定至少一个第二宽带聚焦反射矩阵的步骤，所述第二宽带聚焦反射矩阵在输入处和输出处以聚焦基在与所述第一谱带不同的第二谱带中定义；
-确定在所述聚焦基和所述观察基之间在所述第二谱带中定义的至少一个第二畸变矩阵的步骤，所述第二畸变矩阵是由在所述聚焦基和所述像差校正基之间定义的所述第二宽带反射矩阵与针对模型介质在与上述相同的基之间定义的参考反射矩阵的相位共轭矩阵的逐项矩阵乘积直接获得或在基变换之后获得；
在不同谱带上定义的畸变矩阵可用于确定相同数量的校正后的反射矩阵，然后可以将所述校正后的反射矩阵相加以形成单个宽带反射矩阵。该频谱分析使得可以校正由声速在所述介质中的频率分散而引起的色差。在不同谱带上获得的所述多个畸变矩阵也可用于确定所述声速或所述多重散射率的频谱映射。
根据第二方面，本说明书涉及一种用于对非均匀介质进行非侵入性超声表征的系统，该系统被配置为执行如上所述的所有示例性超声表征方法。根据第二方面的超声表征系统包括：
-至少一个第一换能器阵列，其被配置为在所述非均匀介质的区域中产生一系列入射超声波，并记录由所述区域反向散射的超声波随时间的变化；
-计算单元，其与所述至少一个第一换能器阵列耦合，并被配置为：
ο记录在输入处的发射基和输出处的换能器基之间定义的实验性反射矩阵；
ο基于所述实验性反射矩阵确定至少一个第一宽带聚焦反射矩阵，所述第一宽带聚焦反射矩阵在输入处和输出处以聚焦基在第一谱带中定义；
ο确定在所述聚焦基和观察基之间在所述谱带中定义的第一畸变矩阵，所述第一畸变矩阵是由在所述聚焦基和所述校正基之间确定的所述第一宽带反射矩阵与针对模型介质在与上述相同的基之间定义的第一参考反射矩阵的相位共轭矩阵的逐项矩阵乘积直接获得或在基变换之后获得；
ο基于所述第一畸变矩阵确定所述非均匀介质的物理参数的至少一个映射。
根据本说明书的表征系统可以包括既能发射和又能接收的至少一个换能器阵列，或者包括多个换能器阵列，其中一些专用于超声波的发射而其他的专用于超声波的接收。
附图说明
通过阅读下面的详细描述并参照附图，上述呈现的其他优点和技术特征将变得显而易见。
-图1A至图1C(已描述)示出了用于超声成像和定量化的已知发射/接收机制；
-图2(已描述)示出了根据现有技术基于自适应聚焦来校正超声成像中的像差的技术；
-图3A和图3B(已描述)示意性地示出了等晕域的概念；
-图4示出了用于实施根据本说明书的超声表征方法的示例性超声表征系统；
-图5示出了根据本说明书的方法的一个步骤，旨在根据一个示例性实施例，基于实验性反射矩阵来确定第一宽带聚焦反射矩阵；
-图6示出了为了以聚焦基记录的反射矩阵的投影而进行的基变换矩阵运算；
-图7A、图7B示出了根据本说明书的方法的一个步骤，旨在根据一个示例性实施例，基于宽带聚焦反射矩阵来确定第一畸变矩阵；
-图8示出了在没有像差以及具有像差的情况下体模的超声图像以及相应的聚焦畸变矩阵，该聚焦畸变矩阵在给定深度处确定并且对应于图像中的线。
-图9A和图9B分别示出了在如图8所示的具有像差的情况下与体模相对应的畸变矩阵的奇异值分布(图9A)，以及相关联的第一特征向量在傅立叶平面中的相位(图9B)；
-图10示出了畸变矩阵的示例性用途(体模上的像差水层产生单个等晕域的情况)，更精确而言，示出了未经/经过校正的超声图像、在给定深度处确定并对应于图像中的线的相应聚焦畸变矩阵以及未经/经过校正的扩散函数的表示，该扩散函数是在声穿透区上取平均值并且由聚焦畸变矩阵获得；
图11示出了畸变矩阵的示例性用途(在产生多个等晕域的人体乳房体模上)，并且更精确而言，示出了常规超声图像和使用畸变矩阵的第一向量和畸变矩阵的第二向量校正像差之后的超声图像；
图12A至图12C示出了畸变矩阵用于过滤多重反射的示例性用途，更精确而言，图12A是描述实验装置的示意图(超声探头和体模之间放置有机玻璃板)，图12B是针对均匀速度模型(c＝1800m/s)在从图像中去除镜面反射分量和多重反射前后的超声图像、聚焦畸变矩阵在过滤前后针对超声图像中的一条线的示例、投影到平面波基中的相同畸变矩阵，图12C是说明镜面反射分量和多重反射分量在平面波基中的特征的图。
图13示出了畸变矩阵用于声速的定量测量的示例性用途，并且更精确而言，示出了所研究体模的超声图像，聚焦畸变矩阵针对图像中的线以及均匀速度模型(c＝1540m/s)的示例、在输入处的焦平面和傅立叶平面之间的对应畸变矩阵、通过波长归一化的扩散函数的宽度随波速c变化的函数、畸变矩阵的归一化后的第一奇异值随波速c变化的函数；
图14示出了畸变矩阵用于建立在介质中的声速的映射的示例性用途，并且更精确而言，示出了体模针对均匀速度模型(c＝1540m/s)且穿过有机玻璃板的情况下得到的常规超声图像，以及过滤多重反射并针对图像的每一行优化速度模型后的相同图像、声速在深度范围内的分布；
图15示出了畸变矩阵用于多重散射的局部量化的示例性用途，并且更具体而言，示出了标准的体模超声图像、解释了在畸变矩阵中出现多重散射的情况的图、针对超声图像中的线的聚焦畸变矩阵、从前述矩阵推导的显示与多重散射分量相关联的非相干背景的扩散函数示例、超声图像的每个像素中的多重散射率；
将参照附图描述各个实施例，在这些实施例中，相似或相同的元件采用相同的附图标记。
具体实施方式
在下面的详细描述中，仅详细描述了一些实施例以确保描述的清楚性，但是这些示例并不旨在限制从本说明书中所示出的原理的一般范围。
本说明书中描述的各个实施例和方面可以以多种方式进行组合或简化。特别而言，除非另有说明，否则各种方法中的步骤可以重复、反向或并行执行。
图4示出了根据本说明书的用于实现对非均匀介质20的超声表征方法的示例性超声表征系统40。系统40包括换能器11的至少一个第一阵列10，例如线性或二维阵列。换能器是例如超声压电换能器，其可以采用与介质20接触的刚性棒的常规形式。换能器阵列形成例如探头设备41的一部分。换能器阵列与计算单元42连接，该计算单元42本身可以与显示设备43连接；计算单元向/从每个换能器11传输并记录电信号。然后，超声换能器将这些电信号转换成超声波，反之亦然。计算单元42被配置用于实施计算或处理步骤，尤其是用于实施根据本说明书的方法的步骤。
在图4中，如同在本说明书的其余部分，参考了用于发射和接收的换能器阵列，应当理解，在更一般的情况下，可以同时使用多个换能器阵列。这些换能器阵列可以是既发射也接收，或者仅其中某些用于发射而其他用于接收。
在本说明书中，当参考尤其用于实施方法步骤的计算或处理步骤时，应理解，每个计算或处理步骤可以由软件、硬件、固件、微代码或这些技术的任何适当组合来实现。当使用软件时，每个计算或处理步骤可以通过计算机程序指令或软件代码来实施。这些指令可以存储在或传输至计算机(或计算单元)可读的存储介质和/或可以由计算机(或计算单元)执行，以实施这些计算或处理步骤。
本说明书中所示的各种实验结果(图5至图15)取自两个超声体模，即再现生物组织的机械特性(c＝1540m/s)并产生生物组织的超声散斑特征的合成样品。第一体模使得可以量化成像系统在分辨率和对比度方面的性能。具体而言，该第一体模包括位于各种深度和横向位置的近点状目标以及直径为5mm的非均匀回声内含物(图5至图10以及图12至图15)。第二体模再现了乳腺组织的形状和特性，并显示出明亮的回声，模仿了同一组织内部的微钙化现象(图11)。
本说明书描述了用于对非均匀样品进行非侵入性超声表征的方法和系统。这些方法和系统基于至少一个第一矩阵(在本说明书的其余部分被称为“畸变矩阵”)的确定。根据一些示例，图5至图7示出了用于确定畸变矩阵的方法的步骤。
如图5所示，步骤51，非侵入性超声表征方法包括通过换能器11(图4)的阵列10在非均匀介质的区域中产生一系列入射超声波的步骤以及记录在输入处的发射基i和输出处的换能器基u之间定义的实验性反射矩阵R
所述方法还包括基于所述实验性反射矩阵R
在图5的示例中，离散傅立叶变换操作52使得可以在频域R
图6示出了将反射矩阵R
G′
其中r
其中r
发射和接收时波束成形的运算可以解释为回声向其相关联的散射体的反向传播。在频域中，这理解为传播矩阵的相位共轭。因此，可以通过以下矩阵乘积从实验性反射矩阵推导出聚焦反射矩阵R
其中符号×表示矩阵乘积，指数
下一步是对在谱带Δω
τ是与以聚焦基表示的虚拟换能器阵列相关联的回声时间。宽带聚焦反射矩阵中的每个系数R(r
图7A、图7B示出了根据本说明书的方法的一个步骤，旨在根据一个示例性实施例，基于宽带聚焦反射矩阵R
为了估计波前的像差分量，首先，选择用于校正像差的最佳基。例如，在薄的畸变器位于换能器阵列附近的情况下，换能器基看来是最合适的。在所述表征方法的该说明性示例中所选择的情况下，非均匀介质是平移不变的多层介质，平面波基似乎更合适。这例如是对肝脏采用回波描记术的情况，其中声波必须穿过一系列的肌肉层和脂肪组织层。在下文中考虑该示例，但是本说明书绝不限于从傅立叶基对像差进行校正。
图像71(图7A)示出了在存在像差的情况下记录的宽带聚焦反射矩阵R
R
其中考虑了谱带Δν
其中符号o表示逐项矩阵乘积(哈达玛矩阵乘积)，指数
上面执行的运算由图7B示出。包括从测得的波前减去确定的几何分量，从而隔离其畸变分量。矩阵D
在模型介质不均匀的情况下，参考反射矩阵G′
图7A中的图像73显示了从反射矩阵推导出的畸变矩阵D
也可以在焦平面中研究畸变矩阵。可以通过以下基变换，基于在观察平面中表示的畸变矩阵D
其中
因此，式(8)就矩阵系数而言重写如下：
其中，k
D
作为说明，图8通过回溯图像80的实验图示出了在没有像差的情况下获得的体模的超声图像81。图像82示出了在给定深度处确定的并对应于图像81中的线的聚焦畸变矩阵。图像83是在存在像差(像差水层22，图85)的情况下获得的体模的超声图像，并且图像84示出了在给定深度处确定的并对应于图像83中的线的聚焦畸变矩阵。在没有像差的实验中，由介质反射的所有能量都出现在畸变矩阵的中心行(r
为了进一步分析畸变矩阵D
假定样品引起随机散射性反射，下面描述使用畸变矩阵来确定一个或多个像差定律在观察基中的映射(尤其是建立声穿透区的局部反射的映射)的示例。
为此，寻求畸变矩阵中的不变量(换言之，在声穿透区的等晕域内在空间上不变的像差定律)。本领域技术人员已知各种方法用于搜索这种矩阵中的不变量，例如奇异值分解(或“SVD”)或主成分分析(“PCA”)。
对于提取矩阵的行或列之间的相关性来说，奇异值分解是强大工具。数学上，大小为N×M的矩阵D
U和V是大小为N×M的单位矩阵，其列U
矩阵Σ的对角线元素是矩阵D
因此，矩阵D
其中L＝min(M,N)。
SVD主要将矩阵分解为两个子空间：信号子空间(以其行和/或列之间具有实质相关性为特征的矩阵)和噪声子空间(在其行和列之间不具有相关性的随机矩阵)。信号子空间与最大的奇异值相关联，而噪声子空间与最小的奇异值相关联。一方面，D的SVD因此可以过滤噪声子空间，该噪声子空间既包含实验噪声也包含由多重散射事件引起的反射场的非相干贡献。另一方面，信号子空间的每个奇异状态将使得可以根据输出特征向量U
U
下面的图示出了畸变矩阵的各种应用，无论该畸变矩阵是以焦点基定义还是在焦点基与像差校正基之间定义。
图9A和图9B示出了将SVD应用于图8的图像84中所示的畸变矩阵D
为了校正宽带反射矩阵R
R
然后通过在矩阵R
其中R′
在物理上，相位共轭运算包括重新发射由与所测得的畸变的相位相反的相位调制的波前，以及在输出处施加类似的校正。因此，该操作使得可以完美补偿波在其向外和返回路径上积累的相位畸变。然后可以经由在输入和输出处从瞳孔平面到焦平面的基变换从R′
I′(r
因此，I′(r
图10通过比较像差的超声图像101和根据式(17)描述的校正获得的超声图像103来示出该校正的益处。特别是借助于在校正像差后看起来对比度更强并且分辨率更好的回声目标，可以观察到这两个图像之间的定性增益。
然而，在回波描记术中，大多数时候没有回声目标来判断图像的分辨率。为了量化后者，可以使用聚焦畸变矩阵。从初始的反射矩阵推导出的D
基于聚焦畸变矩阵，还可以确定扩散函数。这对应于在像差校正平面中测量的像差定律的空间傅里叶变换。其在每个等晕域内在空间上都是不变的。该扩散函数也可以从聚焦畸变矩阵的奇异值分解获得。如果矩阵G′
图10的图表105比较了由此获得的扩散函数在校正前后的模。在分辨率方面的增益是明显的，并且在校正像差之后达到了衍射极限。
在进行的实验中，测得的像差不依赖于频率，因为在所研究的频带内，声速在体模和水中的变化很小。实际上，情况并非总是如此，并且畸变矩阵的频率分析可能被证明是明智的，以便确定能最佳地校正像差的空频自适应滤波器(或空时等效)。具体而言，可以在以不同中心频率ω
该像差频率校正在波束成形中等效于由空时滤波器发射和测量的聚焦信号的卷积U
到目前为止，已经利用了平移不变畸变层的情况，因此仅生成一个等晕域。现在描述畸变矩阵的第二示例性用途，其出于非平移不变的体积性畸变器的情况下，从而产生多个等晕域。
在这种情况下，倾向于采用迭代方法(后处理)来校正像差。该过程的步骤#0包括通过畸变矩阵D
其中
步骤#1在于重新计算从
接下来的步骤包括通过交替校正在输出处(偶数迭代)和输入处(奇数迭代)的残留像差来再现相同的过程。然而，在每个步骤中，仍然使用在瞳孔平面中的畸变矩阵的相关矩阵的相位(exp[jargB
图11通过在乳房体模上进行实验来示出根据本说明书的方法的益处，所述乳房体模产生多个等晕域并且在20至25mm之间的深度处包含微钙化。图11示出了常规超声图像111、使用归一化后的相关矩阵exp[jarg{B
对样品中每个点的像差定律的了解，不仅可用于成像，而且可用于将超声波物理性地聚焦在样品的任意点上。例如，在图2的实验设备上，可以经由空时滤波器向换能器阵列施加从畸变矩阵的特征向量推导出的初始聚焦信号的卷积。在存在单个等晕域的情况下，将考虑第一特征向量U
Δt(u)＝-arg{U
在存在多个等晕域的情况下，通过迭代过程确定像差定律(请参见先前段落)。根据向量W，将样品在输入处获得的所有像差校正汇总在一起：
在镜面反射模式下，并且在超声图像包含在单个等晕域中的情况下，可以示出畸变矩阵的第一输出特征向量U
U
其中A(θ
V
通过在输入或输出处仅施加一次特征向量U
R′
在超声图像包含P>1个等晕域的情况下，畸变矩阵的秩为P。输入特征向量V
然而，反射器永远不会是完美的镜面反射，因此其随机散射分量将无法通过式(29)进行最佳校正。此外，镜面反射器可能在其界面之间引起多重反射，这可能会使超声图像模糊。因此，必须能够分离由介质反射的场的镜面反射分量和随机散射分量。根据一个或多个示例性实施例，表征方法因此还包括识别和/或消除反射场的镜面分量和/或在换能器阵列与非均匀介质的各个界面或层之间产生的多重反射。为此，可以在输入处和输出处将畸变矩阵投影到像差校正基上。可以通过以下基变换，基于在观察平面中表示的畸变矩阵D
D
因此，前式就矩阵系数而言重写如下：
如图13C所示，对于精确的反射角和入射角对θ
θ
其中Δθ对应于镜面反射器和换能器阵列之间的角度。如果镜面反射器平行于换能器阵列，则Δθ为零。因此，可以容易地滤除镜面反射分量和多重反射分量，并且仅保留反射场的随机散射分量(散斑)，以便精确确定由有机玻璃板引起的像差并进行校正。然后，对随机散射分量的这种区分使得可以直接获得像差定律，该像差定律要被用在输入处和输出处，以校正所述反射矩阵并获得所述介质的最佳图像[式(17)]，这在镜面分量占主导时是不可能的[式(29)]。另外，在傅立叶平面中过滤所述畸变矩阵使得可以消除多重反射(探头与人体的表面层之间的多次回波)，该多重反射污染例如在浅深度处的超声图像。
图12A至图12C示出了使用畸变矩阵对多重反射进行过滤以及该畸变矩阵后续在像差校正和声速映射中的应用。图12A示出了实验装置的图，即有机玻璃板121放置在超声探头41和体模20之间。相应的超声图像122是使用均匀速度模型(c＝1800m/s)构造的。该图像由于有机玻璃板引起的多重反射和像差而大大降级。图12B示出了针对超声图像的深度(z＝25mm)的畸变矩阵的各种形式。图像124示出了其中扩散函数因有机玻璃板的存在而特别降级的聚焦畸变矩阵的模。图像125示出了傅立叶平面中的畸变矩阵。对于θ
根据一个或多个示例性实施例，所述非均匀介质的物理参数的至少一个映射的确定包括确定在所述声穿透区中的声速。这涉及例如研究所述畸变矩阵的第一奇异值作为速度模型的函数的变化，该速度模型用于对参考矩阵进行建模。最佳速度模型是使畸变矩阵的第一奇异值最大化的模型。一种可选方案是研究聚焦畸变矩阵。具体而言，这使得可以研究扩散函数作为用于计算参考矩阵的速度模型的函数的变化。最佳速度模型是使波长归一化后的扩散函数宽度最小化的模型。
因此，图13示出了畸变矩阵用于声速的定量测量的示例性用途。图13示出了已经在图像81(图8)中示出的体模在没有像差的情况下的超声图像131。图像132示出了聚焦畸变矩阵针对图像中的线以及均匀速度模型(c＝1540m/s)的示例，图像133示出了在输入处的焦平面和傅立叶平面之间的对应的畸变矩阵。图像134示出了由波长λ归一化的扩散函数的宽度δ随波速c变化的函数，图像135示出了畸变矩阵的归一化后的第一奇异值(
图像134中的比率δ/λ清楚地显示为至少大约c＝1539m/s，而体模制作商给出的速度为1540m/s。适当的速度模型意味着超声波在发射和接收时的最佳聚焦，以及相应地有更好的分辨率的扩散函数。因此，聚焦畸变矩阵为定量测量体模中的声速提供了相关标准。类似地，图表135中的畸变矩阵的第一奇异值
除了介质的反射率或样品每个点的像差定律知识外，畸变矩阵还可以对所研究样品的光学指标进行3D层析成像。从“焦平面”畸变矩阵D
图14示出了畸变矩阵用于对媒介中的声速建立映射的示例性用途。图14A示出了在均匀速度模型(c＝1540m/s)且穿过有机玻璃板的情况下，体模的常规超声图像141。由于有机玻璃板引起的像差和多重反射，该超声图像质量较差。相比之下，在超声图像142中成功地示出体模的反射率。后者的构建首先是通过过滤图像中的多重反射(图12)，然后通过针对图像中的每一行优化速度模型(图13)。这些操作使得可以获得仅受衍射限制的最佳分辨率。由于镜面反射分量的过滤，仅目标的低空间频率很难重新组合。图143示出了针对在深度上积分的声速而获得的曲线。该速度在进入有机玻璃层时(z＝14mm)急剧增加，然后随着离开有机玻璃层(z＝17mm)而逐渐降低。然后，一种算法使得可以解决反向问题，该反向问题包括从该积分的速度分布图推断出声速随深度变化的精确分布图。在有机玻璃中测得的声速为2700m/s，如所预计的，在体模中测得的声速为1540m/s。
根据一个或多个示例性实施例，所述非均匀介质的物理参数的至少一个映射的确定包括确定在所述声穿透区中的多重散射率。例如，这涉及研究像差已经被预先校正的聚焦畸变矩阵。超出焦点后获得的非相干强度(多重散射)的平均值与焦点处的强度(所述畸变矩阵的中心行)的比值可以给出局部多重散射率。
图15示出了畸变矩阵对于多重散射的局部量化的优势。图15显示了标准体模超声图像151。图152解释了在畸变矩阵中表现为非相干场的多个散射路径的本质。这些发生在焦平面的上游，因为它们与在焦平面中单一散射的波具有相同的飞行次数。聚焦畸变矩阵在图像153中示出，并且对应于超声图像中的线(z＝75mm)。其中心行(r
聚焦畸变矩阵D
其中符号<…>表示除中心元素以外(|r
图15的图像155示出了在超声图像151的每个点处测得的多重散射率。在较浅的深度，该程度不可忽略，但与电子噪声(信号的非互易性)相关。另一方面，该程度在很深处的增加显然与多重散射现象有关。特别而言，在超过60mm的深度，多重散射率超过50％。
该测量结果已通过先前的体内研究得到证实[A.Aubry et al.,J.Acous.Soc.Am.129,225-233,2011]。本发明的贡献在于聚焦畸变矩阵使得可以局部测量多重散射率，而在此之前，只能获得在每个回波时间或深度处的平均值。[A.Aubry et al.,J.Acous.Soc.Am.129,225-233,2011]。
通过花费比弹道时间[式(4)]更长的时间，扩散函数还可以跟着介质中扩散晕的增长，并从中得到多重散射波的传输参数的局部测量(散射系数或传输平均自由程)。研究短时间内的扩散晕可以获得比现有技术在随机散射层析成像中所获得的更好的空间分辨率[A.Aubry et al.,Appl.Phys.Lett.92,124101,2008]。
尽管通过多个详细的示例性实施例进行了描述，但是用于非侵入性超声表征的方法和系统包括不同的变型、修改和改进，这对于本领域技术人员而言将是显而易见的，应理解，这些不同的变型、修改和改进都落入由所附权利要求书限定的本发明范围内。
机译： 用于使用超声的非侵入性地表征异质介质的方法和系统
机译： 用于使用超声的非侵入性地表征异质介质的方法和系统
机译： 异质介质的非侵入性超声表征方法和系统
