Multi-Scale LG-Mod Analysis for a More Reliable SAR Sea Surface Wind Directions Retrieval

: An improved version of the Local-Gradient-Modiﬁed (LG-Mod) algorithm for Sea Surface Wind (SSW) directions retrieval by means of Synthetic Aperture Radar (SAR) images is presented. A “local” multi-scale analysis of wind-aligned SAR patterns is introduced to improve the LG-Mod sensitivity to SAR backscattering modulations, occurring locally with various spatial wavelengths. The Marginal Error parameter is redeﬁned, and the adoption of the Directional Accuracy Maximization Criterion (DAMC) allows for the novel Multi-Scale (MS) LG-Mod to automatically select the local processing scale that may be regarded as optimal for pattern enhancement, once a discrete set of scales has been already ﬁxed. Hence, this optimal scale successfully gives evidence to guarantee the best achievable local direction estimation. The assessment of the MS LG-Mod is carried on both simulated SAR images and a Sentinel-1 (S-1) dataset, consisting of 350 Interferometric Wide Swath Ground Range Multi-Look Detected High-Resolution images, which cover the region of the Gulf of Maine. In the latter case, the removal of artifacts and non-wind features from SAR amplitudes is mandatory before directional estimations. In situ wind observations gathered by the National Oceanic and Atmospheric Administration National Data Buoy Center (NOAA NDBC) are exploited for validation. The ﬁndings obtained from S-1 data conﬁrm the ones from simulated patterns. The MS LG-Mod analysis performs better than each single-scale one in terms of both percentages of reliable directions and directional Root Mean Square Error (RMSE) values achieved.


Introduction
Sea Surface Wind (SSW) is a very important parameter for studies in many fields of application, such as marine meteorology and oceanography [1], coastal ecosystems monitoring and conservation [2], and wind energy mapping [3].
Among available tools, Synthetic Aperture Radar (SAR) systems are commonly used in most meteorological conditions, and they can provide wind information at high spatial resolution [4][5][6][7][8]. In spite of many methods recently proposed, the approach that exploits the dependency of SAR Normalized Radar Cross-Section (NRCS) on wind parameters (wind speed and direction), and the geometry of acquisitions (radar incidence angle) via empirical Geophysical Model Function (GMF) is still widely used [9][10][11][12][13].
One of the inconveniences of this approach is the requirement of wind direction estimate as an input parameter. In fact, assuming wind direction as a priori information, wind speed can be inferred through the inversion of a GMF. This can lead to errors in SAR wind speed estimation mainly due to those related to wind direction used.
Recently, both the United States National Oceanic and Atmospheric Administration (NOAA) and the European Space Agency (ESA) have implemented an operational service based on the scatterometry-based approach, providing wind field from SAR data exploiting wind direction from Numerical Weather Prediction (NWP) models as input [14,15].
In Reference [31] the use of the Non Subsampled Contourlet Transform (NSCT) algorithm for wind direction estimation was proposed. Authors stated that NSCT transform is suitable for the estimation of wind direction from ocean SAR images, which is a typical multi-scale texture analysis problem.
The Local Gradient (LG) method [22] is still one of the most used techniques for wind direction retrieval since it may estimate local wind direction with smaller resolution cells. For each considered scale, one histogram of local directions is built up and the main direction is evaluated. In Reference [32], different option were tried to choose the optimal among estimated directions: (1) to take the direction closest to an available wind direction provided by a NWP model and (2) to take by visual inspection the direction closest to wind-aligned signatures detected on the SAR image.
To the best of our knowledge, the literature that reports the multi-scale analysis of wind rows suffers from several limitations: (1) some of the methods proposed [28,30] implement a multi-scale analysis that has been designed with the aim of enhancing the modulations belonging to the BLRs and WSs, without any effective selection of the optimal scale; (2) one method [32] makes use of a priori external wind information to select the optimal scale among a number of tested ones, thus making the method not fully automatic; (3) the only method [27] for which a selection of the optimal scale is carried out does not present effective comparisons with single-scales performance.
In the present work, a criterion for the optimal scale automatic selection is introduced to overcome all of these issues, supported by a comparison with single-scales performance as well.
In a previous study, Rana et al. [17] proposed the LG-Mod algorithm and investigated the usefulness of the Single-Scale Marginal Error parameter (SS ME ROI α ), in providing an "optimal" SAR processing pixel-size selection for the matching with wind rows modulation wavelengths. This optimal processing scale was identified, among a set of scales, as the single one which can globally yield in a SAR image the largest percentage value of reliable ROIs, i.e., the ROIs characterized by a SS ME ROI α value less than or equal to an established threshold of acceptance ME TH . Nevertheless, it is worth noting that the optimization adopted was not the best. In fact, this first attempt of multi-scale analysis relied on a global approach. The algorithm was limited to choose only one unique single-scale, among an employed set, for the processing of the whole SAR image. It meant, one scale for one SAR image.
However, considering that a full-frame SAR image may cover hundreds of kilometers (for azimuth and range dimensions) and that the actual wind regime conditions may be very different in wide regions, the limitation in Reference [17] was to be overcome. Indeed, since wind rows wavelengths may vary widely depending on several physical factors, such patterns should be considered intrinsically local and locally analyzed as well.
As consequence, this work introduces a newly defined Multi-Scale Marginal Error parameter, i.e., MS ME ROI α , as an indicator for a local choice of the optimal processing scale through a Directional Accuracy Maximization Criterion (DAMC). Differently from the approach in Reference [17], the novel multi-scale analysis is carried out locally and it really allows the processing of a SAR image using multiple scales. This local approach results in the best wind rows enhancement and wind directions estimation performance.
The novel technique is experimented on a valuable Sentinel-1 dataset and exploiting in situ data measured by anemometers on board stations belonging to the NOAA network. Hence, this study deals with the full validation of the technique on real SAR data, also verifying first results obtained on simulated SAR patterns, as presented in Reference [8].
The peculiarity of the "local" MS LG-Mod also allows the method to establish possible spatial correlations between the chosen optimal scale and environmental parameters, such as wind speed and direction, bathymetry and distance from the coast.
The paper is organized as follows: Section 2 describes the procedure used to produce simulated SAR images, the characteristics of the Sentinel-1 dataset, and the available NOAA in situ wind measurements; Section 2 also presents the methodology adopted for extracting Remote Sens. 2021, 13, 410 4 of 24 wind direction from SAR data and for directional comparison criteria; Section 3 reports results and Section 4 provides their discussion; and Section 5 summarizes the conclusions.

Simulated SAR Images
Both linear and circular patterns were generated to simulate SAR backscattering modulations observable in real cases of wind-forced patterns visibility [8]. These patterns were characterized by spatial wavelengths in the range (500 m-2 km) for wind rows detection. Typical values of SAR amplitudes over the sea were assumed. A speckle noise, modeled with a Rayleigh distribution, was also taken into account for SAR amplitudes. Simulated images were considered of a 3000 by 3000 (azimuth by range) dimension, simulating a pixel size of 10 m by 10 m.
Patterns obtained are shown in Figure 1a-d, both without and with the multiplicative speckle noise, respectively.
Remote Sens. 2021, 13, x FOR PEER REVIEW 4 of 24 extracting wind direction from SAR data and for directional comparison criteria; Section 3 reports results and Section 4 provides their discussion; and Section 5 summarizes the conclusions.

Simulated SAR images
Both linear and circular patterns were generated to simulate SAR backscattering modulations observable in real cases of wind-forced patterns visibility [8]. These patterns were characterized by spatial wavelengths in the range (500 m-2 km) for wind rows detection. Typical values of SAR amplitudes over the sea were assumed. A speckle noise, modeled with a Rayleigh distribution, was also taken into account for SAR amplitudes. Simulated images were considered of a 3000 by 3000 (azimuth by range) dimension, simulating a pixel size of 10 m by 10 m.
Patterns obtained are shown in Figure 1a-d, both without and with the multiplicative speckle noise, respectively.

Sentinel-1 SAR images
To evaluate wind directions estimation through the LG-Mod algorithm by means of real SAR data, a Sentinel-1 dataset of C-band VV-polarized images was built up. Both ascending and descending images, acquired by both Sentinel-1 (S-1) A and B satellites, were downloaded from the Copernicus Open Access Hub. However, only S-1 Interferometric Wide Swath Ground Range Multi-Look Detected High-Resolution (henceforth, IW) images were chosen, since SS LG-Mod directional performance analysis already indicated that IW images (with 10 m by 10 m resolution) are the optimal choice for high-resolution SSW retrieval [17].

Sentinel-1 SAR Images
To evaluate wind directions estimation through the LG-Mod algorithm by means of real SAR data, a Sentinel-1 dataset of C-band VV-polarized images was built up. Both ascending and descending images, acquired by both Sentinel-1 (S-1) A and B satellites, were downloaded from the Copernicus Open Access Hub. However, only S-1 Interferometric Wide Swath Ground Range Multi-Look Detected High-Resolution (henceforth, IW) images were chosen, since SS LG-Mod directional performance analysis already indicated that IW images (with 10 m by 10 m resolution) are the optimal choice for high-resolution SSW retrieval [17].

In Situ Data
Wind observations gathered by the in situ stations of the National Oceanic and Atmospheric Administration National Data Buoy Center (NOAA NDBC) were exploited to assess the performance of the MS LG-Mod processing. Wind measurements were provided at the reference height of 10 m above the mean sea level, with both different time intervals between acquisitions and directional accuracies depending on the stations equipment, as reported in Table 1. This table reports also the geographic location and further information of the 17 stations matching in time and space with the Sentinel-1 acquisitions. Deployment of the NOAA NDBC stations is shown in Figure 2.  . Sentinel-1 (S-1) full-frame coverage is in orange and cyan for descending (i.e., S1B-10391-40-1/2) and ascending (i.e., S1A-21484-62-1/2/5/6, S1B-10515-164-6, S1A-6740-18-1, S1A-15417-120-5, and S1A-7513-91-1) S-1 Interferometric Wide Swath Ground Range Multi-Look Detected High-Resolution (IW) images, respectively. Background image (Data SIO, NOAA, US Navy, NGA, GEBCO; Image Landsat/Copernicus) is from Google Earth.

In Situ Data
Wind observations gathered by the in situ stations of the National Oceanic and Atmospheric Administration National Data Buoy Center (NOAA NDBC) were exploited to assess the performance of the MS LG-Mod processing. Wind measurements were provided at the reference height of 10 m above the mean sea level, with both different time intervals between acquisitions and directional accuracies depending on the stations equipment, as reported in Table 1. This table reports also the geographic location and further information of the 17 stations matching in time and space with the Sentinel-1 acquisitions. Deployment of the NOAA NDBC stations is shown in Figure 2.

Multi-Scale LG-Mod Processing
The MS LG-Mod processing scheme is reported in Figure 3. The input to the algorithm consists in a single SAR image, which could be either a simulated one or a co-polarized (VV or HH) Sentinel-1 Normalized Radar Cross-Section (NRCS) image. The latter must be accompanied by auxiliary products (i.e., land-sea mask and, latitude and longitude maps), e.g., obtained using the Sentinel Application Platform (SNAP) freely available from the European Space Agency (ESA). In particular, land-sea mask is exploited to mask out the land from processing. Latitude and longitude maps are used instead to finally provide geographical location for output wind estimations, since the processing is completely carried out in SAR geometry (i.e., on azimuth by range images).

Multi-Scale LG-Mod Processing
The MS LG-Mod processing scheme is reported in Figure 3. The input to the algorithm consists in a single SAR image, which could be either a simulated one or a co-polarized (VV or HH) Sentinel-1 Normalized Radar Cross-Section (NRCS) image. The latter must be accompanied by auxiliary products (i.e., land-sea mask and, latitude and longitude maps), e.g., obtained using the Sentinel Application Platform (SNAP) freely available from the European Space Agency (ESA). In particular, land-sea mask is exploited to mask out the land from processing. Latitude and longitude maps are used instead to finally provide geographical location for output wind estimations, since the processing is completely carried out in SAR geometry (i.e., on azimuth by range images). The main processing steps of the MS LG-Mod algorithm are reported in the following: • Image Smoothing and Decimation [8,17]: This step consists in a series of smoothing and sub-sampling operations with the aim at both enhancing patterns of interest and reducing the speckle noise of the input SAR image. Smoothing operations are edge preserving in order to preserve as much as possible the directional information of the wind-induced SAR signatures. Sub-sampling (with factor 2) reduces by half both the image dimensions as well as it doubles the (final) scale in both dimensions themselves.
To capture SAR wind rows characterized by wavelengths in the range between 500 m and 2 km, slightly smaller than the one chosen in Reference [30], scales are profitably selected spanning the range from about 100 to 400 m. Note that the use of higher scales would introduce too much smoothing, thus leading to an unpleasant directional information loss, and a large reduction on the number of samples for directional estimations as well. Lower scales instead would cause noisy estimations, although with a large number of samples. As consequence, for the available Sentinel-1 dataset, S = 3 processing scales, i.e., the High Scale (320 m by 320 m), the Medium Scale (160 m by 160 m), and the Low Scale (80 m by 80 m), were selected for wind rows detection. However, this step may be generally repeated with a different number S of scales, i.e., ∀k ∈ {1, . . . , S}.

•
Local Gradients Computation [17,22]: This computation is performed for each previously resized image (one image for each k-th scale applied). Local Gradients (LGs) that correspond to "unusable points" are discarded in the next step. • Unusable Points Masking (UPM): At this step, land pixels from the available landsea mask and those affected by unwanted border effects as a result of the filtering performed at the first step are masked. Furthermore, the masking of those points characterized by extremely high or extremely low Local Gradients is carried out. A thresholding is aimed at excluding most of the sigma nought modulations and LGs that are not related to wind rows. To accomplish this task the condition of 2 × 10 −4 < Amplitude(LG) < 2 × 10 −3 (with experimentally determined thresholds) is required for each pixel within each sub-image (or ROI) to be evaluated in the subsequent steps.

•
Main Wind Directions Estimation: Previous steps are repeated for a number S of selected scales, i.e., ∀k ∈ {1, . . . , S}. Each LG direction image obtained from each of such repetition is divided into ROIs. Those ROIs that are affected by a large percentage of unusable points detected (e.g., more than 30%) are discarded and related estimations are not provided. Otherwise, for each k-th LG direction image as well as for each not-discarded ROI, a confidence interval is provided with a user-defined confidence level (1 − α), as follows: where u α is the upper 1 2 α quintile of the standard normal distribution, and α ROI,k 2 represents the 2nd central trigonometric moment of the doubled LG directions, as expressed by (4). Ranging from 0 to 1, as calculated by (5), R ROI,k is the so-called Mean Resultant Length (MRL), and it represents a measure of the alignment of the LG directions.
Equations (2) and (3) allow the algorithm to estimate, respectively, the mean local direction β ROI,k and the related inverse measure of the directional accuracy ME ROI,k α directly from the set of the LG directions (within each ROI and for each k-th scale). These equations were derived by adapting Directional Statistics findings from circular [33] to axial data [34]. In fact, when quasi-linear patterns occur on SAR amplitude, the distribution of LG directions may be assumed bimodal with an angular distance of π (rad) between the two modes. In other words, these local directions used for estimations should be intrinsically considered as axial data instead of circular data. Then, for a large number of useful points, the estimated direction β ROI,k is given as π-ambiguous, while the parameter ME ROI,k α ranges from 0 to π/4 radians.
Once all single-scale evaluations have been completed, the Directional Accuracy Maximization Criterion (DAMC) is applied. This criterion states that the best directional estimate for each ROI, i.e., β ROI , is provided by the single-scale local directional estimate β ROI,k which corresponds to the minimum value among the all Single-Scale Marginal Errors, ME ROI,k α . The latter value is defined as the Multi-Scale Marginal Error, ME ROI α : The DAMC, which defines the core of the novel "local" multi-scale analysis, consequently allows the definition of a "new" confidence interval for each local direction estimate (and ROI). This interval thus represents the best (smallest) confidence interval among the all single-scale-derived ones: The orthogonal to the best local direction β ROI previously obtained finally provides the π-ambiguous wind direction for each ROI. The related directional accuracy is still supplied by the parameter ME ROI α . • Reliable Wind Directions Selection: This final step is devoted to discard those directional estimates considered not reliable enough. As done for SS LG-Mod outcomes, a suitable threshold of acceptance, i.e., a maximum Marginal Error value ME TH , is defined by the user and then applied to all MS LG-Mod estimated directions: Although derived by the application of the same threshold ME TH , the latter population of directions ROI TH Multi is obviously different compared to each SS one, which is defined as follows:

ROI TH
Single,k = ROI ME ROI,k α ≤ ME TH , ∀k ∈ {1, . . . , S} Remote Sens. 2021, 13, 410 9 of 24 The reliable directions obtained through the above-described algorithm are finally compared with either the theoretical local directions of the simulated patterns or, in the case of S-1 data exploitation, the directional measurements from NOAA in situ stations. In the latter case, the LG-Mod wind directions are derived from both 5 km by 5 km and 12.5 km by 12.5 km ROIs, each one centered on the location of a NOAA station. In situ wind directions are temporally interpolated at SAR acquisition times to obtain co-location in time with LG-Mod wind directions.
The comparisons are conducted for both Single-and Multi-Scale LG-Mod results in order to validate improvements derived from the MS analysis with respect to each SS processing.

Results
The MS and SS LG-Mod processing were carried out performing directional estimations on both simulated and real (i.e., Sentinel-1) SAR images.

Results
The MS and SS LG-Mod processing were carried out performing directional estimations on both simulated and real (i.e., Sentinel-1) SAR images. Figure 4a,b plots the Marginal Errors, evaluated with a 95% confidence level for all the LG-Mod estimated directions, as functions of the concentration parameter MRL [8]. The SS ME , , with k ∈ 1, … , S (where S = 3 is the number of single-scales chosen), and the MS values refer to noisy linear and circular patterns shown in Figure 1c,d, respectively. By definition of the DAMC criterion (6) (Section 2), represents the minimum value in the set ME , |k = 1, … , S , it means the value which maximizes the directional accuracy of each ROI directional estimate.    Figure 5a,b shows that the High Scale (320 m by 320 m, blue) is first chosen, and then the Medium Scale (160 m by 160 m, red) and the Low Scale (80 m by 80 m, yellow), while the spatial wavelength of patterns is decreasing from 2 km to 500 m. Figure  5c-g shows, instead, directions retrieved through each single-scale processing from linear and circular patterns, respectively.  Figure 5 reports, for a visual comparison, multi-and Single-Scale LG-Mod directions obtained applying the same final threshold (i.e., ME TH = 10 • , with α = 0.05) for both linear and circular patterns. In particular, the distribution of colors (one for each single-scale employed) in Figure 5a,b shows that the High Scale (320 m by 320 m, blue) is first chosen, and then the Medium Scale (160 m by 160 m, red) and the Low Scale (80 m by 80 m, yellow), while the spatial wavelength of patterns is decreasing from 2 km to 500 m. Figure 5c-g shows, instead, directions retrieved through each single-scale processing from linear and circular patterns, respectively. obtained applying the same final threshold (i.e., = 10°, with = 0.05) for both linear and circular patterns. In particular, the distribution of colors (one for each single-scale employed) in Figure 5a,b shows that the High Scale (320 m by 320 m, blue) is first chosen, and then the Medium Scale (160 m by 160 m, red) and the Low Scale (80 m by 80 m, yellow), while the spatial wavelength of patterns is decreasing from 2 km to 500 m. Figure  5c-g shows, instead, directions retrieved through each single-scale processing from linear and circular patterns, respectively. The percentage of reliable directions obtained by the MS processing (i.e., P , with reference to (w.r.t.) its own population ROI ) is, by definition, greater than the one derived from each SS processing (i.e., P , w.r.t. the population ROI , ), as shown in Figure 6a,b for linear and circular patterns, respectively. Each of these populations was derived through the application of the same but varying final threshold . The percentage of reliable directions obtained by the MS processing (i.e., P TH Multi , with reference to (w.r.t.) its own population ROI TH Multi ) is, by definition, greater than the one derived from each SS processing (i.e., P TH Single,k w.r.t. the population ROI TH Single,k ), as shown in Figure 6a,b for linear and circular patterns, respectively. Each of these populations was derived through the application of the same but varying final threshold ME TH . Obviously, the percentage of reliable directions from both MS and SS processing reduces while decreasing ME TH values. To evaluate the capability of the MS analysis in the selection of the optimal local processing scale to achieve the best directional estimation, the difference of RMSE values obtained for the MS and each SS processing (i.e., RMSE − RMSE , ) was computed considering the multi-scale population of directions ROI . Theoretical local directions of the simulated SAR images were used for comparison with the ones estimated by the MS and SS LG-Mod processing, and thus for RMSE difference evaluations. As shown in Figure 7, for both the abovementioned test patterns, the RMSE differences (y-axes) are always negative for each applied threshold (x-axes). It means that the MS estimation performance is always better than the one from each SS processing, especially for higher threshold values. The performance of the MS is also better with respect to the SS ones even if the comparisons are made on each single-scale population of directions (i.e., ROI , ), as illustrated in Figure 8. In particular, the performance of both MS and SS processing becomes progressively similar in terms of achieved RMSE values while decreasing threshold val- To evaluate the capability of the MS analysis in the selection of the optimal local processing scale to achieve the best directional estimation, the difference of RMSE values obtained for the MS and each SS processing (i.e., RMSE Multi − RMSE Single,k ) was computed considering the multi-scale population of directions ROI TH Multi . Theoretical local directions of the simulated SAR images were used for comparison with the ones estimated by the MS and SS LG-Mod processing, and thus for RMSE difference evaluations. As shown in Figure 7, for both the abovementioned test patterns, the RMSE differences (y-axes) are always negative for each applied threshold (x-axes). It means that the MS estimation performance is always better than the one from each SS processing, especially for higher threshold values. To evaluate the capability of the MS analysis in the selection of the optimal local processing scale to achieve the best directional estimation, the difference of RMSE values obtained for the MS and each SS processing (i.e., RMSE − RMSE , ) was computed considering the multi-scale population of directions ROI . Theoretical local directions of the simulated SAR images were used for comparison with the ones estimated by the MS and SS LG-Mod processing, and thus for RMSE difference evaluations. As shown in Figure 7, for both the abovementioned test patterns, the RMSE differences (y-axes) are always negative for each applied threshold (x-axes). It means that the MS estimation performance is always better than the one from each SS processing, especially for higher threshold values. The performance of the MS is also better with respect to the SS ones even if the comparisons are made on each single-scale population of directions (i.e., ROI , ), as illustrated in Figure 8. In particular, the performance of both MS and SS processing becomes progressively similar in terms of achieved RMSE values while decreasing threshold val- The performance of the MS is also better with respect to the SS ones even if the comparisons are made on each single-scale population of directions (i.e., ROI TH Single,k ), as illustrated in Figure 8. In particular, the performance of both MS and SS processing becomes progressively similar in terms of achieved RMSE values while decreasing threshold values. As a consequence, a threshold value ME TH can be fixed more or less strict, and thus the MS processing can be less or more advantageous, depending on the level of noise that affects the patterns signal.

LG-Mod Results from Real SAR Images
The aim of this study is to verify whether the DAMC can be an effective criterion for selecting the optimal scale for wind rows detections with the aim to obtain better performances from the MS LG-Mod application rather than the adoption of a single-scale approach. For this reason, the entire S-1 ROIs dataset considered (i.e., 1070 ROIs) was filtered out of those cases not attributable to the presence of wind rows. Specifically, the following actions were adopted: 1. Exclusion of ROIs corresponding to in situ wind speeds less than 2 m/s. This threshold was chosen considering that a weaker wind cannot produce roughness and therefore the sea area appears as flat [35]. As a result, the available dataset reduced to a number of 977 ROIs. 2. Exclusion of ROIs that present a coverage greater than 30% of unusable pixels, masked by the application of the UPM module. The latter allows the masking of bright pixels/areas (Figure 9a), radiofrequency interference (RFI) zones [36] ( Figure  9b) and dark areas (Figure 9c) through the upper and lower bounding of the LG amplitude. In fact, the edge areas and the inner parts of such SAR features cause, respectively, strong and weak LGs, which cannot be associated to wind rows. Land pixels not correctly identified in the land-sea mask are masked by the UPM module as well. After this step, the number of ROIs decreased to 448. 3. Further non-wind features [37,38], such as atmospheric gravity waves (AGWs) (Figure 9d) and atmospheric convective cells (Figure 9e

LG-Mod Results from Real SAR Images
The aim of this study is to verify whether the DAMC can be an effective criterion for selecting the optimal scale for wind rows detections with the aim to obtain better performances from the MS LG-Mod application rather than the adoption of a single-scale approach. For this reason, the entire S-1 ROIs dataset considered (i.e., 1070 ROIs) was filtered out of those cases not attributable to the presence of wind rows. Specifically, the following actions were adopted:

1.
Exclusion of ROIs corresponding to in situ wind speeds less than 2 m/s. This threshold was chosen considering that a weaker wind cannot produce roughness and therefore the sea area appears as flat [35]. As a result, the available dataset reduced to a number of 977 ROIs.

2.
Exclusion of ROIs that present a coverage greater than 30% of unusable pixels, masked by the application of the UPM module. The latter allows the masking of bright pixels/areas (Figure 9a), radiofrequency interference (RFI) zones [36] (Figure 9b) and dark areas (Figure 9c) through the upper and lower bounding of the LG amplitude. In fact, the edge areas and the inner parts of such SAR features cause, respectively, strong and weak LGs, which cannot be associated to wind rows. Land pixels not correctly identified in the land-sea mask are masked by the UPM module as well. After this step, the number of ROIs decreased to 448.

3.
Further non-wind features [37,38], such as atmospheric gravity waves (AGWs) (Figure 9d) and atmospheric convective cells (Figure 9e) cannot be removed by applying the current UPM module. To further filter the dataset, the identification and the discarding of those ROIs showing such features was carried out through visual inspection. This last step reduced the number of available ROIs to 239.
tively, strong and weak LGs, which cannot be associated to wind rows. Land pixels not correctly identified in the land-sea mask are masked by the UPM module as well. After this step, the number of ROIs decreased to 448. 3. Further non-wind features [37,38], such as atmospheric gravity waves (AGWs) (Figure 9d) and atmospheric convective cells (Figure 9e) cannot be removed by applying the current UPM module. To further filter the dataset, the identification and the discarding of those ROIs showing such features was carried out through visual inspection. This last step reduced the number of available ROIs to 239.

Exploitation of Aimed at Selection of the Optimal Local Processing Scale
The assessment of the MS LG-Mod at choosing the optimal local processing scale to reach the best directional estimation was evaluated on a dataset of 239 Sentinel-1 5 km by 5 km ROIs, as previously described. Figure 10 reports the numbers of reliable directions (y-axis), respectively, as function of the applied threshold (x-axis). Moreover, for the Sentinel-1 ROIs dataset, the number of reliable directions obtained by the MS processing (i.e., N w.r.t. population ROI ) was greater than the one derived from each SS processing (i.e., N , w.r.t. population ROI , ). The MS processing provided a percentage improvement increasing from 3.48% to 23.74% w.r.t. the Low Scale, from 4.85% to 33.33% w.r.t. the Medium Scale, and from 6.73% to 129.33% w.r.t. the High Scale, while decreasing the threshold from 44.999 ° to 7.5 °.  ) and for each threshold value . The MS performance was proven to be better than the one from each SS processing, especially for higher threshold values , also exploiting real Sentinel-1 SAR data. From the RMSE values plotted in Figure 11a

Aimed at Selection of the Optimal Local Processing Scale
The assessment of the MS LG-Mod at choosing the optimal local processing scale to reach the best directional estimation was evaluated on a dataset of 239 Sentinel-1 5 km by 5 km ROIs, as previously described. Figure 10 reports the numbers of reliable directions (y-axis), respectively, as function of the applied threshold (x-axis). Moreover, for the Sentinel-1 ROIs dataset, the number of reliable directions obtained by the MS processing (i.e., N TH Multi w.r.t. population ROI TH Multi ) was greater than the one derived from each SS processing (i.e., N TH Single,k w.r.t. population ROI TH Single,k ). The MS processing provided a percentage improvement increasing from 3.48% to 23.74% w.r.t. the Low Scale, from 4.85% to 33.33% w.r.t. the Medium Scale, and from 6.73% to 129.33% w.r.t. the High Scale, while decreasing the threshold from 44.999 • to 7.5 • .

Exploitation of Aimed at Selection of the Optimal Local Processing Scale
The assessment of the MS LG-Mod at choosing the optimal local processing scale to reach the best directional estimation was evaluated on a dataset of 239 Sentinel-1 5 km by 5 km ROIs, as previously described. Figure 10 reports the numbers of reliable directions (y-axis), respectively, as function of the applied threshold (x-axis). Moreover, for the Sentinel-1 ROIs dataset, the number of reliable directions obtained by the MS processing (i.e., N w.r.t. population ROI ) was greater than the one derived from each SS processing (i.e., N , w.r.t. population ROI , ). The MS processing provided a percentage improvement increasing from 3.48% to 23.74% w.r.t. the Low Scale, from 4.85% to 33.33% w.r.t. the Medium Scale, and from 6.73% to 129.33% w.r.t. the High Scale, while decreasing the threshold from 44.999 ° to 7.5 °.  ) and for each threshold value . The MS performance was proven to be better than the one from each SS processing, especially for higher threshold values , Figure 10. Number of reliable directions obtained by MS and SS processing (N TH Multi and N TH Single,k , respectively), evaluated on each own population ( ROI TH Multi and ROI TH Single,k , respectively), plotted as function of the applied threshold ME TH . Dataset: Sentinel-1 5 km × 5 km ROIs selected. Reference wind directions from NOAA NDBC. Figure 11a plots the RMSE values obtained from the MS and each SS processing (i.e., RMSE Multi and RMSE Single,k ) evaluated on the MS population of directions (i.e., ROI TH Multi ) and for each threshold value ME TH . The MS performance was proven to be better than the one from each SS processing, especially for higher threshold values ME TH , also exploiting real Sentinel-1 SAR data. From the RMSE values plotted in Figure 11a, the corresponding difference (i.e., RMSE Multi − RMSE Single,k ), evaluated on the same MS population of directions (i.e., ROI TH Multi ), is obtained and reported in Figure 11b. It can be noticed that, in general, the effect of the application of the threshold is that to reduce gradually RMSE values obtained, especially for lower than 15°. This is true for both MS and SS processing (Figure 11a). Applying the same , the MS procedure performs better than each SS one with an improvement of RMSE ranging from    It can be noticed that, in general, the effect of the application of the threshold ME TH is that to reduce gradually RMSE values obtained, especially for ME TH lower than 15 • . This is true for both MS and SS processing (Figure 11a). Applying the same ME TH , the MS procedure performs better than each SS one with an improvement of RMSE ranging from −2. It can be noticed that, in general, the effect of the application of the threshold is that to reduce gradually RMSE values obtained, especially for lower than 15°. This is true for both MS and SS processing (Figure 11a). Applying the same , the MS procedure performs better than each SS one with an improvement of RMSE ranging from     It is important to mention that MS wind direction improvement also yields to better wind speed retrieval through a GMF inversion. In fact, applying the C-band MODel 5.N (CMOD5.N) backscattering model [10], it was found that wind speed RMSE values are reduced up to -0.3 m/s using the wind direction derived from the MS processing rather than the one from each SS directional estimation.
Furthermore, previous analyses were also carried out on a ROI size of 12.5 km by 12.5 km. However, the number of available ROIs reduces up to 203. This reduction depends on the greater number of larger ROIs containing non-wind features and that are then discarded following the above-described procedure (Section 3.2). Moreover, in the case of NOAA stations located near the S-1 frames boundaries, it was not always allowed to extract 12.5 km by 12.5 km ROIs entirely included within the frames.
As reported in Figure 14, results obtained for the larger ROI size show the following: (1) The effectiveness in the application of the threshold is reduced, In fact, as shown in Figure 14a, the number of selected reliable directions decreases with the threshold less significantly than in the case of the ROI size of 5 km × 5 km (Figure 10). In addition, both RMSE and MBE trends, shown in Figure 14b,c, respectively, appear more flat with respect to the case of the smaller ROI size (Figures 11a and 13, respectively). Thus, the 12.5 km by 12.5 km ROIs appear to have an associated directional content higher than the ones for the 5 km by 5 km case.
(2) Both MS and SS directional estimations provide higher RMSE and (absolute) MBE values than the corresponding ones obtained at 5 km by 5 km ROI size. This can be justified considering that about half of the available in situ wind observations derived from stations with a maximum time distance from SAR acquisitions of 5 min, and the remaining measures does not exceed a 30 min delay. Moreover, the analysis is mainly performed in low-moderate wind speeds (i.e., up to 10 m/s).
(3) The MS estimation performance is still better than the one achieved by each SS processing, but the improvement reduces, especially in terms of achieved RMSE values with respect to the Low Scale. In fact, the improvement of RMSE reduces up to −0.09° w.r.t. the Low Scale, to −2. It is important to mention that MS wind direction improvement also yields to better wind speed retrieval through a GMF inversion. In fact, applying the C-band MODel 5.N (CMOD5.N) backscattering model [10], it was found that wind speed RMSE values are reduced up to −0.3 m/s using the wind direction derived from the MS processing rather than the one from each SS directional estimation.
Furthermore, previous analyses were also carried out on a ROI size of 12.5 km by 12.5 km. However, the number of available ROIs reduces up to 203. This reduction depends on the greater number of larger ROIs containing non-wind features and that are then discarded following the above-described procedure (Section 3.2). Moreover, in the case of NOAA stations located near the S-1 frames boundaries, it was not always allowed to extract 12.5 km by 12.5 km ROIs entirely included within the frames.
As reported in Figure 14, results obtained for the larger ROI size show the following: (1) The effectiveness in the application of the threshold ME TH is reduced, In fact, as shown in Figure 14a, the number of selected reliable directions decreases with the threshold less significantly than in the case of the ROI size of 5 km × 5 km ( Figure 10). In addition, both RMSE and MBE trends, shown in Figure 14b,c, respectively, appear more flat with respect to the case of the smaller ROI size (Figure 11a and Figure 13, respectively). Thus, the 12.5 km by 12.5 km ROIs appear to have an associated directional content higher than the ones for the 5 km by 5 km case. (2) Both MS and SS directional estimations provide higher RMSE and (absolute) MBE values than the corresponding ones obtained at 5 km by 5 km ROI size. This can be justified considering that about half of the available in situ wind observations derived from stations with a maximum time distance from SAR acquisitions of 5 min, and the remaining measures does not exceed a 30 min delay. Moreover, the analysis is mainly performed in low-moderate wind speeds (i.e., up to 10 m/s).
(3) The MS estimation performance is still better than the one achieved by each SS processing, but the improvement reduces, especially in terms of achieved RMSE values with respect to the Low Scale. In fact, the improvement of RMSE reduces up to −0.09 • w.r.t. the Low Scale, to −2.38 • w.r.t. the Medium Scale, and to −7.25 • w.r.t. the High Scale, whereas it reaches, with ME TH = 7.5 • , about −2.23 • , −4.15 • , and −10.52 • , respectively, for the 5 km by 5 km ROI size processing (Figure 11b).

Investigations on the Dependence of the Patterns Modulation Scale on Environmental Parameters
By exploiting wind data from NOAA in situ stations and related information (Table  1), some investigations were carried out in order to understand which physical phenomena may affect wind rows modulation scales, or-more precisely-processing scales lo-

Investigations on the Dependence of the Patterns Modulation Scale on Environmental Parameters
By exploiting wind data from NOAA in situ stations and related information (Table 1), some investigations were carried out in order to understand which physical phenomena may affect wind rows modulation scales, or-more precisely-processing scales locally selected as optimal by the MS LG-Mod.
It is worth noting that the following findings refer to directional estimations obtained by the MS LG-Mod on the Sentinel-1 5 km by 5 km ROIs selected, applying a final thresholding with ME TH = 15 • and α = 0.05 as well. The latter threshold was chosen with the double aim to obtain reliable enough directional estimations and an acceptable number of 226 samples for statistical purpose.
The bar graph in Figure 15 shows the number of occurrences of the optimal processing scales (y-axis) for each one of the 17 NOAA stations (x-axis). The figure evidences the following: The bar graph in Figure 15 shows the number of occurrences of the optimal processing scales (y-axis) for each one of the 17 NOAA stations (x-axis). The figure evidences the following: • The Low Scale (yellow) was the most frequently selected, i.e., 116 of the 226 occurrences (51.3%). The RMSE evaluated on these 116 occurrences was 21.38°.

•
The Medium Scale (red) was selected in 73 of the 226 occurrences (32.3%). The RMSE evaluated on the 73 occurrences was 21.44°.
A further analysis was suggested on the dependence of the wind rows modulation scale on both the distance from the coastline and the bathymetry for the available NOAA stations locations. Figure 16a,b displays, respectively, the bar graphs of the percentage of occurrences of the optimal processing scales (y-axes) grouped into three ranges of (a) distance from the coastline (near, 0≤R<10km; mid, 10km≤R<20km; far, R≥20km) and (b) bathymetry (low, 0≤Z<75m; moderate, 75m≤Z<150m; high, Z>150m) (x-axes). The former ranges are quite uniformly covered, with the 34.9%, the 34.9%, and the 30.2% of samples in the near, mid, and far range, respectively. Most of the in situ stations are located at low (57.9%) and moderate (33.2%) bathymetry, while just a small percentage of stations (8.8%) is characterized by high values of bathymetry. The following points were observed: • The percentage of occurrences of the Low-and the Medium-Scale selection reduces and increases, respectively, while the distance R from the coastline is increasing. On the other hand, the High-Scale selection presents about the same percentage of occurrences in the three range of distance.

•
The Low-and the Medium-Scale selections do not show a well-defined relationship between the percentage of occurrences and the bathymetry Z. Instead, the frequency of the High-Scale selection increases with the bathymetry. A further analysis was suggested on the dependence of the wind rows modulation scale on both the distance from the coastline and the bathymetry for the available NOAA stations locations. Figure 16a,b displays, respectively, the bar graphs of the percentage of occurrences of the optimal processing scales (y-axes) grouped into three ranges of (a) distance from the coastline (near, 0 ≤ R< 10 km; mid, 10 km ≤ R < 20 km; far, R ≥ 20 km) and (b) bathymetry (low, 0 ≤ Z < 75 m; moderate, 75 m ≤ Z < 150 m; high, Z > 150 m) (x-axes). The former ranges are quite uniformly covered, with the 34.9%, the 34.9%, and the 30.2% of samples in the near, mid, and far range, respectively. Most of the in situ stations are located at low (57.9%) and moderate (33.2%) bathymetry, while just a small percentage of stations (8.8%) is characterized by high values of bathymetry. The following points were observed:

•
The percentage of occurrences of the Low-and the Medium-Scale selection reduces and increases, respectively, while the distance R from the coastline is increasing. On the other hand, the High-Scale selection presents about the same percentage of occurrences in the three range of distance. The dependence of the wind rows modulation scale on both wind speed and direction as measured by NOAA in situ stations was also examined. Figure 17a,b shows, respectively, the bar graphs of the percentage of occurrences of the optimal processing scales (y-axes) grouped into three ranges of wind (a) speed (light: 2≤W<5.5 m/s; moderate: 5.5 ≤W<13.8 m/s; Strong: W>13.8 m/s) and (b) direction (on-shore and off-shore notations refer to winds that go from the sea toward the land and vice versa respectively; along-shore means winds that blow almost parallel to the coast) (x-axes). Most of the in situ wind observations show a moderate (the 67.7% of the total amount) wind speed regime, while the remaining observations are in the light (31.4%) and strong (0.9%) range of wind. Wind rows were detected almost exclusively in the light and moderate wind-speed regimes. In particular, the range of moderate wind speeds, from 5.5 to 13.8 m/s, seems to be the most suitable to allow visibility of wind rows on SAR images. Concerning the actual in situ wind directions, available observations are in the off-shore, along-shore, and on-shore ranges with the percentage of occurrences of the 25.7%, the 67.7%, and the 6.6%, respectively. Wind rows were detected mostly in case of along-shore winds rather than off-and on-shore winds. The following points were observed: • From Figure 17a, the percentage of occurrences of the Low-and the Medium-Scale selection reduces and increases, respectively, while the wind speed is increasing from the light to the moderate regimes. The High-Scale selection presents about the same percentage of occurrences in both regimes. The Strong case is poorly sampled, and it may be only stated that strong wind speeds were proven to be a bad condition for wind rows visibility. • From Figure 17b, the Low-, the Medium-, and the High-Scale selections do not show a well-defined relationship between the percentage of occurrences and the wind direction range.
Unfortunately, in situ measurements of air and sea temperature were only partially available. For this reason, it was not possible to extend the same previous analysis to these variables, although these types of information are important for boundary layer roll and wind streak development. The dependence of the wind rows modulation scale on both wind speed and direction as measured by NOAA in situ stations was also examined. Figure 17a,b shows, respectively, the bar graphs of the percentage of occurrences of the optimal processing scales (y-axes) grouped into three ranges of wind (a) speed (light: 2 ≤ W < 5.5 m/s; moderate: 5.5 ≤ W < 13.8 m/s; Strong: W > 13.8 m/s) and (b) direction (on-shore and off-shore notations refer to winds that go from the sea toward the land and vice versa respectively; along-shore means winds that blow almost parallel to the coast) (x-axes). Most of the in situ wind observations show a moderate (the 67.7% of the total amount) wind speed regime, while the remaining observations are in the light (31.4%) and strong (0.9%) range of wind. Wind rows were detected almost exclusively in the light and moderate wind-speed regimes. In particular, the range of moderate wind speeds, from 5.5 to 13.8 m/s, seems to be the most suitable to allow visibility of wind rows on SAR images. Concerning the actual in situ wind directions, available observations are in the off-shore, along-shore, and on-shore ranges with the percentage of occurrences of the 25.7%, the 67.7%, and the 6.6%, respectively. Wind rows were detected mostly in case of along-shore winds rather than off-and on-shore winds. The following points were observed: • From Figure 17a, the percentage of occurrences of the Low-and the Medium-Scale selection reduces and increases, respectively, while the wind speed is increasing from the light to the moderate regimes. The High-Scale selection presents about the same percentage of occurrences in both regimes. The Strong case is poorly sampled, and it may be only stated that strong wind speeds were proven to be a bad condition for wind rows visibility. • From Figure 17b, the Low-, the Medium-, and the High-Scale selections do not show a well-defined relationship between the percentage of occurrences and the wind direction range.  Two Sentinel-1 IW images are shown in Figure 18 as NRCS values with superimposed wind directions retrieved by the MS LG-Mod (with ME TH = 10 • and α = 0.05). SSW directions were estimated by the algorithm over a grid of about 0.05 • by 0.05 • (latitude by longitude). Three processing scales, i.e., the High Scale (320 m by 320 m, blue), the Medium Scale (160 m by 160 m, red), and the Low Scale (80 m by 80 m, yellow) were exploited with the aim at spanning the pixel size range (100-400 m) for wind rows detection. According to the DAMC, each different color represents the scale that is locally regarded as optimal by the MS LG-Mod. SSW directions were estimated by the algorithm over a grid of about 0.05° by 0.05° (latitude by longitude). Three processing scales, i.e., the High Scale (320 m by 320 m, blue), the Medium Scale (160 m by 160 m, red), and the Low Scale (80 m by 80 m, yellow) were exploited with the aim at spanning the pixel size range (100-400 m) for wind rows detection. According to the DAMC, each different color represents the scale that is locally regarded as optimal by the MS LG-Mod. The full frame analysis highlights that the Low Scale (80 m by 80 m, yellow) was the most frequently selected while the other scales were especially found in the coastal areas characterized by moderate bathymetry. These findings partially confirm what was found in the analysis made by using buoys data.

LG-Mod Results from Simulated SAR Images
As pointed out by Figure 5a,b, the adoption of the DAMC criterion allows the MS LG-Mod algorithm to locally choose the processing scale that may be regarded as optimal for patterns enhancement. This scale is clearly correlated on the varying spatial wavelength of the local patterns.
As shown in Figure 6a,b, the same DAMC criterion guarantees, as expected, the maximum number of reliable directions selected by the MS LG-Mod in comparison with those obtained by each SS processing, when the same threshold of acceptance ME TH is applied.
Moreover, the MS LG-Mod optimal scale selection provides the best choice for local directional estimation depending on the set of previously fixed processing scales. This statement is confirmed by the plots in Figures 7 and 8, from which the RMSE value derived from MS LG-Mod directional estimates is always lower than those values obtained from SS LG-Mod ones. This finding is true considering both MS and SS populations of reliable directions obtained from each applied threshold ME TH .
However, it is worth noting that the lower the threshold the better each SS directional estimation, according to the fact that RMSE differences tend progressively to zero for decreasing values of the threshold ME TH (Figure 8). This evidence suggests the idea that SS directional populations become mutually exclusive when the applied threshold ME TH is very low (Figure 6). On the contrary, when the latter constrain is not very strict (i.e., higher threshold), an overlap of SS populations may occur and the MS analysis thus supports the optimal processing scale selection among the different SS employed ones.

LG-Mod Results from Real SAR Images
As detailed in Reference [17], for the SS LG-Mod algorithm, the Marginal Error ME ROI α is a function of the scalar Mean Resultant Length (MRL), which represents a measure of a ROI local directions alignment. It follows that the parameter ME ROI α indirectly provides directional content of a ROI. The smaller the Marginal Error, as well as the better the accuracy of the directional estimate, the higher the ROI directional content. This is also true within the framework of the novel MS LG-Mod.
Hence, ME ROI α is capable of detecting wind rows patterns by means of their directional content estimation only in the case of an effective removal of artifacts and non-wind features onto SAR images. Otherwise, the directional content sensed by the parameter ME ROI α would be obviously not related to the local wind direction to be estimated. As a consequence, such conditions need to be avoided, e.g., as implemented by the above described strategy (Section 3.2).

Exploitation of ME ROI
α Aimed at Selection of the Optimal Local Processing Scale Statistics plotted in Figures 10-12 obtained on the Sentinel-1 ROIs selection clearly confirm results achieved for simulated SAR images. This finding supports the fact that, in the case of artifacts and non-wind features removal from real SAR images, the ME ROI α parameter and the DAMC criterion-defined within the MS LG-Mod-guarantee the automatic selection of the optimal local processing scale for the best achievable local direction estimation. In such case, the novel MS analysis proved to be able to sense local directional content better than each SS processing. The MS directional analysis also leads to an improvement in SAR wind speed retrieval, although the found minimal impact on this retrieval (the maximum RMSE difference is of about −0.3 m/s) may be caused by the actual low to moderate wind conditions of the S-1 dataset.
Even though in Reference [28] and in Reference [30] frequencies of modulations belonging to the BLRs and WSs are enhanced, no effective selection of optimal scale is performed within these two methods.
The two approaches proposed in Reference [32], to choose the optimal estimated direction are based either on the use of a priori information from a NWP model or the visual inspection of patterns visible on SAR amplitude. The first option makes the LG method dependent on external data, which could be not sufficiently accurate. The second one makes the method itself not completely automated.
The criterion proposed in Reference [17], which defines the optimal processing scale as the one which can globally yield in a single SAR the largest percentage of reliable ROIs, gave promising results. However, the different mechanisms that may induce wind aligned streaks on SAR images produce SAR NRCS modulations with wavelengths ranging from few hundred meters to few kilometers [22], which may occur even in the same SAR image. This may hamper the use of a criterion previously described in Reference [17].
In Reference [27] the computation of a directional factor at different scales was proposed to find the optimal scale. The scale with the highest value of that parameter gives the optimal scale of the directional feature. To the best of our knowledge, this is the only attempt to find the local optimal scale by means of an automatic criterion. Nonetheless, no quantitative comparison was made between the performance of the single-and multiple-scale algorithms.
In the present work, a local MS criterion for the optimal pixel size (scale) automatic selection was introduced and it was proved to perform better with respect to the use of a SS one, especially at high-resolution (i.e., 5 km by 5 km ROI size) processing.

Investigations on the Dependence of the Patterns Modulation Scale on Environmental Parameters
The results reported in Section 3.4 revealed that wind rows modulation scales may be determined by several concurring factors, such as distance from the coast, bathymetry and wind characteristics as well. The following points are offered to summarize the results:

•
The Low (80 m × 80 m) and the Medium (160 m × 160 m) Scales represent the more adequate ones for wind rows detection, as shown by both the frequency of occurrences (i.e., 51.3% and 32.3%, respectively) and the accuracy of directional estimations (i.e., RMSE values equal to 21.38 • and 21.44 • , respectively). Furthermore, both scales may occur at different distances from the coast, but it is not clear why the frequency of the Low-and Medium-Scale selection decreases and increases, respectively, while the distance from the cost increases (Figure 17a). A possible interpretation is that windinduced SAR signatures, namely wind rows, may be affected close to the coast by a mixing of different modulations due to both Wind Streaks and Boundary Layer Rolls. The former are typically characterized by spatial wavelengths that are more compatible with the Low Scale, whereas the latter are more detectable with the Medium Scale.

•
The High Scale (320 m × 320 m) resulted in being, instead, the less frequent, as well as the less accurate, scale employed, with a percentage of occurrences and a RMSE equal to 16.4% and 22.33 • , respectively ( Figure 12). This scale seems to be just influenced by the bathymetry, showing an increasing of percentage of occurrences with the bathymetry itself (Figure 16b). To the best of our knowledge, no authors have been investigated on the effects of the bathymetry on SAR backscattering modulations induced by atmospheric phenomena, such as wind rows. However, it is well-known that, in shallow areas, small-scale features of the bathymetry cause strong current gradients, which then in turn modulate the spectrum of the short Bragg waves [39]. Zhang et al., in 2017, found different bathymetric features on SAR imagery under different sea states. In particular, under low to moderate wind speeds (3.1~6.3 m/s), they observed wide bright patterns with an average width of 6 km and quasi-linear features only 1 km wide for high winds (5.4~13.9 m/s) [40]. On the other hand, it is recognized that the approximation of waves to shallow areas results in a decrease in the wavelength, which is linked to water depths [41]. However, although wind rows derived modulations in the SAR image are primarily caused by variations in the nearsurface wind field, the bathymetry may change their spatial wavelengths, as suggested by the increasing of percentage of High-Scale occurrences with the bathymetry.

•
Wind rows are detected mostly in the case of light-moderate wind speed regimes (i.e., W less than 13.8 m/s) and along-shore winds (Figure 17a,b, respectively). This finding confirms the results in References [20,21], where it is reported that wind rows are most commonly observed at wind speeds near 8-9 m/s. Whatever the wind speed and direction, the Low-and the Medium-Scale selections occur more frequently than the High-Scale one. • However, the wind rows modulation scale may be also influenced by any further physical phenomena that cause variations of the sea surface roughness and, in turn, the radar signal backscattered from the sea. As reported in Reference [18], there is a large variability in wavelength and amplitude of the wind rows as a result of several processes acting on the wind rows and related to the inhomogeneities in terrain, roughness, and heating and also to vertical wind shear [42].

Conclusions
A novel MS analysis was introduced and fully developed within the LG-Mod algorithm [8,17] with the aim to improve the reliability of SAR-derived SSW directions retrieval.
The MS LG-Mod was proven to be able to automatically select the optimal local processing scale for wind patterns enhancement in order to provide the best achievable local directional estimation, once a set of processing scales has been fixed.
The assessment of the MS LG-Mod directional estimation was carried out on both simulated SAR data and Sentinel-1 images.
Results achieved from simulated SAR images showed that the MS analysis allows to sense local directional content better than each SS processing. In fact, the MS LG-Mod evidenced to locally choose the processing scale that is optimal for patterns enhancement and clearly correlated on the spatial wavelength of the local patterns.
The MS procedure performed always better than each SS one, with improvements in terms of both directional RMSE values and percentages of reliable output directions achieved.
For the analysis of the S-1 dataset, the removal of artifacts and non-wind features from SAR amplitudes was mandatory before directional estimations. Hence, the application of the MS LG-Mod algorithm always requires taking into account the removal of such features.
Further developments should lead to completely automate this step. Finally, a tentative analysis on the dependence of the patterns modulation scale on environmental parameters was made. The Low (80 m × 80 m) and the Medium (160 m × 160 m) Scales resulted the more adequate ones for wind rows detection on Sentinel-1 data, as evidenced by both their frequency of occurrence and the accuracy of directional estimation. The High Scale (320 m × 320 m) was instead the less frequent and accurate scale employed.
Hence, the analysis revealed that both the Low and the Medium Scales may occur at different distances from the coast, showing a frequency of the Low-and Medium-Scale selection that decreases and increases, respectively, with the distance from the cost itself. Moreover, wind rows were detected mostly in light-moderate wind speed regimes and for along-shore winds, for the examined coastal dataset. However, wind rows visibility and their modulation scales may be determined by several concurring physical parameters, which cause variations of the sea surface roughness. That is the reason to further investigate, e.g., using the full-frame analysis only introduced in this work.