Wednesday, October 8, 2008

FWTools

this here be the post dedicated to FWTools (working)

here be the syntax for helpful batch files, etc. for FWTools commands:



1. convert geodatabase.mdb to shapefiles

ogr2ogr out_shp c:\warehouses\testgdb.mdb

#out_shp represents folder created containing converted shapefiles, rename accordingly.



2. batch reproject folder of shapefiles

mkdir epsg102645

set prj_params=-s_srs epsg:26910 -t_srs epsg:102645

for %%g in (*.shp) do (
ogr2ogr -f "esri shapefile" %prj_params% epsg102645\%%g %%g
)

#epsg102645 represents folder created containing reprojected shapefiles, rename accordingly | reassign epsg codes accordingly.



3. batch merge folder of shapefiles to one shapefile

mkdir merged
for %%f in (*.shp) do (
if not exist merged\fed_wd_merge.shp (
ogr2ogr -f "esri shapefile" merged\fed_wd_merge.shp %%f) else (
ogr2ogr -f "esri shapefile" -update -append merged\fed_wd_merge.shp %%f -nln fed_wd_merge )
)



4. The gdalwarp command can be used to reproject images. Here we reproject
the WGS84 geographic image to the Mercator projection:

gdalwarp -t_srs '+proj=merc +datum=WGS84' geoworld.tif mercator.tif



5. batch reproject raster w/ cubic resampling to 'output' folder

mkdir output

set prj_params=-s_srs epsg:26910 -t_srs epsg:102645

for %%g in (*.tif) do (
gdalwarp %prj_params% -r cubic %%g output\%%g
)

No comments: