Discussion :

[lalie73]

Bonjour,
je voulais savoir quel est le système de projection utilisé ds le geoportail (j'obtiens des points du style 440457.2501299669 5102995.218981583), car j'ai des différences d'estimation de surface entre la commande getArea du polygone et une surface obtenue avec un SIG ou avec CartoExplorer.
Pour info, mes polygones sont stockés en base de donnée mysql au format wkt.
A l'origine, ils ont été tracés en WGS84.
Merci

[dgrichard]

Bonjour,
je voulais savoir quel est le système de projection utilisé ds le geoportail (j'obtiens des points du style 440457.2501299669 5102995.218981583), car j'ai des différences d'estimation de surface entre la commande getArea du polygone et une surface obtenue avec un SIG ou avec CartoExplorer.
Pour info, mes polygones sont stockés en base de donnée mysql au format wkt.
A l'origine, ils ont été tracés en WGS84.
Merci

Le système de coordonnées est GEOPORTALxxx ou xxx représente le territoire d'application. Grosso modo, c'est du WGS84 géographique rendu conforme par un cos(latitude_d_origine) (donc, qui conserve les angles : les ronds restent "rond"). Plus 'dinfo là

Les paramètres sont là.

Pour la mesure de surface, l'utilisation de l'outil "Mesure" permet d'avoir une idée (Cf. cet exemple et cet outil : ).

[lalie73]

Merci :-)
J'ai effectivement ajouté l'outil de mesure "surface" dans la boite de contrôle "Outils" qui me donne une meilleure mesure. Pour exemple, un polygone qui fait 31,4177 ha dans mon logiciel sig : la fonction getArea d'openlayers me le donne à 31,0318 et la mesure manuelle (avec l'outil dédié) du geoportail 31,4614. Comment peut-on alors utiliser l'outil de mesure dédié de manière automatique ? j'entends par là obtenir l'aire d'un polygone directement (au click), sans avoir à le retracer.
Autre question, peut-on se servir de la

[lalie73]

oups ! fausse manip
autre question, peut-on utiliser la géométrie des contours de Commune de la couche administrative ? le but est de comparer la géométrie de mes polygones à celle des communes afin de déterminer s'ils s'intersectent

[dgrichard]

Comment peut-on alors utiliser l'outil de mesure dédié de manière automatique ? j'entends par là obtenir l'aire d'un polygone directement (au click), sans avoir à le retracer.

Si l'objet à mesurer est un objet vectoriel (héritant de OpenLayers.Feature.Vector), il suffit de d'utiliser :

Code :

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

1
2
objet.geometry.getArea();

Sinon (image de fond), pas d'autre solution que de « numériser »

[dgrichard]

autre question, peut-on utiliser la géométrie des contours de Commune de la couche administrative ? le but est de comparer la géométrie de mes polygones à celle des communes afin de déterminer s'ils s'intersectent

Les contours sont des images (Cf. réponse précédente)

De plus, il faudrait comparer avec les bonnes géométries issues de la BDCarto ou de la BDTopo ...

[lalie73]

Code :

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

1
2
objet.geometry.getArea();

oui, c'est ce que je fais, mais je n'obtiens pas la "bonne" valeur (cf message précédent) il y a une erreur d'environ 1% sur un polygone de 30ha. L'idée était d'enregistrer automatiquement ds la bdd la surface issue du geoportail car je ne la récupère pas depuis le fichier shape originel

[dgrichard]

oui, c'est ce que je fais, mais je n'obtiens pas la "bonne" valeur (cf message précédent) il y a une erreur d'environ 1% sur un polygone de 30ha.

En comparant avec quoi (puisque tu ne récupères pas les surfaces du shapefile) ?
En quel système de coordonnées, les surfaces ont-elles été calculées ?

L'idée était d'enregistrer automatiquement ds la bdd la surface issue du geoportail car je ne la récupère pas depuis le fichier shape originel

Il me semble qu'un ogrinfo -al sur le shapefile donne les surfaces ...