一种矿区地表高时空分辨率三维形变估计方法与流程

本发明属于矿区测量技术领域，具体涉及一种基于开采沉陷先验模型的矿区地表高时空分辨率三维形变估计方法。

背景技术：

高时空分辨率矿区地表三维形变监测对于理解矿区开采沉陷动态机理和矿区灾害预防和精确评估具有重要作用。合成孔径雷达干涉测量(interferometricsyntheticapertureradar，insar)技术由于其高空间分辨率、高精度、不受云雨天气干扰等优势在矿区地表形变监测和预计中发挥着巨大的作用。然而，insar技术监测得到的形变是地表真实三维形变在雷达视线向的一维投影，而一维雷达视线向形变很难反应矿区真实地表形变。针对此问题，li等于2015年提出一种名为sip(singleinsarpair)的方法(参见文献1liz,yangz,zhuj,huj,etal.retrievingthree-dimensionaldisplacementfieldsofminingareasfromasingleinsarpair[j].journalofgeodesy,2015,89(1):17–32.)，该方法基于矿区地表形变的物理特性，利用单个insar干涉对反演得到矿区地表三维形变。然而，该方法只能获得组成干涉对的两景合成孔径雷达(syntheticapertureradar，sar)影像期间的矿区地表形变。为了克服该局限，yang等于2017年提出一种sgi(singlegeometryinterferometric)方法将sip方法发展到时序上(参见文献2yangz,liz,zhuj,fengg,etal.derivingtime-seriesthree-dimensionaldisplacementsofminingareasfromasingle-geometryinsardataset[j].journalofgeodesy,2017,92(5):529-544.)，该方法基于一个单几何的insar观测数据集，利用加权最小二乘的方法对多时相干涉对进行联合解算，实现了矿区地表三维时序形变获取方法。

然而，yang的方法存在如下明显的局限：1)对于sgi计算得到的矿区地表时序形变结果，其时间分辨率完全由所选单几何insar数据的影像获取时间间隔决定。目前为止，sar卫星的重访周期通常是固定不变的，但是星载/机载sar卫星的数量却越来越多。而对于sgi方法而言，只要所选sar数据类型固定，则所得到的时序形变结果的时间分辨率将不会高于这一sar传感器的重访周期。对于高度非线性的矿区地表形变而言，这将对矿区地表形变的精确预计和建模，以及对开采导致的地质灾害的精确预计和评估造成严重影响。2)由于矿区大多处于农村以及郊区，其地表覆盖物通常变化较大，加上地表形变速率较快，使得基于单轨insar数据很难获得矿区地表时序形变结果。尤其是当由单轨数据集中存在时间断缺时，利用传统的基于最小二乘的算法将很难获得高度非线性的矿区地表时序形变结果。

技术实现要素：

本发明的目的是提供一种对矿区地表形变能更加精确监测的方法，其利用覆盖待监测矿区并沿轨道飞行的不同合成孔径雷达监测平台获取可干涉sar影像，再基于获取的可干涉sar影像监测时间的差异来增加形变监测的时间分辨率；同时融合多轨道sar影像生成insar数据集，可以对某单轨insar数据集中出现的时间断缺现象进行有效补偿，重构单轨数据集中时间断缺，得到完整的三维时序形变结果，有效地解决单轨数据可能存在时间断缺而导致时序形变解算失败的问题。

本发明提供一种矿区地表高时空分辨率三维形变估计方法，包括如下步骤：

s1：基于每个轨道获取的sar影像分别组成可干涉的insar干涉对集合；

其中，l个轨道的sar影像分别组成l个sar影像集来覆盖待监测矿区，l为大于或等于2的正整数，基于各轨道获得的sar影像集分别生成每个轨道对应的insar干涉对集合，sar影像的总数量为m+1，insar干涉对的总数量为g，如下所示：

g＝g1+…gj…+gl,1≤j≤l

m+1＝m1+…mj…+ml

式中，gj表示第j个轨道生成的insar干涉对集合中insar干涉对数量，mj表示第j个轨道对应的sar影像集中sar影像数量；

s2：基于差分insar技术以及步骤s1中生成的各insar干涉对集合计算得到覆盖待监测矿区地表的一维雷达视线向多时相形变场；

其中，得到的一维雷达视线向多时相形变场的总数量为g；

s3：基于矿区地表水平移动与下沉梯度呈线性比例关系的先验模型分别构建矿区地表沿东西、南北方向的水平移动与下沉梯度之间的函数关系，并基于所述函数关系以及一维雷达视线向多时相形变场，分别计算得到矿区地表沿着垂直、东西和南北方向的三维多时相形变场；

s4：构建矿区地表三维形变速率与三维多时相形变观测场之间的观测方程组；

所述观测方程组为矿区地表分别在垂直、东西和南北方向的三维形变速率矩阵、三维多时相形变观测场矩阵以及干涉对的时间差系数矩阵的关系式方程组；

其中，将所有sar影像按照时间先后排序，并基于排序的sar影像构建干涉对的时间差系数矩阵，以及基于排序的sar影像构建所有相邻sar影像期间地表在垂直、东西和南北方向的三维形变速率矩阵；

所述排序的sar影像的时间t为：t＝[t0,t1,…ti…,tm]，其中，ti为排序的sar影像中第i+1景sar影像的获得时间；

所述干涉对的时间差系数矩阵由每个干涉对的时间差系数向量组成，每个干涉对的时间差系数向量由m个元素组成，按序对应基于排序的m+1景sar影像的时间间隔，其中，第k个insar干涉对的时间差系数向量中第mfk个元素前以及第msk-1个元素以后的所有元素均为0，从第mfk到第msk-1个元素之间的值依次为：

其中，t(k,mfk)、t(k,msk-1)分别表示第mfk和第msk-1元素的值，mfk、msk分别对应第k个insar干涉对的主、辅sar影像的时间索引，分别为生成第k个insar干涉对的主、辅sar影像的时间；分别为基于排序的sar影像的时间t中mfk的下一时间，msk的上一时间；

所述三维形变速率矩阵中的元素为基于排序后的所有相邻sar影像期间在垂直、东西和南北方向的形变速率，所述三维形变速率矩阵包括垂直方向的形变速率矩阵vu、东西方向的形变速率矩阵ve和南北方向的形变速率矩阵vn；

s5：基于排序后的每景sar影像时间构建的干涉对的时间差系数矩阵以及步骤s3获取的矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场求解步骤s4中所述关系方程组得到三维形变速率矩阵；

s6：利用如下公式以及步骤s5计算出的三维形变速率矩阵计算待监测矿区在垂直、东西、南北方向的高时空分辨率三维形变；

du(t0)＝de(t0)＝dn(t0)≡0

其中，du(tk)、de(tk)、dn(tk)分别表示基于排序的sar影像中第k+1景sar影像获取时刻待监测矿区地表在垂直、东西和南北方向的高时空分辨率三维形变累积量；du(t0)、de(t0)、dn(t0)分别表示基于排序的sar影像中第1景sar影像获取时刻待监测矿区地表在垂直、东西和南北方向的三维形变值；

且分别表示步骤s6计算出的三维形变速率矩阵中第i与第i+1景sar影像在垂直、东西、南北方向的形变速率。

进一步优选，步骤s4中的所述观测方程组如下所示：

t·va＝da

且

其中，t为干涉对的时间差系数矩阵，矩阵t的维度为g×m；ut、et、nt分别表示矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场，每个方向的多时相形变观测场的维度均为g×1，三个方向的三维多时相形变观测场共同构成维度为g×3的三维多时相形变观测场矩阵da；vu、ve、vn分别表示垂直、东西、南北方向的形变速率矩阵，三个方向的形变速率矩阵共同构成维度为g×3的三维时序形变速率矩阵va；vui、vei、vni分别表示基于排序的所有sar影像中第i与第i+1景sar影像期间在垂直、东西、南北方向的形变速率。

进一步优选，步骤s5中采用一种广义最小二乘法求解步骤s4中的所述关系方程组得到三维形变速率矩阵，其中，计算公式如下：

其中，(t^t·p·t)⁺表示(t^t·p·t)的moore-penrose伪逆；p为参与解算的各个干涉对所占权重，分别为求解出的三维时序形变速率矩阵，以及垂直、东西、南北方向的形变速率矩阵。

