{"id":188,"date":"2014-11-07T10:58:30","date_gmt":"2014-11-07T09:58:30","guid":{"rendered":"http:\/\/elcep.legtux.org\/?p=188"},"modified":"2014-11-10T14:20:08","modified_gmt":"2014-11-10T13:20:08","slug":"ogr2ogr-et-netlogo","status":"publish","type":"post","link":"https:\/\/elcep.legtux.org\/?p=188","title":{"rendered":"ogr2ogr et netlogo"},"content":{"rendered":"

\"GDALLogoColor.svg\"<\/a>Un tr\u00e8s petit post de blog pour garder sous le coude deux ou trois petites choses qui peuvent servir quand on veut pr\u00e9parer rapidement des donn\u00e9es SIG pour netlogo…<\/p>\n

J’ai eu besoin l\u00e0 tout de suite de d\u00e9couper une couche vectorielle, mais je n’avais pas \u00e0 ma disposition de couche pour faire le clip… Donc j’ai cherch\u00e9 un moyen d’utiliser la bbox (non ce n’est pas la box de chez Bouygues Telecom mais boundary box<\/em> : le rectangle qui contient toutes les entit\u00e9s vectorielles)<\/p>\n

Eh bien j’ai trouv\u00e9 sur un post de Tim Sutton (Linfiniti.com)<\/a> des petites astuces pour me simplifier la vie :<\/p>\n

D\u00e9finir la bbox<\/h2>\n

Tim propose un petit script bash qui fait des merveilles et que je me permets de reprendre ici :<\/p>\n

#!\/bin\/bash\r\nSHPFILE=$1\r\nBASE=`basename $SHPFILE .shp`\r\nEXTENT=`ogrinfo -so $SHPFILE $BASE | grep Extent \\\r\n| sed 's\/Extent: \/\/g' | sed 's\/(\/\/g' | sed 's\/)\/\/g' \\\r\n| sed 's\/ - \/, \/g'`\r\nEXTENT=`echo $EXTENT | awk -F ',' '{print $1 \" \" $4 \" \" $3 \" \" $2}'`\r\necho $EXTENT<\/pre>\n
si l’on consid\u00e8re que mon script s’appelle extent.sh on pourra donc obtenir la bbox avec<\/p>\n
.\/extent.sh monShapeFile.shp\r\n3.058669 42.531248 3.065417 42.526966<\/pre>\n
D\u00e9couper le shapeFile avec l’emprise<\/h2>\n
On peut ensuite bien tranquillement profiter de la commande que je proposais dans le pr\u00e9c\u00e9dent post<\/a> sur le sujet, o\u00f9 en lieu et place du shapeFile nous pouvons passer directement l’emprise!<\/p>\n
ogr2ogr -clipsrc 3.058669 42.531248 3.065417 42.526966 -s_srs epsg:2154 -t_srs epsg:4326 ParcellesVitiSQLITE\/small_parcellesWGS84.shp ParcellesVitiSQLITE\/ParcellesViti_L93.shp\r\n<\/pre>\n
D\u00e9couper un raster avec l’emprise<\/h2>\n
Attention le fonctionnement de gdalwrap est un peut diff\u00e9rent de l’option clipsrc … il faudra r\u00e9organiser les coordonn\u00e9es. On effectuera la d\u00e9coupe comme \u00e7a :<\/p>\n
gdalwarp -te 3.058669 42.526966 3.065417 42.531248 temperature_zone.tif temperature_zone_crop.tif<\/pre>\n
et il ne faudra pas oublier de convertir le fichier en asc avec<\/p>\n
gdal_translate -of \"AAIGrid\" -b 1 -co FORCE_CELLSIZE=TRUE temperature_zone_crop.tif temperature_zone_crop.acs\r\n\r\n<\/pre>\n
Voil\u00e0 on peut vite repasser dans netlogo avec le p’tit post pr\u00e9c\u00e9dant<\/a> (ou gama<\/a> :-p )<\/p>\n","protected":false},"excerpt":{"rendered":"
Un tr\u00e8s petit post de blog pour garder sous le coude deux ou trois petites choses qui peuvent servir quand on veut pr\u00e9parer rapidement des donn\u00e9es SIG pour netlogo… J’ai eu besoin l\u00e0 tout de suite de d\u00e9couper une couche vectorielle, mais je n’avais pas \u00e0 ma disposition de couche pour faire le clip… Donc … Continuer la lecture de ogr2ogr et netlogo<\/span> →<\/span><\/a><\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[2,8],"tags":[28,9,29],"class_list":["post-188","post","type-post","status-publish","format-standard","hentry","category-geomatique","category-sma","tag-gis","tag-netlogo","tag-ogr2ogr"],"_links":{"self":[{"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/posts\/188","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=188"}],"version-history":[{"count":10,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/posts\/188\/revisions"}],"predecessor-version":[{"id":207,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=\/wp\/v2\/posts\/188\/revisions\/207"}],"wp:attachment":[{"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=188"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=188"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/elcep.legtux.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=188"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}