In relation to 3D bathymetric modelling, this article aims to analyze the performance of Kriging approaches in dependence of the location and density of the measured depth points. The experiments were carried out on a multi-beam sonar (MBS) dataset that includes 240,000 soundings covering a sea-bottom area near Giglio Island (Italy). Seven subsets were derived in random way from the initial regular MBS dataset, selecting an increasing number of points uniformly spaced. Seven models were generated for both Ordinary Kriging and Universal Kriging. Each model was submitted to leave-one-out cross-validation to define the exactness of the predictive values and compared with the initial grid to better evaluate the accuracy in dependence of the point number and dissemination. To investigate this relationship, a new index called MVI (Morphological Variation Index) was introduced as a measurement of the level of variation of seabed morphology. The results validate the efficiency of the Kriging methods and remark the influence of the dataset distribution on the 3D model, highlighting MVI as a useful index to represent the seabed variation as a unique value. Finally, in no rugged areas using 1 point every 1000 m², the RMSE of the differences between measured and interpolated values falls below 1 m, while a further increment of soundings is required in the presence of a high level of variation of seabed morphology.

As far as the knowledge of the seabed is concerned, both for safe navigation and for scientific research, 3D models, particularly digital bathymetric models (DBMs), are nowadays of fundamental importance. This work aimed to evaluate the quality of DBMs according to the interpolation methods applied to obtain grid format 3D surfaces from scattered sample points. Other complementary factors affecting DBM vertical accuracy, such as seabed morphological complexity and surveyed points sampling density, were also analyzed by using a factorial ANOVA experimental design. The experiments were performed on a multibeam dataset provided by the Italian Navy Hydrographic Institute (IIM) with an original resolution of 1 m × 1 m grid spacing, covering a surface of 0.24 km². Six different sectors comprising different seabed morphologies were investigated. Eight sampling densities were randomly extracted from every sector, each with four repetitions. Finally, four different interpolation methods were tested, including: radial basis multiquadric function (RBMF), ordinary kriging (OK), universal kriging (UK) and Gaussian Markov random fields (GMRF). The results demonstrated that both RBMF and OK produced very accurate DBM in areas characterized by low levels of seabed ruggedness at sampling densities of only 0.0128 points/m² (equivalent grid spacing of 8.84 m). In contrast, a higher density of 0.1024 points/m² (3.13 m grid spacing) was required to produce accurate DBM in areas with more complex seabed topography. On the other hand, UK and GMRF were strongly influenced by morphology and sampling density, yielding higher vertical random errors and more prone to slightly overestimate seabed depths. In addition, sampling density and morphology were the factors that most influenced the vertical accuracy of the interpolated DBM. In this sense, the highly statistically significant influence of the interaction between sampling density and morphology on the vertical accuracy of the interpolated DBM confirms the need to perform a preliminary analysis of seabed morphological complexity in order to increase, if necessary, the number of surveyed points in cases of complex morphologies.

Coastline extraction techniques from multispectral satellite images are of great interest for protection and monitoring of coastal areas. In this regard, the Sentinel-2 satellites can give a great contribution thanks to their wide coverage of the earth’s surface. These images can be processed by GIS software, so as to detect the sea from all the rest. However, the traditional supervised classification requires the involvement of the operator to create suitable training sites: this approach, in addition to being associated to the operator’s skill, often takes a long time to be completed. This contribution presents a study carried out on Sentinel-2 dataset and proposes the application of an unsupervised classification method, the k-means, on four different classification indices. The coastlines extracted by unsupervised classification are therefore compared with the coastline manually vectorized from the RGB composition. The results demonstrate the effectiveness of k-means for distinguishing, in the images produced by the indices application, two clusters (water /no-water) in a reduced time lapse if compared with the traditional supervised techniques.

Apps available for Smartphone, as well as software for GNSS/GIS devices, permit to easily mapping the localization and shape of an area by acquiring the vertices coordinates of its contour. This option is useful for remote sensing classification, supporting the detection of representative sample sites of a known cover type to use for algorithm training or to test classification results. This article aims to analyse the possibility to produce smart maps from remotely sensed image classification in rapid way: the attention is focalized on different methods that are compared to identify fast and accurate procedure for producing up-to-date and reliable maps. Landsat 8 OLI multispectral images of northern Sicily (Italy) are submitted to various classification algorithms to distinguish water, bare soil and vegetation. The resulting map is useful for many purposes: appropriately inserted in a larger database aimed at representing the situation in a space-time evolutionary scenario, it is suitable whenever you want to capture the variation induced in a scene, e.g. burnt areas identification, vegetated areas definition for tourist-recreational purposes, etc. Particularly, pixel-based classification approaches are preferred, and experiments are carried out using unsupervised (k-means), vegetation index (NDVI, Normalized Difference Vegetation Index), supervised (minimum distance, maximum likelihood) methods. Using test sites, confusion matrix is built for each method, and quality indices are calculated to compare the results. Experiments demonstrate that NDVI submitted to k-means algorithm allows the best performance for distinguishing not only vegetation areas but also water bodies and bare soils. The resulting thematic map is converted for web publishing.

Bathymetric surveys are carried out whenever there is a need to know the exact morphological trend of the seabed. For a correct operation of the echo sounder, which uses the principle of acoustic waves to scan the bottom and determine the depth, it is important to accurately determine the sound velocity in water, as it varies according to specific parameters (Density, Temperature, and Pressure). In this work, we want to analyse the role of sound velocity determination in bathymetric survey and its impact on the accuracy of depth measurement. The experiments are conducted on data set provided by “Istituto Idrografico della Marina Militare Italiana” (IIM), the official Hydrographic Office for Italy, and acquired in the Ligurian sea. In our case, the formulas of Chen & Millero (UNESCO), Medwin, and Mackenzie were applied. The introduction of errors on chemical-physical parameters of the water column (Temperature, Pressure, Salinity, Depth) simulating inaccurate measurements, produces considerable impacts on sound velocity determination and subsequently a decrease of the depth value accuracy. The results remark the need to use precise probes and accurate procedures to obtain reliable depth data.