{"id":36,"date":"2020-03-30T07:00:00","date_gmt":"2020-03-30T07:00:00","guid":{"rendered":"http:\/\/wcwp.rf.gd\/2020\/03\/30\/https-landscapearchaeology-org-2020-direction-viewshed\/"},"modified":"2020-03-30T07:00:00","modified_gmt":"2020-03-30T07:00:00","slug":"direction-viewshed","status":"publish","type":"post","link":"https:\/\/landscapeanalysis.org\/?p=36","title":{"rendered":"Distance and direction filters for QGIS Viewshed analysis"},"content":{"rendered":"

Viewshed analysis produces a model of visible areas from a chosen point, given a 2.5 or 3D elevation model. Most often, we model full range visibility, in all directions from an observer point, but such models may not satisfy all uses. For instance, we are sometimes interested in a scene that an observer can see in front of him\/her, excluding the scenery behind his\/her back. We may also prefer to exclude certain zones for various reasons, for instance when we model different visibility ranges (close, middle, far).<\/p>\n

A single visibility model can easily be tailored by deleting unwanted zones, but the problem is not trivial when dealing with hundreds of viewsheds. Starting with the 1.1 version of the Viewshed analysis plugin for QGIS<\/a> users will be able to define a doughnut & wedge shaped mask for each visibility model by associating three parameters with each observer point. These are entered in the observer points data table, under columns radius_in<\/em> (inside radius), azim_1<\/em> and azim_2<\/em> (azimuth angles that limit the wedge zone). Note that in the circular domain of angular variables both, 300 to 20 and 20 to 300 ranges are valid; they cover two opposite sectors.<\/p>\n

\"20-03-30-shema.png\"<\/p>\n

Directional viewshed analysis: a tutorial<\/h2>\n

Our example will be a section of the beautiful Premu\u017ei\u0107 trail which is following the main ridge of the Velebit mountain, Croatia. The trail was built in the 1930\u2019s specifically for hiking, and is truly a masterpiece (see on NP Velebit<\/a>).<\/p>\n

\"20-03-30-premuzic_path.jpg\"<\/p>\n

Let\u2019s answer the following question: what scenery does a hiker see while walking along the path? Obviously, he\/she cannot see behind his\/her back, so we need to restrict the viewing range, let\u2019s say to 45 degrees on each side. Since we\u2019re interested in the overall scenery, we will exclude the immediate zone up to one kilometre distance. The overall analysis radius will be 25 kilometres, to test for far reaching views.<\/p>\n

We start with a line drawn along the path. Find the \u201cExtract vertices\u201d algorithm and use it to transform the path into points.<\/p>\n

\"20-03-30-extract_vertices.JPG\"<\/p>\n

Take a look at the data table of the new file. We have line direction (angle), distance and an ID number for each vertex (yay, all in one step!).<\/p>\n

\"20-03-30-vertices_table.JPG\"<\/p>\n

Now, create observer points for the viewshed analysis using desired parameters. Use vertex_index<\/em> as feature ID.<\/p>\n

\"20-03-30-create_viewpoints.JPG\"<\/p>\n

The observer file does not contain the azimuth data: we have to create a data join with the previous file.<\/p>\n

\"20-03-30-join.jpg\"<\/p>\n

We will now calculate our azimuths and store the parameter in table fields named azim_1<\/em> and azim_2<\/em>. My formula for azim_1<\/em> is as following, considering a 90\u00b0 view range:<\/p>\n

\n
\n
CASE<\/span> \n  WHEN<\/span> \"Vertices_angle\"<\/span> ><\/span> 45<\/span>\n     THEN<\/span> \"Vertices_angle\"<\/span> -<\/span> 45<\/span>\n  ELSE<\/span>  \"Vertices_angle\"<\/span> -<\/span> 45<\/span> +<\/span> 360<\/span>\nEND<\/span>\n<\/code><\/pre>\n<\/div>\n<\/div>\n
\"20-03-30-azim_calc.JPG\"<\/p>\n
And for azim_2<\/em> field:<\/p>\n
\n
\n
CASE<\/span> \n  WHEN<\/span> \"Vertices_angle\"<\/span> <<\/span> 315<\/span>\n    THEN<\/span> \"Vertices_angle\"<\/span> +<\/span> 45<\/span>\n  ELSE<\/span>  \"Vertices_angle\"<\/span> +<\/span> 45<\/span> -<\/span> 360<\/span>\nEND<\/span>\n<\/code><\/pre>\n<\/div>\n<\/div>\n
If needed, you can also specify the inner radius in the radius_in<\/em> filed. Here, it will be 1000 meters.<\/p>\n
\"20-03-30-radius_in.jpg\"<\/p>\n
Now, we can finally run the viewshed analysis. All vital parameters have been entered into Viewpoints data table, we just choose file input\/output. (25 kilometres is a fairly large radius, the calculation may take a couple of minutes.)<\/p>\n
\"20-03-30-viewshed.JPG\"<\/p>\n
We have, now, a model of directed cumulative viewshed which reveals the areas most often seen by a hiker taking the Premu\u017ei\u0107 trail.<\/p>\n
\"20-03-30-view_go.jpg\"<\/p>\n
But wait, what happens when he\/she goes back? Well, we simply readjust our azim_1<\/em> and azim_2<\/em> parameters to point in the opposite direction. (I\u2019m leaving you the pleasure of figuring out the calculation. Or, you can use \u201cReverse line direction\u201d in QGIS and repeat all the steps, if you\u2019re lazy\u2026 )<\/p>\n
Here is our return viewshed model. Apparently, it is more oriented towards the sea than the previous one.<\/p>\n
\"20-03-30-view_back.jpg\"<\/p>\n
Happy mapping!<\/p>\n","protected":false},"excerpt":{"rendered":"
Viewshed analysis produces a model of visible areas from a chosen point, given a 2.5 or 3D elevation model. Most often, we model full range visibility, in all directions from an observer point, but such models may not satisfy all uses. For instance, we are sometimes interested in a scene that an observer can see […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[6],"tags":[13],"class_list":["post-36","post","type-post","status-publish","format-standard","hentry","category-qgis-plugins","tag-qgis-visibility-plugin"],"_links":{"self":[{"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=\/wp\/v2\/posts\/36","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=36"}],"version-history":[{"count":0,"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=\/wp\/v2\/posts\/36\/revisions"}],"wp:attachment":[{"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=36"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=36"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/landscapeanalysis.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=36"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}