Différences

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

--- informatique:ign_bdortho [22/02/2026 14:06] – [Indre-et-Loire (37)] cyrille
+++ informatique:ign_bdortho [23/02/2026 20:54] (Version actuelle) – [Outils] cyrille
@@ Ligne 12: / Ligne 12: @@
   * Décompressés on obtient 302 images JP2 de 94 Mo chacune soit ~29 Go
     * 1 image découpée 640 pixels fait ~1520 fichiers pour ~70 Mo
-    * au final 459 342 fichiers
+    * au final 459 342 fichiers pour 24 Go
 ===== Outils =====
@@ 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
 import glob
+import shutil
 from rasterio.windows import Window
 from PIL import Image
@@ Ligne 34: / Ligne 34: @@
 from tqdm import tqdm
 from concurrent.futures import ThreadPoolExecutor, as_completed
 def decouper_fichier_jp2(fichier, dossier_sortie, taille=640):
     nom_base = os.path.splitext(os.path.basename(fichier))[0]
     sous_dossier = os.path.join(dossier_sortie, nom_base)
     os.makedirs(sous_dossier, exist_ok=True)
+    # Copie du fichier .tab correspondant
+    fichier_tab = fichier.replace('.jp2', '.tab')
+    if os.path.exists(fichier_tab):
+        shutil.copy(fichier_tab, os.path.join(sous_dossier, f"{nom_base}.tab"))
+    # Découpe le fichier JP2 en tuiles JPG
     with rasterio.open(fichier) as src:
         h, w = src.shape
@@ Ligne 44: / 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 i in range(0, h, taille):
+            for y in range(0, h, taille):
-                for j in range(0, w, taille):
+                for x in range(0, w, taille):
-                    window = Window(j, i, taille, taille)
+                    window = Window(x, y, taille, taille)
-                    if i + taille <= h and j + taille <= w:
+                    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)
-                        img.save(f"{sous_dossier}/{nom_base}_x{i}_y{j}.jpg")
+                        img.save(f"{sous_dossier}/{nom_base}_x{x}_y{y}.jpg")
                         pbar.update(1)
 def decouper_jp2_en_tuiles_parallele(dossier_entree, dossier_sortie, taille=640, nb_threads=4):
     os.makedirs(dossier_sortie, exist_ok=True)
     fichiers = glob.glob(os.path.join(dossier_entree, "*.jp2"))
+    #fichiers = glob.glob(os.path.join(dossier_entree, "37-2025-0505-6715-LA93-0M20-E080.jp2"))
     with ThreadPoolExecutor(max_workers=nb_threads) as executor:
         futures = [executor.submit(decouper_fichier_jp2, fichier, dossier_sortie, taille) for fichier in fichiers]
@@ Ligne 63: / 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/",
+    decouper_jp2_en_tuiles_parallele(
- "tuiles_640/",
+        "BDORTHO_2-0_RVB-0M20_JP2-E080_LAMB93_D003_2025-01-01/ORTHOHR/1_DONNEES_LIVRAISON_2026-01-00061/OHR_RVB_0M20_JP2-E080_LAMB93_D03-2025/",
- nb_threads=4
+        "BDORTHO_tuiles-640_03/",
-)
+        nb_threads=4,
+    )
 </code>