QGISで可視領域分析～北海道百年記念塔を例として～（後編）

[前編](https://qiita.com/nuporopet/items/8e23d92120f03a051411)では、QGISとGRASS GISを使って可視領域分析を実行しました。
後編では得られた可視領域について、実際に塔を見ることができるかを検証します。

# 4. 結果

得られた結果を確認してみます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/72fb1886-192c-4399-286c-f567f35ca67b.png)
この図では、QGISのViewshedの結果（薄い白黒の部分）と、GRASS GISのr.viewshedの結果（薄い黄色の部分）を重ね合わせています。
両者で多少の違いはありますが、ほぼ一致しています。簡単に可視領域分析を行う上では、QGISのViewshedを使い、詳細にパラメタを調整しながら可視領域分析を行うには、GRASS GISのr.viewshedを使うという使い分けをすると良いかも知れません。

結果の確認方法ですが、可視領域分析の結果をもとに、実際に塔が見える、見えないを確認してみます。
図の中で、赤い▼は塔の位置、青い▼は検証地点です。なお、検証地点は、可視領域分析で塔が見えないと判定された地点（①）、見えると判定された地点（②）としました。

まず、塔が見えない地点①からの視界を確認してみます。
残念ながら現地での確認が出来なかったため、ここでは模擬的にWeb上のツールを使います。
使ったのは国土交通省の提供する[PLATEAU VIEW](https://plateauview.mlit.go.jp/)です。
PLATEAU VIEWに北海道札幌市の建物モデルを追加し、検証地点から塔の方向の視界を確認します。

塔が見えないと判定された[地点①](https://plateauview.mlit.go.jp/#start=%7B%22version%22%3A%228.0.0%22%2C%22initSources%22%3A%5B%7B%22stratum%22%3A%22user%22%2C%22models%22%3A%7B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%3A%7B%22highlightColor%22%3A%22%23D54A4A%22%2C%22style%22%3A%7B%22color%22%3A%7B%22conditions%22%3A%5B%5B%22true%22%2C%22color%28%27%23FFFFFF%27%2C1+%29%22%5D%5D%7D%7D%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%5D%2C%22type%22%3A%223d-tiles%22%7D%2C%22%2F%22%3A%7B%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%3A%7B%22knownContainerUniqueIds%22%3A%5B%22%2F%22%5D%2C%22type%22%3A%22group%22%7D%7D%2C%22workbench%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%5D%2C%22timeline%22%3A%5B%5D%2C%22initialCamera%22%3A%7B%22west%22%3A141.48970860451175%2C%22south%22%3A43.014474240183546%2C%22east%22%3A141.49330427546542%2C%22north%22%3A43.01600187821261%2C%22position%22%3A%7B%22x%22%3A-3655148.09216577%2C%22y%22%3A2908320.851790169%2C%22z%22%3A4328600.019294306%7D%2C%22direction%22%3A%7B%22x%22%3A0.533803925249846%2C%22y%22%3A-0.4247360291362878%2C%22z%22%3A0.7311994768470471%7D%2C%22up%22%3A%7B%22x%22%3A-0.5721751920764476%2C%22y%22%3A0.45526720122767633%2C%22z%22%3A0.6821637084004094%7D%7D%2C%22homeCamera%22%3A%7B%22west%22%3A139.63595%2C%22south%22%3A35.53022%2C%22east%22%3A139.95451%2C%22north%22%3A35.81489%7D%2C%22baseMaps%22%3A%7B%22defaultBaseMapId%22%3A%22%2Fbasemap%2F%2F%E5%85%A8%E5%9B%BD%E6%9C%80%E6%96%B0%E5%86%99%E7%9C%9F+%28%E3%82%B7%E3%83%BC%E3%83%A0%E3%83%AC%E3%82%B9%29%22%7D%2C%22viewerMode%22%3A%223d%22%2C%22currentTime%22%3A%7B%22dayNumber%22%3A2459947%2C%22secondsOfDay%22%3A63554.017%7D%2C%22showSplitter%22%3Afalse%2C%22splitPosition%22%3A0.5%2C%22stories%22%3A%5B%5D%7D%5D%7D)からの視界を以下に示します。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/b8ff32f4-7931-f213-fe28-774df5b21cf1.png)
確かに、塔のある方向には前方に小高い丘があり、塔の姿は確認できません。

次に塔が見えると判定された[地点②](https://plateauview.mlit.go.jp/#start=%7B%22version%22%3A%228.0.0%22%2C%22initSources%22%3A%5B%7B%22stratum%22%3A%22user%22%2C%22models%22%3A%7B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%3A%7B%22url%22%3A%22https%3A%2F%2Fd2jfi34fqvxlsc.cloudfront.net%2Fmain%2Fdata%2F3d-tiles%2Fbldg%2F01100_sapporo%2Fnotexture%2Ftileset.json%22%2C%22highlightColor%22%3A%22%23D54A4A%22%2C%22style%22%3A%7B%22color%22%3A%7B%22conditions%22%3A%5B%5B%22true%22%2C%22color%28%27%23FFFFFF%27%2C1+%29%22%5D%5D%7D%7D%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%5D%2C%22type%22%3A%223d-tiles%22%7D%2C%22%2F%22%3A%7B%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%3A%7B%22knownContainerUniqueIds%22%3A%5B%22%2F%22%5D%2C%22type%22%3A%22group%22%7D%7D%2C%22workbench%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%5D%2C%22timeline%22%3A%5B%5D%2C%22initialCamera%22%3A%7B%22west%22%3A141.49267935407292%2C%22south%22%3A43.05393443991406%2C%22east%22%3A141.51512935292203%2C%22north%22%3A43.0634650078641%2C%22position%22%3A%7B%22x%22%3A-3652373.6800049474%2C%22y%22%3A2906750.221433988%2C%22z%22%3A4331952.80313889%7D%2C%22direction%22%3A%7B%22x%22%3A-0.42992798657577663%2C%22y%22%3A-0.8750458631139917%2C%22z%22%3A0.22238854243415548%7D%2C%22up%22%3A%7B%22x%22%3A-0.5721779023790939%2C%22y%22%3A0.4546031788237609%2C%22z%22%3A0.682604129662568%7D%7D%2C%22homeCamera%22%3A%7B%22west%22%3A139.63595%2C%22south%22%3A35.53022%2C%22east%22%3A139.95451%2C%22north%22%3A35.81489%7D%2C%22baseMaps%22%3A%7B%22defaultBaseMapId%22%3A%22%2Fbasemap%2F%2F%E5%85%A8%E5%9B%BD%E6%9C%80%E6%96%B0%E5%86%99%E7%9C%9F+%28%E3%82%B7%E3%83%BC%E3%83%A0%E3%83%AC%E3%82%B9%29%22%7D%2C%22viewerMode%22%3A%223d%22%2C%22currentTime%22%3A%7B%22dayNumber%22%3A2459947%2C%22secondsOfDay%22%3A63554.017%7D%2C%22showSplitter%22%3Afalse%2C%22splitPosition%22%3A0.5%2C%22stories%22%3A%5B%5D%7D%5D%7D)からの視界を確認してみます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/6130d21a-881e-8e9d-80d8-870a38f08a13.png)
確かにここでは正面に塔の姿が見えました。

最後に、塔からの可視領域の距離を15kmに拡大して、いくつかの地点からの視界を確認してみます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/e5f7d1a5-6da8-57ea-b425-52307cd1b7d0.png)

観光地として有名な[羊ヶ丘展望台](https://plateauview.mlit.go.jp/#start=%7B%22version%22%3A%228.0.0%22%2C%22initSources%22%3A%5B%7B%22stratum%22%3A%22user%22%2C%22models%22%3A%7B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%3A%7B%22url%22%3A%22https%3A%2F%2Fd2jfi34fqvxlsc.cloudfront.net%2Fmain%2Fdata%2F3d-tiles%2Fbldg%2F01100_sapporo%2Fnotexture%2Ftileset.json%22%2C%22highlightColor%22%3A%22%23D54A4A%22%2C%22style%22%3A%7B%22color%22%3A%7B%22conditions%22%3A%5B%5B%22true%22%2C%22color%28%27%23FFFFFF%27%2C1+%29%22%5D%5D%7D%7D%2C%22filters%22%3A%5B%7B%22name%22%3A%22%E9%AB%98%E3%81%95%E3%81%A7%E7%B5%9E%E3%82%8A%E8%BE%BC%E3%81%BF%22%2C%22minimumShown%22%3A0%2C%22maximumShown%22%3A200%7D%5D%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%5D%2C%22type%22%3A%223d-tiles%22%7D%2C%22%2F%22%3A%7B%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%3A%7B%22knownContainerUniqueIds%22%3A%5B%22%2F%22%5D%2C%22type%22%3A%22group%22%7D%7D%2C%22workbench%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%5D%2C%22timeline%22%3A%5B%5D%2C%22initialCamera%22%3A%7B%22west%22%3A141.41065949708582%2C%22south%22%3A43.03488328664458%2C%22east%22%3A141.5591939214186%2C%22north%22%3A43.097904314486364%2C%22position%22%3A%7B%22x%22%3A-3651110.1101213503%2C%22y%22%3A2915180.142768218%2C%22z%22%3A4327552.898991197%7D%2C%22direction%22%3A%7B%22x%22%3A-0.05032783627928877%2C%22y%22%3A-0.8563408303810556%2C%22z%22%3A0.5139528102051094%7D%2C%22up%22%3A%7B%22x%22%3A-0.572139971253054%2C%22y%22%3A0.44651604752910085%2C%22z%22%3A0.6879529581254407%7D%7D%2C%22homeCamera%22%3A%7B%22west%22%3A139.63595%2C%22south%22%3A35.53022%2C%22east%22%3A139.95451%2C%22north%22%3A35.81489%7D%2C%22baseMaps%22%3A%7B%22defaultBaseMapId%22%3A%22%2Fbasemap%2F%2F%E5%85%A8%E5%9B%BD%E6%9C%80%E6%96%B0%E5%86%99%E7%9C%9F+%28%E3%82%B7%E3%83%BC%E3%83%A0%E3%83%AC%E3%82%B9%29%22%7D%2C%22viewerMode%22%3A%223d%22%2C%22currentTime%22%3A%7B%22dayNumber%22%3A2459947%2C%22secondsOfDay%22%3A63554.017%7D%2C%22showSplitter%22%3Afalse%2C%22splitPosition%22%3A0.5%2C%22stories%22%3A%5B%5D%7D%5D%7D)からの視界では、遠くに塔の姿を確認できます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/bd6a2533-70cf-1aae-0071-9a8c0bb0fda4.png)
塔までの直線距離は10.4kmで、空気の乾いた晴れた日なら肉眼でも見れるかも知れません。

次に、[藻岩山山頂](https://plateauview.mlit.go.jp/#start=%7B%22version%22%3A%228.0.0%22%2C%22initSources%22%3A%5B%7B%22stratum%22%3A%22user%22%2C%22models%22%3A%7B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%3A%7B%22isOpen%22%3Atrue%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%22%5D%2C%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%3A%7B%22url%22%3A%22https%3A%2F%2Fd2jfi34fqvxlsc.cloudfront.net%2Fmain%2Fdata%2F3d-tiles%2Fbldg%2F01100_sapporo%2Fnotexture%2Ftileset.json%22%2C%22highlightColor%22%3A%22%23D54A4A%22%2C%22style%22%3A%7B%22color%22%3A%7B%22conditions%22%3A%5B%5B%22isNaN%28%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D%29%22%2C%22color%28%27%23FFFFFF%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3D%3D%3D+null%22%2C%22color%28%27%23FFFFFF%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3C%3D+12.0%22%2C%22color%28%27%23382A54%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3C%3D+31.0%22%2C%22color%28%27%235A22C8%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3C%3D+60.0%22%2C%22color%28%27%23A675BE%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3C%3D+120.0%22%2C%22color%28%27%23F0D37B%27%2C1+%29%22%5D%2C%5B%22%24%7B%E8%A8%88%E6%B8%AC%E9%AB%98%E3%81%95%7D+%3C%3D+180.0%22%2C%22color%28%27%23FFCD00%27%2C1+%29%22%5D%2C%5B%22true%22%2C%22color%28%27%23F7FF00%27%2C1+%29%22%5D%5D%7D%7D%2C%22filters%22%3A%5B%7B%22name%22%3A%22%E9%AB%98%E3%81%95%E3%81%A7%E7%B5%9E%E3%82%8A%E8%BE%BC%E3%81%BF%22%2C%22minimumShown%22%3A7%2C%22maximumShown%22%3A200%7D%5D%2C%22knownContainerUniqueIds%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%22%5D%2C%22type%22%3A%223d-tiles%22%7D%2C%22%2F%22%3A%7B%22type%22%3A%22group%22%7D%2C%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%22%3A%7B%22knownContainerUniqueIds%22%3A%5B%22%2F%22%5D%2C%22type%22%3A%22group%22%7D%7D%2C%22workbench%22%3A%5B%22%2F%2FPLATEAU+%E3%83%87%E3%83%BC%E3%82%BF%E3%82%BB%E3%83%83%E3%83%88%2F%E5%8C%97%E6%B5%B7%E9%81%93%2F%E6%9C%AD%E5%B9%8C%E5%B8%82%2F%E5%BB%BA%E7%89%A9%E3%83%A2%E3%83%87%E3%83%AB%EF%BC%88%E6%9C%AD%E5%B9%8C%E5%B8%82%EF%BC%89%22%5D%2C%22timeline%22%3A%5B%5D%2C%22initialCamera%22%3A%7B%22west%22%3A139.63595%2C%22south%22%3A35.53022%2C%22east%22%3A139.95451%2C%22north%22%3A35.81489%7D%2C%22homeCamera%22%3A%7B%22west%22%3A139.63595%2C%22south%22%3A35.53022%2C%22east%22%3A139.95451%2C%22north%22%3A35.81489%7D%2C%22baseMaps%22%3A%7B%22defaultBaseMapId%22%3A%22%2Fbasemap%2F%2F%E5%85%A8%E5%9B%BD%E6%9C%80%E6%96%B0%E5%86%99%E7%9C%9F+%28%E3%82%B7%E3%83%BC%E3%83%A0%E3%83%AC%E3%82%B9%29%22%7D%2C%22viewerMode%22%3A%223d%22%2C%22currentTime%22%3A%7B%22dayNumber%22%3A2459947%2C%22secondsOfDay%22%3A63554.017%7D%2C%22showSplitter%22%3Afalse%2C%22splitPosition%22%3A0.5%2C%22stories%22%3A%5B%5D%7D%5D%7D)からの視界を確認してみます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/393d3620-aed9-3c6f-0dbd-79f8adb55be1.png)
これも何とか確認できます。塔までの直線距離は14.7kmで、目の良い方なら（頑張れば）肉眼で見れるかも知れません。
双眼鏡や望遠鏡があれば、もう少し楽に見ることができるかもしれません。

札幌市の中心部ではどうでしょう？
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/57ecae41-e02b-8a4d-bf4f-ff1e7ec717ec.png)
可視領域分析の結果では、札幌市の中心部でも塔が見える場所があるという結果が出ています。
可視領域がモザイク状になっているのは、使用した数値標高モデルの解像度が5mのものを使用したことによると思います。
なお、ここで使用した数値標高データは地表の標高を表すデータ（DEM; Digital Elevation Model）であり、建物や樹木などの高さを表す数値表層データ（DSM; Digital Surface Model）ではないため、モザイクの原因が建物の高さによるものではないと思います。DEMとDSMの違いについては国土地理院の[ページ](https://www.gsi.go.jp/chubu/minichishiki12.html)をご参照ください。

試しに、道庁の前庭からの視界を確認してみます。
![image.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/2986803/f118483f-620a-7473-698a-e9cfc277a7e5.png)
上の図の左側はPLATEAU VIEWで歩行者モードにした時の視界で、右側は同じ場所で上空に移動した場合の視界です。
歩行者モードでは、周辺のビルに視界を阻まれて、まったく塔を確認できません。
一方、同じ場所で上空から眺めた場合、ビルの間に塔の姿を確認できました。
札幌市中心部では、地上では塔を確認できないものの、高いビルなどからは、塔の姿を確認できるようです。

# 5. 考察～やってみてどうだったか？～

今回、QGISとGRASS GISで可視領域分析を行いました。その結果をPLATEAU VIEWを使い、模擬的に検証しました。
この模擬検証の結果ですが、塔が見えない場所と見える場所について、今の所、計算結果に間違いはなさそうです。
ただし、道庁の前庭のように、周囲に高いビルがある場合には、この計算結果をそのまま利用することはできませんでした。
高いビルなどの地物を考慮した可視領域分析は、別の方法を考える必要がありそうです。

また、今回、現地での確認が出来なかったため、PLATEAU VIEWを使って模擬的に視界を確認しました。
この模擬的な結果と実際とで、食い違いがある可能性があります。
現地での確認はいつかできればと考えています。

# 6. まとめ

QGISとGRASS GISを使った可視領域分析の方法と結果をまとめました。
今回の結果は現地での確認を行っていないため、計算結果との食い違いがある可能性があります。
模擬検証の範囲では、自然の地形による影響を考慮した可視領域では、計算結果が正しいことを確認しました。
また、都市中心部では高いビルにより視界が阻まれ、可視領域分析の結果と食い違いがあることが分かりました。
今回、数値標高モデルを使用して可視領域分析を行いました。
一方、心の中にある北海道百年記念塔は今回の方法では確認出来ませんでした。
心の中に、百年記念塔は見えるでしょうか？

# 謝辞

* 本記事を作成するに当たり参照させて頂いたWeb上の記事の執筆者様に御礼を申し上げます。
* 数値標高データの作成と維持を行っている国土地理院関係者の方々に感謝を申し上げます。
* QGISおよびGRASS GISの開発者コミュニティに御礼を申し上げます。
* [PLATEAU](https://www.mlit.go.jp/plateau/)データの整備と公開、利用を許可して下さっている国土交通省をはじめ、関係者の皆様に感謝を申し上げます。

# 参考文献

* [PLATEAU VIEW ver.1.1](https://www.mlit.go.jp/plateau/plateau-view-app/)（国土交通省）(https://www.mlit.go.jp/plateau/plateau-view-app/)

# 免責事項

* 筆者は可能な限り本記事内容の正確性の確保に努めますが、筆者は本記事内容の正確性を保証しません。
* 筆者は本記事の読者が本記事内容を用いて行う一切の行為について、何らの責任を負うものではありません。
* 筆者は本記事の読者が本記事内容を利用することで発生する損害に対し、責任を負わないものとします。