[Qgis-user] Horizon angles

[Mohammad Omayr]

Hello Everyone!

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.

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 .
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
thanks in advance
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-user/attachments/20121119/43adcc45/attachment.html>