Next Article in Journal
A Multilevel Eigenvector Spatial Filtering Model of House Prices: A Case Study of House Sales in Fairfax County, Virginia
Previous Article in Journal
Analyzing the Spatiotemporal Patterns in Green Spaces for Urban Studies Using Location-Based Social Media Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Case Report

Testing Different Interpolation Methods Based on Single Beam Echosounder River Surveying. Case Study: Siret River

Faculty of Science and Environment, European Center of Excellence for the Environment, “Dunarea de Jos” University of Galati, 111, Domneasca Street, 800201 Galati, Romania
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2019, 8(11), 507; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8110507
Submission received: 12 September 2019 / Revised: 4 November 2019 / Accepted: 6 November 2019 / Published: 10 November 2019

Abstract

:
Bathymetric measurements play an important role in assessing the sedimentation rate, deposition of pollutants, erosion rate, or monitoring of morphological changes in a river, lake, or accumulation basin. In order to create a coherent and continuous digital elevation model (DEM) of a river bed, various data interpolation methods are used, especially when single-beam bathymetric measurements do not cover the entire area and when there are areas which are not measured. Interpolation methods are based on numerical models applied to natural landscapes (e.g., meandering river) by taking into account various morphometric and morphologies and a wide range of scales. Obviously, each interpolation method, used in standard or customised form, yields different results. This study aims at testing four interpolation methods in order to determine the most appropriate method which will give an accurate description of the riverbed, based on single-beam bathymetric measurements. The four interpolation methods selected in the present research are: inverse distance weighting (IDW), radial basis function (RBF) with completely regularized spline (CRS) which uses deterministic interpolation, simple kriging (KRG) which is a geo-statistical method, and Topo to Raster (TopoR), a particular method specifically designed for creating continuous surfaces from various elevation points, contour, or polygon data, suitable for creating surfaces for hydrologic analysis. Digital elevation models (DEM’s) were statistically analyzed and precision and errors were evaluated. The single-beam bathymetric measurements were made on the Siret River, between 0 and 35 km. To check and validate the methods, the experiment was repeated for five randomly selected cross-sections in a 1500 m section of the river. The results were then compared with the data extracted from each elevation model generated with each of the four interpolation methods. Our results show that: 1) TopoR is the most accurate technique, and 2) the two deterministic methods give large errors in bank areas, for the entire river channel and for the particular cross-sections.

1. Introduction

The assessment of water bodies (lakes, rivers, etc.) properties makes use of knowledge regarding water quality, temperature, salinity, volumetric and depth measurements. The last contributes to the accurate maintenance of dredging, and the bathymetric mapping of supraglacial lakes, rivers, and streams, etc. The bathymetric models are the best solution to define the bottom surface of any type of lakes or rivers, and these models are obtained by means of depth measurements. The bathymetric light detection and ranging (LIDAR) scanning technology is currently the most applied surveying method for data collection, and it is applied to capture land and seafloor simultaneously, to create a detailed 3D elevation model along the coastline [1]. Although this is an effective and cost-efficient method, it cannot provide enough detail about the stream or river bathymetry, due to its inability to measure depth accurately where the water is cloudy or turbid [2]. In some cases, additional independent measurements are required to confirm LIDAR bathymetry performance, especially in deep, turbid waters or river cloudy waters, as mentioned e.g., by Kinzel [3] and Saylam [4]. Thus, if no bathymetric LIDAR is available, an accurate DEM (digital elevation model) for river beds can be obtained, when real-time kinematic (RTK) global positioning system (GPS) technologies are combined with multibeam or single-beam echosounders [5,6].
River bathymetry or the so-called “bed topography” plays a critical role in flow dynamics numerical modeling, in sediment transport modeling, geomorphologic modeling or in ecological assessments [7]. The best way to collect precise data for a river is by surveying it on a specific direction, namely on transversal cross-sectional paths [8,9]. The bathymetry measurement method is based on a boat-mounted single-beam or multibeam echosounder device combined with a global positioning system (GPS) with real-time kinematic (RTK) correction. Using this combination of measurements, a series of three-dimensional data (x, y, z) can be collected. The spatial resolution of collected data can be improved either by using a very expensive echosounder or by applying different methods of spatial interpolation in a post-processing step. Spatial interpolation of scattered data points has been an active research topic during the last decade since it has been applied in various fields. For example, Li and Heap [10] present a review of 53 comparative studies which assess the performance of 72 interpolation submethods and combined methods, used in different domains. Special attention is paid to the use of spatial interpolation techniques in hydrology, for interpolating scattered elevation data in order to create digital elevation models (DEM), which include the minor bed of river channel. According to Liffner [11], and Merwade [12] and Merwade [13], the DEM is affected directly by the type of interpolation methods used.
The present study aims at demonstrating the importance of choosing the appropriate spatial interpolation methods for interpolating river channel bathymetry data. The performance of four different interpolation techniques were analyzed by applying different statistical methods of performance assessment, such as: Minim value (MinV), maxim value (MaxV), mean error (ME), mean absolute error (MAE), mean square error (MSE), root mean square error (RMSE), root mean square standardized error (RMSSE) and standard deviation (SD) [13,14,15,16,17,18].

Study Area

The Siret hydrographic basin is located in the eastern part of Romania and it is the drainage basin of River Siret, which is the largest and most important tributary of the Danube [19]. The Siret River basin covers an area of 42,890 km2 in Romania, 28.116 km2 of which are managed by the Siret Water Directorate [20] under the name Siret Hydrographic Space (SHS). The study area covers a 35 km section upstream River Siret, along Galati–Sendreni–Independenta, i.e., from its confluence with the Danube, up to Independenta village. River Siret forms a natural border between Galati and Braila counties. Figure 1 shows an aerial map of the area under scrutiny.
The Siret riverbed material is composed of a combination of fine sand and silt [21]. The river bankfull along the studied area is approximately 100 m wide. Since no accurate bathymetric measurements have been made in this area [22], analyzing the bathymetry of this watercourse is scientifically relevant. Moreover, the measurements and three-dimensional information corresponding to this river section are of utmost contribution for enlarging the existing national databases.
The single-beam echosounder (SBES) is the most common instrument used in ports and in lake and river surveys, because it is an efficient, low cost and accessible instrument. The SBES uses acoustic depth measurements, which are based on measuring the elapsed time that an acoustic pulse takes in order to travel from a transducer to the waterway bottom and back. The general formula of the corrected depth is given by [23]:
d = 1 2 ( v t ) + k + d r .
where: d = corrected depth; v = average velocity of sound in the water; t = time from transducer to bottom and back; k = system index constant; dr = distance from reference water surface to transducer.

2. Materials and Methods

2.1. Materials

The bathymetry data were collected by using a boat-mounted single-beam acoustic depth sounder (SBES) linked to a real-time kinematic (RTK) global positioning system (GPS), which can provide sub-decimeter accuracy for the surveyed points. Figure 2 shows several photos taken during data collection.
Bathymetric depths were measured using an Ohmex SonarMite/BTX single-beam echosounder, with a 235 kHz frequency. The accuracy of the echo sounding equipment used to collect the bathymetric depth data is ±0.025 cm, with a ±4-degree beam spread, able to operate between 0.30 m to 75.00 m depth range (software limited). The sound velocity in water ranges between 1400 to 1600 m/s, but, when a sound velocity profiler (SVP) is not available, the echosounder uses a 1500 m/s average value. Depending on the water depth, the echosounder applies 3 to 6 Hz ultrasonic ping rate, with an output data range of 2 Hz [24].
The measurements were made in 5 days, within the interval 22.03.2017–01.04.2017. Each bathymetry point is characterized by the longitude and latitude obtained from the RTK GPS, and by the elevation (z), which is obtained by subtracting the water depth from the water surface elevation recorded at gauge station. The variation of the daily water level was obtained from the Sendreni gauge station. As Figure 3 shows, the water level during the measurement campaign varied between 435 and 460 cm.
The bathymetry data for the surveyed area (Figure 4 and Figure 5) comprise a set of 216816 (x, y, z) points or approximately 45,000 points/km2 measured along the river channel bed.
Figure 4 and Figure 5 present the 9 sections of the whole surveyed area. The total length of the surveyed path reached 103.22 km with an average boat speed of 1.75 knots (0.9 m/s). The distance between the cross-section ranges from 25 to 100 m. This distance depends on the linearity or sinuosity of the main river channel being larger in areas where the river track is straight, and smaller where the river is meandered. Since the water discharge rate is high in the case of River Siret (about 1000 m3/s for a 4 m water level) [25], it was almost impossible to maintain equal distances between the cross-sections measured.
An ortopthomap with a 0.5 × 0.5 m planimetric resolution, produced in 2016, was used in order to represent the data in terms of geospatial location. This map was provided by the National Agency for Cadastre and Land Registration of Romania (ANCPI). The bathymetric measurements and data processing observed the International Hydrographic Organization (IHO) S-44 regulation since Romania does not have a clearly defined national regulation for river bathymetric measurements. According to this regulation, the accuracy of bathymetric depth needs to comply with the requirements of the Special Order [26] which is used only for those areas where under-keel clearance is critical. The accuracy of depth data by IHO-S44 [26] regulation for this type of area is ±0.25 m, by applying the maximum allowable total vertical uncertainty (TVU), thus the results given by the SBES used for the present research article complies with the IHO S-44 regulations.
The Lat-Long coordinates and elevation (Z) of each point were measured with an RTK GPS – South S82-V equipment. The South S82-V is an RTK GNSS receiver, designed for precision acquisition data. This surveying instrument can receive GPS and satellite signals from GLONASS and GALILEO. The planimetric points were measured in the WGS84 coordinate system and the altimetric point in the Black Sea 1975 Constanta height Datum. By additionally using the RTK method the accuracy increases up to ±8 mm + 1 ppm RMS (Root Mean Square) for horizontal data surveying, and up to ±15 mm + 1 ppm RMS for vertical data surveying [27].
The bank points were measured by using a Trimble M5 surveying total station. The maximum distance between the points which mark the left and right river banks was 30 m. A combination of bathymetric and topographic measurements was used in order to achieve higher accuracy when the four interpolation methods were applied. The accuracy of the DEM depends directly on the number of measured points and it is crucial for flood mapping and hydraulic modeling.
To perform various functions, such as automated mapping functions, data management or spatial analysis function, the ArcGIS version 10.2., provided by ESRI (Environmental Systems Research Institute), was used. A GIS software populated with the correct data can help in the organization, interpretation, and communication of environmental data. For example, this particular GIS software has been used to create digital soil mapping, generated by kriging interpolation, in different urban areas [28]. Also, the GIS software was used [29] for comparing spatial interpolation methods in order to estimate the precipitation distribution. In our study, the GIS software was used for spatial data exploration and surface generation, selecting different interpolation methods and statistical analyses which allow developing a DEM for riverbed channel from various data measured at discrete points.

2.2. Methods

The digital elevation model (DEM) is a widely used product and provides a three-dimensional representation (X, Y, Z) of the studied areas. According to Burrough [30], the term DEM can be defined as “a regular matrix representation of continuous variations of space relief units”. This section describes the four different methods which are commonly used for interpolating scattered point data in order to obtain an accurate DEM for this section of River Siret. Calculations were performed by using the z values resulting from the bathymetric field measurements.

2.2.1. Inverse Distance Weighting (IDW)

IDW represents a deterministic spatial interpolation model, based on an assumption that the value at an unsampled point can be approximated by a weighted average of observed values within a circular search neighborhood [31]. The radius of the search circle is defined by the range of a fixed number of closest points. In most cases, the weights used for averaging the data value are a decreasing function of the distance between the sampled and unsampled points [32,33,34]. The general equation used by IDW in order to estimate a value at an unsampled point is given by:
z * = i = 1 N λ i Z i   .
where Zi (i = 1, 2, …, N) are measured values at N points, and the λi is given by:
λ i = ( 1 d i p ) / (   i = 1 N 1 d i p )
where di is the distance between ith sampled point (zi) and the unsampled point (z*), and p is the exponent variable. A higher power exponent results in less influence from distant points, and vice versa. The standard power value is two.

2.2.2. Simple Kriging (KRG)

Kriging (KRG) is a method developed in the 1960s by the French mathematician Matheron [35] and exemplified in detail by Coburn [36]. This is a geostatistical method of interpolation which is very useful in different fields. Maps generated by this method have a very well-structured visual appearance, all based on irregular spatial data [35]. The general model of this method is based on a constant μ(s) for the data set and random errors ε(s) which have spatial dependence [37,38,39,40]:
Z ( s ) = μ ( s ) + ε ( s ) .  
The prediction of value in a new cell position is done by summing the values of the weight of the multiplied points with the Z value of the data used in the interpolation [41]:
Z ^ ( s 0 ) = i = 1 N λ i Z ( s i ) ,    i = 1 N λ i = 1
Being a very flexible method, it can be used by default or customized to match the data by specifying the nearest semi-variogram model. The semi-variogram consists of two components: The experimental semi-variogram (empirical) and the theoretical model of the semi-variogram [42]. The experimental semi-variogram results from the calculation of the variance of each measured point versus the other points used in spatialization [43]:
λ ^ ( h j ¯ ) = 1 2 N j i = 1 N j [ Z ( s i ) Z ( s i + h ) ] 2
where: Nj is a set of pairs of locations separated by the distance h; h ¯ —represents the average of the distances between distinct pairs Nj.
Data analysis for KRG interpolation was performed using two particular cases of this method, namely: Simple kriging (SKRG) and universal kriging (UKRG).
The SKRG method involves determining the unknown value by estimating the weighted average values of the neighboring points [44]:
Z ( s 0 ) = i = 1 N λ i Z ( s i ) ± ε .
where Z ( s 0 ) represents the estimated value, Z ( s i ) is the value of the neighboring points, λ i is the coefficients of weight, which satisfy the condition i = 1 N λ i = 1 , and ε is the default estimation error.
The UKRG method uses the assumption that the spatial variation of z is dependent on three components: a data set structure, a correlated random component, and a residual error. This is a hybrid method, where the spatial trend is measured by polynomial surfaces of different orders, derived globally or locally [45]:
i = 1 N λ i f k ( s i ) = f k ( s 0 ) .
where Z ( s 0 ) represents the estimated value of the unknown variability, μ ( s 0 ) is the deterministic function, and ε ( s 0 ) is random variation.

2.2.3. Radial Basis Function (RBF)

The radial basis function (RBF) is a interpolation method similar to the kriging family, but it does not benefit from the contribution of space spatial data analysis through the variogram [46]:
Z ^ ( s 0 ) =   i = 1 N ω i φ ( s i s 0 ) + ω n + 1 .
where φ ( r ) is the radial base function, r = s i s 0 is the radial distance between the point for which a new s0 value is calculated and the points with measured values si, and ω symbolizes the weights to be estimated.
Basic radial functions are often used to create neural networks. This method is used to calculate surfaces composed of a very large number of points [47]. According to Dumitrescu [48], the most applied RBF functions are:
  • Multiquadric function (MQ):
    φ ( r ) = ( r 2 + σ 2 ) 1 / 2 .
  • Thin-plate spline (TPS):
    φ ( r ) = ( σ + r ) 2 ln ( σ + r ) .
  • Spline with tension (ST):
    φ ( r ) = ln σ r 2 + K 0 σ r + C e
    where K0(x) is the modified Bessel function.
  • Completely regularized spline (CRS):
    θ ( r ) = n = 1 ( 1 ) n ( σ r ) 2 n n ! n = ln ( σ r 2 ) 2 + E 1 ( σ r 2 ) 2 + C e .
    where E1(x) is the exponential integration function, and Ce is the Euler constant.

2.2.4. Topo to Raster Interpolation (TopoR)

The Topo to Raster interpolation method (TopoR) is specifically designed for creating hydrologically correct digital elevation models (DEMs) [49]. It is based on the ANUDEM program developed by Michael Hutchinson [50,51,52,53]. This method was designed to take advantage of all types of possible input data which typically characterize the landforms (elevation points, contour lines, stream centerline, sink, boundary, lake, cliff, coast) [54].
The Topo to Raster interpolation was developed by the ANU (Australian National University), but this method is often used in many countries and research papers [55,56]. It uses a finite differential interpolation technique and it is optimized to streamline the local interpolation computational method such as IDW interpolation, without losing the continuity of surfaces obtained by the kriging or RBF Spline interpolation methods. The difference is defined as the first and second degree of partial derivation f of the interpolation method described by the following Equations [56]:
J 1 ( f ) = ( f x 2 + f y 2 ) d x d y .  
J 2 ( f ) = ( f x x 2 + f x y 2 + f y y 2 ) d x d y .  
J1 and J2 must be minimized in order to eliminate the maximum and minimum peak effect (excessively smooth or peaked surface) and to obtain a realistic surface terrain. If only J2 is minimized, a very smooth surface is obtained, and vice versa, if only J1 is minimized, maximum and minimum peaks occur [57]. Hutchinson [50] suggests applying a roughness penalty by taking into account the cell resolution so as to reduce this effect:
J ( f ) = 0.5 h 2 J 1 ( f ) + J 2 ( f ) .  
Essentially, TopoR is a combined interpolation method which uses a discrete technique based on spline polynomial functions of degree m and smoothness k, where the roughness can be modified to allow the generation of DEM with steep changes in the field, such as areas influenced by depths of currents, ridges or cliffs. These landforms are also called local maximums or local minimums. In the present research, point elevation features were used.

3. Results and Discussion

3.1. Block Data Performance Analysis

After collecting and validating topo-bathymetric data, digital models were generated (Figure 6), by using the four interpolation methods described in the previous section.
Figure 6a–d presents the four digital elevation models: (a) inverse distance weighting (IDW) interpolation; (b) kriging (KRG) interpolation; (c) radial basis function (RBF) interpolation with completely regularized spline (CRS) method; (d) Topo to Raster interpolation method, for a selected river section of 1500 m. The spatial resolution of 5 × 5 m was used for all DEMs.
All interpolation methods were run by using Spatial Analyst Tools from ArcGIS 10.2. For the simple kriging method, a single spherical variogram model was fitted to the sample variogram of all data by varying different nugget, sill, and range to fit the variogram model. Special attention was paid to matching the slope for the first several reliable lags. Due to this fact a spherical semivariogram, with Nugget = 0.012, Sill = 1.15, Range = 106.44, and Lag Size = 13.30 was adopted.
In the case of the TopoR interpolation method, the differences in elevation ranges were least obvious than those resulting from the IDW and RBF interpolation methods. The TopoR DEM has a much smoother and flatter graphic representation, with lower sinuosity elements.
The basic statistical value which defines the dispersion of the frequency distribution of deviations between the measured and unmeasured interpolated values is the mean-standard deviation (mean SD). Figure 7 shows the Gauss curve representation (a) and boxplot chart (b) for each model, which illustrate how closely the created model predicts the measured values. The Gauss distribution shown in Figure 7a proves that the narrowest line is associated to the TopoR method. A good model will be described by a narrow error distribution and a small, close to zero, SD. Both plots in Figure 7a,b suggest that TopoR is the best model, since the associated error histogram is the narrowest and the mean SD is the smallest.
For a fair comparison of the interpolated points on the chosen section, the same number of points and the 5 m × 5 m spatial resolution of the raster were used. In order to assess the accuracy of each interpolation methods with the SBES measurements, some statistical tests whose results are shown in Table 1 and Table 2, were performed. Table 1 describes the results of obtained mean elevation and variance for all points along the 1500 m river section, and the differences between the depth measurements (SBES) with the four interpolators by using a two-tailed t-test. Values of correlation coefficients and results of the regression analysis are shown in Table 2.
Scatterplots of interpolated (Vi) and measured values (Vm) can be seen in Figure 8a–d, together with the corresponding correlation coefficient. The four methods give reliable results for most elevation values, except for the interval 4–6 m. Correlation coefficients are close to 1 for all methods (see Table 2); however, not all methods give similar results. Table 2 shows that the RBF method has the lowest performance, with the smallest regression coefficient, R2 = 0.947, and a correlation coefficient of 0.973. This is in line with Arseni [18] and Rosu [58], who found a similar result.
The TopoR interpolation method gives the best results, with the highest values for both regression and correlation coefficients (R2 = 0.973 and 0.988). Also, the DEM produced with TopoR has the lowest number of points within the problematic sector of 4–6 m. This is in line with Arseni [22], who found a similar result.

3.2. Local Cross-Section Performance Analysis

In order understand better the qualities of each of the four interpolation methods, the individual behaviour of the models for five cross-sections chosen for validation was also analyzed. Thus, five random cross-sections were chosen for repeated surveying and data validation from the 1500 m studied area (Figure 9).
The SBES profile was obtained from measurements by using the linear method (cross-sectional survey). The four modeled profiles were obtained by extracting each point value from the generated DEMs. Each cross-section surveyed by SBES bathymetric measurements was compared with the values extracted from each DEM characteristic of each interpolation method. The differences between each model and measurements for each cross-section is shown in Figure 10. Detailed information about the five statistical parameters (variables) for models and measurements are presented in Table A1 (Appendix B): minimum value (MIN), maximum value (MAX), mean value (MEAN), median (MEDIAN), and standard deviation (SD), calculated by using the ArcGIS 10.4 software. These parameters were used to assess the success of the data analysis.
The differences between SBES measurements and interpolated elevation value, for each of the five section, are shown in Figure 10 in order to assess, in a direct manner, the performance of each model. Figure 10 shows that all models reproduce very well the river bed except for two regions of about 20–30 m from the river banks, which are basically regions of small depths or positive elevation (as seen in Figure A1/Appendix A).
The deviation increases when approaching both banks for all models. TopoR slightly overestimates the elevation for all sections, while IDW and RBF underestimate the elevation at the shores. The KRG model has an irregular behavior, with values which are either larger or smaller than measurements. Obviously the TopoR model is the best, since the deviations from measurements are the smallest (1 m). On the other hand, IDW has the worst performance, with differences between the model and the field measurements of more than 5 m (5.35 m for cross-Section 4, close to right bank of the river flow). It seems that the IDW, KRG and RBF methods are not reliable as long as the river banks have steep slopes.
Further, all data have been normalized by the value of the elevation (Z), which represents the elevation of the river bed channel. Normalization was performed by recalculating the Z values between 0 and 1, using the following formula:
Z i = X i min ( X ) max ( X ) min ( X )
where X = (X1, X2, …, Xn) and Zi represent the ith normalized data.
This type of normalization is also called “feature scaling”, or “unity-based normalization”. Another type of normalization (like normalization by dividing to maximum value) cannot provide scaled data between 0 and 1, considering the negative elevation values.
The normalized profiles obtained from SBES measurements and the four interpolation methods for each cross-section are shown in Figure 11.
Figure 11 shows that each transversal profile obtained by interpolation methods reasonably describes the measured profile of the selected cross-sections of the river bed channel. The good performance of the TopoR DEM (black line) is illustrated in Figure 11, since the corresponding profile is very close (sometimes even superimposed) to the profile given by SBES bathymetric measurements (red line). The other three models give sometimes erroneous results, especially in the area close to the right and left banks, as shown in Figure 10. This can be also seen in Appendix A, Figure A1, where plots of normalized profiles are compared to real profiles. This is in accordance with the observed higher dispersion seen for low depths (0–5 m) in Figure 8.
Figure 12 shows results of the regression analysis for each of the four methods on SBES measurements (normalized values).
The regression analysis shows that the IDW method has a diffuse data dispersion in all 5 cross-sections. Thus, this is the most inaccurate method. The other two methods, i.e., KRG and RBF, have less regular behavior, with a higher data dispersion for some sections. Again, TopoR is the least dispersive method, whose results are almost perfectly linearly dependent on measurements for all cross-sections.
The statistical analysis was performed by using the Pearson correlation factor (R) and root mean square error (RMSE). The results obtained are shown in Figure 13. The root mean square error (RMSE) or so-called root mean square deviation (RMSD) measures the deviation of output data from the initial input data [59,60]. A small RMSE together with a high R indicate a good fit of the model. Obviously, TopoR has the best performance, closely followed by the KRG method, while IDW is the poorest interpolation method in all cases.
Our results are in line with Curebal’s [61] conclusions for the Keçidere basin, or for models of the Aswan high dam reservoir storage capacity and morphology development [62], which both show that TopoR is the most accurate method for obtaining hydrologically correct DEM. The IDW and KRG methods are reported as being inaccurate when analyzing river bed topography by Goff [63] and Merwade [13]. However, sometimes the kriging (KRG) method estimates the elevation better than IDW [64]. Also, Šiljeg [57] shows that that the results based on TopoR has no significant differences (small errors) in comparison with the aero-photogrammetric measurements (which made up his base model). However, Šiljeg [57] reports that TopoR is not recommended for the micro-level analysis of, e.g., rocks or specific micro-landforms. A complex analysis of many interpolation methods for creating DEMs in the Belgorod region of Podsadneea [65], shows that TopoR is the most accurate method. However, despite the accuracy of interpolation techniques for generating accurate DEM assessed by several studies, there are still no consistent findings about the performances of the spatial interpolators for river bathymetry [66].

4. Conclusions

The creation of the digital elevation model of an area depends on the interpolation method. A good DEM helps calculating the plan and volumetric areas of the river more precisely; thus, DEM can be very useful for evaluating any morphometric changes. This study aimed at identifying the best interpolation method which can give the most accurate representation of a river bed from single-beam river bathymetry measurements.
The bathymetric survey was carried out the on section of River Siret and was based on single- beam echosounder combined with an RTK GPS system. The digital elevation bathymetric models created by single-beam measurement technique depend on the interpolation methods used. In order to assess the accuracy of the interpolated surface, commonly available interpolation methods such as IDW, KRG, RBF and TopoR [61] were used. The Mean SD was one of the statistical indicators which was used to assess the interpolation accuracy. There were no differences between the measured and unmeasured interpolated values in the case of accurate interpolation (e.g., IDW). One can only guess that these differences were determined by cross-validation. The cross-validation of block data showed that the TopoR interpolation has the smallest standard deviation (0.306 m). The coefficients of the determination R2 between the predicted and the measured values were RIDW2 = 0.959, RRBF2 = 0.947, RKRG2 = 0.969 and RTopoR2 = 0.973, confirming that TopoR is the most accurate method. Differences were also observed in the case of the digital elevation model representation. The TopoR interpolation method described an attenuated and smoothed DEM (Appendix A, Figure A2, Figure A3, Figure A4 and Figure A5).
The DEM obtained by using the four interpolation methods were compared with the SBES measurements for the five cross-sections of the river in order to check the accuracy of each method. The TopoR interpolation profile is almost similar to the SBES measurements profile. The profiles obtained by using the IDW, KRG, and RBF methods had sizeable differences as compared to the measured profile, especially in the areas close to the river banks.
The statistical results show that all DEMs are statistically significant, regardless of the method; however, the TopoR interpolation method is the best, with the smallest RMSE and the highest Pearson correlation coefficient, when comparing the modeled and measured depth data. The statistical parameters are: RMSETopoR = 0.276 m and RTopoR = 0.997 m, against the IDW method, which has a mean RMSEIDW = 1.532 m and RIDW = 0.897 m.
Thus, all results show that the TopoR method is the most accurate interpolation method to be used when creating a DEM, for a given number of points and grid.
Ideas of a future study have emerged from the analysis of these tests, and one refers to adding bathymetric longitudinal measurements (along the river channel or perpendicular to cross-sections), thus creating a more regular measurement grid. This may improve the analysis of interpolation methods and will also lead to a higher accuracy. The number of points and pixel size (spatial resolution) is another important parameter worth taking into account. The investigation of the effect of these on digital elevation models is presently on the way.

Author Contributions

Maxim Arseni and Adrian Rosu performed the measurements, establish the methodology, validated the data and analyzed the data; Adrian Rosu designed the figures; Lucian Puiu Georgescu and Catalina Iticescu supervised the progress of teamwork; Maxim Arseni did the statistical analysis and, together with Mirela Voiculescu, prepared the manuscript. All authors commented on the original manuscript.

Funding

This work was partially co-financed by the European Social Fund, through The Human Capital Operational Programme 2014–2020 (POCU), Romania, POCU/380/6/13.

Acknowledgments

Arseni Maxim’s and Rosu Adrian’s research work was supported by the project “ANTREPRENORDOC”, Contract no. 36355/23.05.2019, financed by The Human Capital Operational Programme 2014–2020 (POCU), Romania, and the work of Voiculescu Mirela, Iticescu Catalina and Georgescu Puiu Lucian was supported by the project “EXPERT”, financed by the Romanian Ministry of Research and Innovation, Contract no. 14PFE/17.10.2018.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Profile shape from right to left bank of SBES and extracted from DEM profiles for the Siret River: (a) Cross-Section 1; (b) Cross-Section 2; (c) Cross-Section 3; (d) Cross-Section 4; (e) Cross-Section 5.
Figure A1. Profile shape from right to left bank of SBES and extracted from DEM profiles for the Siret River: (a) Cross-Section 1; (b) Cross-Section 2; (c) Cross-Section 3; (d) Cross-Section 4; (e) Cross-Section 5.
Ijgi 08 00507 g0a1
Figure A2. River Siret depth map obtained by using the IDW interpolation method.
Figure A2. River Siret depth map obtained by using the IDW interpolation method.
Ijgi 08 00507 g0a2
Figure A3. River Siret depth map obtained by using the KRG interpolation method.
Figure A3. River Siret depth map obtained by using the KRG interpolation method.
Ijgi 08 00507 g0a3
Figure A4. River Siret depth map obtained by using the RBF interpolation method.
Figure A4. River Siret depth map obtained by using the RBF interpolation method.
Ijgi 08 00507 g0a4
Figure A5. River Siret depth map obtained by using the TopoR interpolation method.
Figure A5. River Siret depth map obtained by using the TopoR interpolation method.
Ijgi 08 00507 g0a5

Appendix B

Table A1. Statistics of model fits as given by the ArcGIS report.
Table A1. Statistics of model fits as given by the ArcGIS report.
ResultsMINMAXMEANMEDIANSDNumber of Measured PointsWidth [m]
Method
Cross-Section 1SBES−4.6364.272−0.420−0.3962.629237109
IDW−4.5461.512−1.128−1.1291.858
KRG−4.6533.583−0.573−0.5622.445
RBF−4.6682.052−0.880−0.7582.063
TopoR−4.5304.528−0.330−0.2872.702
Cross-Section 2SBES−3.6615.5670.156−0.7782.928190128
IDW−3.6052.163−0.824−0.7061.745
KRG−3.6524.596−0.255−0.8532.585
RBF−3.6492.996−0.611−0.8452.056
TopoR−3.6555.5250.111−0.8442.861
Cross-Section 3SBES−7.9734.569−1.700−1.9894.172174109
IDW−7.9272.304−2.568−2.3663.277
KRG−8.0384.208−2.119−2.3893.833
RBF−8.0003.070−2.347−2.3643.513
TopoR−7.8064.868−1.545−2.0854.267
Cross-Section 4SBES−4.7484.713−0.624−0.8792.755189118
IDW−4.5262.040−1.423−1.5241.894
KRG−4.7514.444−0.797−1.2542.653
RBF−4.7192.901−1.118−1.2522.184
TopoR−4.6885.223−0.481−0.9822.922
Cross-Section 5SBES−2.3944.224−0.326−0.5991.756197137
IDW−2.2942.379−0.890−0.9011.202
KRG−2.3254.180−0.448−0.7371.821
RBF−2.3282.260−0.816−0.7691.209
TopoR−2.2384.962−0.195−0.7131.993

References

  1. Quadros, N.D. Technology in Focus: Bathymetric Lidar. GIM Int. Worldw. Mag. Geomat. 2016, 30, 46–47. [Google Scholar]
  2. Quadros, N.; Collier, P. Integration of bathymetric and topographic LiDAR: A preliminary investigation. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, 36, 1299–1304. [Google Scholar]
  3. Kinzel, P.J.; Legleiter, C.J.; Nelson, J.M. Mapping River Bathymetry With a Small Footprint Green LiDAR: Applications and Challenges. JAWRA J. Am. Water Resour. Assoc. 2013, 49, 183–204. [Google Scholar] [CrossRef]
  4. Saylam, K.; Hupp, J.R.; Averett, A.R.; Gutelius, W.F.; Gelhar, B.W. Airborne lidar bathymetry: Assessing quality assurance and quality control methods with Leica Chiroptera examples. Int. J. Remote Sens. 2018, 39, 2518–2542. [Google Scholar] [CrossRef]
  5. Falkowski, T.; Ostrowski, P.; Siwicki, P.B.M. Channel morphology changes and their relationship to valley bottom geology and human interventions; a case study from the Vistula Valley in Warsaw. Geomorphology 2017, 297, 100–111. [Google Scholar] [CrossRef]
  6. Chénier, R.; Faucher, M.A.; Ahola, R.; Shelat, Y.; Sagram, M. Bathymetric photogrammetry to update CHS charts: Comparing conventional 3D manual and automatic approaches. ISPRS Int. J. Geo-Inf. 2018, 7, 395. [Google Scholar] [CrossRef]
  7. Chaplot, V.; Darboux, F.; Bourennane, H.; Leguédois, S. Accuracy of interpolation techniques for the derivation of digital elevation models in relation to landform types and data density. Geomorphology 2006, 77, 126–141. [Google Scholar] [CrossRef]
  8. Benjankar, R.; Tonina, D.; Mckean, J. One-dimensional and two-dimensional hydrodynamic modeling derived flow properties: Impacts on aquatic habitat quality predictions. Earth Surf. Process. Landf. 2015, 40, 340–356. [Google Scholar] [CrossRef]
  9. Kinsman, N. Single-Beam Bathymetry Data Collected in Shallow-Water Areas near Gambell, Golovin, Hooper Bay, Savoonga, Shishmaref, and Wales, Alaska 2012–2013; Department of Natural Resources. Division of Geological & Geophysical Surveys: Fairbanks, AK, USA, 2015.
  10. Li, J.; Heap, A.D. A review of comparative studies of spatial interpolation methods in environmental sciences: Performance and impact factors. Ecol. Inform. 2011, 6, 228–241. [Google Scholar] [CrossRef]
  11. Liffner, J.; Hewa, G.; Peel, M. The sensitivity of catchment hypsometry and hypsometric properties to DEM resolution and polynomial order. Geomorphology 2018, 309, 112–120. [Google Scholar] [CrossRef]
  12. Merwade, V.; Maidment, D.; Goff, J. Anisotropic considerations while interpolating river channel bathymetry. J. Hydrol. 2006, 331, 731–741. [Google Scholar] [CrossRef]
  13. Merwade, V. Effect of spatial trends on interpolation of river bathymetry. J. Hydrol. 2009, 371, 169–181. [Google Scholar] [CrossRef]
  14. Patel, A.; Katiyar, S.K.; Prasad, V. Performances evaluation of different open source DEM using Differential Global Positioning System (DGPS). Egypt. J. Remote Sens. Space Sci. 2016, 19, 7–16. [Google Scholar] [CrossRef] [Green Version]
  15. Groeneveld, R.A.; Meeden, G. Measuring Skewness and Kurtosis. J. R. Stat. Soc. Ser. D 1984, 33, 391. [Google Scholar] [CrossRef]
  16. Bhunia, G.S.; Shit, P.K.; Maiti, R. Comparison of GIS-based interpolation methods for spatial distribution of soil organic carbon (SOC). J. Saudi Soc. Agric. Sci. 2018, 17, 114–126. [Google Scholar] [CrossRef] [Green Version]
  17. Avila, R.; Horn, B.; Moriarty, E.; Hodson, R.; Moltchanova, E. Evaluating statistical model performance in water quality prediction. J. Environ. Manag. 2018, 206, 910–919. [Google Scholar] [CrossRef]
  18. Arseni, M.; Rosu, A.; Nicolae, A.-F.; Georgescu, P.L.; Constantin, D.E. Comparison of models and volumetric determination for catusa lake, Galati. Tehnomus J. New Technol. Prod. Mach. Manuf. Technol. 2007, 24, 67–71. [Google Scholar]
  19. Romanescu, G.; Stoleriu, C. Causes and effects of the catastrophic flooding on the Siret River (Romania) in July–August 2008. Nat. Hazards 2013, 69, 1351–1367. [Google Scholar] [CrossRef]
  20. Dăscăliţa, D.; Daniela, P.; Olariu, P. Aspects regarding some hydroclimatic phenomena with risk character from Siret hydrographic area. Structural and nonstructural measures of prevention and emergency. Present Environ. Sustain. Dev. 2008, 1, 318–332. [Google Scholar]
  21. Grecu, F.; Zaharia, L.; Ioana-Toroimac, G.; Armaș, I. Floods and Flash-Floods Related to River Channel Dynamics. In Landform Dynamics and Evolution in Romania; Springer: Cham, Switzerland, 2017; pp. 821–844. [Google Scholar] [CrossRef]
  22. Arseni, M. Modern GIS Techniques for Determination of the Territorial Risks Modern GIS techniques for determination of the territorial risks. Ph.D. Thesis, University “Dunarea de Jos” of Galati, Galati, Romania, 2018. [Google Scholar]
  23. US Army Corps of Engineers. Hydrographic Surveying. vol. 5.; Department of the Army: Washington, DC, USA, 2013.
  24. Arseni, M.; Roșu, A.; Georgescu, L.P.; Murariu, G. Single beam acoustic depth measurement techniques and bathymetric mapping for Catusa Lake, Galati. Ann. Univ. Dunarea Jos Galati Fascicle II Math. Phys. Theor. Mech. 2016, 39, 281–287. [Google Scholar]
  25. Murariu, G.; Puscasu, G.; Gogoncea, V.; Angelopoulos, A.; Fildisis, T. Non—Linear Flood Assessment with Neural Network. AIP Conf. Proc. 2010, 1203, 812–819. [Google Scholar] [CrossRef]
  26. International Hydrographyc Organisation (IHO). Standards for Hydrographic Surveys, 5th ed.; Spec. Publ. No. 44; International Hydrographic Bureau: Monaco, 2008. [Google Scholar]
  27. SOUTH S82V Technical Specifications. Available online: https://geo-matching.com/gnss-receivers/s82v (accessed on 19 July 2019).
  28. Panday, D.; Maharjan, B.; Chalise, D.; Shrestha, R.K.; Twanabasu, B. Digital soil mapping in the Bara district of Nepal using kriging tool in ArcGIS. PLoS ONE 2018, 13, e0206350. [Google Scholar] [CrossRef] [PubMed]
  29. De Borges, P.A.; Franke, J.; da Anunciação, Y.M.T.; Weiss, H.; Bernhofer, C. Comparison of spatial interpolation methods for the estimation of precipitation distribution in Distrito Federal, Brazil. Theor. Appl. Climatol. 2016, 123, 335–348. [Google Scholar] [CrossRef]
  30. Burrough, P.A.; MCDonnell, R.A.; Lloyd, C.D. Principles of Geographical Information Systems, 3rd ed.; Oxford University Press: Oxford, UK, 2015. [Google Scholar]
  31. Harman, B.I.; Koseoglu, H.; Yigit, C.O. Performance evaluation of IDW, Kriging and multiquadric interpolation methods in producing noise mapping: A case study at the city of Isparta, Turkey. Appl. Acoust. 2016, 112, 147–157. [Google Scholar] [CrossRef]
  32. He, Q.; Zhang, Z.; Yi, C. 3D fluorescence spectral data interpolation by using IDW. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2008, 71, 743–745. [Google Scholar] [CrossRef]
  33. Li, M.; Yuan, Y.; Wang, N.; Li, Z.; Liu, X.; Zhang, X. Statistical comparison of various interpolation algorithms for reconstructing regional grid ionospheric maps over China. J. Atmos. Sol.-Terr. Phys. 2018, 172, 129–137. [Google Scholar] [CrossRef]
  34. Zhang, H.; Lu, L.; Liu, Y.; Liu, W. Spatial sampling strategies for the effect of interpolation accuracy. ISPRS Int. J. Geo-Inf. 2015, 4, 2742–2768. [Google Scholar] [CrossRef]
  35. Matheron, G. Principles of geostatistics. Econ. Geol. 1963, 58, 1246–1266. [Google Scholar] [CrossRef]
  36. Coburn, T.C.; Yarus, J.M.; Chambers, R.L. Stochastic Modeling and Geostatistics: Principles, Methods, and Case Studies; The American Association of Petrolium Geologists: Tulsa, OK, USA, 2006; Volume 2. [Google Scholar]
  37. Liu, W.; Zhang, H.R.; Yan, D.P.; Wang, S.L. Adaptive surface modeling of soil properties in complex landforms. ISPRS Int. J. Geo-Inf. 2017, 6, 178. [Google Scholar] [CrossRef]
  38. Cressie, N. The origins of kriging. Math. Geol. 1990, 22, 239–252. [Google Scholar] [CrossRef]
  39. Erdogan, S. A comparision of interpolation methods for producing digital elevation models at the field scale. Earth Surf. Process. Landf. 2009, 34, 366–376. [Google Scholar] [CrossRef]
  40. Cheng, M.; Wang, Y.; Engel, B.; Zhang, W.; Peng, H.; Chen, X.; Xia, H. Performance Assessment of Spatial Interpolation of Precipitation for Hydrological Process Simulation in the Three Gorges Basin. Water 2017, 9, 838. [Google Scholar] [CrossRef]
  41. VerH oef, J.M.; Krivoruchko, K. Using ArcGIS Geostatistical Analyst; Esri: Redlands, Australia, 2001. [Google Scholar]
  42. Gundogdu, K.S.; Guney, I. Spatial analyses of groundwater levels using universal kriging. J. Earth Syst. Sci. 2007, 116, 49–55. [Google Scholar] [CrossRef] [Green Version]
  43. Pebesma, E.J. Multivariable geostatistics in S: The gstat package. Comput. Geosci. 2004, 30, 683–691. [Google Scholar] [CrossRef]
  44. Patriche, C.V. Statistic Methods Apllied in Climatology. Terra Nostra: Iasi, Romania, 2009. [Google Scholar]
  45. Tziachris, P.; Metaxa, E.; Papadopoulos, F.; Papadopoulou, M. Spatial Modelling and Prediction Assessment of Soil Iron Using Kriging Interpolation with pH as Auxiliary Information. ISPRS Int. J. Geo-Inf. 2017, 6, 283. [Google Scholar] [CrossRef]
  46. Smith, M.J.; Goodchild, M.F.; Longley, P.A. Geospatial Analysis: A Comprehensive Guide to Principles, Techniques and Software Tools; Troubador Publishing LTD: Leicester, UK, 2007. [Google Scholar]
  47. Buhman, M.D. Radial Basis Functions: Theory and Implementations; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  48. Dumitrescu, A. Spatialisation of Meteorological and Climatic Parameters by SIG Techniques. Ph.D. Thesis, University of Bucharest, Bucharest, Romania, 2012. [Google Scholar]
  49. Childs, C. Interpolating surfaces in ArcGIS spatial analyst. ArcUser July-Sept. 2004, 3235, 569. [Google Scholar]
  50. Hutchinson, M.F. A new procedure for gridding elevation and stream line data with automatic removal of spurious pits. J. Hydrol. 1989, 106, 211–232. [Google Scholar] [CrossRef]
  51. Hutchinson, M.F. A locally adaptive approach to the interpolation of digital elevation models. In Proceedings of the Third International Conference/Workshop on Integrating (GIS) and Environmental Modeling; National Center for Geographic Information and Analysis: Santa Barbara, CA, USA, 1996; pp. 21–26. [Google Scholar]
  52. Hutchinson, M.F. Optimising the degree of data smoothing for locally adaptive finite element bivariate smoothing splines. ANZIAM J. 2000, 42, 774. [Google Scholar] [CrossRef]
  53. Hutchinson, M.F.; Xu, T.; Stein, J.A. Recent Progress in the ANUDEM Elevation Gridding Procedure. Geomorphometry 2011, 19–22. [Google Scholar]
  54. Chen, C.; Li, Y.; Yan, C.; Dai, H.; Liu, G. A Robust Algorithm of Multiquadric Method Based on an Improved Huber Loss Function for Interpolating Remote-Sensing-Derived Elevation Data Sets. Remote Sens. 2015, 7, 3347–3371. [Google Scholar] [CrossRef] [Green Version]
  55. Wang, C.; Yang, Q.; Jupp, D.; Pang, G. Modeling Change of Topographic Spatial Structures with DEM Resolution Using Semi-Variogram Analysis and Filter Bank. ISPRS Int. J. Geo-Inf. 2016, 5, 107. [Google Scholar] [CrossRef]
  56. Salekin, S.; Burgess, J.; Morgenroth, J.; Mason, E.; Meason, D. A Comparative Study of Three Non-Geostatistical Methods for Optimising Digital Elevation Model Interpolation. ISPRS Int. J. Geo-Inf. 2018, 7, 300. [Google Scholar] [CrossRef]
  57. Šiljeg, A.; Lozić, S.; Radoš, D. The effect of interpolation methods on the quality of a digital terrain model for geomorphometric analyses. Teh. Vjesn. Gaz. 2015, 22, 1149–1156. [Google Scholar]
  58. Rosu, A.; Rosu, B.; Constantin, D.E.; Voiculescu, M.; Arseni, M.; Murariu, G.; Georgescu, L.P. Correlations between NO2 distribution maps using GIS and Mobile DOAS measurements in Galati city. Ann. Univ. Dunarea Jos Galati 2018, 41, 23–31. [Google Scholar]
  59. Moharana, S.; Khatua, K.K. Prediction of roughness coefficient of a meandering open channel flow using Neuro-Fuzzy Inference System. Measurement 2014, 51, 112–123. [Google Scholar] [CrossRef]
  60. Xu, Y.; Zhang, J.; Long, Z.; Tang, H.; Zhang, X. Hourly Urban Water Demand Forecasting Using the Continuous Deep Belief Echo State Network. Water 2019, 11, 351. [Google Scholar] [CrossRef]
  61. Curebal, I.; Efe, R.; Ozdemir, H.; Soykan, A.; Sönmez, S. GIS-based approach for flood analysis: Case study of Keçidere flash flood event (Turkey). Geocarto Int. 2016, 31, 355–366. [Google Scholar] [CrossRef]
  62. El-Shazli, A.; Hoermann, G. Development of storage capacity and morphology of the Aswan High Dam Reservoir. Hydrol. Sci. J. 2016, 61, 2639–2648. [Google Scholar] [CrossRef] [Green Version]
  63. Goff, J.A.; Nordfjord, S. Interpolation of Fluvial Morphology Using Channel-Oriented Coordinate Transformation: A Case Study from the New Jersey Shelf. Math. Geol. 2004, 36, 643–658. [Google Scholar] [CrossRef]
  64. Zimmerman, D.; Pavlik, C.; Ruggles, A.; Armstrong, M.P. An Experimental Comparison of Ordinary and Universal Kriging and Inverse Distance Weighting. Math. Geol. 1999, 31, 375–390. [Google Scholar] [CrossRef]
  65. Podsadneaea, E.A.; Pavliuc, B. The interpolation methods for elevation data. City Manag. Theory Pract. 2017, 3, 52–60. [Google Scholar]
  66. Pankalakr, S.S.; Jarag, A.P. Assessment of spatial interpolation techniques for river bathymetry generation of Panchganga River basin using geoinformatic techniques. Asian J. Geoinform. 2016, 15, 9–15. [Google Scholar]
Figure 1. The study area located at the confluence of River Siret with the Danube. The 35 km area measured is positioned between the yellow lines.
Figure 1. The study area located at the confluence of River Siret with the Danube. The 35 km area measured is positioned between the yellow lines.
Ijgi 08 00507 g001
Figure 2. Left: Bathymetric depth measurements with SonarMite single-beam echosounder (SBES) (mounting method on a fiberglass rigid boat) linked to South S82-V RTK GNSS; middle: GNSS and SBES instruments and accessories; right: total station and South S82-V RTK GNSS used for topographic measurements.
Figure 2. Left: Bathymetric depth measurements with SonarMite single-beam echosounder (SBES) (mounting method on a fiberglass rigid boat) linked to South S82-V RTK GNSS; middle: GNSS and SBES instruments and accessories; right: total station and South S82-V RTK GNSS used for topographic measurements.
Ijgi 08 00507 g002
Figure 3. Water level at Sendreni gauge station recorded during the measurements campaign. Flood levels: C.A. (attention level) = 550; C.I. (flood level) = 600; C.P. (risk level) = 650.
Figure 3. Water level at Sendreni gauge station recorded during the measurements campaign. Flood levels: C.A. (attention level) = 550; C.I. (flood level) = 600; C.P. (risk level) = 650.
Ijgi 08 00507 g003
Figure 4. Top and middle: Ortopthomap views of sections S1 (downstream)–S4 (upstream); the blue line shows the trajectory of the boat; the river banks are marked in yellow. Bottom: ortopthomap view of the whole bathymetric survey area, with the first four sections.
Figure 4. Top and middle: Ortopthomap views of sections S1 (downstream)–S4 (upstream); the blue line shows the trajectory of the boat; the river banks are marked in yellow. Bottom: ortopthomap view of the whole bathymetric survey area, with the first four sections.
Ijgi 08 00507 g004
Figure 5. The same as in Figure 4, but for the last five sections of the bathymetric surveying: S5 (downstream)–S9 (upstream).
Figure 5. The same as in Figure 4, but for the last five sections of the bathymetric surveying: S5 (downstream)–S9 (upstream).
Ijgi 08 00507 g005
Figure 6. Digital elevation model (DEM) of a 1500 m river section of the main channel obtained by using the 4 interpolation methods: (a) inverse distance weighting (IDW); (b) kriging (KRG); (c) radial basis function with completely regularized spline (RBF); (d) Topo to Raster.
Figure 6. Digital elevation model (DEM) of a 1500 m river section of the main channel obtained by using the 4 interpolation methods: (a) inverse distance weighting (IDW); (b) kriging (KRG); (c) radial basis function with completely regularized spline (RBF); (d) Topo to Raster.
Ijgi 08 00507 g006
Figure 7. Statistical results of model fit: (a) The frequency of error distribution for each of the four models: IDW (blue), KRG (red), RBF (green), TopoR (magenta); (b) Mean standard deviation between measured and interpolated z-values for each of the four models. IDW, inverse distance weighting; KRG, simple kriging; RBF, radial basis function.
Figure 7. Statistical results of model fit: (a) The frequency of error distribution for each of the four models: IDW (blue), KRG (red), RBF (green), TopoR (magenta); (b) Mean standard deviation between measured and interpolated z-values for each of the four models. IDW, inverse distance weighting; KRG, simple kriging; RBF, radial basis function.
Ijgi 08 00507 g007
Figure 8. The linear regression for the four interpolation methods: (a) IDW; (b) KRG; (c) RBF; (d) TopoR.
Figure 8. The linear regression for the four interpolation methods: (a) IDW; (b) KRG; (c) RBF; (d) TopoR.
Ijgi 08 00507 g008
Figure 9. Locations of the five-cross-section measured with SBES technique.
Figure 9. Locations of the five-cross-section measured with SBES technique.
Ijgi 08 00507 g009
Figure 10. Difference between measurements (by SBES) and interpolated elevation values from each DEM model (SBES–DEM): IDW—red line, KRG—green line, RBF—magenta line, TopoR—black line; negative values correspond to overestimating the elevation.
Figure 10. Difference between measurements (by SBES) and interpolated elevation values from each DEM model (SBES–DEM): IDW—red line, KRG—green line, RBF—magenta line, TopoR—black line; negative values correspond to overestimating the elevation.
Ijgi 08 00507 g010aIjgi 08 00507 g010b
Figure 11. Comparison between the SBES normalized profile (red) and normalized profiles extracted from DEM (orange: KRG, black: TopoR, green: IDW and blue: RBF), for five cross-sections for River Siret: (a) Cross-Section 1; (b) Cross-Section 2; (c) Cross-Section 3; (d) Cross-Section 4; (e) Cross-Section 5.
Figure 11. Comparison between the SBES normalized profile (red) and normalized profiles extracted from DEM (orange: KRG, black: TopoR, green: IDW and blue: RBF), for five cross-sections for River Siret: (a) Cross-Section 1; (b) Cross-Section 2; (c) Cross-Section 3; (d) Cross-Section 4; (e) Cross-Section 5.
Ijgi 08 00507 g011
Figure 12. Simple linear regression of the four DEM-based normalized profile on the normalized SBES profile (left to right: IDW, KRG, RBF, TopoR) for each selected cross-section.
Figure 12. Simple linear regression of the four DEM-based normalized profile on the normalized SBES profile (left to right: IDW, KRG, RBF, TopoR) for each selected cross-section.
Ijgi 08 00507 g012
Figure 13. Root mean square error (RMSE) bar, left axis, and the Pearson correlation coefficient (R) diamonds, right axis, for modeled values for the selected cross-sections: blue: IDW, orange: KRG, grey: RPF, yellow: TopoR. Please note that the best fit is given by high R and small RMSE.
Figure 13. Root mean square error (RMSE) bar, left axis, and the Pearson correlation coefficient (R) diamonds, right axis, for modeled values for the selected cross-sections: blue: IDW, orange: KRG, grey: RPF, yellow: TopoR. Please note that the best fit is given by high R and small RMSE.
Ijgi 08 00507 g013
Table 1. Report on the results from T-test to test the significance of the compared method.
Table 1. Report on the results from T-test to test the significance of the compared method.
ResultsT-Test
Method MeanVarianceCountT-statP(T < = t) Two-Tail
SBES−1.0604.1708709--
IDW−1.0873.72487090.9170.358
KRG−1.0653.96387090.1740.861
RBF−1.0913.70887091.0620.288
TopoR−0.9834.1198709−2.2710.023
Table 2. Results of correlation between SBES and the DEMs obtained by using interpolation methods for the 1500m river section (left column) and regression (right column) for each of the four DEMs.
Table 2. Results of correlation between SBES and the DEMs obtained by using interpolation methods for the 1500m river section (left column) and regression (right column) for each of the four DEMs.
MethodCorrelation TestRegression Analysis
Pearson CoefficientR2SDP-Value
IDW0.9790.9590.4141.53 × 10−38
KRG0.9840.9690.3564.15 × 10−4
RBF0.9730.9470.4681.54 × 10−30
TopoR0.9880.9730.3066.54 × 10−109

Share and Cite

MDPI and ACS Style

Arseni, M.; Voiculescu, M.; Georgescu, L.P.; Iticescu, C.; Rosu, A. Testing Different Interpolation Methods Based on Single Beam Echosounder River Surveying. Case Study: Siret River. ISPRS Int. J. Geo-Inf. 2019, 8, 507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8110507

AMA Style

Arseni M, Voiculescu M, Georgescu LP, Iticescu C, Rosu A. Testing Different Interpolation Methods Based on Single Beam Echosounder River Surveying. Case Study: Siret River. ISPRS International Journal of Geo-Information. 2019; 8(11):507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8110507

Chicago/Turabian Style

Arseni, Maxim, Mirela Voiculescu, Lucian Puiu Georgescu, Catalina Iticescu, and Adrian Rosu. 2019. "Testing Different Interpolation Methods Based on Single Beam Echosounder River Surveying. Case Study: Siret River" ISPRS International Journal of Geo-Information 8, no. 11: 507. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8110507

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop