Lors de l’\u00e9criture d’un raster, surtout en THR (tr\u00e8s haute r\u00e9solution), il est important de compresser au maximum les images pour avoir des fichiers qui soient g\u00e9rables. Il est important d’avoir en t\u00eate que certaines compression se font avec une perte de qualit\u00e9, et d’autres sans perte. Si l’on veut simplement faire de la visualisation, on peut se permettre un peu de perte, par contre pour du traitement ou pour des donn\u00e9es type MNT, on ne voudra aucune perte d’information.<\/p>\n
Toutes ces histoires de compression sont bien expliqu\u00e9es sur ce blog http:\/\/blog.cleverelephant.ca\/2015\/02\/geotiff-compression-for-dummies.html<\/a><\/p>\n Possibilit\u00e9 de tuiler (daller) pour permettre de g\u00e9rer des trop grosses images. C’est parfois indispensable, par exemple si l’on passe au dela de 20 \u00e0 30000 pixels de large pour une image.<\/p>\n On pourra consulter pour plus d’infos cette page https:\/\/sigeo.cerege.fr\/?p=116<\/a><\/p>\n le GeoTIFF permet la gestion de tr\u00e8s grosses images avec un syst\u00e8me de tuilage interne. Tout est donc stock\u00e9 dans 1 fichier TIFF.<\/p>\n Pour les images >= 4Go, il faudra les \u00e9crire en BigTiff.<\/p>\n Parfois cela est r\u00e9alis\u00e9 automatiquement par l’application qui va \u00e9crire l’image. Si l’on utilise GDAL, il est pr\u00e9f\u00e9rable de le signifier explicitement avec l’option “-co BIGTIFF=YES”<\/p>\n pour forcer l’\u00e9criture des tuiles en interne, c’est avec l’option -co TILED=YES<\/p>\n exemple de commande avec ces 2 options: par exemple pour un MNT ou de l’imagerie que l’on veut garder intacte, j’utilisais la compression avec DEFLATE<\/p>\n ex . Maintenant, il existe de nouveaux algos de compression plus efficaces, ZSTED et ESRI LERC. Exemple de compression avec ZSTED:<\/p>\n on peut utiliser la compression JPEG pour r\u00e9duire efficacement la taille des images multi-canaux en 8 bits (ne fonctionne pas pour les images 16 ou 32 bits) :<\/p>\n Le stockage en JPEG\u00a0est encore plus performant dans l’espace de couleurs YCBCR.<\/p>\n Pour des images en 16 ou 32 bits, on utilisera LZW ou PACKBITS pour la compression avec perte de qualit\u00e9.<\/p>\n Les pyramides permettent d\u2019acc\u00e9l\u00e9rer l\u2019affichage en cr\u00e9ant des aper\u00e7us selon les niveaux de zoomes choisis. Ici la perte de qualit\u00e9 n’est pas un probl\u00e8me car c’est simplement pour l’affichage, et la donn\u00e9e originale n’est pas alt\u00e9r\u00e9e.<\/p>\n Les pyramides peuvent \u00eatre stock\u00e9es dans un fichier s\u00e9par\u00e9 (format OVR), ce qui est en g\u00e9n\u00e9ral pr\u00e9f\u00e9rable pour plus de clart\u00e9, ou en interne dans un GeoTiff.<\/p>\n exemple de pyramides pour de l’imagerie 8 bits :<\/p>\n pour un MNT (16 ou 32 bits en g\u00e9n\u00e9ral), on fera plut\u00f4t avec du DEFLATE ou du LZW<\/p>\n Compression pour dummies :\u00a0 http:\/\/blog.cleverelephant.ca\/2015\/02\/geotiff-compression-for-dummies.html<\/a><\/p>\n Format GeoTifff : https:\/\/gdal.gloobe.org\/gdal\/formats\/gtiff.html<\/a><\/p>\n Format Geotiff : https:\/\/www.gdal.org\/frmt_gtiff.html<\/a><\/p>\n gdaladdo pyramides : https:\/\/www.gdal.org\/gdaladdo.html<\/a><\/p>\n","protected":false},"excerpt":{"rendered":" Lors de l’\u00e9criture d’un raster, surtout en THR (tr\u00e8s haute r\u00e9solution), il est important de compresser au maximum les images pour avoir des fichiers qui soient g\u00e9rables. Il est important d’avoir en t\u00eate que certaines compression se font avec une perte de qualit\u00e9, et d’autres sans perte. Si l’on veut simplement faire de la visualisation, … Continue reading “Gros rasters et Compression”<\/span><\/a><\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[12,7],"tags":[],"class_list":["post-213","post","type-post","status-publish","format-standard","hentry","category-gdal-ogr","category-geo-traitements"],"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/213","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=213"}],"version-history":[{"count":6,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/213\/revisions"}],"predecessor-version":[{"id":948,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/213\/revisions\/948"}],"wp:attachment":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=213"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=213"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=213"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}Tuilage<\/h2>\n
Le BigTIFF et les tuiles internes<\/h2>\n
gdal_translate im_in.tif im_out.tif -co BIGTIFF=YES -co TILED=YES<\/code><\/p>\nCompression sans perte de qualit\u00e9<\/h2>\n
gdal_translate dem_in.tif dem_out.tif -co COMPRESS=DEFLATE<\/code><\/p>\ngdal_translate -of GTiff input.tif output_zstd.tif \\\n-co TILED=YES \\\n-co COMPRESS=ZSTD \\\n-co PREDICTOR=3 \\\n-co ZSTD_LEVEL=9 \\\n-co BIGTIFF=IF_SAFER \\\n-co BLOCKXSIZE=512 -co BLOCKYSIZE=512<\/pre>\n
Compression avec perte de qualit\u00e9<\/h2>\n
gdal_translate \\
\n-co COMPRESS=JPEG \\
\n-co PHOTOMETRIC=YCBCR \\
\n-co TILED=YES \\
\ninput.tif output.tif<\/code><\/p>\nCompression des pyramides raster<\/h2>\n
gdaladdo \\
\n--config COMPRESS_OVERVIEW JPEG \\
\n--config PHOTOMETRIC_OVERVIEW YCBCR \\
\n--config INTERLEAVE_OVERVIEW PIXEL \\
\n-r average \\
\n-ro \\
\nim_input.tif \\
\n2 4 8 16 32<\/code><\/p>\ngdaladdo -ro --config COMPRESS_OVERVIEW LZW dem_input.tif 2 4 8 16 32 64<\/code><\/p>\nR\u00e9f\u00e9rences<\/h2>\n