En urbanisme, la densité d'un territoire doit être le reflet de la quantité d'humain qui y vivent par unité de surface. En particulier, elle doit permettre de mettre en lumière l'intensité d'usage du bâti et des services et réseaux qui y sont attachés. Dans un quartier dense, on aura plus d'habitants et donc besoin de plus d'écoles, d'une plus grande fréquence de ramassage des ordures et de plus larges canalisations. C'est un indicateur des politiques urbaines pour justifier la construction de nouveaux logements ou l'implantation de services.

Son calcul se fait généralement à partir de la surface d'une zone administrative (ville, arrondissement, quartier, IRIS) et de la population qui y réside. Le nombre de résidant provient des recensements qui s'intéressent uniquement au lieu de résidence et non au lieu de travail des individus, ce qui introduit immédiatement une imprécision importante pour étudier l'intensité d'usage car elle ommet les travailleurs. Un quartier d'affaire peut être considéré comme dense (en semaine) mais aura une faible population résidente. Ce n'est cependant pas ce problème qui m'intéresse aujourd'hui.

Vue du 18ème arrondissement avec ses zones non-habitées

Ce qui me semble plus gênant est l'usage de la surface de la zone administrative totale, dans le cas de Paris généralement celle de l'arrondissement. En effet, une part non-négligeable de la surface des arrondissements est occupée par des zones non-construites et non-habitées : périphérique, rails, terrains de sport, cimetières, ... Surtout, cette part est très inégalement repartie entre les arrondissements de Paris. Les arrondissements périphériques en possèdent généralement plus que ceux du centre.

Densité selon la surface totale de l'arrondissement, population 2015

Dans l'étude de la densité des arrondissements, celle-ci se retrouve donc sur-évaluée dans les arrondissements avec peu d'espace non-construits (3ème, 11ème) et sous-évaluée dans ceux accueillant de grandes infrastructures (18ème, 15ème) si l'on se base sur la surface totale de l'arrondissement. Voici une carte de la densité des arrondissements de Paris, basée sur ce paramètre. La population est celle du recensement de de 2010, compilées par l'APUR.

Pour obtenir une cartographie plus réaliste, je propose de calculer la densité à partir de l'emprise bâtie totale de l'arrondissement, c'est à dire la somme des surfaces de tous les bâtiments qui la compose. Cela permet d'obtenir de masquer le biais générées par les surfaces non construites inégales. Cela permet également de générer un indicateur représentant réellement « l'intensité d'usage humain du bâti ».

Densité selon l'emprise bâtie totale, population 2015

Quelques notes pour l'utilisation des matériaux PBR (Physically Based Rendering) dans V-Ray Next.

Pavés PBR avec V-Ray Next pour Rhino

Il existe deux modes : Specular et Metalness. On peut deviner auquel on a affaire en regardant le nom des couches (maps) du matériau. Si l'on voit une couche nommée « metal » ou « metallic », ce sera forcément le second mode.

Specular : partir d'un matériau générique (V-Ray BRDF). Les couches se placent comme ceci :

• Base color → Diffuse. C'est l'équivalent de la couleur.
• Specular (niveaux de gris) → Reflection Color et avec Fresnel coché. Réflectivité.
• Glossiness → Reflection Glossiness. Niveau de flou des reflets.
• Normal (RVB) → Bump mode Normal. Texture.
• Height (niveaux de gris) → Displacement mode Normal. Modification de la géométrie de la surface.

Metallic : partir d'un matériau métalique (V-Ray Metallic). Les couches se placent comme ceci :

• Base color → Color. C'est l'équivalent de la couleur.
• Metalness (niveaux de gris) > Metalness (Color Space en Rendering Space). Brillance.
• Roughness → Reflection Glossiness (Color Space en Rendering Space). Niveau de flou des reflets.
• Normal (RVB) → Bump mode Normal (Color Space en Rendering Space). Texture.
• Height (niveaux de gris) → Displacement mode Normal (Color Space en Rendering Space). Modification de la géométrie de la surface.
• IOR (Index of Refraction) : à laisser par défaut ou régler manuellement selon le matériau et la documentation.

La couche Ambient Occlusion doit être mélangée avec Diffuse/Color (ou une autre couche selon le cas). Il faut cliquer sur le pavé de texture puis séléctionner Mix (map) et placer les deux textures.

Les textures PBR étant généralement en grande résolution, il faut enfin s'assurer de ne pas limiter le moteur de rendu, voir Render Params > GPU Textures > Full size.

La loi relative à la solidarité et au renouvellement urbains du 13 décembre 2000, dite loi SRU, impose aux communes de disposer d'un taux minimum de logements sociaux. Ce taux originellement fixé à 20 % a été porté à 25 % en 2013. Les communes carencées sont soumises à une pénalité financière.

Les données compilées par Alexandre Léchenet à partir des balances comptables des communes permettent de générer les cartes du cumul des pénalités versées sur la dernière décennie, de 2010 à 2019.

Cumul des pénalités versées au titre de l'article 55 de la loi SRU par commune sur la période 2010–2019.

Voici le Top-20 :

1. Saint-Maur-des-Fossés (94068), 20 141 112 €
2. Le Cannet (06030), 10 910 556 €
3. Saint-Raphaël (83118), 5 726 747 €
4. Boulogne-Billancourt (92012), 5 685 269 €
5. Saint-Laurent-du-Var (06123), 4 404 465 €
6. Mandelieu-la-Napoule (06079), 4 375 592 €
7. Le Perreux-sur-Marne (94058), 4 245 494 €
8. Nogent-sur-Marne (94052), 3 888 770 €
9. Agde (34003), 3 842 009 €
10. L'Union (31561), 3 364 781 €
11. Courbevoie (92026), 3 349 247 €
12. Villeneuve-Loubet (06161), 3 341 522 €
13. Le Raincy (93062), 3 238 790 €
14. Aix-en-Provence (13001), 3 222 150 €
15. Grasse (06069), 3 202 149 €
16. Roquebrune-sur-Argens (83107), 2 803 872 €
17. Cabriès (13019), 2 637 965 €
18. Allauch (13002), 2 613 121 €
19. Menton (06083), 2 602 351 €
20. Six-Fours-les-Plages (83129), 2 573 161 €

En additionnant les pénalités cumulées des communes à l'intérieur de chaque département on obtient la carte suivante.

Cumul des pénalités versées au titre de l'article 55 de la loi SRU par département sur la période 2010–2019.

Avec un Top-10 :

1. Bouches-du-Rhône (13), 54 639 629 €
2. Alpes-Maritimes (06), 53 621 511 €
3. Val-de-Marne (94), 34 468 691 €
4. Var (83), 33 046 587 €
5. Hérault (34), 28 550 536 €
6. Yvelines (78), 27 668 725 €
7. Rhône (69), 23 778 458 €
8. Essone (91), 21 303 253 €
9. Loire-Atlantique (44), 16 058 016 €
10. Haute-Garonne (31), 14 393 730 €

Pour information, le texte de loi précise que les pénalités sont prélevées par l'Etat sur le produit des taxes et cotisations foncières du territoire de la commune, après déduction des dépenses effectivement exposées par la commune pour la réalisation de logements sociaux. La pénalitée est reversée suivant sa localisation : aux établissements publics de coopération intercommunales, établissements publics foncier (en Alsace, l'EPFA), fonds nationaux et régionaux d'aménagement foncier et urbain. L'argent est ensuite utilisée pour financer des acquisitions foncières et immobilières en vue de la réalisations de logements locatifs sociaux.

Pour supprimer les fichiers de sauvegarde automatique, de prévisualisation, etc. pour Autocad et Rhino.

gci *.bak -Recurse | foreach{rm $_}
gci *.3dmbak -Recurse | foreach{rm $_}
gci plot.log -Recurse | foreach{rm $_}
gci hardcopy.log -Recurse | foreach{rm $_}
gci Drawing1.dwl -Recurse | foreach{rm $_}
gci Drawing1.dwl2 -Recurse | foreach{rm $_}
gci *.cdc -Recurse | foreach{rm $_}
gci *.tmp -Recurse | foreach{rm $_}

A insérer dans un fichier .ps1 dans le dossier à nettoyer, puis clic droit > Exécuter avec Powershell. Le script s'applique au dossier et ses sous-dossiers.

Ce script Python génère un fichier SCR contenant les instructions Autocad pour la publication de tous les DSD du dossier courant. Pour charger ce fichier SCR dans Autocad, il faut utiliser la commande scr, indisponible sur Autocad LT.

#! python3

#creates a SCR file containing the Autocad commands to publish all DSD files in the current folder. The SCR file must be manually loaded in Autocad.

import os

dirpath = "."
files = os.listdir(dirpath)
scrFileName = str(os.getcwd().split(os.sep)[-1])+"_DSD.scr"

out = ""
out += "filedia 0\n"

for file in files:
    fileFull = os.path.abspath(os.path.join(dirpath, file))
    root, extension = os.path.splitext(fileFull)
    if extension == ".dsd":
        out += "-publish "+str(fileFull)+"\n"

out += "filedia 1"

with open(scrFileName, 'w') as outputFile:
    outputFile.write(out)

Pour que les DSD génèrent les fichiers PDF avec le bon nom dans le bon dossier il faut les configurer comme suit depuis la fenêtre de publication d'Autocad (commande publish) :

• Options de publication > Configurer le dossier
• Options de publication > Décocher "Prompt for name" et entrer le nom du fichier PDF désiré
• Bien écraser le DSD après ces règlages.

Même si on a raté le train, il n'est pas déraisonnable de chercher à attraper le suivant. Ainsi, il faut se mettre à Powershell. D'autant que le nouveau Terminal de Windows 10 nous le propose par défaut.

