2010/11/25 sly (sylvain letuffe) <[email protected]> > > > !!ATTENTION!! c'est du Émilie, donc c'est, comment dire, technique, > complexe, > dur à comprendre... mais efficace ! >
Je ne sais pas si je dois le prendre comme un compliment ou comme une insulte :P > > Je te résume l'idée, car tu peux tout aussi bien tenter dans ton coin, sur > le > même principe mais ne pas t'embrouiller avec des requêtes longuent comme un > anaconda : > > Pour la commune X : > 1- tu importes (osm2pgsql) le fichier osm du bati dans des tables à part > (cadastre_**) > > 2- tu importes (osm2pgsql) les batiments de la zone en provenance d'osm > dans > d'autres tables (communeX_**) (attention, il y a building=* mais aussi > place_of_workship, certaines amentity, etc.) > > 3- Tu consitues la requêtes de la mort qui va supprimer (ou marquer comme > non > importable) des tables cadastre_** tous les bâtiments dont il existe une > surface non nulle d'intersection avec ceux de osm. La requette doit faire > usage de > ST_AREA() (obtenir la mesure de la surface) > ST_Intersection(table1.gemometrie1,table2.geometrie2) (obtenir la surface > d'intersections entre deux géométries des deux tables) > et autres joyeusetés > > 4- ensuite tu utilises pgsql2shp pour sortir au format shapefile le contenu > de > ta table nettoyée genre : > $ pgsql2shp -f batiments_qui_superposent_pas.shp "select * from > cadastre_polygone where supperpose=0" > > 5- avec ogr2osm, tu convertis tes fichiers shapefiles en fichiers osm avec > les > bon tags dont tu as besoin, exemple d'utilisation de ogr2osm ici : > > http://wiki.openstreetmap.org/wiki/WikiProject_France/Parcs_nationaux_et_régionaux,_réserves_naturelles/Import_des_données_INPN<http://wiki.openstreetmap.org/wiki/WikiProject_France/Parcs_nationaux_et_r%C3%A9gionaux,_r%C3%A9serves_naturelles/Import_des_donn%C3%A9es_INPN> > tout en bas > Oui cette méthode marche parfaitement. J'ai même commencé à l'implémenter dans Postgis, mais le manque de temps récurrent chez moi m'a empêché de continuer. De plus, j'avais commence à écrire un programme qui fait cela sans avoir besoin de Postgis ou d'une base de donnée et qui travaille à partir de XAPI (au moins sur le papier). L'idée est à terme de faire cela chaque jour pour détecter les changement dans les bâtiments sur le cadastre et de gérer la maintenance. Si besoin est pour cette méthode, je serais plus qu'heureuse d'aider. Emilie Laffray
_______________________________________________ dev-fr mailing list [email protected] http://lists.openstreetmap.org/listinfo/dev-fr
