Blog de la fab des cartes de pollution lumineuse 2016

[frédogoto]

C'est parti.

Ça fait déjà trois jour que je travaille dessus, sans beaucoup avancer.

1 récupérer les données Copernicus

Vous pouvez le voir, c'est pas hyper léger.

Mais alors à maniper c'est carrément l'horreur !

Exemple ci dessous :

Un simple zoom, à la vitesse réelle, alors que j'ai un gros pc avec 24 go de ram.

[bbvideo]

[/bbvideo]

J'ai des soucis de compatibilité des données Corine avec mes logiciels. Je suis en train de devenir chèvre :

pour le moment un logiciel accepte les données : Qgis

Mais c'est pas une fusée ce logiciel, je ne peux rien faire d'autre que les ouvrir...

Je ne peux pas les exporter...

J'ai peut-être une astuce : je vais imprimer les données vectorielles dans un PDF puis importer ce PDF dans illustrator pour maniper les données puis les exporter vers photoshop...

 

Nouvelle difficulté : pour "imprimer" il faut créer un composeur d'impression sur lequel on réimporte les données affichées dans le logiciel. Grrr

[Dieter333]

courage.. on compte sur toi ;-)

[frédogoto]

Merci.

Victoire : l'export' en pdf s'est bien passée.

Reste à transformer tout ça en raster (image matricielle) dans photoshop.

Là, les vrais ennuis commencent (image monumentale) :

[frédogoto]

C'est parti pour la pixellisation.

Image finale 32k*32k pixel soit 1gigapixel.

Et ça, c'est rien car ensuite il faudra ajouter les 29 calques nécessaires, et en 16bit...

htttp://www.avex-asso.org/filerun/wl/?id=dIybPekZgAktK3wcKBwFjSsDXRZGT4Hm

[frédogoto]

Fait !

Mais la résolution finale est trop faible....

htttp://www.avex-asso.org/filerun/wl/?id=nQP6B5DuIT75JhrFZ1wfTClwttzIpVE3

Je vais essayer de la doubler.

[Xavier]

HAL-LU-CI-NANT !!!

Il te faudrait des VM chez un hébergeur. Ils ne font pas ça chez Infomaniak ?

[frédogoto]

ha la vache : je suis limité 32000 pixel ²

[frédogoto]

Compris : dans Photoshop, une image avec calque ne peut pas dépasser cette limite.

Mince... Je vais faire sans.

[frédogoto]

Voila une mini vidéo qui montre la vitesse de rafraîchissement de la carte.

Quand je vous dit que c'est pénible.

Un simple zoom, à la vitesse réelle. Alors que j'ai un gros pc avec 24go de ram.

[frédogoto]

Pour contourner la limite de 32K pixel au carré je vais donc réduire l'emprise de la carte...

Tant pis pour la Turquie, Chypre, l’Islande, et la partie supérieure des pays scandinaves et la partie orientale de la Roumanie et un bout de la Crète...

Accessoirement, pour me faciliter la vie, j'ai changé la couleur des types de sols directement en dégradé de gris selon le type.

[frédogoto]

Et voici la nouvelle emprise que je propose (le carré). Elle n'est pas encore validée, si vous avez des objections c'est le moment...

Cette emprise me permet de doubler pratiquement la résolution des pays bénéficiaire...

[frédogoto]

Voici la résolution finale : c'est mieux.

-asso.org/filerun/wl/?id=HAr5aXkwVls6zaliEDtXStcR8LsPGgxq' alt='?id=HAr5aXkwVls6zaliEDtXStcR8LsPGgxq'>

A comparer.

-asso.org/filerun/wl/?id=nQP6B5DuIT75JhrFZ1wfTClwttzIpVE3' alt='?id=nQP6B5DuIT75JhrFZ1wfTClwttzIpVE3'>

[frédogoto]