Voici une commande pour convertir tous les fichiers d'un dossier contenant des FLAC en MP3, avec ffmpeg.

gci -filter *.flac -file | Foreach-Object {ffmpeg -i $_.FullName -ab 320k -map_metadata 0 -id3v2_version 3 "$($_.BaseName).mp3"}

La commande fonctionne avec tout type de fichier que ffmpeg pourrait réussir à convertir en MP3 et génère simplement une erreur lorsqu'elle recontre un fichier qu'elle ne peut pas traiter.

import rhinoscriptsyntax as rs
import Rhino.DocObjects

def loopBlock(block):
    objectsInBlock = rs.BlockObjects(rs.BlockInstanceName(block))

    for obj in objectsInBlock:
        if rs.ObjectType(obj) != 4096:
            action(obj)
        else:
            loopBlock(obj)

def action(object):
    rhino_object = rs.coercerhinoobject(object)
    rhino_object.Attributes.MaterialSource = Rhino.DocObjects.ObjectMaterialSource.MaterialFromLayer
    rhino_object.CommitChanges()

rs.UnselectAllObjects()

objs = rs.AllObjects(include_lights=True, include_grips=True, include_references=True)
if objs:
    for obj in objs:
        if rs.ObjectType(obj) == 4096: #block
           loopBlock(obj)
        else:
            action(obj)
else:
    print "No objects selectable."

Ce script profite de la routine permettant d'appliquer une action sur tous les objets, même ceux situés à l'intérieur de blocs. Il offre un service bien utile lorsqu'on importe de la géométrie d'un autre modeleur (comme SketchUp par exemple) qui insiste pour configurer tous les matériaux par objet.

On remarquera l'usage de la fonction rs.coercerhinoobject() qui permet de passer du GUID d'un objet à l'objet Rhino lui même. Ce passage est nécessaire pour utiliser les fonctions de l'API RhinoCommon.

Mise à jour : cette fonction est maintenant intégrée à CHL Tools (0.3)

Script Python pour automatiser le téléchargement des images des galeries du site d'architecture Baunetz et Baunetz Architekten. Nécessite Python et les librairies Requests et BeautifulSoup (commandes pip dans le code ci-dessous).

#! python3

import re
import shutil
import requests # "pip install requests" to install package
from bs4 import BeautifulSoup # "pip install beautifulsoup4" to install package
headers = {'User-Agent': 'Mozilla/5.0'}

inputURL = input("Adresse de la page Baunetz à aspirer :")

#Functions
def downloadFile(url):
    url = url.replace("\\","")
    filename = url.split("/")[-1]
    responseImg = requests.get(url, stream=True)
    with open(filename, 'wb') as out_file:
        shutil.copyfileobj(responseImg.raw, out_file)
        del responseImg

def batchDownload(urls):
    print (len(urls),"image(s) trouvée(s).")
    c = 1
    for url in urls:
        downloadFile(url)
        print("Téléchargement", c, "/", (len(urls)))
        c += 1



#Patterns
patternBaunetz = re.compile("^https?://(www.)?baunetz.de/meldungen/([^.]*).html$")
patternBaunetzArchitekten = re.compile("^https?://(www.)?baunetz-architekten.de/([^/]*)/([^/]*)/projekt/([^/]*)$")


# Baunetz
if patternBaunetz.match(inputURL):
    response = requests.get(inputURL, headers=headers)
    soup = BeautifulSoup(response.text,"html.parser")

    scriptTags = soup.find_all("script")

    for scriptTag in scriptTags:
        s = str(scriptTag)
        if s.find("xxlGalerie.xxlimages") >= 0:
           urls = re.findall(r"'url': '(https?://[^']*)'",s)
           batchDownload(urls)

    del response
else:
    print("L'adresse ne correspond pas à une page Baunetz.")


#Baunetz Architekten
if patternBaunetzArchitekten.match(inputURL):
    urls = []
    response = requests.get(inputURL, headers=headers)
    soup = BeautifulSoup(response.text,"html.parser")

    #cover img
    coverImg = soup.select("div.project-detail-image img")
    url = coverImg[0].get("data-src")
    urls.append(url)

    #imgs on landing page
    galleryImgs = soup.select("div.project-detail-gallery__image img")
    for img in galleryImgs:
        url = img.get("data-src")
        urls.append(url)

    #more imgs link (completes the slideshow)
    moreImgs = soup.select("div[data-additional-images]")
    s = str(moreImgs[0].get("data-additional-images"))
    moreImgsUrls = re.findall(r'"src":"(https:[^"]*)',s)
    for moreImgsUrl in moreImgsUrls:
        urls.append(moreImgsUrl)

    batchDownload(urls)
    del response
else:
    print("L'adresse ne correspond pas à une page Baunetz Architekten.")


# End
input("Press Enter to close")