Le portail francophone de la géomatique

ah merci beaucoup, je vais essayer ça, je suis pas expert en QGIS
Je vous mettrez au courant de la suite

Je pense qu'on peut aboutir à quelque chose en utilisant la théorie des graphes et en déterminant certaines caractéristiques inhérentes aux noeuds telles que leur centralité, etc...

Pour déterminer les polygones qui sont au bord, vous pourriez faire une intersection entre l'enveloppe linéaire de vos polygones et vos polygones.

Pour déterminer ceux qui sont au milieu, vous pourriez faire une matrice des distances : QGIS > Outils d'analyse et calculer la distance totale qui sépare un noeud de tous les autres. Et faire ça pour chaque noeud.
Voir http://en.wikipedia.org/wiki/Centrality et les mesures de la centralité. Ici, c'est l'indicateur de centralité de Friedman.

La centralité peut être aussi le nombre de voisins qu'un noeud a.

est il possible de prendre contact avec vous pour que ça soit plus claire pour moi ?
car je pense  que je me suis mal exprimé! les polygones dont je parle ce sont des bâtiments et quand je dis milieu c'est pas genre à l’intérieur du bloc, mais plutôt un polygone qui a un voisin à gauche à droite et par derrière (possible aussi), donc la notion de l'intersection n'est pas une solution à a ma requette car elle me sortira tous les polygones périphériques et dans ce cas là tous les bâtiments car ya pa de bâtiment à l’intérieur des bloc! tous les bâtiments ont une façade, dommage je peux pas joindre une figure pour être vraiment explicite...

J'ai essayé ce que vous m'avez dit mais malheureusement, il fait pas de différence entre polygone milieu et polygone coin ou angle, car tous les polygone touche le polygone enveloppe, j'ai essayé  de calculer la distance minimale entre les nœuds du réseau routier (intersection) et les polygones pour reconnaître ceux formant le coin ou angles, mais malheureusement j'ai pas une couche route complète (mise à jour)..

voila, les polygones représentent des batiments.

Dernière modification par sed.recherche (Mon 24 November 2014 16:30)

Oula, c'est un véritable défi, votre projet.

On peut considérer qu'un bâtiment mitoyen de deux et seulement deux autres sur deux faces orthogonales est dans un coin.
Du coup, je déterminerais les faces qui touchent d'autres voisines pour chaque bâtiment et calculerais leur angle.
Si l'angle est proche de 0, bâtiment enserré au milieu de deux autres.
Si proche de 90, alors dans un coin.

Concrètement, sous QGIS :
- Transformer polygones en lignes
- Eclater les lignes
- voir celles qui touchent d'autres (sans doute nécessité de programmation)
- Déterminer l'angle des faces d'un bâtiment : fonction xat et trigonométrie : arctan

Ou sinon, on peut considérer qu'un bâtiment est dans un coin s'il est distant du milieu du groupe de bâtiments auquel il appartient. Du coup, je proposerais ceci :
- créer les groupes de bâtiments par tampons érodés ou méthode sus-citée sous R
- rasteriser le centroïde du groupe
- créer raster de distance
- calculer pour chaque bâtiment la valeur maximale de distance au centre du groupe : statistiques de zone
- l'écart par rapport à la moyenne, l'utilisation d'un écart type donnerait un score.

Bon, tout ça, c'est du bricolage et cela mériterait d'être mis en oeuvre effectivement afin de vérifier les résultats et, si besoin, optimiser la méthode.

Dans tous les cas, c'est loin d'être simple à mettre en oeuvre..

Quoi qu'il en soit, vous pourriez déjà déterminer si un bâtiment est au bord ou pas selon la méthode que je vous avais décrite précédemment.

En fait il y a un problème à considérer, quand vous parlez  des deux faces orthogonales des deux bâtiments mitoyens, vous considérez que les bâtiments sont des rectangles ! sauf qu'en réalité tous les polygones ont une forme aléatoire!!

J'ai déjà réfléchi en façades consécutives et angles! mais elle est vraiment limitée.. car souvent on peut pas généraliser sur le reste et en faire un algorithme pour toute la problématique...

En fait j'ai une question, dans la quelle je bloque vraiment, je vois que beaucoup de solutions s'offre à nous pour trouver ce que je cherche mais souvent, il s'agit de tracer l'enveloppe extérieur des polygones ayant des faces en commun ! et donc ma question est:

Comment tracer cette enveloppe sous ArcGIS, sachant que le script R dont vous m'avez parler ça n'a pas fonctionné, ni sur Arcgis ni QGIS, y-a t il une autre alternative pour tracer les enveloppes des polygones connectés ensemble ? (concrètement, ces polygones connectés représentent les blocs de bâtiments "agrégat").

Merci beaucoup de votre aide.

ah d'accord je vois,

Est il possible de me dire concrètement ce que  je suis sensé faire sur QGIS pour configurer Sextante (les étapes et les commandes à utiliser)?

J'ai pu avancé dans ma requette, j'ai opté pour l'option de la distance des polygones par rapport aux centroides des enveloppes, quand je les classes par ordre décroissant y a un autre détail qui s'impose! combien faut t'il considérer de polygones dans le classement ? (les plus distant sont les bâtiments coin ou angle)

quand on a une enveloppe de 4 coins, dans ce cas là il faut considérer les 4 premiers!! dans une enveloppe triangle il y a trois coins donc on ne considère que les 3 premiers! comment automatiser ce calcul sachant que j'ai 29 000 enveloppe avec des formes très variables ?

Dernière modification par sed.recherche (Mon 01 December 2014 13:07)