la base étant acquise, il faut maintenant le fond de carte pour le relief (c'est plus jolie et plus facile à utiliser) et aussi qui sert à pondérer la pollution lumineuse en fonction de l'altitude..

ces donnée c'est ce qu'on appelle des DEM (digital elevation model)

les données sont la..

http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Eurasia/

en vrac . un poil chiant de faire le tri dans ce merdier mais on y arrive

après assemblage ça donne ça

les niveaux bas sont sombre, les "hauts" sont blancs

-asso.org/filerun/wl/?id=mIQidcHIdxedjzAfRPkGtU6Redrxu6Oe' alt='?id=mIQidcHIdxedjzAfRPkGtU6Redrxu6Oe'>

on peut remarquer immédiatement que l'image est anamorphosée par rapport aux images de corrine (plus étendue en hauteur)

ci dessous une superposition

-asso.org/filerun/wl/?id=wF4QXnJh8I7F50G9lDqiQk0X7kxfxXZb' alt='?id=wF4QXnJh8I7F50G9lDqiQk0X7kxfxXZb'>

c'est un peu complexe à expliquer en détail, mais pour faire court , il existe des centaines de façons de représenter une sphère sur du "plat". aucune n'est parfaite, et chacune est issue d'un processus de choix

dans le cas des cartes Corinne, centrée sur l'europe Corinne utilise une projection assez sympas qui permet d'avoir une taille entre les pays d’Europe qui soie relativement homogène, mais qui à pour léger désagrément de rendre les états du nord plus grand que ceux du sud

le but est donc de faire correspondre les mode de projection de ces carte l'une a l'autre.

j'utilise donc un logiciel qui permet cette opérations

j'ai lancé l'opération a 20:23 et voici ou nous en sommes à 20:57

c’est chaud

-asso.org/filerun/wl/?id=rOHO8ASouvrbaUEO2VlR6wCtO5HBSgK3' alt='?id=rOHO8ASouvrbaUEO2VlR6wCtO5HBSgK3'>

[frédogoto]

et hop voici le dem dans photoshop

projeté a la mode corrine

-asso.org/filerun/wl/?id=BR8NAOyiYQHng14JgMHa4ud5c1Ah2gYf' alt='?id=BR8NAOyiYQHng14JgMHa4ud5c1Ah2gYf'>

quelque chose d’incroyable s'est produit

:

l'image fait 51336*47034

donc on ne peut pas creer d'image de plus 32k² pixel mais s'il les doc existent : pas de soucis

j'ai ensuite importé le corrine par dessus

-asso.org/filerun/wl/?id=zAwf8jdK0G0NLaKsSPGTXQg8tZoGgjJh' alt='?id=zAwf8jdK0G0NLaKsSPGTXQg8tZoGgjJh'>

[frédogoto]

et ce sera tout jusqu' dimanche, je pars en W& astro avec les copain du club

[frédogoto]

douloureuse reprise des activités

apres plusieurs heure d'effort, la "base de pollution lumineuse" est enfin calée sur les données altimétrique

-asso.org/filerun/wl/?id=giKgI7Sgzy1CFZXVHcdXGHVe6VFQDOxW' alt='?id=giKgI7Sgzy1CFZXVHcdXGHVe6VFQDOxW'>

[frédogoto]

passage en 16bit

-asso.org/filerun/wl/?id=9OlUZAOc8BZCxfuI3xplAXi8Ii3yqpdi' alt='?id=9OlUZAOc8BZCxfuI3xplAXi8Ii3yqpdi'>

[frédogoto]

niveau de la couche alti

-asso.org/filerun/wl/?id=Xzr5rYByR04fp2t1bFkSCRmSBHjRePTP' alt='?id=Xzr5rYByR04fp2t1bFkSCRmSBHjRePTP'>

[frédogoto]

ici j’évide tout ce qui n'est pas PL en vue, plus tard de faire jouer la transparence des textures

-asso.org/filerun/wl/?id=y4aUtDicwXq4JolGelb7PNaIFflv8L8d' alt='?id=y4aUtDicwXq4JolGelb7PNaIFflv8L8d'>

[frédogoto]

transparence définie

-asso.org/filerun/wl/?id=APLG0mdPt016hT1KYg8voFBq3CussFL2' alt='?id=APLG0mdPt016hT1KYg8voFBq3CussFL2'>

[Xavier]

"Pfffffuuiiiii quel boulot comme" dirait Gaston

(un belge ! c'est hélas de circonstance...)

[frédogoto]

duplication des calques de PL en vue de la préparation de l'application de l’algorithme

-asso.org/filerun/wl/?id=7V4Ha3EayWdBnEE5R6flS72hi05hI6zY' alt='?id=7V4Ha3EayWdBnEE5R6flS72hi05hI6zY'>

[frédogoto]

depuis ce matin l'algo tourne

c'est long

ici la passe 3/9

ça devrait etre terminé ce soir sans probleme

le fait d'avoir 20go de ram et deux SSD aide pas mal

-asso.org/filerun/wl/?id=GXkuDPTEXi2Df8bAk2VYdqtvpdPUaDYy' alt='?id=GXkuDPTEXi2Df8bAk2VYdqtvpdPUaDYy'>

[frédogoto]

c'est gros...

a gauche : le poid du fichier d'impression, à droite, la taille du document avec les calques

je crois que mon record est atteints : 133 GO

-asso.org/filerun/wl/?id=S0dOb6gP3meKchYT4nUsUI1bZvN61PxF' alt='?id=S0dOb6gP3meKchYT4nUsUI1bZvN61PxF'>

[frédogoto]

ça commence a s'affiner...

finalement ce sera terminé avant midi

-asso.org/filerun/wl/?id=XRnvyIWgaSVMQyInXp8NcOa44NQjPFab' alt='?id=XRnvyIWgaSVMQyInXp8NcOa44NQjPFab'>

[frédogoto]

ayéeee

-asso.org/filerun/wl/?id=Rwe48lwxIo15hHETaPe6HhPEkkJQTFxd' alt='?id=Rwe48lwxIo15hHETaPe6HhPEkkJQTFxd'>

[frédogoto]

Je suis obligé de supprimer régulierement l'historique de mes actions car mes disque de travail sont saturées

ça crains >:O

[frédogoto]

maintenant, je m'attaque à la pondération de la pollution lumineuse en fonction de l'altitude (pour etre plus exacte, en fonction de la densité et la quantité des aérosol)

le probleme est simple, plus on monte, moins les photons trouvent de matières pour se réfléchir.

j'utilise donc mon fichier altimétrique DEM (sou couche du relief) , don j'inverse les valeurs les blanc deviennet les plaine et les noir les sommets

le calques est passé en incrustations : ceci a pour effet d’amplifier les zones sur leur valeurs, créant le masquage par le relief et surtout permettant un rendu fin de la PL en fonction de l'altitude.

Ne croyez pas cel empirique : j'ai effectué de nombreuse mesure de PL du fait de l'altimétrie et ce travail s'appuie sur ces mesures....

voici une animation qui montre l'effet (attendre une minute que l’animation se charge)

-asso.org/filerun/wl/?id=Lc1wZtHMj9lfNozgFoWlVF7hOd5rSFkP' alt='?id=Lc1wZtHMj9lfNozgFoWlVF7hOd5rSFkP'>

[frédogoto]

pour l'anecdote : j'ai lancé l'enregistrement de sauvegarde du fichier à 22h00.

il est 22h50 et la progression indique 39%

[frédogoto]

23h20 : 46%

je vai me coucher.. ce sera tout pour ce soir

[Fred_76]

Joli travail !

[frédogoto]

plantage de photoshop lors d'une sauvegarde hier soir

la premier pour ce "run" : en cause : plus de place disponible sur les disque de travail..

j'ai fait un ménage colossal se matin : 500go de libéré et je relance les opérations

[frédogoto]

depuis le début d’après midi : travail sur les textures....

c'est trop long

[frédogoto]

maintenant j’attaque la "mise en couleur" je le fait depuis photoshop

mais j'ai vue que cette fonction existe depuis un logiciel de géographie : global mapper 17 64b

je vais d'abord le faire "comme d'habitude" sous photoshop, puis j’essaierais de la faire depuis le logiciel de geo.

l'avantage : dans le logiciel de geo je peu ajuster en quasi temps réel la courbe, en fonction des photos du ciel que les astronomes amateurs mon envoyés

http://www.astrosurf.com/ubb/Forum2/HTML/041618.html

et

http://www.webastro.net/forum/showthread.php?t=131918

ce sera plus pratique et plus juste car sur le logiciel de geo je peux également avoir les coordonnées d'un lieu en temps réel.

dans photoshop non, je dois donc faire des capture d’écran dans google map d'un lieu ou telle photo a été prise, superposer dans Photoshop la capture d’écran à la carte de PL, mesurer la PL de la photo et ajuster la couleur du dégradé de couleur en fonction...

j'ai mieu vous dire que c'est hyper chiant : j'ai 46 photos a traiter de la sorte

[frédogoto]

finalement il n'est pas possible d'appliquer une cartographie des couleur" dans global mapper si les donnée ne sont pas reconue comme etant des donnée altimétrique..

j'ai bien essayé de faire croire au logiciel que les donnée de PL etaient des données d'alti mais rien y fait >

donc retour a la case photoshop avec les interminables manipe qui vont bien pour ajuster et qualibrer

en revanche, un fois fait j'ai pu importer mes cartes dans globalmap et y intégrer un fond de caret altimétrique téléchargé en ligen en temps reel (bien plus precis que mes propres fond actuels

je peux même ajouter des données vectorielle comme les département, des communes. la je cherche a intégrer les nom des villes

sur cette vue vue avez une simu si l'europe passait en un éclairage LED qui remplaceraient à 75% l’éclairage sodium classique