esgn

-- Activation des logs HTTP
SET logging_storage = 'memory';
CALL enable_logging(
    'HTTP',
    storage = 'memory',
    storage_buffer_size = 0
);

-- On affiche le nombre de threads pour contrôle
SELECT current_setting('threads') AS threads;

esgn

Test geoparquet BDTOPO avec DuckDB

Maintenant que le bug sur les headers manquant est corrigé, je reteste à nouveau.
Je teste avec DuckDB v1.5.5 (Variegata) sous Ubuntu 24.04.5 LTS

Version des données

Deux versions sont disponibles sur https://data.geopf.fr/chunk/telechargement/ : 15/03/2026 et 15/06/2026

• https://data.geopf.fr/chunk/telechargement/download/BDTOPO_PQT/BDTOPO_3-5_TOUSTHEMES_GEOPARQUET_WGS84G_FRA_2026-03-15/batiment.parquet

esgn

import sys
import argparse
import geopandas as gpd
import pandas as pd
import json

def parse_args():
    parser = argparse.ArgumentParser(
        "Create JSON PDAL pipelines to download point cloud extract")
    parser.add_argument("--input_footprints", "-i",

esgn

How to push a Docker image to Google Artifact Registry with a GitHub Action without using a Service Account

This how I managed to push a Docker image in a Google Cloud Artifact Registry using a GitHub Actions workflow without using a GCP Service Account.

Before using this tutorial make sure you have :

• A GitHub repository containing a Dockerfile
• A GCP project with sufficient permissions to make changes to it
• The gcloud tool installed, authenticated and configured to use the GCP project. We will not use the --project flag in the following gcloud commands.

Create an Artifact Registry of Docker type

esgn

resource "google_storage_bucket" "osm_dl_planet" {
  project                     = var.project-id
  name                        = var.bucket-name
  uniform_bucket_level_access = true
  location                    = var.project-region
  storage_class               = "REGIONAL"
  force_destroy               = true
}

resource "google_service_account" "scheduler" {

esgn

#!/usr/bin/env bash

# Working area in EPSG:2154
xmin=892000
ymin=6249000
xmax=894000
ymax=6247000

rm -rf footprints
mkdir footprints

esgn

#!/usr/bin/env python3

import os
import laspy
import shutil

file_ext = ".las"
output_file = "output.las"

if os.path.exists(output_file):

esgn

import sys

try:
    from osgeo import ogr, osr, gdal
except:
    sys.exit('ERROR: cannot find GDAL/OGR modules')

# create the spatial reference
srs = osr.SpatialReference()
srs.ImportFromEPSG(2154)

esgn

Extraire une zone du LiDAR HD à distance avec PDAL et PDAL Wrench

À partir d’une emprise exprimée en Lambert-93, ce tutoriel montre comment extraire depuis la Géoplateforme uniquement les points LiDAR HD nécessaires, même lorsque la zone couvre plusieurs dalles ou plusieurs blocs, sans télécharger intégralement les fichiers COPC sources en local.

Trouver les dalles LiDAR HD concernées sur la zone

Passage par un VPC (Virtual Point Cloud)

Le service de téléchargement de la Géoplateforme propose des fichiers VPC pour chaque bloc de LiDAR HD. Exemple : https://data.geopf.fr/chunk/telechargement/download/LiDARHD-NUALID/NUALHD_1-0__LAZ_LAMB93_PQ_2025-03-14/index.vpc

esgn

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/IO/read_points.h>
#include <CGAL/property_map.h>
#include <CGAL/Surface_mesh.h>
#include <CGAL/Shape_detection/Efficient_RANSAC.h>
#include <CGAL/Polygonal_surface_reconstruction.h>
#include <CGAL/Polygon_mesh_processing/triangulate_faces.h>
#include <CGAL/Polygon_mesh_processing/orientation.h>
#include <CGAL/Polygon_mesh_processing/polygon_mesh_to_polygon_soup.h>
#include <CGAL/Polygon_mesh_processing/polygon_soup_to_polygon_mesh.h>