Viewshed analysis


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