Hello Everyone!<div><br></div><div>I have a 30m DEM raster for which I want to calculate Horizon angles (Max angle of sky obstruction) in some specific directions from a point within the DEM. Lets say the point is in middle of DEM and I want to calculate max sky obstruction angle from this point for 30, 60, 90 and 120 degree North and so on.</div>
<div><br></div><div>All I know is to march out in a specific direction and calculate max angle on the basis of height and distance your from point .</div><div>How can I do that in QGis python console After loading the raster through gdal and reading it as numpy arrays? What could be the logic of this calculation? Urgent Help needed</div>
<div>thanks in advance</div>