視野分析 Viewshed analysis

ENGLISH ONLY auto-stripped version of
中英版 Chinese+English version.

７公里內視野 7 km viewshed

points visible within 7 km

這裡（ｘ）所能見的（白），鄉鎮界（黑）Points visible from here (x) marked in white. Grid=230,235km W; 2670,2675 N km TWD67. Black=township lines. 抱歉初次用 Sorry it's my first try with GRASS GIS. Data source: DTM files were obtained from professors I met at GIS conferences.

Further viewshed 遠視野

越遠，能見這些點 Further away, we can see these points. Altitude (m), azimuth (deg), and distance (km)

標高 方位角 距離
2300  63.40 12.2 稍來山南端 Shaolai Shan south edge
3283 140.60 43.9 干卓萬山 Ganzhuowan Shan（八瘦）
3368 146.48 46.2 卓社大山 Zhuosheda Shan（九峨）
1064 147.68  5.5 黑田山 Heitian Shan
 556 152.7  25.9 虎子山（台灣地理中心） Huzi Shan (center of Taiwan)*
3275 153.18 74.7 內嶺爾山 Neilinger Shan（十崇）
2357 153.73 51.5 加年端山 Jianianduan Shan
3255 154.08 73.9 馬路巴拉讓山 Malubalarang Shan（六易）
3619 159.42 65.7 東郡大山 Dongjunda Shan（十崇）
2017 159.81 38.2 大尖山 Dajian Shan
3437 161.34 69.1 櫧山 Zhu Shan（七峭）
3483 162.65 77.5 馬利亞文路山 Maliyawenlu Shan（九平）
3231 163.48 68.3 無雙山 Wushuang Shan（九嶂）
2112 163.99 39.5 水社山（頂） Shuishe Shan
3765 164.34 76.1 馬博拉斯山 Mabolasi Shan（十峻）
2058 164.55 39.1 水社山（樁） Shuishe Shan
2483 164.65 49.3 治茆山北峰 Zhimao Shan north peak
2909 166.34 49.9 治茆山 Zhimao Shan
2863 167.34 50.8 治茆山南峰 Zhimao Shan south Peak
 775 171    36   日月潭 Riyue Tan (Sun Moon Lake)*
3952 173.25 79.4 玉山 Yushan*
1126 186.49  5.8 頭拒山 Touju Shan
250  238.82 81.8 麥寮六輕煙囪 Formosa Plastics sixth plant smokestack

空氣乾燥始得見，有的一年僅一二見。下雪望遠鏡得見。 These points visible only if air dry. Some can only be seen once or twice a year. Snow, if any, is visible in binoculars. * = 地平線之下 below horizon. (Shan means mountain.)

內嶺爾山至馬路巴拉讓山之嶺線位於花蓮縣。The Neilinger Shan to Malubalarang Shan ridge line is in Hualian County. 縣界以白色於此圖 Map: White: Nantou / Hualian County boundary. Green: visible areas. 3 x 3 km square: e=270 w=267 n=2611 s=2608 km TWD67，「十崇」之二。

I remembered to subtract the earth's curvature from the DTMs...

Before I forget, here's how one can compute the above:

$ perl -F, -anwMMath::Trig -e'use strict;
if ( $F[0] eq q(玉山) ) {
    @_ = Math::Trig::cartesian_to_spherical(
        $F[3] - 235480,
        $F[4] - 2675351,
        $F[5] - 777
    );
    @_ = ( $F[5], ( shift @_ ) / 1000, map { 90 - rad2deg($_) } @_ );
    printf "%d %6.2f %4.1f %.2f\n", @_[ 0, 2, 1, 3 ];
    exit;
}' twntwd67v1.csv
3952 173.25 79.4 2.29

目錄（Index）

積丹尼（Dan Jacobson）

Last modified: 2007-01-29 10:22:07 +0800