申请/专利权人：法国电力公司;国家科学研究中心

申请日：2016-05-03

公开（公告）日：2021-01-12

公开（公告）号：CN107710014B

主分类号：G01S15/89(20060101)

分类号：G01S15/89(20060101);G01S7/52(20060101)

优先权：["20150512 FR 1554258"]

专利状态码：有效-授权

法律状态：2021.01.12#授权;2018.03.16#实质审查的生效;2018.02.16#公开

摘要：本发明涉及通过波传播进行探测的方法，其中，向介质1中发射入射波，所述介质包括能够反射所述波的散射体2。捕捉表示由介质对入射波反射的反射波的信号，所述捕捉信号是简单散分量和多重散分量之和。捕捉信号通过将多重散射分量与简单散射分量相分离来处理。

主权项：1.一种利用传感器组件3通过波传播进行探测的方法，所述方法包括：a发射步骤，在该步骤中，传感器组件3向介质1中发射入射波，所述介质散射所述波，b测量步骤，在该步骤中，所述传感器组件3捕捉表示反射波的信号，所述反射波是由介质1对入射波的反射，所述捕获的信号包括：-简单散射分量，表示是由介质的各个散射体单次反射入射波所形成的反射波的波路径，-而且，在适当情况下，多重散射分量，表示是由介质的散射体上多次连续反射入射波在到达传感器组件之前所形成的反射波的波路径，c处理步骤，在该步骤中，处理所述捕获的信号，以便确定介质的特征，通过滤波至少一个频率变换矩阵，在多重散射分量与简单散射分量之中提取所选择的至少一个分量，所述频率变换矩阵表示在传感器组件的传感器之间的响应，并且至少包括下列子步骤：c1确定窗口变换矩阵的子步骤，在该步骤中，确定与元件之间的响应的窗口时间矩阵KT，t＝[kijT，t]相对应的至少一个窗口频率变换矩阵KT，f，所述元件之间响应的窗口时间矩阵在接近时间T和持续时间Δt的时窗内与在传感器组件的传感器之间的时间响应hijt相对应，其中，f是频率，c2数据旋转子步骤，在该步骤中，通过按照第一方向旋转以及从窗口频率变换矩阵提取分量，通过窗口频率变换矩阵KT，f来计算两个矩阵AlT，f和A2T，f，c3滤波子步骤，在该步骤中，将各个矩阵Al、A2中的多重散射分量与简单散射分量相分离，由此得到至少两个滤波矩阵A1F、A2F，所述两个滤波矩阵分别与矩阵Al、A2相对应并且分别表示简单散射分量或者多重散射分量，c4反向数据旋转子步骤，在该步骤中，通过按照与第一方向相反的第二方向旋转以及从滤波矩阵A1F、A2F中提取分量，通过两个滤波矩阵A1F、A2F来计算滤波窗口变换矩阵KFT，f，所述方法的特征在于：-在数据旋转子步骤c2中，窗口频率变换矩阵KT，f的所有分量分布在两个矩阵AlT，f和A2T，f中的一个或另一个矩阵的中心部位，这些矩阵所具有的零点在中心部位之外的部分，以及，-在反向数据旋转子步骤c4中，提取滤波矩阵A1F、A2F的中心部位的所有非零分量并将其放置在滤波窗口变换矩阵KFT，f之中。

全文数据：利用波传播进行探测的方法和设备技术领域[0001]本发明涉及利用波传播进行探测的方法和设备。背景技术[0002]利用波传播进行探测的方法采用传感器组件并且允许例如能够测量介质的特征参数以及或者检测介质的奇异点以及或者创建介质的图像。[0003]这类方法尤其适用于检测和成像系统，例如，声纳、雷达、超声，等等。[0004]在已知的这类方法中，尤其是在超声或雷达成像方法中，采用所捕获信号的简单散射分量:如果各个散射体仅只与波相互作用一次，那么每个回波的到达时间与传感器和产生该回波的散射体的分隔距离之间存在着有效的直接等效。回波在指定时刻的检测就意味着散射体存在于在与回波的到达时间相对应的距离上。然后，如有需要，就可以通过所接收到的信号来构建介质的反射性图像，也就是在介质内的各个散射体的位置。在超声或雷达成像方法中，并不太多地采用多重散射。这些成像方法往往是基于可以忽略多重散射的假设。[0005]在存在着大量多重散射分量的情况下，尤其是在介质中所包含的散射体具有较高散射功率以及或者在介质中非常密集的情况下，普通成像方法会受到严重的干扰并且不再是可靠的。实际上，在这种情况下，在回波的到达时间与传感器和介质散射体的分隔距离的之间不再有等效性性，这就难以构建出介质的图像。[0006]专利申请书W0-2010001027提出了一种能够将多重散射分量与简单散射分量相分离的探测方法。具体而言，该文件的方法包括以下步骤：[0007]a发射步骤，在该步骤中，传感器组件意即组件的某些或全部传感器）向介质中发射入射波，所述介质散射所述波，[0008]⑹测量步骤，在该步骤中，所述传感器组件（意即组件的某些或全部传感器捕捉表示反射波的信号，所述反射波是由介质对入射波的反射，所述捕获的信号包括：[0009]-简单散射分量，表示是由介质的各个散射体单次反射入射波所形成的反射波的波路径，[0010]-而且，在适当情况下，多重散射分量，表示是由介质的散射体上多次连续反射入射波在到达传感器组件之前所形成的反射波的波路径，[0011]c处理步骤，在该步骤中，处理所捕获的信号，以便确定介质的特征所讨论的特征可包括介质的图像，和或介质的参数数值，和或存在或者不存在比如多相性这样的奇异点，等等），通过滤波至少一个频率变换矩阵，在多重散射分量与简单散射分量之中提取所选择的至少一个分量，所述频率变换矩阵表示在传感器组件的传感器之间的响应，并且至少包括下列子步骤：[0012]cl确定窗口变换矩阵的子步骤，在该步骤中，确定与元件之间的响应或者与在组件的传感器之间的响应的窗口时间矩阵KT，t=[kijT，t]相对应的至少一个窗口频率变换矩阵KT，f，所述元件之间响应的窗口时间矩阵在接近时间T和持续时间At的时窗内与在传感器组件的传感器之间的时间响应hut相对应，其中，f是频率，[0013]c2数据旋转子步骤，在该步骤中，通过按照第一方向旋转以及从窗口频率变换矩阵提取分量，通过窗口频率变换矩阵KT，f来计算两个矩阵AlT，f和A2t，f，[0014]c3滤波子步骤，在该步骤中，将各个矩阵A1、A2中的多重散射分量与简单散射分量相分离，由此得到至少两个滤波矩阵A1F、A2F，所述两个滤波矩阵分别与矩阵A1、A2相对应并且分别表示简单散射分量或者多重散射分量，[0015]c4反向数据旋转子步骤，在该步骤中，通过按照与第一方向相反的第二方向旋转以及从滤波矩阵A1F、A2F中提取分量，通过两个滤波矩阵A1F、A2F来计算滤波窗口变换矩阵!^!^。发明内容[0016]本发明的目的旨在精炼和改善上述探测方法。[0017]为此目的，根据本发明的方法的处理步骤使得：[0018]-在数据旋转步骤c2中，窗口频率变换矩阵KT，f的所有分量分布在两个矩阵AlT，f和A2t，f中的一个或另一个矩阵的中心部位，这些矩阵所具有的零点在中心部位之外的部分，以及，[0019]-在反向数据旋转步骤c4中，提取滤波矩阵A1F、A2F的中心部位的所有非零分量并将其放置在滤波窗口变换矩阵KfT，f之中。[0020]通过根据本发明的方法的这些特征，可以在处理步骤中考虑来自介质的更大量信息。因此，利用波传播来探测介质也就更为有效，并且使之能够探测到介质的更大深度。[0021]此外，由此所产生的介质图像的空间分辨率可以更高。因此，有可能使用具有较少数量传感器的传感器组件，使得该方法的成本较低。[0022]此外，例如，分量的分离使之能够：[0023]-仅仅只利用简单的散射分量来构建介质的图像，或者，[0024]-通过多重散射分量来确定介质的特征参数，或者，[0025]-确定很少有或者没有多重散射，以此表述介质的特征。[0026]在根据本发明的方法的各个实施例中，有可能采用下列设置中的一项或多项。[0027]根据方法的一方面，窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且在数据旋转步骤c2中，两个矩阵AlT，f=[aluvT，f]和A2=[a2uvT，f]可以通过窗口频率变换矩阵KT，f=[1¾T，f]按照下列关系来计算：[0028]-如果u+v-1-N-I2和v-u+N+32-1都在1与N之间且包含1和N，[0029]那么[0030]否则aluvT，f=0,并且，[0031]-如果u+v-N-I2和v-u+N+32-1都在1与N-I之间且包含1和N-I，[0032]那么[0033]否则a2uvT，f=0。[0034]根据方法的一方面，数据旋转子步骤c2包括：[0035]-第一个子步骤c21，在该子步骤中，将窗口频率变换矩阵KT，f引入尺寸更大的广义矩阵KeT，f中，除了来自窗口频率变换矩阵的元素之外，广义矩阵的所有元素都为零，以及，[0036]-第二个子步骤c22，在该子步骤中，两个矩阵AlT，f和A2T，f可以通过广义矩阵KeT，f来计算。[0037]根据方法的一方面，窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且：[0038]-在第一个子步骤c21中，广义矩阵KeT，f二化叫⑴^定义如下：[0039]如果ie[1+N-D2;N+N-I2]，以及，[0041]那么[0042]否则keijT，f=〇,以及，[0043]-在第二个子步骤C22中，两个矩阵AlT，f=[aluvT，f]，A2T，f=[a2uvT，f]可通过下列关系来计算：[0044]，以及，[0046]根据方法的一方面，窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且在反向数据旋转步骤c4中，滤波窗口变换矩阵可通过滤波矩阵按照下列关系来计算：[0047]如果i-j是偶数，[0048]那么[0049]否则[0050]根据方法的一方面，在确定窗口变换矩阵的子步骤cl中，通过窗口时间矩阵KT，t的小波变换来确定各个窗口频率变换矩阵KT，f，所述窗口时间矩阵KT，t在接近时间T和持续时间Δt的所述时窗内与在传感器组件的传感器之间的时间响应t相对应。[0051]根据方法的一方面，各个窗口变换矩阵KT，f的时窗至少成对地重叠。[0052]根据方法的一方面，在滤波子步骤c3中：[0053]-从矩阵KT，f的各个系数kijT，f中提取简单散射分量ksijT，f，以便简单散射分量Sk\T，f在所述矩阵KT，f的各个反对角相互相关，而且，[0054]-从矩阵KT，f的各个系数kijT，f中提取多重散射分量kMijT，f。[0055]根据方法的一方面，在发射步骤中，发射进入介质中的入射波以预定的角度朝所述介质的外表面倾斜。[0056]本发明还涉及实施上述探测方法的设备，包括传感器组件，适合向散射介质中发射入射波并且捕获表示反射波的信号，所述反射波是由介质反射入射波所形成的，所述捕获的信号包括：[0057]-简单散射分量，表示反射波是由介质中的各个散射体单次反射入射波所形成的波路径，[0058]-以及，在适当情况下，多重散射分量，表示反射波是由介质的散射体上多次连续反射入射波在到达传感器组件之前所形成反射波的波路径，[0059]所述设备进一步包括处理装置，适合处理所述捕获的信号，以便确定介质的特征，[0060]处理装置适合通过至少一个频率变换矩阵的滤波，从多重散射分量和简单散射分量中提取所选择的至少一个分量，所述频率变换矩阵表示在传感器组件的传感器之间的响应，并且处理装置至少执行以下子步骤：[0061]cl确定窗口变换矩阵的子步骤，在该步骤中，确定与元件之间响应的窗口时间矩阵KT，t=[kijT，t]相对应的至少一个窗口频率变换矩阵KT，f，所述元件之间响应的窗口时间矩阵在接近时间T和持续时间Δt的时窗内与在传感器组件的传感器之间的时间响应hij⑴相对应，其中，f是频率，[0062]c2数据旋转子步骤，在该步骤中，通过按照第一方向旋转并且从窗口频率变换矩阵中提取分量，通过窗口频率变换矩阵KT，f来计算两个矩阵AlT，f和A2t，f，[0063]c3滤波子步骤，在该步骤中，将各个矩阵A1、A2中的多重散射分量与简单散射分量相分离，由此得到至少两个滤波矩阵A1F、A2F，所述两个滤波矩阵分别与矩阵A1、A2相对应并且分别表示简单散射分量或者多重散射分量，[0064]c4反向数据旋转子步骤，在该步骤中，通过按照与第一方向相反的第二方向旋转并且从滤波矩阵A1F、A2F中提取分量，通过两个滤波矩阵A1F、A2F来计算滤波窗口变换矩阵!^!^。[0065]本发明的设备使得处理装置适合于：[0066]-在数据旋转步骤c2中，窗口频率变换矩阵KT，f的所有分量都分布在两个矩阵AlT，f和A2t，f中的一个或另一个矩阵的中心部位，这些矩阵所具有的零点在中心部位之外的部分，以及，[0067]-在反向数据旋转步骤c4中，提取滤波矩阵A1F、A2F的中心部位的所有非零分量并将其放置在滤波窗口变换矩阵KfT，f之中。附图说明[0068]通过以非限制性实例的方式列出的本发明实施例的以下说明，参考附图，本发明的其它特征和优点显而易见。[0069]在附图中：[0070]-图1是实施根据本发明实施例的方法的设备的示意图，[0071]-图2是展示了两个简单散射路径的图表，在图1所示的实例中，所述简单散射路径有助于形成探测介质在时间T=2Rc反向散射的场，[0072]-图3是展示了两个多重散射路径的图表，在图1所示的实例中，所述多重散射路径有助于形成探测介质在时间T=2Rc反向散射的场，[0073]-图4和图5阐释了数据旋转子步骤，所述数据旋转子步骤是根据现有技术方法的实施例的一部分，其中，从各个窗口频率变换矩阵K中提取两个矩阵Al和A2，[0074]-图6阐释了反向数据旋转子步骤，在该子步骤中，对所述矩阵滤波之后，通过上述矩阵Al和A2来确定滤波变换矩阵Kf，[0075]-图7阐释了滤波阵列Al和A2通过奇异值分解的迭代过程，[0076]-图8阐释了根据本发明第一个实施例的示例性广义矩阵Ke，[0077]-图9和图10显示了通过图8的第一个实施例的广义矩阵所获得的或者根据本发明第二个实施例从窗口频率变换矩阵K所直接获得的两个矩阵Al和A2，[0078]-图11a、图Ilb和图lie阐释了在本发明的方法中可用于形成窗口频率变换矩阵K的母小波函数的三个实例，以及[0079]-图12是设备的变体的示意图，具有相对于介质外表面倾斜并且适合实施根据本发明方法的实施例的多个传感器。具体实施方式[0080]图1显示了适合通过发射和接收波来探测介质1的设备的一个实例。在下文中，本发明针对所讨论的波为超声波例如，频率为2MHZ至4MHz的情况作更具体的说明，当然，本发明也适合任何类型的波，例如，除超声压缩波以外的机械波、电磁波或其它波。[0081]介质1是所讨论的波的散射介质，意即所述介质是均质的以及或者包含随机分布的散射体2，并且所述散射体能够反射发射入介质1中的波。[0082]所讨论的介质1可以是，例如人体的一部分，而散射体2可以是包含在介质1中的未溶解的小颗粒这种散射体会在超声中产生“散斑”图像）。当然，待探测的介质1可以是其它介质，例如，工业项目中的一部分，有待于在非破坏性试验的环境下对其结构进行测试。[0083]图1所示的探测设备包括传感器4的组件3,例如，超声传感器的线性阵列，其典型的形式通常为与介质1相接触的刚性阵列。该组件3包括N个传感器，例如，N介于100至500之间。例如，可以采用一百个传感器4的阵列，每个传感器的宽度约为0.39毫米。值得注意的是，传感器4在此为超声传感器，它能够将电气信号转换为超声压缩波，反之亦然，但是就本专利的意义而言，传感器更普遍地是指能够将电气信号转换为某个类型的波机械波、电磁波等并且反之亦然的任何设备。[0084]组件3的各个传感器4可以由中央处理器5UC单独控制，例如，所述中央处理器包括数字信号处理装置，该中央处理器5能够例如在屏幕6上显示介质1的图像。[0085]为了探测介质1，中央处理器5向传感器4发送电气信号，通过所述传感器将所述电气信号转换为发射入介质2中的波，在这种情况下，为超声压缩波并且这些波由介质中包含的散射体2局部反射。一部分散射波或回波）因此返回传感器4,所述传感器捕获并将其转换为电气接收信号，随后由中央处理器5对其进行处理。[0086]这些波返回传感器4:[0087]-或者，经散射体2反射一次之后（例如，沿着图1的虚线按照图示所示的路径），换言之，在简单散射之后，[0088]-或者，经过多个散射体2连续反射之后（例如，沿着图1的实线按照图示所示的路径），换言之，在多重散射之后。[0089]经介质1散射并返回传感器4的整体波因此包括两个作用：[0090]简单散射”作用，与入射波在回到传感器4之前仅经过介质1各个散射体1的单次反射的情况相对应，[0091]多重散射”作用，与入射波在到达传感器4之前经过在介质散射体2上的多次连续反射的情况相对应。[0092]这种方法能够通过滤波将这两个作用相分离，所以只采用其中一个，或者分别对其进行处理。例如：[0093]-简单散射作用有助于构建介质的图像，尤其是在此所述的情况下通过超声来构建介质的图像但是，在地震波的情况下，可以是震波的图像，或者在电磁波的情况下，可以是雷达的图像，或者为其它的图像）：事实上，因为各个散射体2仅仅只与入射波相互作用一次，所以每个回波的到达时间与在传感器和产生该回波的散射体之间的距离之间有效地直接等效。以指定时间来检测回波信号可以表示在与回波的到达时间相对应的距离存在着散射体2。然后，可以通过记录的信号来构建介质反射性的图像，意即介质内各个散射中心的位置；[0094]-多重散射作用还有助于表述无序介质的特征，因为它能够测量统计参数，例如平均畅通路径或者通过随机介质传播的波的散射系数;这种多重散射作用的相对重要性使之能够表述例如通过超声所获得的介质1图像的可靠性的特征。[0095]为了把简单散射作用与多重散射作用相分离，首先记录组件3的各对传感器4的元件之间的响应。[0096]为此目的，如图2所示，在发射步骤a中，由组件3的第i个传感器横坐#Xl发送与脉冲信号相对应的入射波，并且在测量步骤b中，由第j个传感器横坐标记录在发射入射波之后由介质背向散射反射的波。在这个操作过程中，各个传感器j捕获在传感器i发射之后的时间信号t。该操作适合于组件3的所有可能存在的发射器接收器对，意即实际上适合于传感器4的所有对进行该操作。[0097]N2个响应的集构建元件之间响应的时间矩阵Ht=[h^t]，即尺寸为NxN的矩形矩阵，这是介质1的整体响应。注意，例如，通过文件W0-A-2004086557所描述的程序，有可能更迅速地获取元件之间响应的时间矩阵Ht，而无需组件3的各个传感器i连续发射脉冲信号。[0098]信号⑴通过中央处理器5进行数字化例如，按照9位或相类似），采样例如，例如针对超声波以20MHz的频率进行采样），以及储存。[0099]中央处理器5然后执行处理步骤c，包括确定窗口变换矩阵的初始子步骤cl。[0100]在该子步骤cl中，首先是将各个脉冲响应hijt缩短窗口）为持续时间Δt的连续时窗。[0101]由此得到尺寸为NxN的一系列窗口时间矩阵KT，t=[1^1'4]，其中，1^1'4是与时窗[T-Δt2窗口时;T+Δt2窗口]相对应的hij⑴的贡献，为：[0103]使得：[0104]-若，则Wr⑴=1，否则Wr⑴=0，[0105]-R=c.T2是沿着z轴从传感器组件3测量的距离（见图2，所述z轴垂直于所述传感器组件的纵向X在此所考虑的实例中构成线性阵列），[0106]-C是波在介质1中的速度在超声的具体情况下，为在水中或在人体的软组织中，c=1·5mmys，[0107]-Δt是预定的持续时间，对应于在上述时间间隔[T-Δt2时;T+Δt2]内波所行进的距离AR=C.At，有利的是，这是连续的，意即其完全包括所确定的时间段换言之，所进行的处理包括视野深度，意即也是连续的R的一系列值）。[0108]在子步骤cl中，进行矩阵KT，t系数的离散傅里叶变换，从而获得尺寸为NxN的变换矩阵的T的各个数值，所述变换矩阵也称之为窗口频率变换矩阵KT，f=[kijT，f]，其中，kijT，f是kijT，t的离散傅里叶变换以及f是其频率。[0109]根据这些窗口频率变换矩阵KT，f，在随后的滤波子步骤c3中，可以通过滤波将简单散射作用和多重散射的作用相分离，所述随后的滤波子步骤c3是处理步骤c的一部分。[0110]具体的说，在该随后的滤波子步骤C3中，可以基于所述窗口频率变换矩阵KT，f的各个反对角上的窗口频域变换矩阵KT，f的系数KijT，f的相干性将各个窗口频率变换矩阵KT，f中的多重散射分量和简单散射分量相分离反对角常称之为所述矩阵的系数KijT，f队列，以便i+j是常量窗口。[0111]实际上，简单散射波沿着矩阵KT，f的反对角具有明显的相干性，而多重散射波具有随机性且在所述矩阵KT，f不呈现出优选的相干方向。通过仔细滤波这些反对角，可以把两个作用相分离。[0112]该特性可作如下解释。[0113]各个脉冲响应t可分解为以下形式：[0115]其中，和.分别与来自简单散射⑸和多重散射⑽的信号相对应。[0116]同样，窗口频率变换矩阵KT，f的系数T，f可分别分解为的形式，其中，是简单散射的作用，是多重散射的作用。[0117]可将各个作用视为与多个简单散射路径相关的分波的总和，其两个实例路径Cl1和路径d2如图2所示。在时窗[T-Δt2;T+At2_]内，可以到达传感器阵列的简单散射波与初始波对在位于厚度为AR=c.Δt、深度为R=cT2的介质区域内的散射体的反射相对应。[0118]我们用术语“等体积”表示一组点，这组点在指定时刻T有助于在阵列中所捕获的信号。实际上，等体积并不全是与阵列表面平行的部分，而是由椭圆叠加产生的，所述椭圆的焦点是发射器元素⑴和接收器元素（j。在远场中，意即在R足够大的情况下，等体积相似于平行于阵列且与其距离为R=cT2的厚度为ΔR的部分。[0119]在元素i和j之间的响应可以分解为通过在等体积的Nd个散射体上反射产生的分波的总和。采用两维方式，在近轴近似的情况下，我们可以把I写为：[0121]其中，整数d表示有助于在时间T所接收到的信号的第dth个简单散射路径，Xd是散射体d沿着轴X的横向位置，k是周围介质中的波数k=2ja，A是波长），以及Ad是表述散射体d反射性特征的振幅。[0122]值得注意的是，在上述等式⑶中以及在当前专利申请书的其它等式中，j是虚数，所以在不作为下标的情况下，j2=-l，但是在作为下标的情况下，则表示矩阵元素的位置。[0123]多重散射的作用还可以分解为与多重散射路径相对应的分波，所述多重散射路径的长度在区间[R-AR2;R+ΔR2]内，如图3针对指数P和q的两个路径所示。[0124]我们可以把写为：[0126]其中，整数p表示所涉及的多重散射路径的指数。在图3所示的实例中，对子和分别表示路径P的第一个和最后一个散射体的坐标。仏相当于与路径P相关的分波的复振幅，从第一个散射事件到最后一个散射事件[0127]尽管散射体2在介质1中的分布是完全随机的而且在散射体之间没有关联，但是与简单散射相关的信号却有特殊的相干性，这与多重散射的作用不同。实际上，等式⑶可以写为：[0129]在等式5之和的前面所出现的项与散射体的精确分布无关，并因此起着决定性的作用，表述了简单散射的特征。右边的项是随机的，因为它明确地取决于散射体的位置。[0130]相比之下，与多重散射相关的信号等式4就不能这样来写成。[0131]源自简单散射的信号的该特性可以表达为沿着各个矩阵KT，f的反对角的特殊相干性，如图4所示。实际上，沿着各个反对角，意即针对发射器i接收器j对，其中，i+j=常数，的随机作用是常数。对于指定的散射介质而言，在简单散射中，因此位于相同反对角上的各个矩阵KT，f的系数之间存在确定性的相位关系。[0132]然而，在多重散射中，则不再满足这个特性，而且矩阵KT，f不具有特殊的相干性:元素不受彼此影响。[0133]本方法利用了这个特性，以便在用实验方法测量的信号的基础上，利用各个矩阵KT，f内的简单散射作用的特殊对称性，通过滤波就能将简单散射作用与多重散射作用相分离。因此，可以过滤掉：[0134]-矩阵KT，f的各个系数KljT，f的简单散射分量，使得简单散射分量Sks1」T，f在所述矩阵KT，f的各个反对角相互相干，[0135]-以及矩阵KT，f的各个系数KijT，f的多重散射分量，所以分量是随机的且彼此不关联。[0136]适用于分离两种作用的两种示例性滤波技术如下文所述。[0137]按照这两种技术，在确定窗口变换矩阵的子步骤cl之后，处理步骤C包括以下两个子步骤：[0138]-数据旋转子步骤c2，通过旋转各个矩阵KT，f的数据并构建两个子阵^T，f和A2T，f，[0139]-上文提及的滤波子步骤c3，在该子步骤中，对矩阵^和如进行滤波，尤其是通过投影技术1或者通过奇异值分解技术2，由此可以得到两个滤波矩阵，表示为和[0140]-反向数据旋转子步骤c4，由此使之能够通过和重新构建包含简单散射的作用或者多重散射的作用的滤波变换矩，这取决于所应用的滤波的类型。[0141]下面详细说明了这些子步骤c2至c4。[0142]子步骤c2:数据旋转[0143]在该子步骤c2中，中央处理器5通过各个矩阵KT，f来计算两个矩阵AlT，f=[aluvT，f]和A2=[a2uvT，f]，[0144]其中：[0148]选择N，使得M为整数:例如，N=125和M=32。如果传感器3的总数使得M不是整数，我们通过减少的数量N个传感器进行计算，使得M=N+34为整数在此所考虑的具体实例中，例如，我们可以利用128个传感器的阵列且仅仅只采用N=125进行计算）。[0149]矩阵Al和A2是由按照逆时针方向旋转45°的矩阵KT，f的子集所构建的矩形矩阵。图4和图5分别阐释了这些矩阵Al和A2，可将其理解如下：[0150]-Al是尺寸为2M-1X2M-1的矩形矩阵，其中，各行由矩阵KT，f中的某些对角线海隔一个所构成，包括在图4中圈出的矩阵KT，f的系数KT，f的系数在图4和图5中用点表示），[0151]-A2是尺寸为2M-2X2M-2的矩形矩阵，其中，各行由矩阵KT，f中的其它对角线所构成，包括在图5中圈出的矩阵KT，f的系数。[0152]在下文中，我们可以用Ar=[ar^]来表示矩阵Al和A2r=l或2，而且我们用L表示矩阵Ar的尺寸关于矩阵Al，我们因此得到L=2M-1，关于矩阵A2，得到L=2M-2。[0153]由于空间的相互作用，矩阵K相对于其主对角线DGuj=I^1对称。矩阵Ar因此也具有对称性:其上部每一行与其下部的一行相同，相对于与主对角线D相对应的水平中心线对称。因此，通过下部就可以直接推断出矩阵Ar的上部。于是，矩阵Al的各列仅仅只包括M个独立系数，只要其尺寸LM，以及矩阵A2的各列包括M-I个独立系数。更普遍而言，矩阵Ar的独立系数的数量因此为数量Mr，使得如果r=1，则Mr=M;如果r=2，则Mr=M-I。[0154]子步骤c3:滤波[0155]在滤波子步骤c3中，中央处理器5将各个矩阵Ar中的多重散射分量与简单散射分量相分离，r是等于1或2的指数，因此我们得到至少两个滤波矩阵A，分别与两个矩阵Ar相对应并且分别表示简单散射分量或多重散射分量。[0156]尤其是，该滤波可以根据上文所提及的技术1或者根据技术2进行。[0157]1.技术1:投影滤波[0158]按照这第一种技术，在滤波子步骤c3中，中央处理器5计算分别表示简单散射的两个滤波矩阵A。[0159]通过下列公式来计算各个所述滤波矩阵：[0161]其中：[0162]-S=[su]是列向量，tS*是向量S的共辄向量的转置阵列，[0163]-向量S的分量Su是复数，等于：[0164]-S是常数，实际上等于1因此下文未提及），[0165],其中，如果r=l，则U=i-j2+M，如果r=2,则U=i-j-D2+M,[0166]-XdPxj是沿着轴X的指数i和j的传感器的横坐标，传感器组件至少沿着所述轴X延伸，[0167]-若r=1则L=2M-1，若r=2则L=2M-2。[0168]该公式调整如下。[0169]各个矩阵Ar为两项之和，即Ars和A，分别表示简单散射和多重散射的作用：[0171]数据旋转，换言之，从KT，f变换到Ar，在数学上表达为坐标变换Xi，Xj—yu，yv：[0173]然后，根据该新基础把等式⑸写为：[0175]其中[0176]因此，对于指定介质1而言，矩阵Ars的各列都完全确定地取决于各列的指数U。[0177]然而，不能对多重散射的作用等式4进行如此简单的因式分解。即使在旋转矩阵之后，散射体的位置的随机性仍存在于矩阵A的各列和各行。[0178]因此，通过将整个矩阵Ar的各列投影到由坐标向量S所产生的“简单散射的特征空间”中，就可以对简单散射信号进行滤波：[0180]分母中存在确保向量S的归一化。通过如此投影所产生的行向量P可写为：[0181]P=tS^Ar9,[0182]其中，V是向量S的共辄转置阵列。[0183]矢量P的坐标可由下式给出：[0185]剩余项与向量S上的多重散射信号的投影相对应。[0186]然后，用列向量S乘以行向量P得到滤波矩阵A:[0187]Arp=SP=S-V-A11[0188]于是，把矩阵A的坐标写为：[0190]第一项直接等于简单散射分量等式7。因此，我们得到：[0192]就矩阵而言，等式（13写为：[0194]矩阵ArF实际上并不包含与简单散射As相关的作用，还包含与多重散射StSlw的存在相关的剩余项。该项的存留是因为多重散射信号并非直接与通过向量S所产生的简单散射的特征空间相正交。[0195]所获得的滤波因此并不完美;然而可以评估残留噪声的程度。[0196]实际上，如涉及数据旋转的段落所述，矩阵Al中的各列都仅仅只有M个独立的系数，矩阵A2中仅仅只有M-I个独立系数;因此，多重散射的作用经滤波之后减少了因子由于简单散射的作用保持不变，所以，信噪比，或更具体而言，“简单散射多重散射”比与相似。[0197]上述滤波技术技术1尤其适用于想要提取被多重散射所遮盖的简单散射作用的情况，意即在与多重散射的信号相比，简单散射信号的介质的振幅非常低的情况。这特别适合于检测埋藏在散射介质中的目标的情况。[0198]2.技术2:通过奇异值分解的滤波[0199]这第二项技术包括通过对旋转之后所得到的矩阵AjPAdA行奇异值分解SVD将简单散射和多重散射相分离。SVD的特性为将矩阵分解为两个子空间：“信号空间”（矩阵具有在其各行和或各列之间存在着显著相关性的特征）以及“噪声空间”（在元素之间没有相关性的随机矩阵）。通过将SVD应用于在旋转之后所得到的矩阵Ar，信号空间与矩阵Ars相对应简单散射的作用，特征在于沿着其各列的高度相关性），噪声空间与矩阵ArM多重散射的作用有关，其中，Ar=Ars+ArM在关于技术1的段落已经描述了等式6。[0200]矩阵Ar的SVD写为：[0202]其中，U和V是尺寸为L的单一矩形矩阵。其各列仏和1与本征向量相对应，所述本征向量与有关奇异值有关。Λ是尺寸为L的矩形对角矩阵，其对角元素与按照降序排列的奇异值sM,r相对应。在关于数据旋转的段落中，强调了矩阵Ar的特殊对称性:该矩阵仅仅只包括Mr个独立行，因此阶为Mr〈L。因此，矩阵Ar仅仅只具有Mr个非零奇异值，等式（15写为：[0204]因为简单散射的特征在于，在数据旋转之后，沿着矩阵Ar各列具有很高的相干性，所以SVD揭示了信号空间中的这个作用简单散射作用与最高奇异值相关），与此同时，多重散射的作用与最低奇异值相关。在此，与第一种滤波技术不同，在简单散射的情况下，没有关于在矩阵KT，f反对角上的现有相干性的形式的先验假设:只是简单地假设存在这种相干性。[0205]问题是确定奇异值的哪一阶对应于将“信号”空间（与简单散射相关和“噪声”空间（与多重散射相关相分离的阈值。如果等式⑸是严格准确的，则只有第一个奇异值可能与信号空间相对应。如果假设等式5是不严格准确的，则简单散射的作用不是1阶的，而且多个奇异值有承载该作用的痕迹。必须构建关于分离简单散射的作用信号空间）与多重散射的作用噪声空间）的标准。[0206]为此，我们利用通过随机矩阵理论所得出的结果。按照惯例和简单起见，奇异值Xi,r通过其均方根作归一化：[0208]对于大的矩阵，其中的系数是完全随机的而且其间没有相关性，第一个奇异值$从不超过值（在矩形矩阵的情况下，）。[0209]实验上，多重散射作用并非准确地与完全随机的矩阵相对应，因为传感器之间存在剩余相关尤其是由于组件3的相邻传感器之间的机械耦合或者电耦合），由此修改根据[由A.M.Sengupta和P.P.Mitra编著的《关于某些随机矩阵的奇异值分布》详见1999年《物理评论》E版第60⑶卷，第3389-3392页]，可以构建关于这种矩阵的奇异值的新概率论并且用其来确定的数值，由此使之能够定义关于分离信号空间与噪声空间的客观标准。[0210]在旋转实验数据之后，待处理的矩阵Ar见等式6因此是与简单散射相关的p〈M阶矩阵Ars以及与想要滤波的多重散射相关的M阶矩阵A之和。[0211]所提出的技术如下：在实施SVD之后，中央处理器5考虑在归一化之后的第一个奇异值。如果该值大于，则意味着第一本征空间与简单散射相关。[0212]然后，在2阶并且如有必要在更高阶重复该过程。[0213]如图7所示，中央处理器5所执行的程序因此是迭代程序，包括以下子步骤：[0214]c31将q初始化为1，[0215]c32利用Aq,r计算归一化的奇异值：[0217]c33如果.至少等于预定的阈值，则归因于简单散射，并在q增1之后重复子步骤c32，[0218]c34如果小于阈值，则.以及随后的所有奇异值都归因于多重散射。[0219]如果的阶表示为P+1，则由此得到：[0222]于是，矩阵Ars包含简单散射的作用（附加剩余多重散射的作用），矩阵A与多重散射相关。[0223]值得注意的是，技术2假设第一个归一化奇异值超出阈值s它不能用于高的散射介质，意即对所述介质而言，多重散射的作用相对于简单散射而言是主要的。在这种情况下，利用将反对角投影到简单散射空间上的滤波技术来代替技术1，从而提取简单散射的作用。相反，如果简单散射的作用是主要的，或者与多重散射的数量级相同，则可以采用矩阵A滤波的SVD技术技术2并由此提取多重散射的作用。[0224]子步骤c4:反向数据旋转[0225]在反向数据旋转子步骤c4中，中央处理器5执行子步骤cl中所述变换的逆向变换，并由此计算滤波窗口变换矩阵_，其中：[0226]-在i-j是偶数的情况下：[0227]-在i-j是奇数的情况下：[0228]矩阵KFT，f是尺寸为2M-1X2M-1的矩形矩阵，包含源自简单散射或多重散射的信号，根据所应用的滤波类型。图6阐释了反向数据旋转程序，图6显示了矩阵KfT，f包括矩阵AlT，f和A2T，f的反对角的某些系数。[0229]然后，可以通过各种方式利用滤波矩阵KfT，f:[0230]-如果KfT，f与简单散射分量相对应，可将其具体用于检测介质1的奇点或者构建介质1的图像。为此目的，例如，可以采用两种成像方法：[0231]按照第一种成像方法，例如，中央处理器5可以通过反傅立叶变换由各个矩阵KfT，f来计算元件之间响应的滤波窗口时间矩阵HfT，t，然后求不同时间矩阵HfT，t之和，以便得到元件之间响应的滤波时间矩阵Hft，表示在传感器组件的传感器之间的响应。然后，中央处理器通过形成基于所述元件之间响应的滤波时间矩阵Hft的路径来构建介质图像。[0232]按照第二种成像方法，称为分解时间反转算子的方法（“第二种方法（D0RTmethod”，由Prada等人所定义的，详见的爱思唯尔出版商的《波动》201994151-163中的《时间反转算子的本征模式：多目标介质中选择性聚集的方案》），中央处理器5可以执行以下步骤：[0233]dl针对各个滤波变换矩，确定时间反转算子0RTT，f=KF\T，fKFT，f，它与所有或各个变换矩阵KfT，f相对应，KF$T，f是#T，f的复共辄矩阵，[0234]d2确定所述时间反转算子的至少某些本征向量和或本征值，[0235]d3至少根据所述本征向量和或所述本征值来检测介质中的至少一个奇点。在该步骤d3中，可以确定时间反转算子的本征向量Vif，然后在介质1的数字模型中用数字传播这些本征向量，从而确定与各个本征向量相对应的在介质1中反射物体的位置，[0236]-如果矩阵#T，f与多重散射分量相对应，则可将其具体用于通过所述多重散射分量来计算表示介质中多重散射的重要性的指数。在这种情况下，例如，可以通过已知方式或者上述方法来构建介质的图像超声或其它），并且根据表示多重散射重要性的所述指数来量化所述图像的可靠性。有利的是，可以在多个点（尤其是以不同深度R来计算该指数，并且根据表示多重散射的重要性的所述指数来量化与所述多个点相对应的图像的多个部分的可靠性。[0237]本发明的方法改进了上述方法。[0238]尤其是，避免在数据旋转c2和反向数据旋转c4步骤中丢失来自矩阵的信息。[0239]实际上，如图4至图5所示，现有技术的数据旋转步骤c2在构建两个矩阵AlT，f和A2T，f时没有利用窗口频率变换矩阵KT，f的角，所以对于方法的其余计算而言，丢失了窗口频率变换矩阵KT，f的很大一部分数据。该窗口频率变换矩阵的四个角Kl、K2、K3、Κ4没有用于构建矩阵Al见图4或者构建矩阵Α2见图5。因此，该数据旋转步骤大体上排除了窗口频率变换矩阵的一半数据。[0240]根据本发明方法的第一个实施例，数据旋转子步骤c2包括：[0241]-第一个子步骤c21，在该子步骤中，将窗口频率变换矩阵KT，f引入最好在中心尺寸较大的广义矩阵KeT，f中，除了来自窗口频率变换矩阵KT，f的元素外，所述广义矩阵Ke的所有元素都为零，以及，[0242]-第二个子步骤c22与现有技术的数据旋转子步骤相似），在该子步骤中，通过从按照第一方向旋转45°的子阵的广义矩阵KeT，f中的提取，由广义矩阵KeT，f来计算两个矩阵和[0243]调整广义矩阵Ke的尺寸，以便来自窗口频率变换矩阵的元素都包含在倾斜45°的正方形Ker之中，而且其本身也包含在广义矩阵Ke中。尤其是，广义矩阵Ke的尺寸大于或等于2xN-l〇[0244]图8显示了一个实例，在该实例中，窗口频率变换矩阵KT，f的尺寸为9x9N=9，而且将该矩阵引入较大的广义矩阵Ke的中心。于是，广义矩阵的尺寸为17x17。[0245]广义矩阵Ke因此包括角的部分1^1、1^2、163、164，侧面部分165、166、、167、168以及中心部位Kec。侧面部分也是正方形Ker的角。除了中心部位Kec以外的广义矩阵Ke的所有元素包含窗口频率变换矩阵KT，f的元素都为零。换言之，角的部分Kel、Ke2、Ke3、Ke4*侧面部分1^5、1^6、167、168的所有元素都为零。[0246]正方形Ker的作用是通过采用广义矩阵Ke的所有其它元素构成第二个子步骤的矩阵Al和A2的基础，如下文所述。[0247]例如，如果N是奇数，则在第一个子步骤c21中，广义矩阵的尺寸为2xN-l，并利用下列关系式或者任何其它等效公式，根据矩阵KT，f将广义矩阵Ke定义为KeT，f=[ke„T，f]:[0248]如果，以及，[0250]贝Ij[0251]否则keijT，f=0。[0252]如果N是偶数，则可以创建等效的公式。[0253]于是，在第二个步骤c22中，通过广义矩阵KeT，f=[kei」T，f]来计算两个矩阵和[0254]在图8所示的情况下，矩阵Al的尺寸为9x9以及矩阵A2的尺寸为8x8。[0255]实际上，矩阵Al的尺寸等于窗口频率变换矩阵KT，f的尺寸，意即尺寸等于NxN，矩阵A2的尺寸等于窗口频率变换矩阵KT，f的尺寸减1，意即尺寸等于N-IXN-I。[0256]图9显示了从图8所示的广义矩阵Ke中提取的矩阵A1。该矩阵Al采用通过由大的点所表示的元素，使其按照第一方向（逆时针方向）旋转45°并将其放置在矩阵Al的中心部位Ale。在图8所示的广义矩阵Ke用al表示的元素，在图9的矩阵Al中也用al来表示，以便显示广义矩阵Ke与矩阵Al之间的元素的位移。矩阵Al的角的部分A11、A12、A13和A14的元素分别来自侧面部分Ke5、Ke6、Ke7和Ke8并且全部为零。[0257]图10显示了从图8所示的广义矩阵Ke中提取的矩阵A2。该矩阵A2采用通过由小的点所表示的元素，使其按照第一方向（逆时针方向）旋转45°并将其放置在矩阵A2的中心部位A2c。在图8所示的广义矩阵Ke用a2表示的元素8,在图10的矩阵A2中也用a2来表示，以便显示广义矩阵Ke与矩阵A2之间的元素的位移。矩阵A2的角的部分A21、A22、A23和A24的元素分别来自侧面部分Ke5、Ke6、Ke7和Ke8并且全部为零。[0258]更具体而言，如果N是奇数，在第二个子步骤c22中，按照下列关系式或者任何其它等效公式，通过广义矩阵KeT，f=[keijT，f]来计算矩阵AlT，f和A2T，f:[0259],以及，[0261]如果N是偶数，则可以创建等效的公式。[0262]根据本发明方法的第二个实施例，直接执行数据旋转子步骤c2无需构建广义矩阵Ke。[0263]例如，如果N是奇整数，则利用下列关系式或者任何其它等效公式，通过窗口频率变换矩阵KT，f=[kijT，f]来计算两个矩阵AlT，f=[aluvT，f]和A2=[a2uvT，f]:[0264]-如果u+v-1-N-I2和v-u+N+32-1都在1与N之间，包含1和N，[0265]贝Ij[0266]否则aluvT，f=0,以及，[0267]-如果u+v-N-I2和v-u+N+32-1都在1与N-I之间，包含1和N-I，[0268]贝Ij[0269]否则a2uvT，f=0。[0270]如果N是偶数，则可以创建等效的公式。[0271]按照本发明的方法，滤波子步骤c3是与例如现有技术的滤波子步骤相同的。尤其是，通过矩阵Al和A2来计算两个滤波矩阵A14PA2f，例如，或者利用上文所述的投影滤波技术，或者利用上文所述的奇异值分解滤波技术，或者利用任何其它技术。[0272]滤波矩阵A14PA2f所具有的尺寸与通过数据旋转步骤c2所确定的矩阵Al和A2的尺寸相同。[0273]尤其是，如果N是奇数，那么有利的是，滤波矩阵尺寸为NxN以及滤波矩阵A2f的尺寸为N_lXN_l。[0274]此外，在反向数据旋转子步骤c4中，执行针对数据旋转子步骤c2所述变换的大体逆向变换。在该子步骤中，通过在之前的滤波子步骤C3中所计算的滤波矩阵和来计算滤波变换矩阵[0275]该子步骤与上文所述的现有技术的反向数据旋转子步骤c4非常相似或者完全一致。在该子步骤中，通过按照与第一方向（按照顺时针方向旋转45°相反的第二方向旋转，使得滤波矩阵A14PA2f的中心部位的元素都包含并设置在滤波变换矩阵#中。滤波矩阵A1F、A2F的角的部分的元素未包含在滤波变换矩阵#中。这些元素无论如何都是零，只要这些元素都是矩阵Al和A2的角的部分的元素。[0276]例如，如果N是偶整数，则通过下列关系式或者任何等效公式，根据两个矩阵AIf和来计算滤波变换矩阵[0277]如果i-j是偶数，[0280]与现有技术的方法不同，滤波窗口频率变换矩阵KfT，f的尺寸因此为NxN，该尺寸与初始窗口频率变换矩阵KT，f的尺寸相同。[0281]按照现有技术的方法见图6，滤波窗口频率变换矩阵^^幻所具有的尺寸为N+12xN+12。其所包含的元素的数量要少得多:大约比现有技术方法的数据旋转子步骤c2中要丢失的元素少四倍。[0282]通过本发明的方法，滤波窗口频率变换矩阵KfT，f包含大量的元素并且有可能将大量的多重散射考虑在内。因此，可以通过波传播到更深的深度来探测介质。此外，想要生成的介质的图像具有更高的空间分辨率。有可能使用具有少量传感器的传感器组件3,因此该方法和设备的成本较低。[0283]此外，在有利的实施例中，在确定窗口变换矩阵的子步骤cl中，通过窗口时间矩阵KT，t的小波变换来确定各个窗口频率变换矩阵KT，f。[0284]时间函数f的小波变换定义如下：[0286]其中：[0287]Φ是母小波函数，[0288]礼,s是小波系的小波，[0289]*表示复共辄，[0290]u是平移因子，以及，[0291]S是膨胀因子。[0292]因此，小波系与滤波器组相对应。该滤波器组必须包含待研究的频率。于是，小波变换的计算等效于各个滤波器的脉冲响应与各个时间信号的卷积，脉冲响应t。窗口时间矩阵KT，t是在接近时间T和持续时间Δt的时窗内传感器组件的传感器之间的所述脉冲响应hijt。窗口频率矩阵KT，f可以通过小波变换来计算，例如通过上述定义的卷积。[0293]可以在下列各项中选择母小波函数Φ:[0294]-Morlet小波，例如，如图IIa所示：[0295]以及，[0296]-Maar小波二阶高斯导数），例如，如图IIb所示：[0297]其中，.以及，[0298]-渐进小波，例如，如图IIc所示：[0300]因此，可以选择小波系，使之与传感器带宽以及待处理的脉冲信号类型（因此与介质材质最佳地相对应。[0301]通过利用所述小波变换，本发明的方法使之能够通过波传播更深的深度来探测介质且具有更高的精度和分辨率。[0302]此外，根据有利的实施例，在发射步骤a中，传感器组件3发射入射波，所述入射波按照预定角度α倾斜于所述介质1的外表面la，如图12所示。[0303]为此，将传感器的组件3置位于所示的耦合介质7之中。该耦合介质7可置位于组件3与待探测的介质的外表面1之间，或者使该耦合介质7与探针集成为一体，所述探针至少包含组件和所述耦合介质。[0304]例如，耦合介质7所具有的折射指数m不同于待探测的介质1的折射指数m。因此，在该介质中的传播速度分别为Vl和V7。[0305]于是，按照滤波子步骤c3的投影滤波技术，通过下列修改的公式来确定列向量S=[su]的分量：[0307]其中：[0308]df是焦距，它是通过两个介质之间的指数变化进行校准的在耦合介质7中行进的第一距离d7与在介质中行进的第二距离Cl1之和;换言之：[0310]由于采用倾斜入射波的发射，本发明的方法还得以改进并且使之能够通过更大深度的波传播来探测介质且具有更高的精度和分辨率。尤其是，可以得到更多关于探测介质1中任何缺陷的几何图形平面或立体的信息。

权利要求：1.利用传感器组件⑶通过波传播进行探测的方法，所述方法包括：a发射步骤，在该步骤中，传感器组件3向介质（1中发射入射波，所述介质散射所述波，⑹测量步骤，在该步骤中，所述传感器组件3捕捉表示反射波的信号，所述反射波是由介质（1对入射波的反射，所述捕获的信号包括：-简单散射分量，表示是由介质的各个散射体单次反射入射波所形成的反射波的波路径，-而且，在适当情况下，多重散射分量，表示是由介质的散射体上多次连续反射入射波在到达传感器组件之前所形成的反射波的波路径，c处理步骤，在该步骤中，处理所述捕获的信号，以便确定介质的特征，通过滤波至少一个频率变换矩阵，在多重散射分量与简单散射分量之中提取所选择的至少一个分量，所述频率变换矩阵表示在传感器组件的传感器之间的响应，并且至少包括下列子步骤：cl确定窗口变换矩阵的子步骤，在该步骤中，确定与元件之间的响应的窗口时间矩阵KT，t=[kijT，t]相对应的至少一个窗口频率变换矩阵KT，f，所述元件之间响应的窗口时间矩阵在接近时间T和持续时间Δt的时窗内与在传感器组件的传感器之间的时间响应hi」⑴相对应，其中，f是频率，c2数据旋转子步骤，在该步骤中，通过按照第一方向旋转以及从窗口频率变换矩阵提取分量，通过窗口频率变换矩阵KT，f来计算两个矩阵AlT，f和A2t，f，c3滤波子步骤，在该步骤中，将各个矩阵Al、A2中的多重散射分量与简单散射分量相分离，由此得到至少两个滤波矩阵A1F、A2F，所述两个滤波矩阵分别与矩阵A1、A2相对应并且分别表示简单散射分量或者多重散射分量，c4反向数据旋转子步骤，在该步骤中，通过按照与第一方向相反的第二方向旋转以及从滤波矩阵A1f、A2f中提取分量，通过两个滤波矩阵A1f、A2f来计算滤波窗口变换矩阵KfT，f，所述方法的特征在于：-在数据旋转步骤c2中，窗口频率变换矩阵KT，f的所有分量分布在两个矩阵AlT，f和A2T，f中的一个或另一个矩阵的中心部位，这些矩阵所具有的零点在中心部位之外的部分，以及，-在反向数据旋转步骤c4中，提取滤波矩阵A1F、A2F的中心部位的所有非零分量并将其放置在滤波窗口变换矩阵KfT，f之中。2.根据权利要求1所述的方法，其特征在于，所述窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且在数据旋转步骤c2中，两个矩阵AlT，f=[aluvT，f]和A2=[a2uvT，f]可以通过窗口频率变换矩阵KT，f=[kijT，f]按照下列关系来计算：-如果u+v-1-N-I2和v-u+N+32-1都在1与N之间且包含1和N，男|3么aluvT，f—ku+v-1-N-12,v-u+N+32-1T，f，否则aluvT，f=0,并且，-如果u+v-N-I2和v-u+N+32-1都在1与N-I之间且包含1和N-I，那么a2uvT，f—ku+v-N-I2,v-u+N+32-1Tjfj否则a2uvT，f=0。3.根据权利要求1所述的方法，其特征在于，所述数据旋转子步骤c2包括：-第一个子步骤c21，在该子步骤中，将窗口频率变换矩阵KT，f引入尺寸更大的广义矩阵KeT，f中，除了来自窗口频率变换矩阵的元素之外，广义矩阵的所有元素都为零，以及，-第二个子步骤c22，在该子步骤中，两个矩阵AlT，f和A2T，f可以通过广义矩阵KeT，f来计算。4.根据权利要求3所述的方法，其特征在于，所述窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且：-在第一个子步骤c21，广义矩阵KeT，f=[1«^1'，幻]定义如下：如果ie[1+N-I2;N+N-I2]，以及，je[1+N-I2；N+N-I2],那么keijT，f=Ki-N-12,j-N-12Cr，f，否则keijT，f=O，以及，-在第二个子步骤c22中，两个矩阵AlT，f=[aluvT，f]，A2T，f=[a2uvT，f]通过下列关系来计算：ElluvT，f—keu+v-1，v-u+N+32-1+N-12T，f，以cl2uvTjf—keu+v,v-u+N+32-1+N-12T，f〇5.根据权利要求I至4中任一项所述的方法，其特征在于，所述窗口变换矩阵KT，f是尺寸为NxN的矩阵，其中，N是奇数，而且在反向数据旋转步骤c4中，滤波窗口变换矩阵KfT，f=[kFijT，f]可通过滤波矩阵AIfT，f=[alFuvT，f]，A2fT，f=[a2FuvT，f]按照下列关系来计算：如果i_j是偶数，那么kFijT，f=alFi—j2+N+l2,i+j2，否贝1JkFijCT，f=a2Fi—j—i2+N+i2,i+j-12。6.根据权利要求I至5中任一项所述的方法，其特征在于，在确定窗口变换矩阵的子步骤cl中，通过窗口时间矩阵KT，t的小波变换来确定各个窗口频率变换矩阵KT，f，所述窗口时间矩阵KT，t在接近时间T和持续时间Δt的所述时窗内与在传感器组件的传感器之间的时间响应hij⑴相对应。7.根据权利要求1至6中任一项所述的方法，其特征在于，各个窗口变换矩阵KT，f的时窗至少成对地重叠。8.根据权利要求1所述的方法，其特征在于，在滤波子步骤c3中：-从矩阵KT，f的各个系数kijT，f中提取简单散射分量ksijT，f，以便简单散射分量Sks1」T，f在所述矩阵KT，f的各个反对角相互相干，而且，-从矩阵KT，f的各个系数kijT，f中提取多重散射分量kMijT，f。9.根据权利要求1所述的方法，其特征在于，在发射步骤中，发射进入介质（1中的入射波以预定的角度朝所述介质的外表面倾斜。10.适用于实施根据前述权利要求任一所述探测方法的设备，包括传感器组件（3，适合向散射介质（1中发射入射波并且捕获表示反射波的信号，所述反射波是由介质反射入射波所形成的，所述捕获的信号包括：-简单散射分量，表示反射波是由介质各个散射体单次反射入射波所形成的波路径，-以及，在适当情况下，多重散射分量，表示反射波是由介质的散射体上多次连续反射入射波在到达传感器组件之前所形成的反射波的波路径，所述设备进一步包括处理装置5，适合处理所述捕获的信号，以便确定介质（1的特征，处理装置5适合通过至少一个频率变换矩阵的滤波，从多重散射分量和简单散射分量提取所选择的至少一个分量，所述频率变换矩阵表示在传感器组件的传感器之间的响应，并且处理装置5至少执行以下子步骤：cl确定窗口变换矩阵的子步骤，在该步骤中，确定与元件之间响应的窗口时间矩阵KT，t=DujT，t]相对应的至少一个窗口频率变换矩阵KT，f，所述元件之间响应的窗口时间矩阵在接近时间T和持续时间Δt的时窗内与在传感器组件的传感器之间的时间响应hij⑴相对应，其中，f是频率，c2数据旋转子步骤，在该步骤中，通过按照第一方向旋转并且从窗口频率变换矩阵中提取分量，通过窗口频率变换矩阵KT，f来计算两个矩阵AlT，f和A2T，f，c3滤波子步骤，在该步骤中，将各个矩阵Al、A2中的多重散射分量与简单散射分量相分离，由此得到至少两个滤波矩阵A1F、A2F，所述两个滤波矩阵分别与矩阵A1、A2相对应并且分别表示简单散射分量或者多重散射分量，c4反向数据旋转子步骤，在该步骤中，通过按照与第一方向相反的第二方向旋转并且从滤波矩阵A1f、A2f中提取分量，通过两个滤波矩阵A1f、A2f来计算滤波窗口变换矩阵KfT，f，所述设备特征在于所述处理装置⑸适合于：-在数据旋转步骤c2中，窗口频率变换矩阵KT，f的所有分量分布在两个矩阵AlT，f和A2t，f中的一个或另一个矩阵的中心部位，这些矩阵所具有的零点在中心部位之外的部分，以及，-在反向数据旋转步骤c4中，提取滤波矩阵A1F、A2F的中心部位的所有非零分量并将其放置在滤波窗口变换矩阵KfT，f之中。

百度查询： 法国电力公司;国家科学研究中心 利用波传播进行探测的方法和设备

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD