Extraction depuis le S3 LIDAR HD

Question pratique: Comment extraire un morceau de LIDAR HD depuis les dalles COPC mises à disposition par l'IGN en utilisant 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 nuage de points directement sans télécharger toute la donnée en local au préalable.

Ce type de pipeline PDAL fonctionne de la manière souhaitée (pas de téléchargement complet préalable):

[
  {
    "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 pipeline PDAL ne fonctionne pas de la manière souhaitée (ici, toute la donnée est chargée en local au préalable):

[
  {
    "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"
]

A noter qu'à moins de disposer d'un VPC en entrée il faut préciser à PDAL la liste des dalles à partir desquelles extraire la donnée finale souhaitée. Ici, les flux WFS IGN suivant peuvent s'avérer utiles :

Ensemble des blocs LIDAR HD disponibles : https://data.geopf.fr/private/wfs/?SERVICE=WFS&VERSION=2.0.0&APIKEY=interface_catalogue&REQUEST=GetFeature&TYPENAMES=ta_lidar-hd:bloc&OUTPUTFORMAT=application/json
Ensemble des dalles LIDAR HD disponibles (exemple de filtrage par emprise) : https://data.geopf.fr/private/wfs/?SERVICE=WFS&VERSION=2.0.0&APIKEY=interface_catalogue&REQUEST=GetFeature&TYPENAMES=ta_lidar-hd:dalle&OUTPUTFORMAT=application/json&BBOX=881000,6220000,936000,6283000,EPSG:2154

esgn/poc-flux-lidar-copc.md