Submarine topography construction method based on multi-source water depth data integration

摘要

The invention discloses a submarine topography construction method based on multi-source water depth data integration, which achieves reasonable splicing of various data such as multi-beam sounding, single-beam sounding, historical topography map and global topography and submarine topography construction. The invention comprises the following steps of: converting a historical submarine topography map into grid data, converting water depths of other sources into grids, using an overlapping and contrasting method to evaluate the accuracy of the water depths of different sources; using a multi-map-layer method to carry out integration, cutting and splicing on multi-source water depth data, constructing an error distribution model to reasonably revise and integrate various-source water depth data; and constructing multi-source water depth submarine topography and grid. The submarine topography construction method has the advantages that in marine mapping and charting process, the multi-source data can be effectively used, the historical water depth data accuracy and the comprehensive utilization efficiency are promoted, and the submarine topography construction method has important practical application values in marine mapping, submarine topography charting and submarine scientific studies.

主权项

1. A computer implemented submarine topography construction method executed on a processor based on multi-source water depth data integration, characterized by comprising the steps as follows and in the following order:
step 1: converting a historical submarine topography map into grid data; (1) vectorizing an original submarine topography map: scanning a paper pattern submarine topography map through a flatbed scanner, performing electronic vectorization on a plane topography map, and reserving depth attribute of each counter, wherein the vectorized electronic topography map needs to be perfectly matched with the plane topography map, the vectorized electronic topography map forms an electronic topography data set D1={t1i}, and each counter is composed of a water depth set {x1j, y1j, z1j}, and t1i={x1j, y1j, z1j}, wherein i and j are natural number; (2) adding feature auxiliary lines: checking the vectorized electronic topography map and adding auxiliary lines for regions having saddle topography and positive and negative value topography, wherein the auxiliary lines added are consistent with the trend of a vector line in the electronic topography map, and a data set D2={t2i} is formed after the auxiliary lines are added, wherein each counter t2i is composed of a water depth set {x2j, y2j, z2j}, and t2i={x2j, y2j, z2j}; (3) converting the topography map into a grid: employing a Kriging gridding method to perform data conversion, wherein grid data Ghis(I,J) is formed after the conversion, wherein I, J represent water depth value; (4) re-charting a submarine topography map: for the grid after the conversion, re-charting a submarine topography map according to the same charting parameters as that of the plane topography map in step (1) to form a re-charted topography data set D3={t3i}, wherein each counter t3i is composed of a water depth set {x3j, y3j, z3j}, and t3i={x3j, y3j, z3j}, and wherein the charting parameters refer to projective mode, coordinate system and counter spacing; or, method of overlapping, contrasting and evaluating conversion accuracy: contrasting the data sets formed in step (1) and step (4) through a map layer overlapping and contrasting method; (a) calculating the spatial distance between the old and new counters (t1i={x1j, y1j, z1j} and t3i={x3j, y3j, z3j}: dj=√{square root over ((x1j−x3j)2+(y1j−y3j)2)} point by point, (b) calculating the mean error value of the distance between the old and new counters:d_=∑j=1n⁢dj, (c) when d>d, getting back to step (2) to re-add auxiliary lines, wherein d refers to a pre-determined error value of the distance between the old and new counters, when d≦d, outputting a converted grid document Ghis(I,J), or, employing a method of extracting data of the same point and evaluating the conversion accuracy: in the grid data Ghis(I,J), extracting water depth value z4j in sequence based on the coordinates (x1j, y1j) of each counter t1i={x1j, y1j, z1j} in the topography data set and in accordance with an inverse distance weighting method to form a new data set D4={t4i}, wherein each new water depth data t4i is composed of a water depth point set, and t4i={x4j, y4j, z4j}; (a) a method of extracting water depth data: reading t1i={x1j, y1j, z1j} in sequence, querying the position in the grid document Ghis(I,J) according to the coordinates (x1j, y1j), and calculating the water depth value z4j of each point among four proximal points according to the inverse distance weighting method:z⁢⁢4j=∑I=kk+1⁢∑J=ll+1⁢w(I,J)×dep(I,J)∑I=kk+1⁢∑J=ll+1⁢w(I,J)w(I,J)=1(x⁢⁢1i-x(I,J))2+(y⁢⁢1i-y(I,J))2 where, x(I,J) and y(I,J) are a coordinate value of the grid data; w(I,J) is a weighted value; and dep(I,J) is a water depth value in the grid document Ghis (I,J); (b) calculating the mean error value of the water depth between the old and new counters:z_=∑j=1n⁢z⁢⁢4j-z⁢⁢1jn; (c) when z>z, getting back to step (2) to re-add auxiliary lines, wherein d refers to a pre-determined error value of the water depth between the old and new counters, when z≦z, outputting the converted grid document Ghis (I,J), step 2: converting water depths of other sources into grids; (1) measuring the submarine water depth data by using multi-beam sensors, and gridding the measured submarine water depth data through the inverse distance weighting method to form grid data Gmb(I,J); (2) measuring the submarine water depth data by using single-beam sensors, and gridding the measured submarine water depth data through a Kriging gridding method to form grid data Gsb(I,J); (3) measuring the submarine water depth data by using a satellite altimetry and global topography system, and gridding the measured submarine water depth data through Gaussian spline interpolation or inverse distance weighting method to form grid data Goth(I,J); gridding the submarine water depth data measured through the different methods above by the same parameters, including interpolation method, interpolation parameters, coordinates, projection and water depth data measurement datum plane; step 3: evaluating the accuracy of the submarine water depth data measured through the different methods, wherein the following steps are employed to evaluate: (1) using a map layer overlapping and contrasting method to carry out overlapping and contrasting on the submarine water depth data measured through the different methods, and determining the accuracy thereof through the water depth difference of different grids in the water depth point of the overlapped region; (2) contrasting the original submarine water depth data if any, sorting according to the high-low sequence of the water depth accuracy, and reserving the water depth data having the highest accuracy; in case of no original submarine water depth data, taking precedence of the submarine water depth data measured through multi-beam or single-beam method, followed by the submarine water depth data reflected by the historical topography map and finally the submarine water depth data measured through satellite altimetry and global topography library method; step 4: integration, cutting and splicing on multi-source water depth data, wherein the following steps may be employed: (1) taking the water depth data having the highest accuracy determined in step 3 in the overlapped region through a map layer overlapping and contracting method, and correcting the data of the non-overlapped region through the mean error value Δd of the water depth data of the overlapped region; or, employing a gridding method to correct the error water depth data, i.e. constructing a water depth error grid Gerr(I,J) through a spline or Bezier curve according to the water depth difference data set Derr={t1} of the data having the highest accuracy determined in step 3 and other-source data in the same position of the overlapped region, and then using the error grid to correct the data of the non-overlapped region point by point; (2) reserving the data having the highest accuracy in the overlapped region and cutting other data to respectively form data grids Ghis(I,J), Gmb(I,J), Gsb(I,J) and Goth(I,J) after cutting; step 5: constructing multi-source water depth submarine topography and the grid, wherein the submarine topography is constructed through the following steps: (1) converting the data grids Ghis(I,J), Gmb(I,J), Gsb(I,J) and Goth(I,J) formed in step 4(2) into discrete water depth data Dhis={ti}, Dmb={ti}, Dsb={ti} and Doth={ti}; (2) re-constructing the submarine topography grids for the discrete water depth data in step 5(1) by employing a uniform gridding method to form a uniform data grid Gfus(I,J); or, reading various water depth data grids formed in step 4(2) and integrating the grids to form a uniform data grid Gfus(I,J); (3) re-charting a submarine topography map based on the integrated grid Gfus (I,J), overlapping and contrasting the submarine topography map with the topography map in step 1, and evaluating the integrated data accuracy according to the method of method of overlapping, contrasting and evaluating conversion accuracy or the method of data extracting and evaluating conversion accuracy in step 1(4).