Extraction depuis le S3 LIDAR HD

poc-flux-lidar-copc.md

Extraction depuis le S3 LIDAR HD

Question pratique: Comment extraire un morceau de LIDAR HD depuis les dalles COPC avec PDAL ?

Extraction sur une emprise incluse dans une seule dalle

On streame la donnée pour n'en récupérer qu'un morceau sur une emprise (bounds)

[
  {
    "type":"readers.copc",
    "filename":"https://storage.sbg.cloud.ovh.net/v1/AUTH_63234f509d6048bca3c9fd7928720ca1/ppk-lidar/PQ/LHD_FXX_0890_6263_PTS_C_LAMB93_IGN69.copc.laz",
    "bounds":"([890250,890750],[6262250,6262750])"
  },
  "extract_on_single_tile.laz"
]

Extraction sur une emprise à cheval sur plusieurs dalles

La question ici est de savoir comment récupérer le morceau de nuages de points directement sans télécharger toute la donnée au préalable.

La solution suivante est fonctionnelle:

[
  {
    "type":"readers.copc",
    "filename":"https://storage.sbg.cloud.ovh.net/v1/AUTH_63234f509d6048bca3c9fd7928720ca1/ppk-lidar/PQ/LHD_FXX_0890_6263_PTS_C_LAMB93_IGN69.copc.laz",
    "bounds":"([890250,890750],[6262500,6263500])"
  },
  {
    "type":"readers.copc",
    "filename":"https://storage.sbg.cloud.ovh.net/v1/AUTH_63234f509d6048bca3c9fd7928720ca1/ppk-lidar/PQ/LHD_FXX_0890_6264_PTS_C_LAMB93_IGN69.copc.laz",
    "bounds":"([890250,890750],[6262500,6263500])"
  },
  "extract_on_two_tiles.laz"
]

Ce type de solution ne semble pas fonctionnel. Téléchargement des données avant le crop a priori.

[
  {
    "type":"readers.copc",
    "filename":"https://storage.sbg.cloud.ovh.net/v1/AUTH_63234f509d6048bca3c9fd7928720ca1/ppk-lidar/PQ/LHD_FXX_0890_6263_PTS_C_LAMB93_IGN69.copc.laz"
  },
  {
    "type":"readers.copc",
    "filename":"https://storage.sbg.cloud.ovh.net/v1/AUTH_63234f509d6048bca3c9fd7928720ca1/ppk-lidar/PQ/LHD_FXX_0890_6264_PTS_C_LAMB93_IGN69.copc.laz"
  },  
  {
    "type":"filters.crop",
    "bounds":"([890250,890750],[6262500,6263500])"
  },
  "extract_on_two_tiles_slow.laz"
]