{"id":16040,"date":"2025-12-01T10:00:00","date_gmt":"2025-12-01T09:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=16040"},"modified":"2025-11-28T21:45:47","modified_gmt":"2025-11-28T20:45:47","slug":"when-deep-learning-dives-beneath-the-surface-mapping-corals-with-qgis-and-pytorch","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/en\/when-deep-learning-dives-beneath-the-surface-mapping-corals-with-qgis-and-pytorch\/","title":{"rendered":"When Deep Learning Dives Beneath the Surface: Mapping Corals with QGIS and PyTorch"},"content":{"rendered":"\n

Deep learning is revolutionizing satellite image analysis.<\/p>\n\n\n\n

Long reserved for large laboratories or proprietary software, it is now becoming available to the wider world thanks to PyTorch and QGIS.<\/p>\n\n\n\n

This article explores the principles of Deep Learning applied to geomatics, compares ESRI models with those that can be used in QGIS, and concludes with a concrete example: the automatic detection of coral reefs using Sentinel-2 images.<\/p>\n\n\n\n

\n\n\n

\n

\n

Contenu <\/p>\nToggle<\/span>

Introduction to Deep Learning applied to geomatics<\/a>

Practical example: coral segmentation using Sentinel-2 images in QGIS<\/a>

Conclusion<\/a><\/li><\/ul><\/nav><\/div>\n\n\n\n\n

<\/span>Introduction to Deep Learning applied to geomatics<\/span><\/h2>\n\n\n\n

Deep Learning is a branch of artificial intelligence inspired by the functioning of the human brain. It is based on artificial neural networks capable of learning from examples, without being explicitly told all the rules.<\/p>\n\n\n\n

Unlike traditional classification methods, where the indicators and thresholds are chosen by the user, Deep Learning automatically discovers relevant structures and patterns in the data.<\/p>\n\n\n\n

In the field of geomatics, this approach opens up impressive possibilities:<\/p>\n\n\n\n

recognition of urban, agricultural, or forest areas from satellite images,<\/li>\n\n\n\n
detection of changes or natural disasters,<\/li>\n\n\n\n
identification of specific elements (roads, roofs, corals, ships, etc.),<\/li>\n\n\n\n

detailed segmentation of landscapes from Sentinel, PlanetScope, or drone images.<\/li>\n<\/ul>\n\n\n\n

The principle is simple: a model is trained on a large set of annotated images until it learns to reproduce the desired task (for example, distinguishing between water, vegetation, and sand). Once trained, the model can be applied to new areas to automatically generate high-resolution thematic maps.<\/p>\n\n\n\n

<\/span>Why \u201cdeep\u201d?<\/span><\/h3>\n\n\n\n

The term \u201cdeep\u201d comes from the fact that these networks have many successive layers\u2014sometimes dozens.<\/p>\n\n\n\n

Each layer learns to recognize increasingly complex patterns:<\/p>\n\n\n\n

the first layers detect edges and textures,<\/li>\n\n\n\n
the next identify shapes or structures,<\/li>\n\n\n\n

and the last understand entire objects or contexts.<\/li>\n<\/ul>\n\n\n\n

It is this hierarchy of representations that gives deep learning its power, but also its appetite for data and computing power.<\/p>\n\n\n\n

<\/span>Deep learning and remote sensing<\/span><\/h3>\n\n\n\n

In remote sensing, the most commonly used models are image segmentation models, which are capable of assigning a class to each pixel.<\/p>\n\n\n\n

Architectures such as U-Net, DeepLab, and Mask R-CNN have become benchmarks for automatic mapping based on multispectral imagery.<\/p>\n\n\n\n

These models are generally trained with frameworks such as PyTorch or TensorFlow, then deployed in GIS environments.<\/p>\n\n\n\n

Both major worlds have taken an interest in this:<\/p>\n\n\n\n

ESRI<\/strong>, with its .dlpk (Deep Learning Package) model format integrated into ArcGIS Pro and ArcGIS Online;<\/li>\n\n\n\n

QGIS<\/strong>, which allows PyTorch or TensorFlow models to be used via the Processing Toolbox or custom Python scripts.<\/li>\n<\/ul>\n\n\n\n
\n\n\n\n
<\/span>ESRI’s DLPK format<\/span><\/h3>\n\n\n\n
ESRI was one of the first GIS players to integrate deep learning directly into its ArcGIS ecosystem.<\/p>\n\n\n\n
To facilitate the exchange and deployment of models, the company created a standard format: DLPK (Deep Learning Package).<\/p>\n\n\n\n
A .dlpk file is not just a model: it is a complete and portable package containing all the elements necessary for its execution.<\/p>\n\n\n\n
It generally contains:<\/p>\n\n\n\n
\n
The trained model (often in PyTorch .pth or TensorFlow .h5 format)<\/li>\n\n\n\n
A JSON definition file describing the model architecture, expected parameters, and class names<\/li>\n\n\n\n
Metadata: input type (raster, tile, image), patch size, number of bands, normalization, etc.<\/li>\n\n\n\n
Training samples (optional) to document or retrain the model.<\/li>\n<\/ul>\n\n\n\n
Thanks to this organization, ArcGIS Pro or ArcGIS Online can automatically interpret the model without the need to write code.<\/p>\n\n\n\n
Tools such as \u201cClassify Pixels Using Deep Learning\u201d or \u201cDetect Objects Using Deep Learning\u201d directly read the .dlpk, load the model, and perform inference on a raster or image set.<\/p>\n\n\n\n
This turnkey approach has two major advantages:<\/p>\n\n\n\n
\n
Interoperability: a model trained elsewhere can be used by any ArcGIS user, without complex dependencies.<\/li>\n\n\n\n
Replicability: metadata ensures that the model is applied under the same conditions as when it was trained.<\/li>\n<\/ol>\n\n\n\n
The downside, of course, is the closed format: .dlpk remains tied to the ArcGIS ecosystem and is not always easy to use elsewhere.<\/p>\n\n\n\n
\n\n\n\n
<\/span>The open source counterpart: PyTorch and QGIS<\/span><\/h3>\n\n\n\n
On the open source side, QGIS does not impose a proprietary format.<\/p>\n\n\n\n
Models are simply saved in their native format (often .pth for PyTorch or .pt) and executed via Python scripts integrated into the Processing Toolbox.<\/p>\n\n\n\n
The idea is the same as with ESRI:<\/p>\n\n\n\n
a multispectral image is loaded, a trained model is applied, and a class or probability map is generated.<\/p>\n\n\n\n
But instead of relying on a packaged format such as .dlpk, QGIS gives developers complete freedom:<\/p>\n\n\n\n
\n
The model is read with torch.load().<\/li>\n\n\n\n
The input bands can be selected dynamically (for example, B4-B3-B2 for RGB, or B8-B4-B3 in false color).<\/li>\n\n\n\n
The Python script controls the entire processing flow: normalization, water masking (NDWI), block cutting, merging of results, etc.<\/li>\n<\/ul>\n\n\n\n
This approach allows for maximum flexibility, which is particularly useful for research or experimentation.
For example, a U-Net model trained in PyTorch can be applied directly in QGIS via a custom script\u2014without relying on ArcGIS Pro.<\/p>\n\n\n\n
QGIS then becomes a true AI spatial analysis laboratory, where users can:<\/p>\n\n\n\n
\n
test several models (.pth) from the community,<\/li>\n\n\n\n
adapt pre-processing according to the area (coastal strip, forest, urban, etc.),<\/li>\n\n\n\n
automate an entire workflow using Processing algorithms,<\/li>\n\n\n\n
and combine the results with other GIS layers (vegetation, bathymetry, land cover, etc.).<\/li>\n<\/ul>\n\n\n\n
\n\n\n\n
<\/span>Summary<\/span><\/h2>\n\n\n\n
Aspect<\/th> ESRI (DLPK)<\/th> QGIS (PyTorch)<\/th><\/tr><\/thead>
Format<\/strong><\/td> .dlpk<\/code> (complete package)<\/td> .pth<\/code> or .pt<\/code> (model only)<\/td><\/tr>
Interoperability<\/strong><\/td> Simple, but proprietary<\/td> Free and modifiable<\/td><\/tr>
Use<\/strong><\/td> Integrated ArcGIS tools (\u201cClassify Pixels,\u201d \u201cDetect Objects\u201d)<\/td> Custom Python Processing scripts<\/td><\/tr>
Customization<\/strong><\/td> Limited<\/td> Total<\/td><\/tr>
Learning curve<\/strong><\/td> Simpler for the end user<\/td> More flexible for the developer or researcher<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
\n\n\n\n
<\/span>Practical example: coral segmentation using Sentinel-2 images in QGIS<\/span><\/h2>\n\n\n\n
Now that we have explored the logic behind deep models and the formats used by ESRI and QGIS, let’s move on to a concrete example: analyzing coral reefs using Sentinel-2 satellite images.
The goal is to distinguish marine areas (sandy bottoms, seagrass beds, corals) from terrestrial or turbid areas using a pre-trained U-Net model in PyTorch.<\/p>\n\n\n\n
\n\n\n\n
<\/span>Data used<\/span><\/h3>\n\n\n\n
The input image comes from Sentinel-2, a mission of the European Copernicus program.
These free multispectral images offer a resolution of 10 to 20 m and cover several bands in the visible and infrared:<\/p>\n\n\n\n
BandNameWavelength (nm)Main useB2Blue490Water, turbidityB3Green560Vegetation, marine environmentB4Red665Soil, vegetation, coralsB8NIR842Land\/water differentiationB11SWIR11610Humidity, sand, clouds<\/p>\n\n\n\n
Band<\/th> Name<\/th> Wavelength (nm)<\/th> Main use<\/th><\/tr><\/thead>
B2<\/td> Blue<\/td> 490<\/td> Water, turbidity<\/td><\/tr>
B3<\/td> Green<\/td> 560<\/td> Vegetation, marine environment<\/td><\/tr>
B4<\/td> Red<\/td> 665<\/td> Soil, vegetation, corals<\/td><\/tr>
B8<\/td> NIR<\/td> 842<\/td> Land\/water differentiation<\/td><\/tr>
B11<\/td> SWIR1<\/td> 1610<\/td> Humidity, sand, clouds<\/td><\/tr><\/tbody><\/table><\/figure>\n\n\n\n
The raster is derived from S2DR3 processing – Sentinel-2 super-resolution at 1 m, which improves the resolution of all Sentinel 2 image bands to approximately 1 m (see \u201cUsing S2DR3 in Google Colab to study corals in Mauritius\u201d). The image used is the one that served as an example in the article cited.
The user simply loads the multispectral .tif file, for example:
Palmar_MS.tif.<\/p>\n\n\n\n
\n\n\n\n
<\/span>Processing in QGIS<\/span><\/h3>\n\n\n\n
The analysis is performed using a Python algorithm integrated into the Processing toolbox.<\/p>\n\n\n\n
This script loads a PyTorch model (unet_coraux.pth) and applies segmentation to the image in several steps:<\/p>\n\n\n\n
\n
Reading the raster and normalizing the bands
The values are scaled to a level compatible with the model training.<\/li>\n\n\n\n
Optional masking of land areas
A Normalized Difference Water Index (NDWI) is calculated to isolate the sea.
Pixels with high NDWI values are considered marine and processed by the model; the others are masked.<\/li>\n\n\n\n
Cutting into blocks (patches)
The image is processed in portions to avoid memory overload.
Each block is analyzed independently, then the results are merged.<\/li>\n\n\n\n
Application of the U-Net model<\/li>\n\n\n\n
The model performs pixel-by-pixel segmentation: it assigns each pixel a probability of belonging to the \u201ccoral\u201d or \u201cnon-coral\u201d class.
The result is an output raster containing probability or class values.<\/li>\n\n\n\n
Saving the output raster
The result is saved in GeoTIFF format (palmar_model_9.tif), ready to be overlaid with other GIS layers.<\/li>\n<\/ol>\n\n\n\n
\n\n\n\n
<\/span>Procedure<\/span><\/h3>\n\n\n\n
<\/span>Preliminary operations<\/span><\/h4>\n\n\n\n
We will create a QGis processing script. As a prerequisite, you must install<\/p>\n\n\n\n
<\/div>
OsGeo4W shell<\/div>
<\/div><\/div>
<\/p>\n\n\n\n
python -m pip install torchgeo\npython -m pip install torch torchvision\npython -m pip install segmentation-models-pytorch\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n Then we need to download the desired model. To do this, enter the following script in the QGIS Python console:<\/p>\n\n\n\n <\/div>Python Console<\/div><\/div><\/div><\/p>\n\n\n\nimport segmentation_models_pytorch as smp\nimport torch\n\n# U-Net RGB pr\u00e9-entra\u00een\u00e9\nmodel = smp.Unet(\n encoder_name=\"resnet34\",\n encoder_weights=\"imagenet\",\n in_channels=3,\n classes=1 # corail\/non-corail\n)\n\n# Sauvegarde pour QGIS\ntorch.save(model, \"c:\/models\/unet_coraux.pth<\/em><\/strong>\")\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n In this example, we save the downloaded models in a directory c:\/models<\/p>\n\n\n\n \n\n\n\n<\/span>Setting up the processing script<\/span><\/h4>\n\n\n\nIn QGIS \u2192 Processing Toolbox \u2192 Scripts \u2192 New Script,<\/p>\n\n\n\n <\/div>traitement<\/div><\/div><\/div><\/p>\n\n\n\n# -*- coding: utf-8 -*-\n\"\"\"\nAlgorithme QGIS : Segmentation Coraux (U-Net) avec option masquage terrestre\nCompatible QGIS 3.44+\n\"\"\"\nfrom qgis.core import (\n QgsProcessing,\n QgsProcessingAlgorithm,\n QgsProcessingParameterRasterLayer,\n QgsProcessingParameterFile,\n QgsProcessingParameterEnum,\n QgsProcessingParameterRasterDestination,\n QgsProcessingParameterBoolean,\n)\nfrom qgis.PyQt.QtCore import QCoreApplication\nimport torch\nimport torch.serialization\nimport numpy as np\nfrom osgeo import gdal\nfrom scipy.signal import windows\n\n\nclass SegmentationCoraux(QgsProcessingAlgorithm):\n \"\"\"Segmentation des images Sentinel-2 \u00e0 1m avec un mod\u00e8le U-Net\"\"\"\n\n BAND_OPTIONS = [\"RGB (4-3-2)\", \"Fausse couleur (8-4-3)\"]\n\n def initAlgorithm(self, config=None):\n self.addParameter(\n QgsProcessingParameterRasterLayer(\n \"raster_input\",\n self.tr(\"Image Sentinel-2 (GeoTIFF multibandes)\")\n )\n )\n\n self.addParameter(\n QgsProcessingParameterEnum(\n \"band_set\",\n self.tr(\"S\u00e9lection manuelle des bandes\"),\n options=self.BAND_OPTIONS,\n defaultValue=0\n )\n )\n\n self.addParameter(\n QgsProcessingParameterFile(\n \"model_path\",\n self.tr(\"Mod\u00e8le PyTorch (.pth)\"),\n extension=\"pth\"\n )\n )\n\n self.addParameter(\n QgsProcessingParameterRasterDestination(\n \"output_raster\",\n self.tr(\"Raster de sortie (masque)\")\n )\n )\n\n # --- Nouvelle option : masquage terrestre ---\n self.addParameter(\n QgsProcessingParameterBoolean(\n \"mask_land\",\n self.tr(\"Masquer les zones terrestres (NDWI)\"),\n defaultValue=True\n )\n )\n\n def processAlgorithm(self, parameters, context, feedback):\n model_path = self.parameterAsFile(parameters, \"model_path\", context)\n mask_land = self.parameterAsBool(parameters, \"mask_land\", context)\n band_set = self.parameterAsEnum(parameters, \"band_set\", context)\n\n feedback.pushInfo(\"Chargement du mod\u00e8le PyTorch...\")\n try:\n from segmentation_models_pytorch.decoders.unet.model import Unet\n torch.serialization.add_safe_globals([Unet])\n except Exception as e:\n feedback.pushInfo(f\"Avertissement : impossible d\u2019ajouter Unet aux classes s\u00fbres : {e}\")\n\n model = torch.load(model_path, map_location=torch.device('cpu'), weights_only=False)\n model.eval()\n\n raster_input = self.parameterAsRasterLayer(parameters, \"raster_input\", context)\n ds = gdal.Open(raster_input.source())\n feedback.pushInfo(f\"Ouverture du raster : {raster_input.source()}\")\n nrows, ncols = ds.RasterYSize, ds.RasterXSize\n\n # --- D\u00e9tecter automatiquement les canaux d'entr\u00e9e ---\n conv1_weights = model.encoder.conv1.weight.data\n n_channels = conv1_weights.shape[1]\n feedback.pushInfo(f\"Le mod\u00e8le attend {n_channels} canaux d'entr\u00e9e.\")\n\n if n_channels == 3:\n # Analyser l\u2019importance des canaux\n channel_mean = conv1_weights.abs().mean(dim=(0, 2, 3)).numpy()\n sorted_idx = list(np.argsort(-channel_mean))\n feedback.pushInfo(f\"Classement des canaux par importance : {sorted_idx}\")\n\n band_selection = [4, 3, 2] # RGB classique\n feedback.pushInfo(f\"Utilisation automatique des bandes : {band_selection} (RGB)\")\n else:\n band_selection = [4, 3, 2] if band_set == 0 else [8, 4, 3]\n feedback.pushInfo(f\"Utilisation manuelle des bandes : {band_selection}\")\n\n # --- Lecture des bandes s\u00e9lectionn\u00e9es ---\n img_list = []\n for b in band_selection:\n band = ds.GetRasterBand(b)\n arr = band.ReadAsArray().astype(np.float32)\n img_list.append(arr)\n img = np.stack(img_list, axis=0)\n\n # Normalisation\n feedback.pushInfo(\"Normalisation des valeurs...\")\n img = (img - img.min()) \/ (img.max() - img.min() + 1e-6)\n\n # --- Calcul NDWI si option activ\u00e9e ---\n if mask_land and ds.RasterCount >= 8:\n try:\n B3 = ds.GetRasterBand(3).ReadAsArray().astype(np.float32)\n B8 = ds.GetRasterBand(8).ReadAsArray().astype(np.float32)\n ndwi = (B3 - B8) \/ (B3 + B8 + 1e-6)\n water_mask = ndwi > 0\n feedback.pushInfo(\"NDWI calcul\u00e9 : les zones terrestres seront exclues.\")\n except Exception as e:\n water_mask = np.ones((nrows, ncols), dtype=bool)\n feedback.pushInfo(f\"NDWI non calculable ({e}), toutes les zones seront trait\u00e9es.\")\n else:\n water_mask = np.ones((nrows, ncols), dtype=bool)\n feedback.pushInfo(\"Masquage terrestre d\u00e9sactiv\u00e9.\")\n\n # --- Traitement par blocs ---\n block_size = 2048\n overlap = 512\n output_mask = np.zeros((nrows, ncols), dtype=np.float32)\n weight = np.zeros((nrows, ncols), dtype=np.float32)\n\n total_blocks = ((nrows \/\/ (block_size - overlap) + 1) *\n (ncols \/\/ (block_size - overlap) + 1))\n done = 0\n feedback.pushInfo(\"D\u00e9but du traitement par blocs...\")\n\n for y in range(0, nrows, block_size - overlap):\n for x in range(0, ncols, block_size - overlap):\n if feedback.isCanceled():\n break\n\n block = img[:, y:y+block_size, x:x+block_size]\n if block.size == 0:\n continue\n\n mask_block_water = water_mask[y:y+block_size, x:x+block_size]\n if not mask_block_water.any():\n continue\n\n # ======= BLOC ORIGINAL PR\u00c9SERV\u00c9 =======\n with torch.no_grad():\n block_tensor = torch.from_numpy(block).unsqueeze(0)\n pred = model(block_tensor)\n mask_block = pred.squeeze().numpy()\n\n h, w = mask_block.shape\n mask_block *= mask_block_water[:h, :w]\n\n # Fen\u00eatre Hanning\n win_y = windows.hann(h)[:, None]\n win_x = windows.hann(w)[None, :]\n weight_block = win_y * win_x\n\n output_mask[y:y+h, x:x+w] += mask_block * weight_block\n weight[y:y+h, x:x+w] += weight_block\n\n done += 1\n progress = int(100 * done \/ total_blocks)\n feedback.setProgress(progress)\n # ======= FIN DU BLOC ORIGINAL =======\n\n # Fusion finale\n output_mask \/= np.maximum(weight, 1e-6)\n output_mask = np.clip(output_mask, 0, 1)\n\n # --- Sauvegarde GeoTIFF ---\n feedback.pushInfo(\"Sauvegarde du masque final...\")\n driver = gdal.GetDriverByName(\"GTiff\")\n out_path = self.parameterAsOutputLayer(parameters, \"output_raster\", context)\n out_ds = driver.Create(out_path, ncols, nrows, 1, gdal.GDT_Float32)\n out_ds.SetGeoTransform(ds.GetGeoTransform())\n out_ds.SetProjection(ds.GetProjection())\n out_ds.GetRasterBand(1).WriteArray(output_mask)\n out_ds.FlushCache()\n\n feedback.pushInfo(\"\u2705 Segmentation termin\u00e9e avec succ\u00e8s.\")\n return {\"output_raster\": out_path}\n\n # --- M\u00e9tadonn\u00e9es ---\n def name(self):\n return \"segmentation_coraux\"\n\n def displayName(self):\n return self.tr(\"Segmentation des images Sentinel 2 \u00e0 1m (U-Net)\")\n\n def group(self):\n return self.tr(\"Deep Learning\")\n\n def groupId(self):\n return \"deeplearning\"\n\n def tr(self, string):\n return QCoreApplication.translate(\"SegmentationCoraux\", string)\n\n def createInstance(self):\n return SegmentationCoraux()\n\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n You will find explanations of the different parts of the code in a future article.<\/p>\n\n\n\n The QGIS script presented here allows you to apply U\u2011Net-based image segmentation models saved in PyTorch (.pth) format. It is designed to process multispectral images, such as those from Sentinel\u20112 satellites, and automatically adapts the bands used according to the number of input channels expected by the model. For example, if the model is trained on RGB images, the script will select the Red, Green, and Blue bands; if the model expects more channels, it will offer a manual selection or use all available bands. The code also includes the ability to mask terrestrial areas using NDWI calculations, allowing segmentation to focus on water areas, for example to identify corals. In practice, any properly saved U-Net PyTorch model can be loaded and applied, provided that the architecture and weights are included in the .pth file.<\/p>\n\n\n\n \n\n\n\n<\/span>Usage<\/span><\/h4>\n\n\n\nRun the script. The following window will open:<\/p>\n\n\n\n <\/a><\/figure>\n\n\n\nManual band selection is only used if the script is unable to determine which bands to use based on the template.<\/p>\n\n\n\n \n\n\n\n<\/span>Interpreting the results<\/span><\/h3>\n\n\n\nThe raster generated from the model represents a probability map: values close to 1 indicate a high presence of corals, while values close to 0 correspond to non-coral areas (sand, algae, depth, etc.).<\/p>\n\n\n\n Appropriate symbology (from light blue to red) makes it easy to visualize the spatial distribution of probable coral areas. By combining this map with other data (bathymetry, substrate, turbidity), it becomes possible to estimate the vulnerability or degradation of reefs over time. <\/p>\n\n\n\n<\/a><\/figure>\n\n\n\n \n\n\n\nAdvantages of this approach<\/p>\n\n\n\n The integration of PyTorch into QGIS opens up new possibilities for AI-assisted environmental mapping:<\/p>\n\n\n\n \nOpen source and reproducible: the entire process can be shared, modified, or adapted to other coastal areas.<\/li>\n\n\n\n Local autonomy: no need for ArcGIS Pro or expensive licenses to test or apply deep learning models.<\/li>\n\n\n\n Flexible experimentation: other architectures (SegNet, DeepLabV3, etc.) can be tested, or preprocessing can be adapted to the specific characteristics of each area.<\/li>\n<\/ul>\n\n\n\n \n\n\n\n<\/span>Towards an ecosystem of open models<\/span><\/h3>\n\n\n\nUltimately, we could imagine a shared library of open source environmental models\u2014the free equivalent of the DLPK format\u2014where each .pth file would be accompanied by its description file (bands, normalization, classes).<\/p>\n\n\n\n QGIS could then offer an interface for importing, testing, and documenting these models, facilitating their reuse in tropical, coastal, or forest contexts.<\/p>\n\n\n\n \n\n\n\n<\/span>Conclusion<\/span><\/h2>\n\n\n\nThe rise of deep learning marks a new stage in the evolution of geomatics.<\/p>\n\n\n\n While traditional image processing relied on thresholds and spectral indices, neural models now learn to recognize complex shapes, textures, and signatures directly in pixels.<\/p>\n\n\n\n Thanks to tools such as ArcGIS Pro (with its DLPK) and QGIS (via PyTorch and custom scripts), this power is now accessible to everyone: researchers, technicians, and environmental mapping enthusiasts.<\/p>\n\n\n\n The example presented here\u2014the segmentation of corals from Sentinel-2 images\u2014illustrates the potential of these approaches for detailed analysis of coastal environments and the preservation of marine ecosystems.<\/p>\n\n\n\n The challenge is no longer just technical, but also collective: pooling models, documenting their inputs, sharing methods, and making deep learning more transparent, reproducible, and open.<\/p>\n\n\n\n The future of remote sensing will likely be built at the crossroads of these two worlds\u2014software engineering and field knowledge\u2014to transform satellite data into true indicators of ecological status.<\/p>\n\n\n\n \n","protected":false},"excerpt":{"rendered":"Deep learning is revolutionizing satellite image analysis. Long reserved for large laboratories or proprietary software, it is now becoming available to the wider world thanks to PyTorch and QGIS. This article explores the principles of Deep…<\/p>\n","protected":false},"author":1,"featured_media":15925,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[3857],"tags":[3873,3875,3877],"class_list":["post-16040","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-iaen","tag-deep-learning","tag-ia-en","tag-pytorch-en"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2025\/11\/deeplearningavant.jpeg?fit=272%2C198&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-4aI","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/16040","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=16040"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/16040\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/15925"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=16040"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=16040"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=16040"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}