diff --git a/.DS_Store b/.DS_Store deleted file mode 100644 index 8c8e9e3..0000000 Binary files a/.DS_Store and /dev/null differ diff --git a/.gitignore b/.gitignore new file mode 100644 index 0000000..496ee2c --- /dev/null +++ b/.gitignore @@ -0,0 +1 @@ +.DS_Store \ No newline at end of file diff --git a/README.md b/README.md index 31e7d36..e5fc853 100644 --- a/README.md +++ b/README.md @@ -7,18 +7,22 @@ As described in Daniel Huffman's [Creating Shaded Relief in Blender](https://som # How to Use 1) Clone/download the repo. -2) Make sure [GDAL/OGR](https://gdal.org/) is installed. (gdalinfo --version) -3) Create a vector specifying your area of interest and add to the project folder. Default points to "box.geojson"- to update, open "blenderize.sh" with your text editor of choice and update the "bounding_box" variable. -4) Find and download DEM(s) which cover at least your area of interest. Err on the side of including more than you think you might need. Add your DEM's to the project folder. -5) Specify your target projection by updating the "target_projection" variable in "blenderize.sh". EPSG codes or proj4 strings are acceptable. -6) Open a command line window, navigate to the project folder, and run: - sh blenderize.sh [space-separated list of DEM's] +1) Make sure [GDAL/OGR](https://gdal.org/) is installed. (gdalinfo --version) +1) Make sure the GDAL/OGR python module is installed. +1) Create a vector specifying your area of interest and add to the project folder. Default points to "box.geojson"- to update, open "blenderize.sh" with your text editor of choice and update the "bounding_box" variable. + - [geojson.io]() is a useful tool for building a simple rectangular bounding box in the correct geojson format. +1) Find and download DEM(s) which cover at least your area of interest. Err on the side of including more than you think you might need. Add your DEM's to the project folder. +1) Specify your target projection by updating the "target_projection" variable in "blenderize.sh". EPSG codes or proj4 strings are acceptable. +1) Open a command line window, navigate to the project folder, and run: +`sh blenderize.sh [space-separated list of DEM's]` For example: - sh blenderize.sh my_dem.tif //just one DEM - sh blenderize.sh *.tif //all DEM's in project folder + `sh blenderize.sh my_dem.tif` for just one DEM + `sh blenderize.sh *.tif` for all DEM's in project folder -#Notes -By default, the script will mosaic(if there are multiple dems), reproject to your target projection, clip to your bounding box, and output a TIF with 16-bit unsigned integer values, scaled from 0 to 65,535. The script *does not* resample the DEM, instead retaining the full resolution of the input data. If I want to quickly resample my elevation data, I just pull the output into photoshop and make sure to resize with resampling set to cubic or bilinear. Here's a couple good resources for more on generalizing elevation data for hillshading: -https://cartographicperspectives.org/index.php/journal/article/view/cp67-leonowicz-et-al/pdf -http://shadedrelief.com/tutorials.html -https://somethingaboutmaps.wordpress.com/2011/10/18/on-generalization-blending-for-shaded-relief/ +# Notes + +By default, the script will mosaic (if there are multiple dems), reproject to your target projection, clip to your bounding box, and output a TIF with 16-bit unsigned integer values, scaled from 0 to 65,535. The script *does not* resample the DEM, instead retaining the full resolution of the input data. If I want to quickly resample my elevation data, I just pull the output into photoshop and make sure to resize with resampling set to cubic or bilinear. Here's a couple good resources for more on generalizing elevation data for hillshading: + +- https://cartographicperspectives.org/index.php/journal/article/view/cp67-leonowicz-et-al/pdf +- http://shadedrelief.com/tutorials.html +- https://somethingaboutmaps.wordpress.com/2011/10/18/on-generalization-blending-for-shaded-relief/