datagistips

load spatial;

-- point (adresse)
-- Ex. 34+bis+avenue+philippe+solari
-- Utilisation de la BAN
create or replace table point as (
	with a as (select unnest(features, recursive:=true) from read_json_auto(format('https://api-adresse.data.gouv.fr/search/?q={}', '34+bis+avenue+philippe+solari')))
	select st_point(coordinates[1], coordinates[2]) as geom from a
);

datagistips

start <- "2.337306%2C48.849319"
end <- "2.367776%2C48.852891"
library(jsonlite)

get_route <- function(start, end) {
  url <- glue("https://wxs.ign.fr/calcul/geoportail/itineraire/rest/1.0.0/route?resource=bdtopo-osrm&start={start}&end={end}&profile=car&optimization=fastest&constraints=%7B%22constraintType%22%3A%22banned%22%2C%22key%22%3A%22wayType%22%2C%22operator%22%3A%22%3D%22%2C%22value%22%3A%22autoroute%22%7D&getSteps=true&getBbox=true&distanceUnit=kilometer&timeUnit=hour&crs=EPSG%3A4326")
  j <- read_json(url)
  coords <- j$geometry$coordinates
  out <- lapply(coords, function(elt) {
    long <- elt[[1]]

datagistips

# Auparavant, on trouvait facilement les emplacements des mairies depuis Le Répertoire Géographique des Communes
# https://geo.data.gouv.fr/fr/datasets/2d67a8d0da4ba99004a88cbe93fc0e0fceb43bb8
# La donnée n'est plus disponible
# Mais l'API geo.api.gouv.fr permet de retrouver l'emplacement des mairies.
# Cependant, il y a des soucis avec la requête https://geo.api.gouv.fr/communes?geometry=mairie&format=geojson
# qui retourne 'Bad request'
# On utilise ici une itération le long des différents départements
# On exporte le résultat en GéoPackage

library(sf)

datagistips

library(tidyverse)
library(jsonlite)
library(glue)
library(sf)
library(leaflet)

make_url <- function(cityname) {
  cityname <- URLencode(cityname)
  glue("https://nominatim.openstreetmap.org/search?city={cityname}&format=geojson&polygon_geojson=1")
}

datagistips

library(shiny)
library(banR) # Pour le géococage et l'appel à l'API BANO
library(leaflet) # Pour la carte interactive

# ui est la partie affichage
ui <- fluidPage(
  
  fluidRow(column(12, 
                  h1("Géocodage d'adresse BANO"),
                  textInput("adresse", label = NULL, "", placeholder = "Tapez une adresse"),

datagistips

library(magick)
library(rjson)
library(glue)

birdseye <- function(coords, key) {
  
  # This script is an R adaptation from python # https://gist.github.com/WvanWaas/b0e60d88f22d13a9d431a53208c8ed7c
  # You must first have a Bing Maps account : https://learn.microsoft.com/en-us/bingmaps/getting-started/bing-maps-dev-center-help/getting-a-bing-maps-key
  
  # Get latlon as text

datagistips

// 1. Go to Google Drive
// 2. New > Google Apps Script > Name it 'Auto Google Form - Schema'
// 3. Copy-Paste below code
// 4. Select 'process' function next to 'Debugging' button
// 4. Save
// 5. Press 'Execute'

function getItem(field) {

  var s1 = field.name + ' ('+field.title+')';

datagistips

Bing Maps Birds Eye

• Aller sur https://www.bingmapsportal.com/Account/Register
• Créer le compte my account
• Créer la clé <maClé>
• La copier <maClé>

Tester la couverture

https://www.bing.com/api/maps/sdk/mapcontrol/isdk/birdseyev2

datagistips

library(jsonlite)
library(glue)
library(leaflet)
library(magrittr)
library(htmltools)

# Teste si la chaîne est une référence cadastrale
isParcelle <- function(s) {
  grepl("^([013-9]\\d|2[AB1-9])\\d{3}(0|[A-Z])[A-Z][0-9]{4}[a-z]?$", s)
}

datagistips

library(jsonlite)
library(sf)
library(leaflet)
library(glue)
library(magrittr)

getContourCommune <- function(code_insee) {
  # geo.api.gouv.fr arguments
  # geometry = contour pour avoir le contour
  # fields = geometry car seule la géométrie nous intéresse