Le portail francophone de la géomatique

Bonjour,

A présent que vous avez tous les champs dans une seule table, la table carreaux, je présume, utilisez la calculatrice de champ pour calculer la colonne IMP

bien cordialement

Bonjour,

J'avais, il y a quelques temps, calculé un indice de densité bâtie sur la ville de paris (visible en ligne) - probablement avec QGIS 1.8.

D'un point de vue méthodologique, ca ressemble en tout point à ton IMP. J'avais du convertir mes données shapefiles en PostGRESQL pour effectuer mes calculs.

Je doute que la version 2.2 de QGIS soit assez puissante pour effectuer le traitement. La syntaxe native de QGIS étant encore un peu faiblarde.

Bon finalement j'ai trouvé comment faire (même si je ne comprend toujours pas pourquoi ce message d'erreru s'affiche). Je poste ici la démarche complète pour les intéressés :

CALCUL DE L’IMP SOUS QGIS

DONNEES :

BATI (BDPARCELLE)
PARCELLE (BDPARCELLE)

METHODO :

Créer les centroïdes pour la couche BATI.

Avec le plugin QMarxan, créer une grille avec des carreaux de 200m de coté (4ha). Puis, avec le plugin QMarxan, joindre le nombre de bâtiments de la couche BATI sur la couche GRILLE :
-    NBR_BATI

Réaliser une requête pour ne sélectionner que les carreaux avec 5 bâtiments ou plus.

Avec l’outil requête spatiale, sélectionner les parcelles bâties se situant dans la grille :
-    Créer nouvelle sélection des parcelles se situant dans l'emprise de la GRILLE : BD PARCELLE intersecte GRILLE.
-    Retirer les parcelles non bâtie de la sélection : BD PARCELLE disjoint BATI
-    Créer la couche : PARCELLE_BATIE

Inverser la sélection (PLUGIN SELECTPLUS) et avec l’outil requête spatiale, sélectionner les parcelles non bâties se situant dans la grille :
-    Retirer de la sélection : BD PARCELLE disjoint GRILLE
-    Créer la couche : PARCELLE NON BATIE

Réaliser l’intersection des couches PARCELLE BATIE (découpée) et BATI (découpage).

Avec le plugin QMarxan, joindre les surfaces des 3 couches crées précédemment sur la couche grille :
-    PARC_BATI
-    PARC_N_BATI
-    AIRE_BATI

Rendre la couche GRILLE éditable et créer une nouvelle colonne IMP complétée avec la formule suivante :

(surface emprise des bâtiments / surface parcelles bâties) + (1 - (surface des parcelles non bâties / surface du carreau))

Soit ici :

(AIRE_BATI / PARC_BATI) + (1 - (PARC_N_BATI/ 40000))

Réaliser une requête pour ne sélectionner que les carreaux avec 5 bâtiments ou plus.

Légende, hiérarchisée, 4 classes, mode quantiles

Et voilà !

Dernière modification par vince100 (Tue 08 April 2014 14:54)

J'ai corrigé un peu la méthodo ce week-end. Evidemment, la démarche est plus précise avec les données BDPARCELLE mais ici je n'avais pas l'ensemble des données

Maintenant j'aimerai savoir si certains d'entre vous ont déjà travaillé avec le carroyage de l'INSEE ? J'aimerai adapter cette méthode en croisant les données IGN et INSEE notamment pour calculer la densité résidentielle.

@ Dof : oui tout à fait. Une proposition/estimation pour effectuer le calcul que tu proposes ? Combien de personnes par mètre cube dans les bâtiments résidentiels ?

Du coup est ce que quelqu'un a le carroyage de l'INSEE en shape ? J'ai essayé de récupérer les fichiers depuis le site de l'INSEE mais l'ordi du boulot rame lorsque je veux les exploiter. L'idéal serait effectivement de travailler sur les fichiers MAJIC (même si leur exploitation ne me semble pas très aisée) mais à l'échelle d'un SCoT je pense qu'une certaine approximation peut être tolérée.

Bonjour,

Je ne suis pas très familier de la thématique traitée mais pourquoi ce choix de 5 bâtiments, sachant que notamment dans les lieux où on a des immeubles, on peut avoir une seule géométrie pour le bâti mais qui corresponde à plusieurs portes d'entrée (et avec un son grand jardin).
Et si on reste dans la restriction que vous avez retenue, il n'est pas bizarre le carreau blanc en plein milieu des bleus, à l'est?

Puis, avec le plugin QMarxan, joindre le nombre de bâtiments de la couche BATI sur la couche GRILLE :
-    NBR_BATI



Comment faites-vous cela ? Je suis bloquée à cette étape.....