J'ai eu l'idée de ce blog en finissant mon énième CDD et en rangeant mon bureau plein de post-its, de notes dans les coins de cahier, sur les feuilles de brouillon ...


Dommage de perdre toutes ces informations glanées au fur et à mesure des années et puis finalement pourquoi ne pas en faire profiter la sphère des heureux sigistes et fières de l'être?


J'ai travaillé sous GRASS, Arcview, ArcGis, MapInfo,QGIS et je suis devenue accroc à l'open source.


Alors voila, il était une fois ...

jeudi 20 août 2015

Le Méridien de Greenwich ne se trouve pas là où nous le croyons... il est 102 mètres trop à l'ouest

Source: http://www.atlantico.fr/decryptage/meridien-greenwich-ne-se-trouve-pas-ou-croyons-est-102-metres-trop-ouest-2289658.html

Le Méridien de Greenwich ne se trouve pas là où nous le croyons... il est 102 mètres trop à l'ouest

Le décalage entre le point de marquage du méridien et l'endroit où il se trouve vraiment a été découvert dans les années 1960.
En cherchant le célèbre méridien de Greenwich (ou premier méridien) avec un GPS, il est possible d'observer que cette ligne nord-sud se trouvant à une longitude zéro est marquée de manière inexacte. Le GPS ne conduit pas l'utilisateur sur le monument qui est censé marquer le premier méridien, mais indique un autre point qui se situe 102 mètres plus loin, à l'est. 
Le décalage entre le point de marquage du méridien et l'endroit où il se trouve vraiment a été découvert dans les années 1960 lorsque les satellites ont commencés à être utilisés pour la localisation précise. Il y a plusieurs décennies, la position du méridien a été déterminée avec du mercure. La surface de ce liquide était connue pour être parallèle à la Terre, avec un niveau de précision alors considéré comme très satisfaisant. Mais les GPS dont nous nous servons massivement aujourd'hui s'avèrent bien plus précis que ce procédé, et cette faille conséquente n'a pu être expliquée que très récemment. 
Des scientifiques de l'Université de Virginie ont découvert que l'aplatissement de la Terre au niveau des pôles est la cause de ce décalage de coordonnées.
Le GPS prend en considération la légère bosse présente à la surface de la Terre, ainsi que sa gravité irrégulière, et les lignes de longitude ne s'étendent pas sur notre planète aussi harmonieusement que ne l'indiquent les données exactes de l'astronomie. Les résultats de ces travaux ont été publiés dans la revue Geodesy. 
Sur ce cliché satellite, se trouve le monument censé marquer le méridien, et sa localisation réelle, 102 mètres plus à l'est.
image: http://www.atlantico.fr/sites/atlantico.fr/files/u61243/2015/08/greenwich_2.jpg


"Je pense qu'il serait opportun de placer un marqueur sur la position réelle du méridien, cela pourrait apporter une explication précise aux personnes qui utilisent un GPS et ne comprennent pas cette anomalie. Cela indiquerait où se trouve le 'GPS zéro', une sorte de nouvelle attraction pour les touristes" commente Ken Seidelmann, chercheur à l'Université de Virginie, pour le site IFLScience. 
Selon le docteur Marek Kukula, les personnels présents sur le site demande régulièrement des explications sur cette anomalie. Mais la situation du méridien de Greenwich n'est pas une exception. Les scientifiques ont découvert que les marqueurs de nombreuses autres lignes de longitude à travers le monde ne sont pas exacts. Les travaux des chercheurs révèlent également que ces décalages ne sont pas identiques selon les méridiens, et varient d'un endroit à l'autre du globe, ce qui confirme les théories des scientifiques sur la gravité. Le plus gros écart relevé a été observé en Amérique du Sud, au Chili, où les équipements astronomique se situent à haute altitude dans les montagnes. Si cette récente étude n' pas provoqué de chamboulement dans le monde scientifique, ces nouvelles découvertes constituent un nouveau défi aux allures de puzzle, que les puristes et les spécialistes se réjouiront de reconstituer. 

Read more at http://www.atlantico.fr/decryptage/meridien-greenwich-ne-se-trouve-pas-ou-croyons-est-102-metres-trop-ouest-2289658.html#UtbvxKIYa3ktsX7u.99

lundi 10 août 2015

Se libérer du format propriétaire ECW et passer au Geotiff en gardant un super rapport qualité/poids

J'ai eu un soucis lorsque j'ai installé un serveur QGIS, associé au web client Lizmap : toutes mes orthophotos livrées par l'IGN étaient au format propriétaire ECW. Or mon serveur ne supporte pas ce format.

J'ai  voulu convertir mes ECWen un format libre tout en gardant l'avantage du poids de l'ECW.

