High Resolution Satellite Images for Instantaneous Shoreline Extraction Using New Enhancement Algorithms

: Knowledge of a territory is an essential element in any future planning action and in appropriate territorial and environmental requalification action planning. The current large-scale availability of satellite data, thanks to very high resolution images, provides professional users in the environmental, urban planning, engineering, and territorial government sectors, in general, with large amounts of useful data with which to monitor the territory and cultural heritage. Italy is experiencing environmental emergencies, and coastal erosion is one of the greatest threats, not only to the Italian heritage and economy, but also to human life. The aim of this paper is to find a rapid way of identifying the instantaneous shoreline. This possibility could help government institutions such as regions, civil protection, etc., to analyze large areas of land quickly. The focus is on instantaneous shoreline extraction in Ortona (CH,


State-of-the-Art of Shoreline Detection
There is no single method for modeling the shoreline, although over the last ten years, great progress has been made in mapping technology, ranging from the development of new and increasingly-precise measurement techniques such as GNSS (Global Navigation Satellite System) to high and very high resolution satellite images and the development of particularly efficient devices for coastal topography assessment. The most common methods currently used for the identification of the shoreline include: 1. Traditional topographical surveys (leveling, surveys with total stations) and GPS [8][9][10]: These allow high precision measurements in individual points, without any possibility to repeat the survey; 2. Video systems [8,9]: The technique is used mainly with a network of fixed terrestrial cameras (e.g., ARGUS [11]) and Siren) from a few units up to tens, installed at prominent points in the landscape or on specially-positioned supports. They acquire at intervals of up to 10 hz, encompassing 180° views and allowing total coverage of about 4-6 km of beach. The acquired images are oblique and require orthorectification operations, as well as georeferencing with algorithms derived from classical photogrammetry; 3. Aerophotogrammetric/UAV survey [12][13][14][15][16][17][18][19][20]: This does not provide a detailed relief, but represents the entire study area at the time of acquisition. Achievable precisions are at the centimeter-subcentimeter level, but the costs are high. "Ad hoc" flights must be planned, and above all, it is significantly influenced by weather conditions (sunny and minimal wind conditions are needed); this leads to limitations on flight seasons; 4. Satellite remote sensing [21,22]: This refers to the latest-generation remote sensing satellites, which are becoming more effective as compared to those available previously. Traditionally, medium resolution satellites (e.g., Landsat and Sentinel-2) have been used advantageously for coastline studies that did not require very high accuracy. Satellites with high and very high resolution are required because they are particularly advantageous compared to traditional photogrammetric aerial acquisitions. High spatial resolution remote sensing satellites allow data to be acquired and processed more quickly, with comparable precision, while offering a very high level of detail. Furthermore, the fundamental ability to capture the scene in several spectral bands allows more information to be extracted than is extractable from images covering only the visible part of the electromagnetic spectrum and thus allows thematic maps of the territory to be created through multispectral classification. Another advantage, which makes this technique more attractive than the others, is that satellite revisit times are very short (a few days), which allows the area to be studied on images taken at different times.
From the literature, it is clear that the use of satellite images is the most diffused for the coastal environment by now, and this paper will deal with the study of innovative filters for the extraction of the shoreline. This work is a first step for a comparison between ACM filters (never used in the engineering fields) and filters well known in the world of the software used.

Case Study
The following experimentation concerns the shoreline in Ortona (CH, Italy) using WorldView-2 satellites. Common filters were initially used. Then, ACM filters were tested and compared with previous ones. The images were visually analyzed in order to find rapid ways to extract the instantaneous shoreline.

WorldView-2
WorldView-2 satellites acquire eight-band multispectral images with a resolution of 2 m to provide the best natural representation of the Earth from space, and panchromatic images with a resolution of 50 cm [23]. In particular, WorldView-2 is the first commercial high resolution satellite capable of acquiring, in addition to the four typical blue, green, red and near-infrared bands, four additional multispectral bands, previously unavailable on a very high resolution satellite: coastal, yellow, red edge, and Near Infrared 2 bands [23]. The coastal band is designed to study the shoreline, because of its ability to penetrate under the water's surface. WorldView-2 operates at an altitude of 770 km, and its advanced onboard system is capable of capturing images above the Earth [23]. These images provide unprecedented detail, enabling very accurate detection.

Data Pre-Processing
The "raw" images of Ortona's coastline, taken by the WorldView-2 satellite ( Figure 1) by DigitalGlobe on 29 June 2010, contain geometric and radiometric distortions, so specific software based on rigorous algorithms or based on rational polynomials must be used to correct them. The panchromatic and multispectral images are georeferenced to the WGS84 reference system (WGS84/UTM Zone 33 N, EPSG: 32633). The image correction steps are [21,24] For the Ortona images, 23 GCPs (Ground Control Points) were acquired during a previous research [21] by means of double-frequency GNSS receivers and processed in post-processing mode compared to the network of permanent stations of the Abruzzo region [6]. Eight of these were used as CPs (Check Points), to evaluate the accuracy of the measurements. To take into account land elevations, the basis for the relief displacements due to terrain elevations, DEM in software library was used. To correct geometric distortion of WorldView-2 satellite images, Rational Polynomial Coefficients (RPCs) were used. After that, we proceeded with a resampling, using the nearest neighbor technique. Then, the "HCS (Hyperspherical Color Space) Resolution Merge" technique was adopted [25] to perform pansharpening.
After these steps were complete, the images were ready to be processed using the so-called "enhancement techniques".

Data Processing with Common Filters
A series of filters (enhancement technique) [26] was tested on the image (in this experimentation, ERDAS IMAGINE 2015 software was used) to find the best ones that showed a clear and defined shoreline.

Decorrelation Stretch
The first step is PCA (Principal Component Analysis), which is used in statistics for the simplification of source data. It is often used as a data compression method: it allows redundant data to be compressed into fewer bands using linear transformation. Four bands were chosen (coastal, red, green, and blue bands) [24]. Then, the decorrelation stretch filter was applied. This allows the whole range of DN (Digital Numbers) varying from 0-255 to be used and, therefore, a greater amount of information to be extrapolated. The decorrelation stretch is a sort of "stretching" that is applied only to the chosen four bands (which are called "main components") and not to the original image. The potential of this technique is exploited to highlight the elements that are almost invisible to the naked eye, to improve and facilitate visualization. The resulting image enabled clear shoreline distinction due to the radiometric response of the elements on the Earth's surface ( Figure 2). In Figure 2, the darker color of the water indicates the presence of shallow waters, which correspond to artificial barriers (used to mitigate the erosive effect of the sea). The difference between wet and dry sand is revealed in the change of color from dark red to light red, respectively.

Normalized Difference Vegetation Index
The NDVIs [27] are mainly used for forests or related to agriculture applications [28][29][30], but they also give excellent results in shoreline definition. As is well known, the main objective of this index is to analyze the differences between the characteristic components of the vegetation's spectral signature. This denominator assumes the role of normalization factor, while the relevant characteristic of the index is the numerator. NDVI combines the red and the near-infrared (Equation (1)) pixel by pixel: According to this definition, NDVI ranges mathematically between −1 and +1, but, in practice, between −0.4 (water bodies) and +0.7 (very dense and broad-leaved vegetation). In general, still water has low reflectance. In addition, NDVI is able to compress the size of the data by a factor of two because it replaces two spectral bands with a single band. The result computed for the shoreline is shown in Figure 3. The shoreline between water and beach is visually well defined: the contrast between them is much accentuated by the abrupt transition between black and white, respectively.

WorldView Water Index
The WorldView Water Index (WV-WI) proved to be the most powerful algorithm. Above all, it is specific to images taken by WorldView satellites and, therefore, suitable for this case study. It is good for detecting water or shadows [24]. Traditionally, the blue band is normalized with the NIR2 band, because the radiometric answer is clearly visible. The WorldView Water Index considers coastal and NIR2 bands [24] (Equation (2)) because the difference between these bands, in terms of wavelengths, is greater (coastal band 400-450 nm, blue band 450-510 nm, NIR2 band 860-900 nm [27,31]), so it provides an even more discreet threshold to identify water rather than land or vegetation areas: The results of applying it to the WorldView-2 Images are shown in Figure 4: Water has a much higher value on the coastal blue band than on the NIR2 band: this is why it is white in color. The contrast between the water and the ground is very high. The differences between deeper and less deep water and even between dry and wet sand are not evident, but this algorithm is definitely the best of those tested from the point of view of shoreline definition because it shows a clear contrast between water and land.

Unsupervised and Supervised Classification
"Classification" refers to the grouping of pixels and similar values into groups (called "clusters") that identify a specific category within the data [4]. Pixels are divided into categories based on their spectral properties. The most delicate problem is defining the type and number of classes. As is well known, there are two types of classifications: unsupervised and supervised. K-means and ISODATA algorithms were used for the unsupervised classification [24], taking 30 clusters into account. In the supervised classification, 12 clusters were identified, and the maximum likelihood parametric method was used [4].
The shoreline derived from both methods is well defined ( Figure 5), but the main information that can be extrapolated in the results of unsupervised classification is the difference in depth in the various areas of water (from deep to less deep water, color degrades from blue/grey, pink, blue, orange, to light blue at the shore). The same can be said for the difference between dry and wet sand where the drier sand is green, while wet sand varies from light blue to red. Comparing the results of supervised classification into 12 classes with unsupervised classification into 30 classes, the supervised classification reveals less information about sea depth.

Data Processing with the ACM System
In this section, some experimental algorithms are tested on WorldView-2 images to determine their effectiveness for the case study. These algorithms have been developing by Professor Paolo Massimo Buscema, mathematician, Director of the Semeion Research Centre of Science of Communication of Rome and Full Professor Adjoint at the University of Colorado (USA).
This patent concerns Active Connections Matrix systems (ACM), according to which each image is considered as an active matrix (network) of connected elements (pixels) that develops over time. The main idea upon which this theory is based states that each digital image stores the maximum amount of information within the pixel values and their relationships. Furthermore, it is possible to obtain important information by analyzing the reciprocal positions occupied by pixels. For a complete presentation of ACM algorithms, see [32,33].
Any digital image is a matrix made of as many rows as the pixels that determine the width and as many columns as the pixel related to the height. Each pixel is identified by its coordinates i ∈ 1, ..., R, j ∈ 1, ..., C and

Analysis of Results
For the case study, software "ACM Batch" was used [34]. To test possible similarities or differences, the algorithms were tested on: • a four-band image (coastal, red, green, and blue); • an image containing only the coastal blue band (which is specific for coastal studies); • an image with the coastal blue band plus NIR2; • an image with only the NIR2 band.
These choices were carried out because these bands best respond to shoreline detection. In this paper, for the sake of brevity, the results are shown in Appendix B. Here, only the results obtained by applying the contractive maps on the NIR2 image are shown in Figure 7a, because it had similarities with two of the common filters: NDVI and WVWI (Figure 7b,c). In every image, the high reflectance between land and water shows a defined shoreline, which is metrically detected taking three profiles (diagrams of the pixel values as a function of a distance) over three sample distances. As a reference, a pansharpened image is considered (Figure 8). In this image, every color corresponds to a profile, and in particular, Profile 1 is blue, Profile 2 green, and Profile 3 red. These colors will be the same as the profiles plotted in Figure 9. The letter "A" refers to the shoreline extracted from the pansharpened image, used as a reference. It is also worth noting the peaks at the rocks (Lines B and C) because they show similar trends in almost all the graphics. Line D in Figure 9 refers to the shorelines extracted from, respectively, NDVI, WVWI, and CM over NIR2 images. Line D was determined in correspondence with the maximum slope on the jump for the CM over the NIR2 image and in correspondence with the middle point of the jump for the NDVI and WVWI. The differences between the reference image (pansharpened, Line A) and the other filter (Line D) was studied. It is shown in Table 2. For Profile Number 2 on the NDVI image ( Figure 9b) and Profile 3 on the WVWI image (Figure 9f), there was not a clear jump, which emphasizes the shoreline, as in the other cases. For this reason, Line D on these two graphs is not drawn. As shown in the table, the differences were the same for Profile 1, but were closer to the value taken from the pansharpened image in the case of Profiles 2 and 3 for the CM filter. Figure 9. Profiles over three sample distances on: NDVI (a-c), WVWI (d-f) and CM over NIR2 image (g-i). A, B and C, D are respectively referred to the shoreline extracted from the pansharpened image, used as reference, the artificial rocks and the shoreline extracted from NDVI, WVWI and CM algorithms.
A similar trend can be seen in the graphs. Peaks and hollows, which correspond to the passage from the mainland to the sea (and vice versa), are displayed at the same distance in all the graphs. In other words, in Profile 1, the first peak/hollow can be seen around 70 m and the second one around 150 m. This behavior is the same in all three for Profile 1 (Figure 9a,d,g). The same can be said for Profiles 2 (Figure 9b,e,h) and 3 (Figure 9c,f,i). This means that the algorithms worked in a very similar way. The original elaborations are shown in Appendix B.

Discussion of Results and Conclusions
In the present work, the potential of the ACM models for the automatic identification of the coastline from very high resolution multispectral satellite images has been tested. The interest for the definition of a certain and robust procedure in this regard derives from the increase in the phenomenon of coastal erosion that affects low coasts, as well as the Italian ones. The ACM models already successfully applied in other (mainly health) applications have been adapted and tested here for the first time on satellite images.
To verify the accuracy and robustness of this approach, the results of the experiments on a sample area were compared with those obtained on the same area using all the algorithms and filters most commonly used for these same applications [35].
By combining the various types of images with the various models, it was possible to obtain a wide range of tests in order to validate the model exhaustively on the specific test site. For space reasons, the complete results of all the tests have been reported in the Appendix B. From the above-mentioned experiments, we can give some general considerations on the ACM algorithms: • their ability to define the contours of real images was excellent; • their adaptation to satellite imagery made them a powerful tool for digital image processing, and they can be used for automatic territory analysis; • they were capable of distinguishing noise and salient forms within the image.
More in detail, we can say that by comparing the profiles obtained automatically by the various algorithms, the ACM models always showed an accuracy greater than or equal to that of the most common methodologies. Moreover, the measured deviation from what can be observed on the screen was always the result of the order of the pixel and often subpixel, making the method even more accurate than the much more expensive manual acquisition. This is potentially very useful for future change detection between two images taken in different time periods, generating a concrete possibility that an appropriate intervention plan can be developed, indispensable if government institutions, such as civil protection agencies, need to carry out rapid monitoring.

Author Contributions:
The authors contributed equally to this work.
Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: In the following, the most useful algorithms for this case study will be illustrated.

Appendix A.1. Automata Rule
The Automata Rule (AR) [32] serves to initialize the values of the connections in the ACM with fixed connections (first order of complexity). This algorithm defines the connections between each pixel-unit as a non-linear transformation of their different brightnesses: where: c The The resulting image will highlight the edges of the pictures in the image. In fact, the more pixels are similar, the more the connection values will look like one (white), so the mean rate will be high. If the pixels have different values, the weights will be closer to −1 (black), as will the average.

Appendix A.2. New IAC
As with all the ACM algorithms, New Interactive Activation and Competition [32] (New IAC) considers each pixel as an agent connected to the others in its neighborhood. The connections linking uij to all the uxy ∈ Iuij are initialized by means of the automata rule, then they remain fixed until the end of the process; therefore, New IAC belongs to the first order family. Each uij receives excitations and inhibitions from the surroundings and changes its state accordingly at each iteration. Equations (A6) and (A7) allow the computation of the excitation and inhibition values released by the neighborhood of the pixel uij at the iteration n. Equation (A8) determines the input arriving at it as a weighed sum of those quantities. The final value of the activation of the pixel uij is then computed by means of Equation (A9): [n]  [32] (New CS) is an ACM system with fixed connections based on the constraint satisfaction network model [1]. The main idea behind New CS is to consider each pixel-agent of the matrix as a "hypothesis", which may be true or false in a fuzzy way, in relation to the brightness of the pixel itself. New CS tries to make every hypothesis of the matrix "true" by maximizing the brightness of each pixel-agent, having as a constraint the value of fixed connections previously calculated with AR (Equations (A1)-(A4)). The maximization problem to solve is then:  Table A1 shows some of them.

Pij
At each iteration, the values of pixels of the image at the previous step are replaced by those computed by means of one of the methods of Table A1. Then, the values of connections are updated considering the original value of pixels. The process ends when the dynamic system reaches its attractor and