Raster extent to polygon<\/p>\n\n\n\n
remarque g\u00e9n\u00e9rale: les fichiers rasters interm\u00e9diaires produits peuvent \u00eatre \u00e9crits en format VRT (raster virtuel). Cela permet d’avoir de tous petits fichiers.<\/p>\n\n\n\n
exemple avec un fichier in.tif en entr\u00e9e, et \u00e9criture d’un fichier out.vrt en sortie<\/p>\n\n\n\n
l’option -co NBITS=1 permet de passer l’image en binaire<\/p>\n\n\n\n pour \u00e9crire un fichier en sortie au format VRT, utiliser l’option -of VRT<\/p>\n\n\n\n Le canal alpha est cr\u00e9\u00e9 avec toutes les donn\u00e9es en nodata gr\u00e2ce \u00e0 l’option -dstalpha<\/p>\n\n\n\n on peut aussi prendre en compte une valeur repr\u00e9sentant le nodata dans le fichier en entr\u00e9e gr\u00e2ce \u00e0 l’option -srcnodata<\/b> value<\/em><\/p>\n\n\n\n Bien v\u00e9rifier le raster out1.vrt en sortie pour voir dans quel canal se trouve la bande alpha \u00e0 utiliser par la suite.<\/p>\n\n\n\n l’option -b 1 permet de prendre ce qui est dans la bande 1<\/p>\n\n\n\n on calcule un raster avec une expression logique l\u00e0 o\u00f9 il y a des valeurs<\/p>\n\n\n\n on convertit en polygone le raster:<\/p>\n\n\n\n Raster extent to polygon Conversion en vecteur de l’emprise valide (sans nodata) d’un raster remarque g\u00e9n\u00e9rale: les fichiers rasters interm\u00e9diaires produits peuvent \u00eatre \u00e9crits en format VRT (raster virtuel). Cela permet d’avoir de tous petits fichiers. Etape 1 : Conversion en 1 bit exemple avec un fichier in.tif en entr\u00e9e, et \u00e9criture d’un fichier out.vrt … Continue reading “Emprise d’un raster”<\/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":[37,38,39],"class_list":["post-191","post","type-post","status-publish","format-standard","hentry","category-gdal-ogr","category-sig","tag-extent","tag-limit","tag-raster"],"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/191","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=191"}],"version-history":[{"count":6,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/191\/revisions"}],"predecessor-version":[{"id":930,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/191\/revisions\/930"}],"wp:attachment":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=191"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=191"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=191"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}gdal_translate in.tif out.vrt -ot byte -of VRT -co NBITS=1<\/code><\/p>\n\n\n\nEtape 2 : Cr\u00e9ation d’un canal alpha<\/h3>\n\n\n\n
gdalwarp -dstalpha out.vrt out1.vrt -co NBITS=1<\/code><\/p>\n\n\n\nEtape 3: Vectorisation du raster<\/h3>\n\n\n\n
gdal_polygonize.py out1.vrt -b 1 -f \"ESRI Shapefile\" out2.shp<\/code><\/p>\n\n\n\nMethode Alternative : avec gdal_calc<\/h3>\n\n\n\n
gdal_calc.py -A input --outfile=tmp.tif --type=Byte --calc=\"isnan(A)\" <\/code><\/pre>\n\n\n\ngdal_polygonize.py tmp.tif -f \"ESRI Shapefile\" output<\/code><\/pre>\n","protected":false},"excerpt":{"rendered":"