| 主权项 |
一种探地雷达数据比对算法，该算法建立在概率论和信息论的理论基础上，特征在于，依次采用道间配准算法和道内配准算法对两组探地雷达数据进行配准，采用相关系数作为道间配准算法的评价参数，采用条件熵作为道内配准算法的评价参数，对配准后的两组探地雷达数据，采用加窗相关系数比较两组探地雷达数据的差异；其中，所述的道间配准算法包括以下步骤：步骤(A1)：待配准探地雷达数据的配准起始道Y<sub>c</sub>(c<sub>1</sub>),c<sub>1</sub>∈Z<sup>+</sup>和探地雷达参照数据的配准起始道Y<sub>r</sub>(r<sub>1</sub>),r<sub>1</sub>∈Z<sup>+</sup>满足|Y<sub>c</sub>(c<sub>1</sub>)‑Y<sub>r</sub>(r<sub>1</sub>)|＜K,K∈Z<sup>+</sup>，且c<sub>1</sub>＜2K,r<sub>1</sub>＜2K，在待配准探地雷达数据和探地雷达参照数据的Y<sub>c</sub>(1)～Y<sub>c</sub>(2K),Y<sub>r</sub>(1)～Y<sub>r</sub>(2K)范围内搜索配准起始道；步骤(A2)：将2K分为M′等份，求<img file="FDA0000953765120000011.GIF" wi="1163" he="119" />记作Y<sub>c</sub>(X<sub>1</sub>),Y<sub>c</sub>(X<sub>2</sub>),...,Y<sub>c</sub>(X<sub>M′</sub>)和<img file="FDA0000953765120000012.GIF" wi="1163" he="119" />记作Y<sub>r</sub>(X<sub>1</sub>),Y<sub>r</sub>(X<sub>2</sub>),...,Y<sub>r</sub>(X<sub>M′</sub>)之间的最大相关系数<img file="FDA0000953765120000013.GIF" wi="979" he="79" />其中p,q∈{1,2,...,M′}，配准起始道Y<sub>c</sub>(c<sub>1</sub>),Y<sub>r</sub>(r<sub>1</sub>)在Y<sub>c</sub>(X<sub>p</sub>),Y<sub>r</sub>(X<sub>q</sub>)范围内；步骤(A3)：重复步骤(A2)，逐步缩小范围，找到配准起始道Y<sub>c</sub>(c<sub>1</sub>),Y<sub>r</sub>(r<sub>1</sub>)；步骤(A4)：搜索道数为α,α∈Z<sup>+</sup>，若某道待配准探地雷达数据Y<sub>c</sub>(c<sub>a</sub>),c<sub>a</sub>＝1,2,...,N′与探地雷达参照数据Y<sub>r</sub>(r<sub>b+1</sub>),Y<sub>r</sub>(r<sub>b+2</sub>),...,Y<sub>r</sub>(r<sub>b+α</sub>),r<sub>b+1</sub>,r<sub>b+2</sub>,...,r<sub>b+α</sub>∈Z<sup>+</sup>的相关系数最大值满足<img file="FDA0000953765120000014.GIF" wi="1510" he="87" />ε＞0是阈值，则第r<sub>l</sub>道是探地雷达参照数据的最相关道，在待配准探地雷达数据的第c<sub>a‑1</sub>道和第c<sub>a</sub>道之间线性插值l‑a道探地雷达数据，若<img file="FDA0000953765120000015.GIF" wi="651" he="79" />将Y<sub>c</sub>(c<sub>a</sub>)删除；步骤(A5)：通过<img file="FDA0000953765120000016.GIF" wi="1395" he="87" />找到最不相关待配准探地雷达数据和探地雷达参照数据范围是Y<sub>c</sub>(X<sub>p′</sub>),Y<sub>r</sub>(X<sub>q′</sub>)，阈值是最不相关探地雷达数据相关系数的上确界<img file="FDA0000953765120000017.GIF" wi="1323" he="103" /> |
