Thursday, 5 October 2017

raster - How to iterate GRASS region definition?


I am using GRASS r.skyview command (v.7, Windows) to determine "urban canyon effect" on pollution due to the presence of buildings. The command nicely works, unfortunately I need to process a very large raster (representing buildings in an entire municipality) with a high resolution (10260x6762 pixels) and the command, which process every single pixel of the raster, run very slowly (about 24h to run on the entire area).


Since in the end I need the value of the skyview factor only for some specific points in space, I am wondering if it is possible to automatically do de following:


1) import in GRASS a vector layer representing point locations


2) select a single point


3) set GRASS region to a 100m x 100m square around selected point


4) launch r.skyview command


and loop through all points in vector layer.





No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...