Différences

Ci-dessous, les différences entre deux révisions de la page.

--- informatique:ign_bdortho [23/02/2026 11:55] – [Outils] cyrille
+++ informatique:ign_bdortho [08/03/2026 13:05] (Version actuelle) – [Outils] cyrille
@@ Ligne 19: / Ligne 19: @@
 <code python>
-'''
+"""
-Script généré par https://chat.mistral.ai
 Installation:
 pip install rasterio pillow numpy tqdm
+"""
-'''
 import rasterio
 import os
@@ Ligne 35: / Ligne 32: @@
 from tqdm import tqdm
 from concurrent.futures import ThreadPoolExecutor, as_completed
+import argparse
+import sys
 def decouper_fichier_jp2(fichier, dossier_sortie, taille=640):
@@ Ligne 40: / Ligne 39: @@
     sous_dossier = os.path.join(dossier_sortie, nom_base)
     os.makedirs(sous_dossier, exist_ok=True)
-    # Copie du fichier .tab correspondant
+    # Copie du fichier .tab
     fichier_tab = fichier.replace('.jp2', '.tab')
     if os.path.exists(fichier_tab):
@@ Ligne 50: / Ligne 49: @@
         nb_tuiles_w = w // taille
         total_tuiles = nb_tuiles_h * nb_tuiles_w
+        # réserve un buffer pour limiter une allocation à chaque tuile.
+        buffer = np.empty((src.count, taille, taille), dtype=np.uint8)
         with tqdm(total=total_tuiles, desc=f"Traitement de {nom_base}") as pbar:
             for y in range(0, h, taille):
@@ Ligne 55: / Ligne 56: @@
                     window = Window(x, y, taille, taille)
                     if y + taille <= h and x + taille <= w:
-                        tile = src.read(window=window)
+                        tile = src.read(window=window, out=buffer)
                         tile_rgb = np.moveaxis(tile, 0, -1).astype(np.uint8)
                         img = Image.fromarray(tile_rgb)
@@ Ligne 70: / Ligne 71: @@
             future.result()  # Attend la fin de chaque tâche
-# Exemple d'utilisation
+if __name__ == "__main__":
-decouper_jp2_en_tuiles_parallele(
-  "/BDORTHO_2-0_RVB-0M20_JP2-E080_LAMB93_D037_2025-01-01/ORTHOHR/1_DONNEES_LIVRAISON_2025-12-00193/OHR_RVB_0M20_JP2-E080_LAMB93_D37-2025/",
+    parser = argparse.ArgumentParser(description="Découp les images BDOrtho IGN en tuiles de 640 pixels.")
-  "/BDORTHO_tuiles_640/",
+    parser.add_argument("folder_in", help="Dossier des images à traiter (obligatoire)")
-  nb_threads=4
+    parser.add_argument("folder_out", help="Dossier de sorties pour les tuiles découpées (obligatoire)")
-)
+    parser.add_argument("--threads", type=int, default=4, help="Combien de threads pour la découpe(défaut: 4).")
+    args = parser.parse_args()
+    if not os.path.isdir(args.folder_in):
+        print(f"Error: Folder '{args.folder_in}' not found.")
+        sys.exit(1)
+    decouper_jp2_en_tuiles_parallele(
+        #"/media/cyrille/Seagate Backup Plus Drive/BDORTHO_2-0_RVB-0M20_JP2-E080_LAMB93_D032_2025-01-01/ORTHOHR/1_DONNEES_LIVRAISON_2025-10-00030/OHR_RVB_0M20_JP2-E080_LAMB93_D32-2025/",
+        #"/mnt/nas/IGN_BDOrtho/BDORTHO_2-0_RVB-0M20_JP2-E080_LAMB93_D005_2025-01-01/ORTHOHR/1_DONNEES_LIVRAISON_2026-01-00129/OHR_RVB_0M20_JP2-E080_LAMB93_D05-2025/",
+        args.folder_in,
+        #"/mnt/nas/IGN_BDOrtho/BDORTHO_tuiles-640_05/",
+        args.folder_out,
+        nb_threads=args.threads,
+    )
 </code>