This paper presents an effective equal-area triangulated partition method (EATP) for parallel viewshed computation based on Xdraw algorithm that is a viewshed analysis algorithm with lower overhead. Our approach is to first partition the target domain into eight 45 • octants along a coordinate-axis direction and diagonal line direction around a given observer. Each octant remains mutually independent on data dependency based on the principle of Xdraw algorithm.Then, each octant is repartitioned into equal-area sub-regions with the shape of the triangle according to the requirements of system platform such as the number of computing nodes or processors and the size of memory of each node. Our method is beneficial to group these triangle sub-regions for determination of their dependent relationships during the computation of visibility and distribution of each triangle sub-region to a suitable node. The unified expression of the sub-regions is simple for implementation of parallel Xdraw viewshed algorithm. The experimental results are evaluated to demonstrate a distinct improvement in computation performance compared with other partition methods such as equal-angle and simple equal-area methods. KEYWORDS equal-area triangulated partition, parallel computing, terrain visibility analysis, Xdraw viewshed algorithm
INTRODUCTIONTerrain visibility between two cells with an elevation value, respectively, denotes whether they are visible each other on a terrain, which is also called inter-visibility analysis. A viewshed in terrain is a particular region of the terrain that is visible to a specified observer on the terrain. The viewshed analysis denotes the process whether the part of the terrain is determined to be visible to a specified observer at a given location on the terrain, which is widely applied spatial analysis to determine the visible regions of the terrain. 1,2 There are many applications in GISs (geographic information systems) including the optimal placement of observers according to visibility requirements, line-of-sight (LOS) communication problems, and terrain path planning with certain visibility properties. 3 However, the viewshed computation for a large-scale terrain viewshed analysis will be very time-consuming because a mass of interpolation operations had to be done in many viewshed analysis algorithms such as R2, R3, Xdraw, 4 reference plane algorithms, and so on. It has been a significant challenge of reducing a large amount of numeric calculation in the existing viewshed analysis algorithms. 5 In particular, the current DEM (digital elevation model) permits the applications to process large terrain with high resolution. 6 Viewshed calculation is usually implemented by emanating a LOS ray from the viewpoint or observer to a target point given an elevation value, respectively. 4,7 Considering the conventional sequential computation approaches are unable to provide effective support, the viewshed computation for large-scale terrain had to entail a massive number of LOS-based calculations. Hence, in ...