Ce document d\u00e9crit diff\u00e9rents scripts bas\u00e9s sur GDAL en Shell pour des dalles raster.
\nL\u2019exemple d\u2019utilisation est la r\u00e9cup\u00e9ration de donn\u00e9es aupr\u00e8s de l\u2019IGN, au format TIFF ou ECW<\/p>\n
<\/p>\n
gdalinfo src_dataset<\/code><\/div>\nAffiche l\u2019information sur un jeu de donn\u00e9es raster. Permet notamment d\u2019afficher la r\u00e9f\u00e9rence spatiale et l\u2019\u00e9tendue des donn\u00e9es, le nombre de pixel, la taille du pixel. Documentation officielle\u00a0gdalinfo<\/a><\/p>\nAffecter une r\u00e9f\u00e9rence spatiale \u00e0 un ensemble de dalles<\/h3>\ngdal_translate -a_srs srs_def src_dataset dst_dataset<\/code><\/div>\nAffecte le syst\u00e8me de coordonn\u00e9es \u00ab\u00a0srs_def\u00a0\u00bb \u00e0 l\u2019image en sortie. On utilisera la syntaxe EPSG:n (ex. EPSG\u00a0:2154 pour du Lambert93 RGF93, 27572 pour du NTF Lambert II Etendu). ex. pour du NTF- Lambert II Etendu\u00a0:<\/p>\n
gdal_translate -a_srs EPSG:27572 im_src.tif im_dst.tif<\/code><\/div>\nScript pour un lot de dalles\u00a0:<\/p>\n
#!\/bin\/bash
\nfor FILE in *.tif
\ndo
\nBASE=`basename $FILE .tif`
\nNEWFILE=${BASE}_c.tif
\ngdal_translate -of GTiff -a_srs EPSG:2154 $FILE $NEWFILE
\ndone<\/code><\/div>\noptions de gdal_translate<\/p>\n
\n- compresser une image\u00a0: l\u2019option -co COMPRESS=DEFLATE permet de compresser l\u2019image TIFF sans perte de qualit\u00e9. L\u2019option -co TILED=YES permet de sp\u00e9cifier que l\u2019on fait un tuilage.<\/li>\n
- codage de l\u2019image\u00a0: l\u2019option -co NBITS=1 permet de coder l\u2019image en sortie \u00e0 1 bit. Documentation officielle\u00a0gdal_translate<\/a>
\nCr\u00e9ation d\u2019une mosaique au format \u00ab\u00a0Virtual Raster\u00a0\u00bb VRT<\/h3>\ngdalbuildvrt mosaic.vrt *.tif<\/code><\/div>\n<\/li>\n<\/ul>\nCr\u00e9e un mosa\u00efque \u00ab\u00a0mosaic.vrt\u00a0\u00bb \u00e0 partir d\u2019un ensemble d\u2019une collection d\u2019images tif. Ce format vrt est tr\u00e8s l\u00e9ger et semblable aux algorithmes d\u2019Er-Mapper, \u00e0 savoir que les donn\u00e9es d\u2019origine sont conserv\u00e9es et inchang\u00e9es. Le vrt est simplement un catalogue d\u2019images. Toutes les images de la collection doivent avoir le m\u00eame syst\u00e8me de coordonn\u00e9es. Pour \u00e9viter les zones noires l\u00e0 o\u00f9 il n\u2019y a pas de donn\u00e9es on rajoutera une option \u00ab\u00a0-addalpha\u00a0\u00bb. De cette mani\u00e8re, avec un \u00e9diteur d\u2019image g\u00e9rant le canal alpha (ex. ArcMap), on aura de la transparence dans les zones sans donn\u00e9es source, et de l\u2019opacit\u00e9 dans les zones avec donn\u00e9es source. Ex.\u00a0:<\/p>\n
gdalbuildvrt -addalpha mosaic.vrt *.tif<\/code><\/div>\nDocumentation officielle\u00a0gdalbuildvrt<\/a><\/p>\nConstruire des pyramides<\/h3>\n
Cette fonction permet d\u2019acc\u00e9l\u00e9rer l\u2019affichage en cr\u00e9ant des aper\u00e7us selon les niveaux de zoomes choisis. C\u2019est l\u2019\u00e9quivalent de la cr\u00e9ation de pyramides dans ArcGIS.<\/p>\n
gdaladdo -ro data_src 2 4 8 16 32<\/code><\/div>\nL\u2019option \u00ab\u00a0-ro\u00a0\u00bb permet de cr\u00e9er un aper\u00e7u en tant qu\u2019image externe et non inclus dans l\u2019image tiff elle-m\u00eame. Attention, si vous n\u2019utilisez pas l\u2019option -ro, les pyramides seront construites dans l\u2019image d\u2019origine et elle pourra \u00eatre corrompue.
\nPour compresser le plus possible en JPEG, on utilisera ce type d\u2019options\u00a0:<\/p>\n
gdaladdo -ro --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config JPEG_QUALITY_OVERVIEW 80 data_src 2 4 8 16<\/code><\/div>\nDocumentation\u00a0gdaladdo<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"Ce document d\u00e9crit diff\u00e9rents scripts bas\u00e9s sur GDAL en Shell pour des dalles raster. L\u2019exemple d\u2019utilisation est la r\u00e9cup\u00e9ration de donn\u00e9es aupr\u00e8s de l\u2019IGN, au format TIFF ou ECW Afficher l\u2019information sur une image raster gdalinfo src_dataset Affiche l\u2019information sur un jeu de donn\u00e9es raster. Permet notamment d\u2019afficher la r\u00e9f\u00e9rence spatiale et l\u2019\u00e9tendue des … Continue reading “Gestion des dalles raster avec GDAL”<\/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,10],"tags":[],"class_list":["post-114","post","type-post","status-publish","format-standard","hentry","category-gdal-ogr","category-sig"],"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/114","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=114"}],"version-history":[{"count":1,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/114\/revisions"}],"predecessor-version":[{"id":115,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/114\/revisions\/115"}],"wp:attachment":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=114"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=114"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=114"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}