Le portail francophone de la géomatique

Salut,

si votre version de postgis le permet vous pouvez
utiliser la fonction ST_MinimumBoundingCircle, au nom
assez explicite.

Ensuite avec la fonction st_Buffer, vous pourrez faire "gonfler"
vos cercles.

Oui, et comme d'habitude trop vite ,
j'avais compris que c'était le tronçon qui devait être entouré...

Dernière modification par tumasgiu (Thu 29 August 2019 11:13)

Si ce sont bien les adresses qui doivent etre entourées, voici une fonction plpython, appelable en SQL, qui geocode une serie d'adresses et renvoie un ensemble de docs JSON avec champs remplis par le geocodeur (meme structure que le doc renvoyé par le geocodeur):

• La fonction prend une requete SQL qui doit renvoyer les champs ids, voie, ville, cp: ids est un tableau des id à géocoder, histoire de pouvoir faire un update sur la table (tableau car dans mon cas, j'avais une table avec plusieurs adresses identiques et je en voulais pas appeler le geocodeur plusieurs fois pour la meme adresse)

• Elle construit le csv en mémoire et l'envoie au géocodeur, puis retourne un objet jsonb par adresse géocodée, cet objet contient les lat/long de l'adresse:

drop function geocode_add(query text);
create or replace function geocode_add(query text)
  RETURNS setof jsonb as
$$
import csv
import io
from collections import OrderedDict
import requests
import json

with io.BytesIO(b"") as stream:
    writer = csv.writer(stream, delimiter=',', quotechar='"', quoting=csv.QUOTE_MINIMAL)
    writer.writerow(["ids", "voie", "ville", "cp"])

    for row in plpy.cursor(query):
        # writes row in csv to csv stream
        writer.writerow([row['ids'], row['voie'], row['ville'], row['cp']])

    requests_session = requests.Session()
    kwargs = {
        'data': OrderedDict([
            ('columns', ['voie', 'ville', 'cp']),
            ('postcode', 'cp')
        ]),
        'method': 'post',
        'files': OrderedDict([
            ('data', ('addpg.csv',stream.getvalue()))
        ]),
        'stream': True,
        'url': 'https://api-adresse.data.gouv.fr/search/csv/'
    }

    response = requests_session.request(**kwargs)
    for row in csv.DictReader(response.text.encode('utf-8').splitlines()):
        yield json.dumps(row)

  $$ language plpythonu VOLATILE;

La fonction s'utilise comme ceci pour extraire les infos des documents retournés:

select t->>'ids' as ids, t->>'result_label' as adr,
       st_astext(st_transform(
    st_setSRID(
        st_makePoint((t->>'longitude')::float, (t->>'latitude')::float), 4326), 2154)) as geom
from geocode_add($$select '{1}'::text[] as ids, '2, place du capitole' as voie,  'TOulouse' as ville, '31000' as cp$$) as t;

ids    adr                                               geom
['1']    Place du Capitole 31000 Toulouse    POINT(574291.035317137 6279607.04961392)

Nico

Lorsque j'interroge un point adresse voici ce qui est renseigné dans le champ géométrie


SQL> select geometrie from num_voirie where cle_oracle='4680877';

GEOMETRIE(SDO_GTYPE, SDO_SRID, SDO_POINT(X, Y, Z), SDO_ELEM_INFO, SDO_ORDINATES)

--------------------------------------------------------------------------------

SDO_GEOMETRY(2001, 2154, SDO_POINT_TYPE(1024570,4, 6282014,4, NULL), NULL, NULL)


SQL
Faut il géocoder? il ne me semble pas
Merci

Non c'est bon, la table num_voirie contient bien une geom Oracle point.

Si vous voulez limiter les transferts de données, la conversion d'une bd à l'autre, etc., il existe un Foreign Data Wrapper Oracle (https://github.com/laurenz/oracle_fdw).

Il vous permettra de définir dans PG une table virtuelle pointant sur la table Oracle, vous pourrez alors faire des requetes sur cette table depuis Postgis.
Par exemple, pour créer dans Postgis une table des cercles des adresses:

create table voirie_cercle as 
select st_buffer(geometrie, 50) from num_voirie_fdw where cle_oracle='4680877';

ici, num_voirie_fdw est la FOREIGN TABLE pointant sur la table Oracle, qu'il faudra créer avant (cf doc FDW)

Nicolas