Skip to content

Instantly share code, notes, and snippets.

@esgn
Last active December 20, 2023 13:59
Show Gist options
  • Save esgn/4bbf298ad76f4d72e9f3c133cbc96cf1 to your computer and use it in GitHub Desktop.
Save esgn/4bbf298ad76f4d72e9f3c133cbc96cf1 to your computer and use it in GitHub Desktop.
Extraction depuis le S3 LIDAR HD

Extraction depuis le S3 LIDAR HD

Question pratique: Comment extraire un morceau de LIDAR HD depuis ces 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"
]
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment