User Tools

Site Tools



# preparing raster data for later modelling...
# make environmental date equal resolution an extent based on a shapefile
for file in /home/user/LVM_share/styria/ras/*.asc; do filename=$(basename $file)
gdalwarp -tr 10 10 -r bilinear -cutline ~/LVM_share/styria/vec/stmk.shp -crop_to_cutline -co COMPRESS=DEFLATE -co ZLEVEL=9 -overwrite $file /home/user/LVM_share/styria/output/crop_$filename.tif 
# write a list of files in the output directory
ls ~/LVM_share/styria/output > ~/LVM_share/styria/filenames.txt
# creates a vrt to stack rasters
gdalbuildvrt -separate -input_file_list /home/user/LVM_share/styria/output/filenames.txt /home/user/LVM_share/styria/output/inraster.vrt
# extract raster values to a point shapefile
pkextract -i ~/LVM_share/styria/output/inraster.vrt -s ~/LVM_share/styria/vec/campull.shp -f "ESRI Shapefile" -o ~/LVM_share/styria/vec/extracted.shp
wikistud/mat2017_patrick_shwager.txt · Last modified: 2017/12/05 22:53 (external edit)