权重矩阵p是由各个轨道sar传感器发射的电磁波的波长、获得的sar影像所组成的干涉对的相干性和精度以及矿区开采导致的地表形变量级和速率等进行确定，其通过多次的实验来得到最优值设置。moore-penrose伪逆是求秩亏方程组的“广义逆”的一种方法。

进一步优选，步骤s3中矿区地表东西、南北方向的水平移动与下沉梯度之间的函数关系如下所示：

κe＝b·r(i,j),κn＝b·r(i,j)(i＝1:rmax，j＝1:cmax，h＝1:g)

式中，和分别表示第h个多时相形变观测场中(i,j)像元位置地表沿东西、南北和垂直方向的位移，κe和κn分别为东西、南北两个方向上的水平移动与下沉梯度之间的比例系数，b和r(i,j)分别表示(i,j)位置处地表形变的水平移动系数和主要影响半径，其由矿区实地采矿条件给定；和表示第h个多时相形变观测场中(i,j)像元位置处的垂直下沉在东西和南北方向上的沉降梯度，分别表示每个像元沿东西和南北方向的地面分辨率，rmax和cmax分别表示影像的最大行列数；

所述一维雷达视线向多时相形变监测值与矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场的关系方程如下所示：

且：

式中，表示第h个一维雷达视线向多时相形变观测场数据，θ^h和α^h分别表示第h个多时相insar干涉对对应的sar传感器的入射角和飞行方位角；分别表示由第h个多时相一维雷达视线向形变场解算得到的垂直、东西、南北方向的三维形变场。

将上述矿区地表水平移动与下沉梯度的比例关系函数带入上述一维雷达视线向形变与三维形变之间的关系方程中，并对所有g个多时相一维雷达视线向形变场依次进行逐行/列迭代求解，即可得到矿区地表全盆地三维多时相形变监测结果ut、et和nt。

进一步优选，步骤s2中采用的差分insar技术为“二轨法”差分insar技术。

有益效果

本发明提供了一种矿区地表高时空分辨率三维形变估计方法，该方法首先将各轨道sar数据分别组成可干涉的insar干涉对集合；再获得各insar数据集联合监测得到的矿区地表一维雷达视线向多时相形变监测值；将一维雷达视线向多时相形变监测与矿区开采沉陷导致的水平移动与下沉梯度呈线性比例关系的先验模型进行联合解算，从而得到矿区地表沿着垂直、东西和南北方向的三维多时相形变场；以最早获得的一景sar影像的时间为参考，在时间域对多轨道sar影像数据集中的所有sar影像进行排序，之后建立排序后时间相邻的sar影像期间的三维形变速率与估计的三维多时相形变观测场之间的观测方程；最后解算观测方程并估计矿区地表高时空分辨率三维形变。该方法计算简单，实现方便；本发明的技术效果主要体现在以下几点：

第一、相对于传统基于不同轨道数据的观测几何之间的差异进行三维形变解算的技术，本方法利用获取时间存在差异的各轨道数据来增加观测的时间分辨率，以获取时间最早的一景影像为形变时间参考点，统一时间基准后能够获得融合后的数据集内的所有sar影像期间的形变结果，解决了现有技术中时间分辨率低的问题。

第二、基于单轨insar数据反演矿区地表三维时序形变的方法，譬如yang的方法，其存在如下缺陷：由于矿区地表时序形变的高度非线性特性、矿区地表覆盖物复杂多变以及insar观测数据集中可能存在的获取时间断缺现象等，单轨insar数据集很难获得准确、完整的矿区地表三维时序形变结果。尤其是当单轨insar数据集中存在时间断缺时，高度非线性的矿区时序形变很难通过最小二乘算法等通常算子进行解算，从而导致矿区地表三维时序形变解算失败。而对于单轨insar数据解算矿区地表三维时序形变时可能的由于数据集内存在时间断缺现象而导致时序形变解算失败的情况，本发明提出的多轨道数据(l≥2)融合解算方法可以通过多个轨道的数据进行融合解算，从而利用多个轨道生成的insar数据集对在某单轨insar数据集中出现的时间断缺进行有效补偿，重构单轨数据集中时间断缺，得到完整的三维时序形变结果，有效地解决单轨数据可能存在时间断缺而导致时序形变解算失败的问题；

第三、在yang等提出的单轨insar数据反演矿区地表三维时序形变的方法中，其首先将一维雷达视线向形变转换到垂直方向，在计算出垂直向时序形变之后，再利用矿区地表形变先验模型将时序垂直向形变转化到水平方向。这种方法在计算中两次应用到先验模型，不可避免的会将先验模型中的不确定性以及估计的时序下沉值中的误差直接甚至加大传递到水平移动中。然而，本发明所提方法对多轨道数据进行处理后，直接利用水平移动与沉降梯度之间的先验模型获得对应的多时相三维形变分量，之后对得到的多时相三维形变分量进行三维时序形变求解。该方法减少了一个转换步骤的同时，计算方案更加优越，在一定程度上可以增强三维时序解算结果抵抗误差的能力，增加解算结果的稳定性。

第四，本发明引入矿区开采沉陷先验模型解决由一维雷达视线向形变解算地表三维形变时方程组病态的问题；根据这种巧妙的思路，我们可以克服利用insar数据进行三维时序形变解算时，尤其是针对多轨insar数据进行三维时序形变解算，传统方法各观测几何的数据获取时间难以同步，时间基准难以统一的问题(尤其是对于矿区高度非线性的形变特性)，从而大大降低了对不同insar数据集观测几何和数据获取时间的苛刻要求，并且提高了预估结果的精度和可靠性；

第五、相比于单轨数据，本发明提供的多轨数据解算方法融合更多的观测数据，增加多余观测量，具有使时序形变结果的精度和可靠性得到提高的潜能。此外，对于矿区高度非线性的形变特性，融合更多观测数据从而得到更高时间分辨率的形变解算结果将更有利于反映矿区变形的细节信息，为矿区安全生产、揭示矿山开采沉陷规律以及灾害防治与评估提供更加精确可靠的技术参考。

附图说明

图1本发明所述方法的流程示意图；

图2基于logistic模型模拟的两个观测数据集以及由它们组成的干涉对集合构成的网络；

图3基于矿区单点沉降模型的模拟实验，其中a图表示模拟的存在一个时间断缺的单轨数据集t1解算得到的时序形变解算结果，b图为用本发明方法融合t1和t2数据集数据得到的时序形变解算结果；

具体实施方式

下面将结合实施例对本发明做进一步的说明。

如图1所示，本发明实施例提供一种对矿区地表形变能更加精确预计的方法，包括如下步骤：

s1：基于每个轨道获取的sar影像分别组成可干涉的insar干涉对集合。

其中，l个轨道的sar影像分别组成l个sar影像集来覆盖待监测矿区，l为大于或等于2的正整数。每个sar影像集的数量分别表示为：m1,m2，…，ml，sar影像的总数量为m+1，即，m+1＝m1+…mj…+ml。

基于每个轨道的sar影像集分别生成每个轨道对应的insar干涉对集合，每个insar干涉对集合中insar干涉对的数量分别为：g1,g2,…,gl，insar干涉对的总数来那个为g，即，g＝g1+…gj…+gl,1≤j≤l。

需要说明的是，本实施例步骤s1中每个轨道对应的insar干涉对数据集是根据不同insar数据源特征和待监测矿区形变量级构建的，具体是根据不同insar数据源的影像质量和雷达波长等特征以及待监测矿区开采引起的地表形变速率共同确定的干涉对时空基线阈值生成的，但阈值的选取和确定具有一定的经验性，即以上述内容为参考设定的经验值。而根据干涉对时空基线阈值如何构建insar干涉对数据集为现有技术内容，因此本发明对此不进行具体的说明。

s2：依靠外部数字高程模型(digitalelevationmodel，dem)去除地形相位，利用“二轨法”差分insar(differentialinsar，dinsar)技术分别处理各生成的insar干涉对，从而获得多轨insar数据联合监测得到的待监测矿区地表一维雷达视线向多时相形变监测值。

其中，得到的一维雷达视线向多时相形变监测值的总数量为g，其中，dinsar技术是现有成熟技术，使用dinsar技术处理insar干涉对也是本学科经典数据处理技术手段之一，因此，本发明对此不进行具体说明。

在实际应用中，由于不同轨道获取的sar影像的成像几何参数不同，而同一轨道的sar影像集的成像几何参数大致相同，故本发明所有sar平台的成像几何参数可分为l组。基于得到的不同sar成像几何的一维雷达视线向多时相形变场进行解算前，是需要对其进行坐标转换的。本实施例选择用于去除insar干涉对中地形相位的外部数字高程模型(digitalelevationmodel,dem)的坐标框架几何参数为参考，并选用gamma软件的diff&geo模块进行处理，将监测得到的各成像几何不同的insar一维雷达视线向多时相形变场进行地理编码和数据重采样，得到统一地理坐标框架后的一维雷达视线向多时相形变场，再执行步骤s3。其中，外部dem的坐标系统通常为wgs1984标准大地坐标系。

s3：基于矿区地表水平移动与下沉梯度呈线性比例关系的先验模型分别构建矿区地表沿东西、南北方向的水平移动与下沉梯度之间的函数关系，并基于所述函数关系以及一维雷达视线向多时相形变监测值与矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场的关系方程分别计算出矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场。

其中，矿区地表东西、南北方向的水平移动与下沉梯度之间的函数关系如下所示：

κe＝b·r(i,j),κn＝b·r(i,j)(i＝1:rmax，j＝1:cmax，h＝1:g)

式中，和分别表示第h个多时相形变观测场中(i,j)像元位置地表沿东西、南北和垂直方向的位移，κe和κn分别为东西、南北两个方向上的水平移动与下沉梯度之间的比例系数，b和r(i,j)分别表示(i,j)位置处地表形变的水平移动系数和主要影响半径，东西、南北方向上的b和r(i,j)分别相同；和表示第h个多时相形变观测场中(i,j)像元位置处的垂直下沉在东西和南北方向上的沉降梯度，分别表示每个像元沿东西和南北方向的地面分辨率，rmax和cmax分别表示影像的最大行列数；

一维雷达视线向多时相形变监测值与矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场的关系方程如下所示：

且：

式中，dlos为一维雷达视线向多时相形变观测场，θ和α分别表示各多时相insar干涉对对应的sar传感器的入射角和飞行方位角；ut、et、nt分别表示矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场，每个维度的多时相形变观测场的数量均为g，分别为g个多时相形变观测场中的第h个多时相形变观测场。

将上述矿区地表水平移动与下沉梯度的比例关系函数带入上述一维雷达视线向形变与三维形变之间的关系方程中，并在insar监测得到的矿区全盆地多时域形变影像中分别进行逐行/列迭代求解，即可得到矿区地表全盆地三维多时相形变监测结果。

s4：构建矿区地表三维形变速率与三维多时相形变观测场之间的观测方程组；

其中，将每景sar影像按照时间先后排序，排序的sar影像的时间t为：t＝[t0,t1,…ti…,tm]，其中，ti为排序的sar影像中第i+1景sar影像时间；

再基于排序的sar影像构建干涉对的时间差系数矩阵t。干涉对的时间差系数矩阵t由每个干涉对的时间差系数向量组成，每个干涉对的时间差系数向量由m个元素组成，按序对应所述m+1景sar影像的时间间隔。因此，本实施例中，时间差系数矩阵t的维度为g×m，每行对应为一个干涉对的时间差系数向量。其中，第k个insar干涉对的时间差系数向量(本实施例中第k行)中第mfk个元素前以及第msk-1个元素以后的所有元素均为0，从第mfk到第msk-1个元素之间的值依次为：

其中，t(k,mfk)、t(k,msk-1)分别表示第mfk和第msk-1元素的值，mfk、msk分别对应第k个insar干涉对的主、辅sar影像的时间索引，分别为生成第k个insar干涉对的主、辅sar影像的时间；分别为基于排序的sar影像的时间t中mfk的下一时间，msk的上一时间。

此外，假设排序后的相邻sar影像期间，矿区地表沿着垂直、东西和南北方向的三维形变速率为恒定值，即相邻sar影像期间在垂直、东西和南北方向的三维形变速率用三维形变速率矩阵表示，其具体为如下：

vu＝[vu1,…vui…,vum]^t，ve＝[ve1,…vei…,vem]^t和vn＝[vn1,…vni…,vnm]^t；

其中，vui、vei、vni分别表示基于排序的sar影像中第i与第i+1景sar影像在垂直、东西、南北方向的形变速率。

基于上述干涉对的时间差系数矩阵t以及三维形变速率矩阵的形式，构建的观测方程组如下所示：

s5：基于步骤s3获取的矿区地表沿着垂直、东西、南北方向的三维多时相形变观测场以及基于每个sar影像时间构建的干涉对的时间差系数矩阵采用广义最小二乘法求解步骤s4中的所述关系方程组得到三维形变速率矩阵。

其中，计算公式如下：

其中，分别为求解出的三维时序形变速率矩阵，以及垂直、东西、南北方向的形变速率矩阵，(t^t·p·t)⁺表示(t^t·p·t)的moore-penrose伪逆；p为参与解算的各个干涉对所占权重。

s6：利用如下公式以及步骤s5计算出的三维形变速率矩阵计算出待监测矿区在垂直、东西、南北方向的高时空分辨率三维形变；

du(t0)＝de(t0)＝dn(t0)≡0

其中，du(tk)、de(tk)、dn(tk)分别表示基于排序的sar影像中第k+1景sar影像获取时刻待监测矿区地表在垂直、东西和南北方向的高时空分辨率三维形变累积量；du(t0)、de(t0)、dn(t0)分别表示基于排序的sar影像中第1景sar影像获取时刻待监测矿区地表在垂直、东西和南北方向的三维形变值；

为了更加清楚的说明本发明及其优势，下述将结合部分图来进一步解释本发明提供的所述方法。在基于矿区开采沉陷先验模型解算得到三维多时相形变观测量后，三维时序形变是基于相同的逐点解算系统完成的。鉴于本发明的特征，本节以矿区单点沉降为例，对本发明进行具体实施例演示与分析。本实施例利用logistic模型(一种矿区地表单点沉降预计数学模型)模拟了地下开采导致的地表单点沉降，如图2a中灰色曲线所示，作为本实施例的原始观测数据模型。设置alospalsar卫星的重访周期：46天作为相邻观测的时间间隔，分别模拟得到两个观测数据集(用t1和t2表示)，如图2a中实心圆和十字星号所示。以某矿区的实际监测数据的时空基线为参考，通过给定的时空基线阈值模拟生成由t1和t2组成的insar干涉对集合。由这些insar干涉对组成以及时空基线情况绘制出如图2b所示的基线集网络，其中黑色网络由t1数据集数据组成，灰色网络由t2数据集数据组成，浅灰色阴影部分表示了在t1数据集网络中存在的时间断缺部位。利用logistic模型将生成的干涉对集合内的各干涉对涵盖时间间隔内包含的沉降量依次计算出来，并分别加入均值为0，均方根误差为3cm的高斯白噪声，从而模拟得到矿区多轨道多时相沉降监测数据集。

对于单轨数据集t1，我们利用经典最小二乘算法进行解算，并和利用本方案得到的多轨道数据融合解算结果一并与相应的模拟真值进行对比，如图3所示。我们可以发现：当单轨数据中存在时间断缺时(如图3a)，由其解算得到的时序沉降结果在时间断缺附近将存在较大误差，并在此之后的时序形变解算结果中进行累积，从而导致整个时序形变解算结果出现较大错误；与模拟真值相比差异非常大，均方根误差达到了25cm，占形变量最大值的25％，如此大的误差在矿区地表形变监测中是不可以被接受的，从而使得时序形变解算失败。然而，当加入多轨数据进行融合解算时(如图3b)，所得到的结果能有效地重建由于某单轨数据中存在的时间断缺部分，消除由于时间断缺导致的解算错误，与模拟真值相比，均方根误差也得到大幅度减小，降低至1.5cm，只占最大形变值的1.5％；且相对于单轨数据，融合更对轨道的insar数据将能够有效提高监测结果的时间分辨率，在本实施例中，相对于单轨数据解算结果，基于两个轨道数据集融合解算的时序形变解算结果的时间分辨率提高了大约43.7％；得到了高精度的时序形变解算结果并增强了形变结果抵抗粗差的能力。

需要强调的是，本发明所述的实例是说明性的，而不是限定性的，因此本发明不限于具体实施方式中所述的实例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式，不脱离本发明宗旨和范围的，不论是修改还是替换，同样属于本发明的保护范围。