## 2 Materials and Methods
### Extraction of Land Parcels
#### 2.4.1 Mosaicking and Color Correction  
Only the scenes acquired on 29 May and 7 July captured the entire study area; all the other dates included two or more images, meaning that more than one image was acquired on the same date. The aim of the mosaicking step was to make the dataset homogeneous, modeling it through the rule \"one image per date\", to obtain a comprehensive set of twenty-five images, each referring to only one acquisition date. In instances in which the  
\begin{table}
\begin{tabular}{c c} \hline
**Spectral Bands** & **Wavelength Range** \\ \hline Coastal Blue & 431–452 \\ Blue & 465–515 \\ Green I & 513–549 \\ Green II & 547–583 \\ Yellow & 600–620 \\ Red & 650–680 \\ Red-Edge & 697–713 \\ NIR & 845–885 \\ \hline \end{tabular}
\end{table}
Table 1: PlanetScope’s spectral bands and their respective wavelength ranges.  
\begin{table}
\begin{tabular}{c c} \hline
**Month** & **Dates** \\ \hline March & 7; 16; 23; 30 \\ April & 7; 27 \\ May & 5; 29 \\ June & 2; 21; 27 \\ July & 7; 12; 17; 22 \\ August & 1; 8; 14; 18; 27 \\ September & 1; 4; 27 \\ October & 1; 7 \\ \hline \end{tabular}
\end{table}
Table 2: PlanetScope acquisition dates.  
mosaicking of two images was necessary to ensure comprehensive coverage of the test area, the raster pairs were mosaicked using the workflow proposed by Seamless Mosaic, a tool available within the software ENVI (6.1 version). For each pair, the reference and the adjusted scenes were selected on the basis of their dimensions: the largest scene was used as the reference scene, while the smallest scene was used as the adjusted scene (Figure 3).  
Consequently, the reference scene was layered on top of the mosaic. Then, the image footprints of both the reference and the adjusted images were clearly delineated in order to compute which portion of the adjusted scene would fill the empty portion of the reference scene. The color-correction process was then applied to each mosaic. This algorithm maps discrete greyscale levels from the histogram of the adjusted scene to the corresponding greyscale levels in the reference scenes, considering only the area where the two images overlap. The procedure subsequently interpolates the statistics of the two image histograms, and the resulting matching algorithm is applied to the area of the adjusted image where the reference does not overlap. This approach minimized any radiometric artefacts resulting from the mosaicking. The resultant dataset consisted of twenty-five images, with each image corresponding to a specific date (for a complete list of the resulting image IDs, refer to Table 3).  
\begin{table}
\begin{tabular}{c c c} \hline
**ID** & **Acquisition Date** & **File Name** \\ \hline
1 & 03/07/2023 & 20230307\_090404 & 47 24b2\_3B. AnalyticMS\_SR\_8b\_clip \\
2 & 03/16/2023 & 20230316\_085941 & 53 2430\_3B. AnalyticMS\_SR\_8b\_clip \\
3 & 03/23/2023 & 20230323 & 093937\_17.227b.3B. AnalyticMS\_SR\_8b\_clip \\
4 & 03/30/2023 & 20230330\_090739 & 19 24 bbe.3B. AnalyticMS\_SR\_8b\_clip \\
5 & 04/07/2023 & 20230407\_090308 & 37 2465\_3B. AnalyticMS\_SR\_8b\_clip \\
6 & 04/27/2023 & 20230427\_093907 & 91 2482\_3B. AnalyticMS\_SR\_8b\_clip \\
7 & 05/05/2023 & 20230505 & 093834 78\_2479.3B. AnalyticMS\_SR\_8b\_clip \\ \hline \end{tabular}
\end{table}
Table 3: Scene ID, acquisition date and file name of each of the 25 scenes composing the dataset. The scenes captured by single images are highlighted in green, while the scenes generated by the mosaicking of two images are highlighted in orange.  
Figure 3: The figure shows the case of the scene of 23 March, acquired by the platform 227b: (**a**) this panel contains the bigger image (the reference scene), which was acquired at 9:39:37; (**b**) this panel contains the smaller image (the adjusted scene), which was acquired at 9:39:34.