Blog de la fab des cartes de pollution lumineuse 2016

[frédogoto 385]

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 3]

courage.. on compte sur toi ;-)

[frédogoto 385]

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 385]

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 385]

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 39]

HAL-LU-CI-NANT !!!

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

[frédogoto 385]

ha la vache : je suis limité 32000 pixel ²

[frédogoto 385]

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

Mince... Je vais faire sans.

[frédogoto 385]

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 385]

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 385]

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 385]

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 385]

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 385]

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 385]

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

[frédogoto 385]

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 385]

passage en 16bit

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

[frédogoto 385]

niveau de la couche alti

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

[frédogoto 385]

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 385]

transparence définie

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

[Xavier 39]

"Pfffffuuiiiii quel boulot comme" dirait Gaston

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

[frédogoto 385]

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 385]

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 385]

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 385]

ça commence a s'affiner...

finalement ce sera terminé avant midi

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