2. Background
This section presents a survey of existing methods for LiDAR data
processing. The state of the art in both LiDAR ground filtering and
DTM interpolation is presented.
2.1. Ground filtering algorithms
Ground filtering algorithms fall into four major classes: (1)
slope based, (2) interpolation based, (3) mathematical morphology
based, and (4) hybrid ground filters. A brief overview of these algorithms,
which from this point on may also be simply referred to as
filters, is given below.
2.1.1. Slope-based filters
Slope based filters work by comparing the slope between two
adjacent points to a pre-determined threshold value. The assumption
they make is that the slope (and hence the height) between
two adjacent ground points is smaller than the slope between a
ground point and a non-ground point. Hence, if the slope between
two points is larger than the threshold value, the point with a higher
altitude value is considered to be a non-ground point. An exam ple of a slope based filter is that proposed by Vosselman (2000).
The algorithm compares the height difference between two points
as a function of their distance apart against the acceptable height
difference. A point is regarded as a ground point only if there is
no other point such that the height difference between the two
points, computed as a function of their distance apart, is larger
than the allowed maximum.
Several studies have compared the performance of ground filtering
algorithms (e.g., Sithole and Vosselman, 2004; Liu, 2008;
Meng et al., 2010). These studies show that slope-based algorithms
work well with LiDAR data from relatively flat terrain but perform
poorly with data from complex and steep terrain. One reason for
the poor performance is the fact that in steep terrain the height
or slope between two ground points can be higher than the predetermined
threshold value, which causes many ground points to
be misclassified as non-ground. Moreover, determining a good
slope threshold may prove to be a very difficult task. To overcome
these difficulties, several adaptive filters have been proposed. For
example, to overcome the problem of determining a threshold, Sithole
(2001) proposed a modified slope-based filter in which the
threshold varied with the slope of the terrain. Using this approach
fewer ground points were misclassified.
2.1.2. Interpolation-based filters
While most slope-based algorithms work with two data points
at a time, interpolation-based filters operate on all points in the
data set at once. These filters are iterative algorithms which work
by fitting a surface on the data and iteratively filtering points based
on their residuals from the fitted surface. In the first iteration, an
averaging surface is fit into the data and the residuals of the data
points relative to the surface are computed. Ground points are
more likely to have negative residuals, so they are given more
weight in subsequent iterations. Kraus and Pfeifer (1998) were
the first to propose this approach. Their algorithm is based on linear
prediction, whereby in the first iteration all points are given
equal weights. A special function is used to compute weights for
each point based on their residuals. In each iteration, points with
large negative residuals are given maximum weights and thus they
attract the computed surface towards themselves. Points with
medium residuals get smaller weights and therefore have little
influence on the computed surface.
Like slope-based filters, interpolation-based filters perform
poorly in steep terrain. A similar method based on thin plate spline
surface interpolation has, however, been found to yield good results
in steep terrain. Mongus and Z˘alik (2012) have recently proposed
a ground filtering algorithm in which a thin plate spline
interpolated surface is used iteratively to approximate the ground
surface. The approach uses residuals of points from the surface in
each iteration with a gradually decreasing window size to filter
ground points.i While test results showed that the algorithm can
cope well with steep terrain, the results also showed that the algorithm
did not perform well with the low point densities commonly
found in forest data sets.
2.1.3. Mathematical-morphology-based filters
Mathematical-morphology-based filters are based on the concept
of mathematical morphology used in digital image analysis.
Unlike slope-based and interpolation-based filters, which work
on the original LiDAR data, filters based on mathematical morphology
transform the LiDAR data into digital images before performing
actual filtering. Filtering is done using morphological operations
such as opening and closing (Chen et al., 2007; Zhang et al., 2003).
These filters are simple and easy to implement and perform relatively
well in filtering surface objects, such as buildings, provided
that the data contain a high density of ground points. However, the
filters perform poorly with low point density in the data and steep terrains (Sithole and Vosselman, 2004; Liu, 2008; Meng et al.,
2010).
2.1.4. Hybrid filters
Some filtering algorithms do not directly fall into any of the
above classes. They either use a different approach to filtering or
make use of a combination of techniques used in other filters.
For example, Kobler et al. (2006) proposed a two-stage algorithm
whereby in the first stage negative outliers are removed using
some existing filtering algorithm (e.g., morphological filter) and
most non-ground points are then removed using a slope-based filter.
This stage produces partially filtered data which enters the second
stage of the algorithm. In the second stage, elevations at
unsampled locations are estimated from repeated independent
estimates of elevation samples taken from sampled points. By
combining different filtering methods, the algorithm performed
better than the individual filters.
2.2. DTM interpolation algorithms
With the exception of morphology-based filters, ground filtering
algorithms produce irregularly distributed ground points.
Interpolation is used with these points to estimate ground elevation
at places where measurements are missing. Interpolation
algorithms can be categorized based on the following criteria
(Ali, 2004): (1) exact or inexact (2) deterministic or stochastic,
and (3) local or global.
Exact interpolation algorithms produce DTM surfaces which pass
through every original data point, while approximative interpolation
methods produce surfaces which do not pass through every
original data point but rather follow the general trend in the data.
Examples of exact interpolation algorithms are inverse distance
weighting (IDW) and spline interpolation. An example of an inexact
interpolation algorithm is trend functions (surfaces). Deterministic
interpolation algorithms assume that every data point has influence
at an unsampled location that is inversely proportional to
its distance from the location. Stochastic algorithms, on the other
hand, make use of geostatistics to estimate a surface with some level
of uncertainty. Local interpolation algorithms use only neighboring
data points to estimate a value at an unsampled location,
while global algorithms use the entire data set to estimate values
at unsampled locations. Examples of commonly used interpolation
algorithms in DTM generation include inverse distance weighting
(IDW) (exact, deterministic, local), Kriging (approximative, global,
stochastic), and cubic splines (exact, local, deterministic) (Ali,
2004; Liu, 2008).
The choice of an interpolation algorithm to be implemented depends
on several factors, including the quality of data and required
level of accuracy. For example, IDW offers high accuracy but only
when data density is high.
2 . พื้นหลัง
ส่วนนี้นำเสนอวิธีการสำหรับการประมวลผลข้อมูลสำรวจของที่มีอยู่
ไลด้า รัฐของศิลปะทั้งในดินและในการกรอง LIDAR DTM า
.
2.1 . ขั้นตอนวิธีขั้นตอนวิธีการกรองกรองพื้นดิน
ตกสี่ประเภทหลัก : ( 1 )
ความชันตาม ( 2 ) ( 3 ) การใช้โครงสร้างทางคณิตศาสตร์
ตาม ( 4 ) และกรองดินผสมภาพรวมโดยย่อของขั้นตอนวิธีการเหล่านี้
ซึ่งจากจุดนี้อาจเป็นเพียงแค่เรียกว่า
ตัวกรองจะได้รับด้านล่าง .
2.1.1 . ใช้ตัวกรองตัวกรองตามความลาดชัน
งานโดยเปรียบเทียบความลาดชันระหว่างที่อยู่ติดกันสอง
จุดที่กำหนดไว้ล่วงหน้าของค่า อัสสัมชัญ
พวกเขาให้เป็นที่ลาดชัน ( และดังนั้น ความสูงระหว่าง
)สองจุดพื้นดินที่อยู่ติดกันมีขนาดเล็กกว่าความลาดชันระหว่างพื้นดินและพื้นดิน
จุดที่ไม่ใช่จุด ดังนั้น ถ้าความชันระหว่าง
สองจุดมีขนาดใหญ่กว่าค่าเกณฑ์จุด ด้วยมูลค่าสูงกว่าระดับความสูง
ถือว่าไม่ใช่พื้นดินจุด กรองตามสอบ ple ของความลาดชันที่เสนอโดย vosselman ( 2000 ) .
ขั้นตอนวิธีการเปรียบเทียบความสูงระหว่างสองจุด
เป็นฟังก์ชันของระยะทางของพวกเขาแยกกับความสูงต่างกัน
ที่ยอมรับได้ จุดที่ถือว่าเป็นจุดเดียว หากมีพื้น
ไม่จุดอื่น ๆเช่นที่ความสูงระหว่างสอง
จุด คำนวณเป็นฟังก์ชันของระยะทางของพวกเขาออกจากกัน มีขนาดใหญ่กว่าได้รับอนุญาตสูงสุด
.
หลายการศึกษาเปรียบเทียบประสิทธิภาพของขั้นตอนวิธีการกรอง
พื้น ( เช่น sithole vosselman และ ,2004 , หลิว , 2008 ;
เมิง et al . , 2010 ) การศึกษาเหล่านี้แสดงให้เห็นว่าความชันตามขั้นตอนวิธี
ทำงานได้ดีกับข้อมูล LIDAR จากภูมิประเทศที่ค่อนข้างแบน แต่แสดง
ไม่ดีด้วยข้อมูลจากที่ซับซ้อน และ ภูมิประเทศที่สูงชัน เหตุผลหนึ่งสำหรับ
ประสิทธิภาพที่ดีคือความจริงที่ว่าในภูมิประเทศ ความสูงหรือความลาดชัน
ระหว่างสองพื้นดินจุดสามารถสูงกว่าเกณฑ์ที่กำหนดไว้ล่วงหน้า
ค่าซึ่งสาเหตุพื้นดินหลายจุด
จะ misclassified เป็นบนพื้นดิน นอกจากนี้ การกำหนดเกณฑ์ลาดดี
อาจพิสูจน์ให้เป็นงานที่ยาก เพื่อเอาชนะปัญหาเหล่านี้
, ตัวกรองแบบต่าง ๆ ได้ถูกเสนอขึ้นมา สำหรับ
ตัวอย่าง เพื่อเอาชนะปัญหาการกำหนดเกณฑ์ , sithole
( 2001 ) ได้เสนอการใช้ตัวกรองที่
ความลาดชันเกณฑ์ตามความลาดของภูมิประเทศ การใช้วิธีการนี้ มีจุด misclassified พื้นดิน
.
2.1.2 . การใช้ตัวกรอง
ในขณะที่ส่วนใหญ่ความลาดชันขั้นตอนวิธีการทํางานตามด้วยสองจุดข้อมูล
ครั้ง ตัวกรองการใช้ผ่าตัดทุกจุดใน
ชุดข้อมูลในครั้งเดียว ตัวกรองเหล่านี้จะใช้ซ้ำซึ่ง
โดยการปรับผิวบนและกรองข้อมูลซ้ำจุดตามค่าของพวกเขาจากการติดตั้ง
บนพื้นผิว ในรูปแรกเป็น
เฉลี่ยพื้นผิวพอดีกับข้อมูลและความคลาดเคลื่อนของข้อมูล
จุดเทียบกับพื้นผิวจะคํานวณ จุดดิน
น่าจะมีความคลาดเคลื่อนเชิงลบดังนั้นพวกเขาจะได้รับน้ำหนักมากขึ้น
ในตามมาซ้ำ .และ เคราส์ ไฟเฟอร์ ( 1998 ) ได้เสนอวิธีการ
ครั้งแรกนี้ ขั้นตอนวิธีการของพวกเขาจะขึ้นอยู่กับตรงคำทำนาย
ซึ่งในรูปแรก ทุกจุดจะได้รับ
น้ำหนักที่เท่ากัน เป็นฟังก์ชันพิเศษที่ใช้เพื่อคำนวณน้ำหนัก
แต่ละจุดขึ้นอยู่กับค่าของพวกเขา ในแต่ละซ้ำ จุดที่มีค่าเป็นลบจะได้รับน้ำหนักมาก
สูงสุด และดังนั้นจึงดึงดูดคำนวณพื้นผิวต่อตนเอง คะแนนกับค่ารับน้ำหนักปานกลาง
มีขนาดเล็กและดังนั้นจึงมีอิทธิพลน้อยในการคำนวณพื้นผิว เช่นตัวกรองตัวกรอง
ลาดตามการใช้แสดง
ไม่ดีในภูมิประเทศที่สูงชัน วิธี ที่คล้ายกันอยู่ในบางจานสลัก
พื้นผิวการแก้ไขได้ อย่างไรก็ตาม พบว่าผลผลิต
ผลลัพธ์ที่ดีในภูมิประเทศที่สูงชันmongus และ Z ˘ลิก ( 2012 ) เพิ่งได้เสนอขั้นตอนวิธีการกรอง
พื้นที่บางจานสลัก
ขัดผิวใช้ซ้ำประมาณผิวดิน
การใช้ค่าคะแนนจากพื้นผิวในแต่ละซ้ำด้วย
ค่อยๆลดขนาดหน้าต่างกรอง
จุดพื้น ฉันในขณะที่ผลการทดสอบพบว่าอัลกอริทึม
รับมือกับภูมิประเทศที่สูงชัน ผลการศึกษายังพบว่าขั้นตอนวิธี
ทำไม่ดีกับจุดมีความหนาแน่นต่ำโดยทั่วไปที่พบในชุดข้อมูลป่าไม้
.
ทาง . โครงสร้างทางคณิตศาสตร์พื้นฐานทางคณิตศาสตร์ที่ใช้ตัวกรองตัวกรอง
จะขึ้นอยู่กับแนวคิดของคณิตศาสตร์สัณฐานวิทยาที่ใช้ในการวิเคราะห์ภาพดิจิตอล ซึ่งแตกต่างจากการใช้ใช้ความลาดชันและ
ตัวกรองทำงานซึ่งในข้อมูล LIDAR เดิมกรองตามโครงสร้างทางคณิตศาสตร์การแปลงข้อมูล LIDAR
เป็นภาพดิจิตอลก่อนการปฏิบัติจริงกรอง การกรองโดยใช้ลักษณะทางสัณฐานวิทยา
เช่นการเปิดและปิด ( Chen et al . , 2007 ; Zhang et al . , 2003 ) .
ตัวกรองเหล่านี้จะง่ายและง่ายต่อการใช้และดำเนินการค่อนข้าง
ดีในการกรองพื้นผิววัตถุ เช่น อาคารโดย
ที่ข้อมูลมีความหนาแน่นสูงของจุด G อย่างไรก็ตาม
กรองทํางานกับความหนาแน่นของจุดต่ำในข้อมูล และภูมิประเทศที่สูงชัน ( sithole และ vosselman , 2004 ; Liu , 2008 ; เมิง et al . ,
( )
2.1.4 . ตัวกรองบางขั้นตอนวิธีลูกผสม
ไม่ตรงอยู่ในใด ๆของ
เหนือชั้น พวกเขาใช้วิธีการที่แตกต่างกันเพื่อกรองหรือ
ใช้ประโยชน์จากการรวมกันของเทคนิคที่ใช้ในตัวกรองอื่น ๆ .
ตัวอย่างเช่น kobler et al . ( 2549 ) ได้เสนอขั้นตอนวิธีในการแจกแจงแบบ
โดยในระยะแรกลบจะถูกลบออกโดยใช้ขั้นตอนวิธีการกรอง
บางอย่างที่มีอยู่ ( เช่น ลักษณะทางสัณฐานวิทยาและจุดพื้นดินไม่กรอง )
ส่วนใหญ่จะลบออกแล้วใช้ความลาดชันกรอง .
เวทีนี้ผลิตบางส่วนกรองข้อมูล ซึ่งจะเข้าสู่
2ขั้นตอนของขั้นตอนวิธี ในขั้นตอนที่สอง , เอนไซม์ที่
สถานที่ unsampled ประมาณการจากแบบประเมินซ้ำ
ตัวอย่างการถ่ายจากตัวอย่างคะแนน โดย
รวมที่แตกต่างกันวิธีการกรองขั้นตอนวิธีปฏิบัติ
ดีกว่าตัวกรองแต่ละ .
2.2 . DTM ขั้นตอนวิธีการแก้ไข
ยกเว้นตัวกรองของกรอง
ตามพื้นดินขั้นตอนวิธีการผลิตการกระจายไม่สม่ำเสมอพื้นดินจุด การใช้กับจุดนี้
ความสูงประมาณพื้นที่ ที่วัดหายไป ขั้นตอนวิธีการแก้ไข
สามารถแบ่งตามเกณฑ์
ต่อไปนี้ ( Ali , 2004 ) : ( 1 ) แน่นอนหรือไม่ละเอียด ( 2 ) ติดตั้งใช้งานหรือ Stochastic
( 3 ) ท้องถิ่นหรือระดับโลก .
ขั้นตอนวิธีการแก้ไขที่แน่นอนผลิต DTM พื้นผิวซึ่งผ่าน
ผ่านข้อมูลเดิมทุกจุด ในขณะที่วิธีการประมาณพื้นผิวที่ไม่ผลิต
ต้นฉบับข้อมูลผ่านทุกจุด แต่ตามแนวโน้มทั่วไปในข้อมูล .
ตัวอย่างของขั้นตอนวิธีการแน่นอนตรงกันข้ามระยะทาง
( ( idw ) และเส้นโค้งการประมาณค่าในช่วง ตัวอย่างของการไม่ละเอียด
การประมาณค่าฟังก์ชันขั้นตอนวิธีคือแนวโน้ม ( พื้นผิว ) ขั้นตอนวิธีการติดตั้งใช้
สมมติว่าทุกจุดข้อมูลมีอิทธิพล
ที่ unsampled ที่ตั้งที่เป็นปฏิภาคผกผันกับ
ห่างไกลจากสถานที่ สโมสรฟุตบอลไกเซอร์สเลาเทิร์น บนมืออื่น ๆ
, ให้ใช้ธรณีสถิติประมาณการพื้นผิวที่มีระดับบาง
ของความไม่แน่นอนขั้นตอนวิธีการเพื่อนบ้าน
ท้องถิ่นใช้เฉพาะจุดข้อมูลประมาณการค่าสถานที่ unsampled
, ในขณะที่ขั้นตอนวิธีทั่วโลกใช้ทั้งชุดข้อมูลเพื่อประมาณค่า
unsampled สถานที่ . ตัวอย่างที่ใช้กันทั่วไปของขั้นตอนวิธีในการสร้างเครื่องกำเนิดไฟฟ้า
( idw น้ำหนักรวมระยะทางผกผัน ) ( แน่นอน ติดตั้งเครื่องท้องถิ่น ) คริกกิ้ง ( ประมาณ , ทั่วโลก ,
Stochastic )และลูกบาศก์ 500 ( แน่นอน , ท้องถิ่น , deterministic ) ( Ali ,
2004 , หลิว , 2008 ) .
ทางเลือกของการสอดแทรกวิธีการที่จะใช้ขึ้นอยู่กับ
ปัจจัยหลายประการรวมถึงคุณภาพของข้อมูลและระดับที่ต้องการ
ของความถูกต้อง ตัวอย่างเช่น idw มีความแม่นยำสูงแต่เพียง
เมื่อความหนาแน่นของข้อมูลสูง
การแปล กรุณารอสักครู่..
