Discussion :

[ben53]

Bonjour à tous,
Je récupère des données vecteur de bâtiments à travers les géoservices. Dans les propriétés des features renvoyés, on trouve une donnée nommée "hauteur", qui semble correspondre à la hauteur des bâtiments.
Première question: Même si les intitulés semblent évidents, existe t-il un document/url/requête qui explique à quoi correspond chaque propriété d'une feature?
Ensuite, on constate que certains bâtiments ont une hauteur nulle. Est-ce normal?
Enfin, certains bâtiments récupérés se trouvent aux abords de la BBOX, mais en dehors. Comment faire pour n'avoir que des bâtiments strictement situés à l'intérieur de la BBOX?

[azarzelli]

Bonjour,

Vous trouverez la documentation de la BD TOPO ici : https://geoservices.ign.fr/ressource...BDTOPO_3-0.pdf
L'attribut qui vous intéresse est décrit page 67. Concernant les bâtiments ayant une hauteur nulle, la valeur "0" signifie que la hauteur n'est pas connue ou est non significative. (cela est écrit page 67).

Concernant votre problème, cela est dû au fait que les bâtiments récupérés sont ceux qui intersectent (même partiellement) la bbox. Pour obtenir ceux qui sont strictement contenus à l'intérieur, je vous conseille de faire un post-traitement utilisant par exemple la méthode Within() de GDAL/OGR : https://gdal.org/api/vector_c_api.ht...12OGRGeometryH

Cordialement,

[ben53]

Merci pour votre réponse, et désolé de ne pas avoir trouvé moi-même la documentation, ce n'était pourtant pas difficile... ^^'
Concernant les bâtiments, je comprends et je trouve normal que ceux qui intersectent la BBOX soient récupérés, mais mon problème concerne des bâtiments qui sont entièrement en dehors de la bbox.
Par exemple, avec la requête suivante:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

https://wxs.ign.fr/CLEF/geoportail/wfs?SERVICE=WFS&VERSION=2.0.0&request=GetFeature&typename=BDTOPO_BDD_WLD_WGS84G:bati_indifferencie&outputFormat=application/json&srsname=EPSG:2154&BBOX=1191000,6081000,1192000,6082000,EPSG:2154

La feature n° 33 contient des coordonnées qui se situent à presque 40 mètres en dehors de la BBOX.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
33	
type	"Feature"
id	"bati_indifferencie.34151853"
geometry	
type	"MultiPolygon"
coordinates	
0	
0	
0	
0	1191219.0993895
1	6080958.60429742
2	7.1
1	
0	1191215.79990629
1	6080962.90451934
2	7.1
2	
0	1191210.59976272
1	6080959.30469827
2	7.1
3	
0	1191213.49970165
1	6080954.80490138
2	7.1
4	
0	1191219.0993895
1	6080958.60429742
2	7.1

Ce n'est pas un cas isolé: Pour cette requête, il y en a une vingtaine.

[elias couppe]

Bonjour,

le problème a été transmis aux équipes compétentes par azarzelli. Nous sommes toujours en attente de leur retour.

Nous ne manquerons pas de mettre à jour ce thread quand nous aurons de nouveaux éléments de réponse.

[ben53]

Merci!

[elias couppe]

Bonjour,

Les données sont stockées en BDD en WGS84G, donc les étapes qui sont faites lors de votre requête sont les suivantes :
1 - reprojection de la BBOX Lambert 93 en WGS84G
2 - intersection entre la BBOX et l'objet en WGS84G
3 - renvoi de la géométrie projetée en L93


En pratique :
Prenons une requête de BBOX en L93 : BBOX=1191000,6081000,1192000,6082000,EPSG:2154

Reprojection de la BBOX en WGS84G : POLYGON((8.88608625122707 41.6654208708387,8.88697834843709 41.6743768133159,8.89892360668584 41.6737071903522,8.89802971279575 41.6647513565942,8.88608625122707 41.6654208708387))

Enveloppe en WGS84G : BBOX(8.88608625122707 41.6647513565942,8.89892360668584 41.6743768133159) => requête en BDD

Du coup, si l'on prend un objet stocké en BDD en WGS84G, par exemple : MULTIPOLYGON(((8.88866614 41.66490352,8.88863057 41.66494424,8.88856525 41.66491548,8.88859587 41.66487324,8.88866614 41.66490352))) => la requête d'intersection (intersects) renvoie true !

Le résultat de la requête WFS est finalement un objet en L93

--> C'est un défaut du WFS : la BBOX de recherche est augmentée quand il y a une reprojection, cela permet d'optimiser les recherches en utilisant la projection native

La requête :
https://wxs.ign.fr/MACLEF/geoportail...3159,EPSG:4326 renvoie bien le polygone

La requête :
https://wxs.ign.fr/MACLEF/geoportail...3159,EPSG:4326 ne renvoie pas le polygone

Pour ne pas faire de reprojection de la BBOX, la requête suivante fonctionne aussi : https://wxs.ign.fr/MACLEF/geoportail...000 6081000)))