Après avoir posté mon problème sur le forumSIG (http://www.forumsig.org/showthread.php/40755-Support-de-l-ECW),  en suivant les explications de Mdouchin, et en feintant un peu, je pense avoir trovué THE solution: passer en Geotiff !

Tout d'abord cela dépend :
- du nombre de de dalles à convertir
- de la taille des dalles à convertir
- de la présence ou non d'une bande de transparence dans vos orthos

Dans mon cas j'avais plus de 1000 dalles en ECW pour couvrir mon territoire qui étaient assemblées en une seule ortho.

1. ECW sans bande de transparence                                 

  • Conversion de l'ECW en Geotiff
gdal_translate -co COMPRESS=JPEG -co JPEG_QUALITY=90 -co PHOTOMETRIC=YCBCR -co TILED=YES -co BIGTIFF=YES -of GTiff C:/repertoire1/orthophoto.ecw C:/repertoire2/orthophoto.tiff
  • Construction des pyramides pour un affichage optimisé:
gdaladdo -r AVERAGE -ro --config JPEG_QUALITY_OVERVIEW 90 --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config BIGTIFF_OVERVIEW YES -clean C:/repertoire2/orthophoto.tiff 2 4 8 16 32 64

Ajouter --config GDAL_CACHEMAX 1000 --config GDAL_MAX_DATASET_POOL_SIZE 1014 lors de la conversion et --config GDAL_CACHEMAX 1000 lors de la construction des pyramides.

Les options -co PHOTOMETRIC=YCBCR et --config PHOTOMETRIC_OVERVIEW YCBCR ne fonctionnent que sur des dalles sans bande transparente. 

Le résultat est surprenant, on peut descendre à une qualité JPEG de 40 en gardant un résultat tout à fait acceptable!

Exemple: Geotiff + pyramides = 10 Go (Dont 3.21 Go de pyramides) au lieu de 17 pour l'ecw.

    2. ECW avec bande de transparence


Ça se corse si on a une bande transparente que l'on voudrait bien garder. J'ai voulu convertir ma dalle unique mais le résultat, surtout des pyramides, était énorme en poids, donc pas acceptable dans mon cas.

Exemple: conversion d'un ECW de 10 Go en Geotiff, je me retrouve avec un fichier de 32 Go + 18 Go de pyramides.

J'ai tout testé et au final je me suis retrouvée avec des poids variables mais toujours supérieurs à l'ECW : de 8 à 40 Go pour un ECW de 2 Go au départ.

Je dispose donc de plusieurs dalles ortho, assemblées en un seul ECW de 2 Go.

J'ai tenté de convertir la dalle en entier en Geotiff puis de faire des pyramides: poids énorme!

J'ai tenté de convertir la dalle unique en VRT puis de faire les pyramides : énorme!

J'ai tenté de convertir chaque dalle en tiff puis de faire un VRT puis les pyramides: énorme!
et au final j'ai tenté de ...

  • Convertir chaque dalle ECW en Geotiff compressé à 75 (avec un bash pour la conversion en lot, cf. article sur les bash) :
gdal_translate -co COMPRESS=JPEG -co JPEG_QUALITY=75 -co  -co TILED=YES -co BIGTIFF=YES -of GTiff C:/repertoire/orthophoto.ecw C:/repertoire2/orthophoto.tiff
  • Créer pour chaque Geotiff les pyramides avec aussi un bash) :
gdaladdo -r AVERAGE -ro --config JPEG_QUALITY_OVERVIEW 75 --config COMPRESS_OVERVIEW JPEG   --config BIGTIFF_OVERVIEW YES -clean C:/repertoire2/orthophoto.tiff 2 4 8 16 32 64
  • assembler en VRT dans QGIS ou via GDAL (buildvrt)

J'obtiens en tout : pyramides + Geotiff + VRT 2.38 Go 

PS: configuration Windows 7 (je sais, je devrais être sous LINUX, c'est prévu :) )













mercredi 5 août 2015

Paramétrer les couleurs d'un symbole svg créé sous Inkscape:

C'est sympa de créer ce fameux rond dégradé sous Inkscape, ce dégradé qui crée l'illusion de la 3D et de l'utiliser comme symbole dans Qgis, mais recréer et ré-importer ce symbole chaque fois que l'on veut changer sa couleur, ça devient vite lourd, heureusement il y a ce tuto sur le super site du "portail SIG" :

http://www.portailsig.org/content/ajouter-et-personnaliser-des-symboles-svg-depuis-qgis

Il suffit donc d'éditer, via notepad ++, par exemple, le fichier SVG et de remplacer certains paramètres:


  • fill="param(fill)#010101 que l'on remplace par fill="param(fill)#FFF pour modifier le remplissage
  • stroke="param(outline)"#000" modifie la couleur de bordure
  • sroke-width="param(outline-width) 1" modifie l'épaisseur du contour
  • etc,

Pour aller plus loin aller visiter la page citée ci-dessus :)

Quelle est cette projection

Parceque cela nous arrive à tous .... On nous file, on récupère des données sans système de coordonnées et maintenant ô miracle, on va voir sur ce site génial:

http://dogeo.fr/_apps/projection/