Skip to content

Instantly share code, notes, and snippets.

@dayne

dayne/watana-extract-notes

Last active August 29, 2015 14:23
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save dayne/c6387171220b28f7f728 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save dayne/c6387171220b28f7f728 to your computer and use it in GitHub Desktop.
Download ZIP
watana extract
Raw
watana-extract-notes
# Found some old notes (2012) of the steps I used to create some contour lines
# derived from the SDMI IFSAR DEM
####
# clipping out
gdal_translate -co TILED=YES -co COMPRESS=LZW -a_nodata -99999 -projwin 252983.504358 1469878.71763 340772.495973 1421881.9325 -of GTiff ../SDMI.IFSAR.2010.FUGRO/dirty-org/DSM.vrt SDMI-IFSAR-2010-DSM-watana_extract.tif
add_overviews.rb SDMI-IFSAR-2010-DSM-watana_extract.tif
gdal_translate -co TILED=YES -co COMPRESS=LZW -a_nodata -99999 -projwin 252983.504358 1469878.71763 340772.495973 1421881.9325 -of GTiff ../SDMI.IFSAR.2010.FUGRO/dirty-org/DTM.vrt SDMI-IFSAR-2010-DTM-watana_extract.tif
add_overviews.rb SDMI-IFSAR-2010-DTM-watana_extract.tif
#gdaladdo -r average SDMI-IFSAR-2010-DTM-watana_extract.tif 2 4 8 16 32 64 128 256 512 1024 2048 4096 8192
gdaldem slope SDMI-IFSAR-2010-DSM-watana_extract.tif SDMI-IFSAR-2010-DSM-slope-watana_extract.tif -s 1.0 -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-DSM-slope-watana_extract.tif
gdaldem aspect SDMI-IFSAR-2010-DSM-watana_extract.tif SDMI-IFSAR-2010-DSM-aspect-watana_extract.tif -zero_for_flat -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-DSM-aspect-watana_extract.tif
gdaldem color-relief SDMI-IFSAR-2010-DSM-watana_extract.tif /home/dbroders/apps/share/qgis/python/plugins/GdalTools/tools/terrain.txt SDMI-IFSAR-2010-DSM-color_relief-watana_extract.tif -co TILED=YES -co COMPRESS=LZW -of GTiff
add_overviews.rb SDMI-IFSAR-2010-DSM-color_relief-watana_extract.tif
gdaldem hillshade SDMI-IFSAR-2010-DSM-watana_extract.tif SDMI-IFSAR-2010-DSM-hillshade-watana_extract.tif -z 1.0 -s 1.0 -az 315.0 -alt 45.0 -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-DSM-hillshade-watana_extract.tif
#
gdal_contour -a ELEV -i 22.0 /home/dbroders/scratch/watana-extract/SDMI-IFSAR-2010-DTM-watana_extract.tif /home/dbroders/scratch/watana-extract/watana_extract-22ft-contours/watana_extract-contours-22ft.shp
Raw
watana-extract-shades
for I in DSM DTM; do
echo "$I slope"
gdaldem slope SDMI-IFSAR-2010-$I-watana_extract.tif SDMI-IFSAR-2010-$I-slope-watana_extract.tif -s 1.0 -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-$I-slope-watana_extract.tif
echo "$I aspect"
gdaldem aspect SDMI-IFSAR-2010-$I-watana_extract.tif SDMI-IFSAR-2010-$I-aspect-watana_extract.tif -zero_for_flat -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-$I-aspect-watana_extract.tif
echo "$I color-relief"
gdaldem color-relief SDMI-IFSAR-2010-$I-watana_extract.tif /home/dbroders/apps/share/qgis/python/plugins/GdalTools/tools/terrain.txt SDMI-IFSAR-2010-$I-color_relief-watana_extract.tif -co TILED=YES -co COMPRESS=LZW -of GTiff
add_overviews.rb SDMI-IFSAR-2010-$I-color_relief-watana_extract.tif
echo "$I hillshade"
gdaldem hillshade SDMI-IFSAR-2010-$I-watana_extract.tif SDMI-IFSAR-2010-$I-hillshade-watana_extract.tif -z 1.0 -s 1.0 -az 315.0 -alt 45.0 -alg ZevenbergenThorne -of GTiff -co TILED=YES -co COMPRESS=LZW
add_overviews.rb SDMI-IFSAR-2010-$I-hillshade-watana_extract.tif
done
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment