Analysis of Very High Spatial Resolution Images for Automatic Shoreline Extraction and Satellite-Derived Bathymetry Mapping

: The amount of Earth observation images available to the public has been the main source of information, helping governments and decision-makers tackling the current world’s most pressing global challenge. However, a number of highly skilled and qualiﬁed personnel are still needed to ﬁll the gap and help turn these data into intelligence. In addition, the accuracy of this intelligence relies on the quality of these images in times of temporal, spatial, and spectral resolution. For the purpose of contributing to the global e ﬀ ort aiming at monitoring natural and anthropic processes a ﬀ ecting coastal areas, we proposed a framework for image processing to extract the shoreline and the shallow water depth on GeoEye-1 satellite image and orthomosaic image acquired by an unmanned aerial vehicle (UAV) on the coast of San Vito Lo Capo, with image preprocessing steps involving orthorectiﬁcation, atmospheric correction, pan sharpening, and binary imaging for water and non-water pixels analysis. Binary imaging analysis step was followed by automatic instantaneous shoreline extraction on a digital image and satellite-derived bathymetry (SDB) mapping on GeoEye-1 water pixels. The extraction of instantaneous shoreline was conducted automatically in ENVI software using a raster to vector (R2V) algorithm, whereas the SDB was computed in ArcGIS software using a log-band ratio method applied on the satellite image and available ﬁeld data for calibration and vertical referencing. The results obtained from these very high spatial resolution images demonstrated the ability of remote sensing techniques in providing information where techniques using traditional methods present some limitations, especially due to their inability to map hard-to-reach areas and very dynamic near shoreline waters. We noticed that for the period of 5 years, the shoreline of San Vito Lo Capo sand beach migrated about 15 m inland, indicating the high dynamism of this coastal area. The bathymetric information obtained on the GeoEye-1 satellite image provided water depth until 10 m deep with R 2 = 0.753. In this paper, we presented cost-e ﬀ ective and practical methods for automatic shoreline extraction and bathymetric mapping of shallow water, which can be adopted for the management and the monitoring of coastal areas.


Introduction
Very high resolution (VHR) images acquired by means of Earth observation satellites or unmanned aerial vehicles (UAV) contain an important amount of Earth surface's information, which can be used to tackle the world's most pressing global challenges, such as climate change, population growth, and natural resource depletion. Satellite images with different spatial resolutions have been particularly used to evaluate coastline morphological change evolution [1], as well as morphological modeling of coastline [2]. Multi-temporal shorelines extracted on VHR satellite images [3] have allowed the reconstruction of shoreline evolution and, eventually, the assessment of coastal erosion process. Advancement both in sensors and algorithms capacity has facilitated onshore [4], as well as offshore [5], studies performed using satellite images. Proxy-remote sensing using UAV allows obtaining VHR data with very detailed information, especially used in agriculture [6] or for other environment studies [7][8][9].
In coastal areas, heavily affected by the impact of climate change [10], VHR images are used for accurate mapping of shoreline evolution [11] and for predicting and determining the coastal storm impacts [12]. Although these VHR images are the main source of information for our past, present, and future Earth system, their acquisition requires sufficient financial means. Their analysis and processing [13] also require highly skilled and qualified personnel in order to extract end-user intelligence. In the current paper, we presented a framework of image analysis for automatic extraction of shoreline both from GeoEye-1 satellite image and orthomosaic image acquired by an unmanned aerial vehicle (UAV), as well as shallow water depth extracted on GeoEye-1 satellite image. The methodology, applied on the coast of San Vito Lo Capo-a large Pocket Beach in the north-west Sicily (Italy), included binary images creation of water and non-water pixels on which shorelines were automatically extracted using ENVI software via raster to vector (R2V) algorithm. The binary image obtained from GeoEye-1 was later used for water masking of blue and green bands, which were used for satellite-derived bathymetry (SDB) mapping in conjunction with available in-situ data. Whereas the past remote sensing studies conducted using R2V approach were land-based studies [14] or for identification of shoreline changes using moderate spatial resolution images [15], the current study proposed the application of R2V analysis on VHR images by using only one input near infrared (NIR) band for GeoEye-1 image and red for UAV orthomosaic). Analysis of the GeoEye-1 image was completed by conducting SDB mapping in order to take the full advantage of VHR images application for on land, as well as on offshore, studies.
Mapping shallow water bathymetry in coastal and near-shore areas is of greatest interest to explorers, hydrographers, and scientists who want to have a better understanding of this highly dynamic environment [16][17][18]. Water depth information is essential for many applications, especially for coastal environment impact assessment and protection, nautical charting, construction planning and leisure, hydrodynamic modeling, and environmental exploration [19,20]. With ongoing climate change, temporal monitoring of shallow water depth can reveal geomorphological and land cover changes occurring in coastal areas, which are important parameters for hydrodynamic and wave modeling [19]. For example, the position of the shoreline is influenced by the level of the sea, and knowledge of water depth can help identify the location of the sand beach and sediment deposits and also improve coastal erosion studies [21]. In underwater archeology, water depth information is used for mission planning and, therefore, facilitates the discovery of lost objects and artifacts [22]. In seismic survey exploration, water depth information constitutes an important component for seabed determination and improves natural resources' mapping.
Usually, ocean, sea, and lake bathymetry studies are conducted using traditional methods, such as multi-beam echo sounder and single beam echo sounder (SBES), fixed on ships and vessels [23,24], and Lidar fixed on aircraft [25,26]. These methods need onsite deployment, exploration license permits, and their equipment are very expensive [27]. They are highly precise but time-consuming when studying large areas. In addition, they present some challenges when they are used in very near-shore due to area inaccessibility and hidden seabed rocks, which can halt the passage of ships and vessels [28]. For Lidar, the main challenge is in its high cost, which makes it unsuitable for small projects; in addition, administrative procedures needed for flight authorization can also complicate the planned missions.
To overcome these issues, SDB is presented as an answer for shallow water bathymetry mapping in order to complete the gap left by traditional methods, as well as to update previously available bathymetric data [29]. Since the 1970s, different SDB approaches have been used to extract depth water information using multispectral satellite images [30]; however, advances in satellite sensors and the development of new algorithms are making this technology more attractive, and it is now being used in many sectors and industries, such as oil and gas, coastal engineering [31], aquaculture, and ports and infrastructures.
Recently, SDB has been adopted as a monitoring technique by different companies, such as EOMAP and TCARTA, offering commercial services related to shallow water bathymetry and seabed surveys [32,33]. These companies provide bathymetry information to different organizations but also to different national hydrographic offices, especially for updating and improving nautical charts. As an example, in October 2015, the BA 2066 chart of Southern Antigua was published by the UK Hydrographic office [34].
With intense activities and high levels of interaction occurring in coastal areas, the shoreline and SDB mappings are regarded as important tools for the monitoring and surveillance of anthropic and natural changes. Application of these tools on the pocket beach of San Vito Lo Capo-presented as the case study-constituted a typical remote surveillance methodology, which can be adopted for small and naturally protected beaches in order to access the impact of changing climate at local as well as regional scales. We intended to explore their ability to identify shoreline reference features, used for shoreline extraction, and eventually coastline analysis. On the other hand, we wanted to determine if SDB could provide very near shoreline sediment positions involved in coastal erosion process occurring in this highly dynamic environment.

Area of Study
San Vito Lo Capo is located in northwestern Sicily in the northernmost sector of the karstic peninsula ( Figure 1). This area has important archeological and touristic sites and is considered as the second touristic destination, in Sicily, after Taormina. San Vito Lo Capo hosts a large sand beach with gradual deep clear water and a harbor facility. The area was selected because of its constant sand shifting of the beach, as testified by the local population, with drifting sand from East to West, depositing into the harbor; this process affects both maritime transport by reducing the space for boat maneuvering and other stakeholders because of the loss of valuable space.
On the pocket beach of San Vito Lo Capo, significant changes in the beach morphology are restricted to severe weather episodes when wave-related processes are enhanced [35]. Moreover, both the biogenic carbonate sediments are provided by the communities living in the seagrass meadows (i.e., Posidonia oceanica). The relict sediments also play an important role in the natural supply of sediments to the beaches. This supply is of even greater importance in areas where there is a lack of sediment from rivers; hence, relict sediments may maintain the beach in a state of equilibrium [36].
The seafloor of the San Vito bay is colonized by extended meadows of the seagrass, which grow on rocky outcrops from depths of 10-20 m and continue along the whole coast of the peninsula (Figure 1).
During the year, prevailing winds are mainly from the first quadrant, sometimes in the form of severe storms, especially during winter; this generates a sediment drift from East to West, which cannot be equilibrated by spring waves, which, in the past, before the construction of the harbor, maintained the equilibrium of the shoreline, pushing the sediments eastward. Based on shoreline and water depth evolutions, it is possible to determine regions occupied by accumulated and lost sediments [37] and help improve the functioning of the harbor and possibly contribute to regional economic development.
Geosciences 2020, 10, x FOR PEER REVIEW 4 of 20 cannot be equilibrated by spring waves, which, in the past, before the construction of the harbor, maintained the equilibrium of the shoreline, pushing the sediments eastward. Based on shoreline and water depth evolutions, it is possible to determine regions occupied by accumulated and lost sediments [37] and help improve the functioning of the harbor and possibly contribute to regional economic development.

Data Used
The satellite image analyzed in this study was acquired on 18 October 2014 on the coast of San Vito Lo Capo by GeoEye-1, a satellite of DigitalGlobe Company. The main characteristics of this image are presented in Table 1. The image had four multispectral bands (blue, green, red, and NIR) with a spatial resolution of 2 m. It also had a panchromatic band acquired with 50 cm. Automatic shoreline extraction was performed using the NIR band, while SDB mapping was conducted using blue and green bands. For SDB mapping, these bands were analyzed in conjunction with in-situ bathymetry acquired in December 2013 (Table 2).

Data Used
The satellite image analyzed in this study was acquired on 18 October 2014 on the coast of San Vito Lo Capo by GeoEye-1, a satellite of DigitalGlobe Company. The main characteristics of this image are presented in Table 1. The image had four multispectral bands (blue, green, red, and NIR) with a spatial resolution of 2 m. It also had a panchromatic band acquired with 50 cm. Automatic shoreline extraction was performed using the NIR band, while SDB mapping was conducted using blue and green bands. For SDB mapping, these bands were analyzed in conjunction with in-situ bathymetry acquired in December 2013 (Table 2).
UAV orthomosaic image used in the current study was a result of preprocessing images acquired using DJI Mavic pro 2 aircraft during the fieldwork conducted on 28 May 2019 on the coast of San Vito Lo Capo. Table 3 presents the main characteristics and aircraft setting parameters for UAV images acquisition. After fieldwork acquisition, these images were preprocessed in Pix4Dmapper software. During preprocessing, a dataset of 1691 images was calibrated with a median of 8862.13 matches per calibrated image. In addition, 54 ground control points (GCPs) were acquired using a global positioning system (GPS) real-time kinematic (RTK) rover. These GCPs were used to enhance UAV orthomosaic georeferencing precision, which resulted in a root mean square error (RMSE) of 0.098 m.

Methodology
In this study, two VHR images were analyzed: (1) GeoEye-1 satellite image was preprocessed in order to convert top-of-atmosphere (TOA) spectral radiance into water reflectance via atmospheric correction process. The atmospherically corrected image was used to extract shoreline position and later to compute SBD, and (2) UAV orthomosaic image was only analyzed in order to extract the shoreline position, which was later used to analyze coastline parameters. Figure 2 shows the steps adopted in order to extend the preliminary results obtained by analyzing the GeoEye-1satellite image for SDB mapping [38] and by integrating UAV orthomosaic for shoreline position mapping.

Methodology
In this study, two VHR images were analyzed: 1) GeoEye-1 satellite image was preprocessed in order to convert top-of-atmosphere (TOA) spectral radiance into water reflectance via atmospheric correction process. The atmospherically corrected image was used to extract shoreline position and later to compute SBD, and 2) UAV orthomosaic image was only analyzed in order to extract the shoreline position, which was later used to analyze coastline parameters. Figure 2 shows the steps adopted in order to extend the preliminary results obtained by analyzing the GeoEye-1satellite image for SDB mapping [38] and by integrating UAV orthomosaic for shoreline position mapping. During the acquisition of satellite image, top-of-atmosphere (TOA) reflectance (Equation 1) saved over water varies with absorption (α), backscattering (β), and water reflectance (ρ).
where ρ (TOA) is the top-of-atmosphere reflectance, ESUNλ is mean solar extra-atmospheric irradiance expressed in W × m -2 × µm-1 , d is Earth-sun distance obtained from Julian day acquisition time, Lλ represents top-of-atmosphere spectral radiance, and θs is solar zenith angle equal to 90 degrees minus the sun elevation angle at the time of image acquisition in the integration manager database (IMD) files.
This reflectance is a function of water depth (Z) and is expressed as follows: However, this reflectance ρ (TOA) is affected by atmospheric composition, and it has been noted that only less than 1% of ρ (TOA) contains information about the water [39], thus an atmospheric correction is needed to remove atmospheric noise, which can affect the quality of the image used for water depth estimation. This atmospheric correction was computed via fast line-ofsight atmospheric analysis of hypercube (FLAASH) [40], a tool that is integrated into ENVI and IDL software (v. 5.5). FLAASH is a physics-based atmospheric correction that corrects wavelengths from During the acquisition of satellite image, top-of-atmosphere (TOA) reflectance (Equation (1)) saved over water varies with absorption (α), backscattering (β), and water reflectance (ρ).
where ρ (TOA) is the top-of-atmosphere reflectance, ESUNλ is mean solar extra-atmospheric irradiance expressed in W × m −2 × µm −1 , d is Earth-sun distance obtained from Julian day acquisition time, Lλ represents top-of-atmosphere spectral radiance, and θs is solar zenith angle equal to 90 degrees minus the sun elevation angle at the time of image acquisition in the integration manager database (IMD) files. This reflectance is a function of water depth (Z) and is expressed as follows: However, this reflectance ρ (TOA) is affected by atmospheric composition, and it has been noted that only less than 1% of ρ (TOA) contains information about the water [39], thus an atmospheric correction is needed to remove atmospheric noise, which can affect the quality of the image used for water depth estimation. This atmospheric correction was computed via fast line-of-sight atmospheric analysis of hypercube (FLAASH) [40], a tool that is integrated into ENVI and IDL software (v. 5.5). FLAASH is a physics-based atmospheric correction that corrects wavelengths from the visible to near-infrared and shortwave infrared regions by removing absorption (α) and backscattering (β) to obtain the water reflectance. The atmospheric correction was followed by the orthorectification process designed for topographic variations correction. In addition, the pan-sharpening process was performed in order to enhance spatial resolution from 2 m to 50 cm. All these processes were computed in ENVI software (v. 5.5).

Automatic Shoreline Extraction
In order to automatically extract shoreline on VHR images, binary images were firstly created both from the GeoEye-1 NIR band and UAV orthomosaic red band. The creation of these images was aiming at distinguishing water and non-water pixels in VHR images. Then, each pixel was coded to 1 for water or 0 for non-water. The attribution of these digital numbers was based on threshold reflectance values extracted on the GeoEye-1 NIR band ( Figure 3) and UAV orthomosaic red band (Figure 4). Table 4 shows the maximum and minimum reflectance values used for water and land masking.
Geosciences 2020, 10, x FOR PEER REVIEW 7 of 20 the visible to near-infrared and shortwave infrared regions by removing absorption (α) and backscattering (β) to obtain the water reflectance. The atmospheric correction was followed by the orthorectification process designed for topographic variations correction. In addition, the pansharpening process was performed in order to enhance spatial resolution from 2 m to 50 cm. All these processes were computed in ENVI software (v. 5.5).

Automatic Shoreline Extraction
In order to automatically extract shoreline on VHR images, binary images were firstly created both from the GeoEye-1 NIR band and UAV orthomosaic red band. The creation of these images was aiming at distinguishing water and non-water pixels in VHR images. Then, each pixel was coded to 1 for water or 0 for non-water. The attribution of these digital numbers was based on threshold reflectance values extracted on the GeoEye-1 NIR band ( Figure 3) and UAV orthomosaic red band ( Figure 4). Table 4 shows the maximum and minimum reflectance values used for water and land masking.  The spectral band of the VHR image could be represented as a matrix with each pixel represented by a reflectance value. This band matrix representation allowed the transformation of the original band into a binary image by applying the masking threshold reflectance values. Figure 5 illustrates an example of the transformation of the original GeoEye-1 NIR band into a binary image by applying threshold reflectance values between land and the water. The image obtained was a binary image composed of two digital numbers-1 for water pixels and 0 for non-water pixels.   Shallow water with low values The spectral band of the VHR image could be represented as a matrix with each pixel represented by a reflectance value. This band matrix representation allowed the transformation of the original band into a binary image by applying the masking threshold reflectance values. Figure 5 illustrates an example of the transformation of the original GeoEye-1 NIR band into a binary image by applying threshold reflectance values between land and the water. The image obtained was a binary image composed of two digital numbers-1 for water pixels and 0 for non-water pixels.  Secondly, the binary images were transformed into vectors via the R2V algorithm installed in ENVI software. This algorithm was capable of detecting the edge pixels between 1 and 0 (water and land) ( Figure 5A') and transform this edge into polyline (shoreline). The process was automatic and could be summarized as follows: (1) detected the edge between water and non-water pixels, (2) established the connectivity between consecutive edges, (3) identified the location of consecutive edges and traced their continuity, (4) exported the polylines with spatial georeference information.  Secondly, the binary images were transformed into vectors via the R2V algorithm installed in ENVI software. This algorithm was capable of detecting the edge pixels between 1 and 0 (water and land) ( Figure 5A') and transform this edge into polyline (shoreline). The process was automatic and could be summarized as follows: (1) detected the edge between water and non-water pixels, (2) established the connectivity between consecutive edges, (3) identified the location of consecutive edges and traced their continuity, (4) exported the polylines with spatial georeference information. Figures 6 and 7 show examples of shorelines extracted on the GeoEye-1 NIR band and UAV orthomosaic red band. Coastline geomorphological analysis conducted, based on the positions of these shorelines, provided information on changes occurred on the coast of San Vito Lo Capo.

Satellite-derived Bathymetry
In this paper, we used an optical-based approach, a log-band ratio method using empirical calibration parameters suitable for bathymetry estimation, focusing on two spectral bands-green and blue [41]. The third band, near-infrared, was used for land/water masking, while field data were used for calibration and vertical referencing to the local datum.
The log band ratio method was applied only on the GeoEye-1 satellite image, as shown in Figure 2. As the image was acquired in good conditions (Table 2), no sun glint correction [42] was needed to remove sea surface reflection.
Thus, after atmospheric corrections, the shallow water reflectance equation could be defined as it is shown in Equation (3).
where Rw is water reflectance, ρ b is bottom reflectance, Z is the depth of the seabed, e is the logarithmic base, and 2K is the (up-welling and down-welling) operational two-way diffuse attenuation coefficient.
Assuming that the area of study has a uniform bottom reflectance and water attenuation, we applied Stumpf et al. [41] equation (Equation (4)) for SDB calculation using two spectral bands in blue and green. The application of this equation minimized the influence of bottom reflectance and water attenuation parameters on the area of study, with respect to water depth variation. It took into account known water depth for calibration and vertical referencing. The method adopted was a log-band ratio method using an empirical iterative solution and attenuation coefficient assumption for bathymetry estimation.
where Z is absolute depth (SDB), m1 (gradient of the line) is a tunable constant to scale the ratio to the depth, Rw is water reflectance, n is fixed value for all areas to ensure the logarithm will be positive and the relationship will be linear, m0 is the offset for the depth of 0 m (z = 0), the interception with the Y axis, λi = band i, λj = band j. The log-band ratio method was computed in ArcGIS software (v. 10.2.2, ESRI) in two main steps: first by extracting relative bathymetry from two bands (blue and green), after applying a low-pass 3-by-3 (kernel size 3 × 3) filter, and, finally, by conducting vertical referencing of this relative bathymetry to local datum. This second step began by determining extinction depth on relative bathymetry and the constants (m1 and m0) used for vertical referencing to local datum. The determination of these constants was performed using Microsoft Excel (v. 2010, Microsoft Corporation).
These constants (m1 = 87,868 and m0 = 73,281) were introduced in Equation (5), in order to transform relative bathymetry into absolute bathymetry. The calibration process required information obtained from available in-situ water depth ( Table 2).

Results
Automatic shoreline extraction and SDB methods applied to VHR images, as described in the current study, provided results with reliable intelligence, demonstrating the impact of natural phenomena in the studied area.

Shoreline Variability
The R2V method applied on satellite image acquired on 18 October 2014 and UAV orthomosaic image acquired on 28 May 2019 showed that for the period of 5 years, different areas of the pocket beach of San Vito Lo Capo have dramatically changed. This change was mainly observed on shoreline reference features observed (Figure 8) during the fieldwork conducted on 28 May 2019. We noticed that this coastal area is dominated by very fine sand (in the center), hard rocks forming the headland, and the harbor located in the western central part. The coastal erosion was noticed on the sand beach indicated by 15 m inland migration of the shoreline. The accretion process was observed in the harbor with 3 to 6 m offshore alternative migration of the shoreline (Figure 9). Table 5 provides information on the changes occurred on the pocket beach of San Vito Lo Capo from 2014 to 2019. It showed an expansion length of 103 m for the sand beach shoreline and 10 m for beach distance. This might be explained by the migration of sand, from East to West, which was deposed in the harbor. For this period of 5 years, coastal erosion led to the submersion of 15 m distance inland with 17,446 m 2 of sand beach eroded. However, a gain of 695 m 2 of the sand was observed, especially in the harbor. These coastline changes indicated the vulnerability of coastal areas, even for the ones with headland believed to be their natural protection. and the harbor located in the western central part. The coastal erosion was noticed on the sand beach indicated by 15 m inland migration of the shoreline. The accretion process was observed in the harbor with 3 to 6 m offshore alternative migration of the shoreline (Figure 9).   Table 5 provides information on the changes occurred on the pocket beach of San Vito Lo Capo from 2014 to 2019. It showed an expansion length of 103 m for the sand beach shoreline and 10 m for beach distance. This might be explained by the migration of sand, from East to West, which was deposed in the harbor. For this period of 5 years, coastal erosion led to the submersion of 15 m distance inland with 17,446 m 2 of sand beach eroded. However, a gain of 695 m 2 of the sand was observed, especially in the harbor. These coastline changes indicated the vulnerability of coastal areas, even for the ones with headland believed to be their natural protection.

Satellite-derived Bathymetry
The application of the log-band ratio method in a highly dynamic area of San Vito Lo Capo using blue and green spectral bands provided, in the first step, the relative bathymetry. This bathymetry was later calibrated and vertically referenced to the local datum using available in-situ data in order to obtain absolute bathymetry. Figure 10 shows SBD obtained from the GeoEye-1 satellite image; we noticed the variation of water depth from 0 m to 9 m with high values observed in the northwestern part and low values near the shoreline of the sand beach. The western near-shoreline part showed relatively high values, suggesting the impact of the wind transporting sand from East to West as testified by the local population during our visit on 28 May 2019 to this area.

Satellite-derived Bathymetry
The application of the log-band ratio method in a highly dynamic area of San Vito Lo Capo using blue and green spectral bands provided, in the first step, the relative bathymetry. This bathymetry was later calibrated and vertically referenced to the local datum using available in-situ data in order to obtain absolute bathymetry. Figure 10 shows SBD obtained from the GeoEye-1 satellite image; we noticed the variation of water depth from 0 m to 9 m with high values observed in the northwestern part and low values near the shoreline of the sand beach. The western near-shoreline part showed relatively high values, suggesting the impact of the wind transporting sand from East to West as testified by the local population during our visit on 28 May 2019 to this area. In comparison with traditional method results obtained using SBES, Figure 11 shows that SDB provided additional information, especially in near-shore waters. In addition to water depth information, the use of satellite images provided information on the position of the shoreline, on the extent of sand beach, on the sand beach dynamics, as well as on built-up areas. These are useful information for scientists, hydrographers, and policymakers operating in coastal areas, especially in coastal management planning. In comparison with traditional method results obtained using SBES, Figure 11 shows that SDB provided additional information, especially in near-shore waters. In addition to water depth information, the use of satellite images provided information on the position of the shoreline, on the extent of sand beach, on the sand beach dynamics, as well as on built-up areas. These are useful information for scientists, hydrographers, and policymakers operating in coastal areas, especially in coastal management planning.
The comparison of SDB information with in-situ bathymetry showed a good relationship with a coefficient of determination R 2 = 0.7534 ( Figure 12). This high correlation between SBD and in-situ bathymetry demonstrated the ability of our SDB methodology to provide shallow water bathymetry, even in a highly dynamic environment.
Qualitative assessment of the results (Table 6) showed values for different depth ranges on which errors were calculated in function with available in-situ bathymetry in order to determine the overestimated and underestimated values on SDB. While some of these vertical uncertainties were in the range of maximum allowable total vertical uncertainty (TVU) values in ideal environmental conditions, the outrange ones might result in the presence of suspended sediments in this area. The comparison of SDB information with in-situ bathymetry showed a good relationship with a coefficient of determination R 2 = 0.7534 ( Figure 12). This high correlation between SBD and in-situ bathymetry demonstrated the ability of our SDB methodology to provide shallow water bathymetry, even in a highly dynamic environment. Qualitative assessment of the results (Table 6) showed values for different depth ranges on which errors were calculated in function with available in-situ bathymetry in order to determine the overestimated and underestimated values on SDB. While some of these vertical uncertainties were in the range of maximum allowable total vertical uncertainty (TVU) values in ideal environmental conditions, the outrange ones might result in the presence of suspended sediments in this area.  The comparison of SDB information with in-situ bathymetry showed a good relationship with a coefficient of determination R 2 = 0.7534 ( Figure 12). This high correlation between SBD and in-situ bathymetry demonstrated the ability of our SDB methodology to provide shallow water bathymetry, even in a highly dynamic environment. Qualitative assessment of the results (Table 6) showed values for different depth ranges on which errors were calculated in function with available in-situ bathymetry in order to determine the overestimated and underestimated values on SDB. While some of these vertical uncertainties were in the range of maximum allowable total vertical uncertainty (TVU) values in ideal environmental conditions, the outrange ones might result in the presence of suspended sediments in this area.

Discussion
The current study explored the potential of very high spatial resolution images for shoreline and shallow water bathymetry mappings based on raster to vector (R2V) and log-band ratio methods, respectively. The presented methodology applied the R2V method for automatic shoreline extraction on GeoEye-1 satellite image and unmanned aerial vehicle (UAV) orthomosaic image, whereas the log-band ratio method was applied on GeoEye-1 satellite image for SDB. By analyzing spectral reflectance of water and non-water areas, very high resolution (VHR) images open up opportunities for applications once reserved for in-situ fieldworks. Based on threshold reflectance values of water and non-water pixels, we could extract the shoreline position on satellite, as well as on UAV images. Such shoreline position extracted on multi-temporal remote sensing data allows the assessment of the impact of the natural phenomenon occurring in coastal areas. In addition, based on the reflectance of water pixels, we demonstrated that water depth information could be extracted on the satellite image and calibrated with in-situ bathymetry, eventually allowing the mapping of shallow water bathymetry.

Automatic Shoreline Extraction on GeoEye-1 Satellite and UAV Orthomosaic Images
The results indicated that the coastal area of San Vito Lo Capo has distinct remotely sensed coastline features (Figure 8), useful for identification of the shoreline position on very high spatial resolution images acquired by satellite, as well as by UAV. Thus, the application of the R2V algorithm on binary images contributes to an automated shoreline mapping, a method much faster than traditional manually digitalization shoreline mapping. The ability to create a binary image from one single band (NIR or Red) and the availability of R2V in software, such as ENVI or ERDAS image, make the adopted mapping process a cost-effective tool, useful for monitoring shoreline position variability and eventually allowing coastal erosion impact assessment. Using the NIR band with 0.5 m spatial resolution and UAV orthomosaic image with 1.58 cm spatial resolution, the results obtained showed that from 2014 to 2019, the shoreline has migrated 15 m inland equals to approximately 17% loss of the sand beach. Such a level of detail is only achieved by analyzing very high spatial resolution images. While this mapping process provides accurate shoreline position due to the small size of a pixel of images used for the current study, it can also be adopted for the same studies conducted using freely-available open access moderate spatial resolution satellite images, such as Copernicus Sentinel-2 or Landsat 8. This would enhance the same area revisit monitoring (5 to 3 days) and help identify the impact of sea-level rise affecting the low-lying areas. However, the integration of both very high and moderate spatial resolution images would largely improve spatial and temporal shoreline monitoring.

Satellite-derived Bathymetry Mapping
The log-band ratio method was applied to extract shallow water bathymetry from the GeoEye-1 satellite image in conjunction with in-situ bathymetry data. The log-band ratio method, an empirical satellite-derived bathymetry approach, allows the extraction of the depth of the seabed from water reflectance and training in-situ data. Initial preprocessing steps involve radiometric calibration and atmospheric correction, using the FLAASH algorithm, and water and non-water pixels distinction, based on threshold NIR reflectance values. The satellite image was acquired in good conditions, with a calm sea and 0 cloud cover, and thus, there was no need for sunlight correction or cloud masking.
The application of the log-band ratio requires two spectral bands-blue and green-which allow the mapping of relatively shallow water depth, and later empirically calibrated to local datum based on a minimum number of in-situ training data. This in-situ bathymetry is used in part for the estimation of empirical calibration parameters and for accuracy assessment. The applied method is based on the model developed by Stumpf et al. [41], which allows the determination of water depth with multispectral images acquired over areas with different bottom sea types. This is the case for the pocket beach of San Vito Lo Capo and many other Sicilian pocket beaches where the bottom sea is largely occupied by seagrasses, such as Posidonia oceanica, sand, or rocks. The bathymetry extracted in this area and regression analysis conducted for shallow water less than 11 m of the depth produced a good correlation between the extracted satellite-derived bathymetry and in-situ bathymetry with correlation coefficient R 2 equals to 0.753. Since the method is applicable for multispectral images with at least blue, green, and NIR bands, it could be more useful to conduct shallow water bathymetry using 10 m spatial resolution freely and open access Copernicus Sentinel-2 satellite image and an available limited number of in-situ bathymetry for regular shallow water monitoring as well as for shallow water bathymetry updating.
The methods presented in this paper allow: (i) to conduct coastal erosion assessment based on the automatically extracted shoreline. It is a cost-effective tool as the number of spectral bands used is limited in addition to the reduced time for the processing. (ii) To conduct shallow water monitoring and shallow water bathymetry updating in areas where traditional survey techniques using vessels and ships could not operate properly. Such monitoring can be improved by integrating freely-available and open-access satellite images combined with a limited number of in-situ bathymetry. (iii) To explore available tools and technologies, which can be used to fight the effect of climate change and help decision-makers adopt data-driven policies intended for coastal areas management.

Conclusions
The current study constituted the first coastal erosion assessment conducted on the coast of San Vito Lo Capo using VHR images. In order to obtain information on the current and past dynamic conditions of this coastal area, two VHR images were analyzed to automatically extract the shoreline position in 2014 and 2019, as well as the extraction of shallow water depth on satellite images acquired in 2014.
The R2V approach adopted allowed extracting automatically the position of shoreline with only one input band. This showed that intelligence on coastline dynamic situations could be obtained on a very limited volume of data and also with small computation capability, leading to more sustained coastal erosion monitoring. For the period of 5 years, we noticed that approximately 17% of the sand beach was submerged, confirmed by 15 m inland shoreline migration, and 15,450 m 2 of the sand beach was lost. On the other hand, SDB mapping demonstrated the application of satellite images in shallow water bathymetry studies. The Stumpf et al. equation applied on the GeoEye-1 satellite image in conjunction with available in-situ data provided useful water depth information for the first 9 m. This depth was most challenging for the traditional bathymetry survey, and the methodology used could be adopted as a regular monitoring technique in vulnerable and rapidly evolving coastal areas.
Advances, both in algorithms and VHR images availability, improve the capability of monitoring vulnerable areas; we demonstrated that these images could be analyzed with simple methods, and obtain precision variable information in coastal areas. In addition, the ability to cover inaccessible, large, and high dynamic areas with very high spatial resolution images makes R2V and SDB technique alternative and cost-effective methods for bathymetric and coastal erosion studies.