NEW HORIZONTAL ACCURACY ASSESSMENT TOOLS AND TECHNIQUES FOR LIDAR DATA INTRODUCTION
|
|
- Magdalene Margaret Stevenson
- 5 years ago
- Views:
Transcription
1 NEW HORIZONTAL ACCURACY ASSESSMENT TOOLS AND TECHNIQUES FOR LIDAR DATA John A. Ray Ohio Department of Transportation Columbus, Ohio Lewis Graham GeoCue Corporation Madison, Alabama ABSTRACT Much of the focus to date in LiDAR geometric data assessment has been on vertical accuracy. However, with the increasing use of LiDAR derived images such as LiDAR orthos and LiDAR synthetic stereo models, knowledge of the horizontal accuracy of LiDAR data is essential to project error modeling. We present new assessment tools and techniques for both assessing the horizontal geometric accuracy characteristics of LiDAR data and for effecting global data adjustments. We provide practical examples in highway project analysis as well as ideas for future work. Application of these and similar tools show great promise for deploying LiDAR as an engineering tool of choice. INTRODUCTION The Ohio Department of Transportation (ODOT) performs a portion of the remote sensing and digital mapping for the highway improvement projects delivered annually by ODOT in an effort to fulfill its mission. In 2004, ODOT implemented airborne LiDAR technology as a means to reduce the project development duration and improve the vertical accuracy of the spatial data generated from remote sensing operations. Within a few months, the processes and procedures established were capable of yielding spatial data with a vertical RMSE (Root Mean Square Error) approaching 0.10 feet on hard surfaces. The procedures were the result of research and testing of various flight operation parameters, GPS attributes, LiDAR control methodologies, and data processing techniques. Some of the procedures adopted were somewhat unique, such as adjustment of the LiDAR dataset to points collected using special targets as control (see Figure 1). Figure 1. LiDAR control target being positioned with GPS.
2 When delivering any geospatial product, it is important to be able to report the accuracy of the final product(s) being delivered, but it is critical that the accuracy be known when using the data for engineering purposes. In an effort to express the quality of the LiDAR derived data, ODOT began providing a Quality Control Report for each project. The Quality Control Report details the points used to adjust and determine the vertical accuracy of the LiDAR generated TIN (Triangulated Irregular Network). Reporting the vertical accuracy has also been beneficial in helping educate designers about the potential uses for LiDAR data. After gaining some additional experience, ODOT transitioned to using pavement surfaces for LiDAR control instead of the specially designed circular targets. The targets functioned well, but were labor intensive to install and required placement just prior to the flight, which became problematic. As ODOT s LiDAR process continued to evolve, multiple quality control procedures were implemented. One such quality control procedure is to overlay the planimetric features (including pavement markings) collected using photogrammetric methods on top of the LiDAR intensity image also known as a LiDAR ortho ( Figure 2). This simple check helps identify any potential problems with the GPS/INS system, scale factor and coordinate transformation issues. It soon became apparent that the visual agreement between the LiDAR data and the planimetric features were very good where a clear difference in the intensity of the LiDAR data was evident. Survey control points were also used to check for horizontal positioning. Note that the density of the LiDAR data (typically 12 pts. /m 2 ) certainly has a positive impact on the ability to identify features such as pavement markings. Figure 2. Survey control points on a LiDAR ortho. With ODOT s desire to report the horizontal accuracy in addition to the vertical accuracy of the final products being created, ODOT discussed the issue with The Ohio State University s Center for Mapping. The discussions resulted in The Ohio State University (OSU) proposing a research and development project entitled Airborne LiDAR Reflective Linear Feature Extraction for Strip Adjustment and Horizontal Accuracy Determination. The goal of this project was to automatically extract features within the vicinity of control points in order to determine the horizontal accuracy of the LiDAR dataset. Even though OSU is very good at research and development, they are not in the commercial software and support business. With the ultimate goal being the integration of the horizontal accuracy determination and adjustment capability into an existing workflow, the new functionality would ideally be placed in a product that was currently used within the process. After consideration of many factors, we discussed the idea with GeoCue Corporation of Madison Alabama. The GeoCue software was a logical place to house the new functionality for
3 several reasons. One of the primary reasons was the existing use of the GeoCue software as a LiDAR processing and management tool. The GeoCue software was also used by many other LiDAR data providers which would then create an opportunity for reporting the horizontal accuracy component on an industry-wide basis. When GeoCue was approached regarding the issue, they expressed an interest in developing the functionality within their suite of software products. A project soon followed to create a module known as FMAS (Feature Matching and Adjustment System). While the primary objective of the GeoCue effort was to integrate the OSU developed algorithms seamlessly into their software suite, GeoCue developed other interactive techniques for measuring and adjusting data based on measuring point and linear features, thereby increasing system functionality and the end users options. Of particular utility, and the subject of this paper, is a novel algorithm developed by GeoCue that allows measurement and correction of data on non-conjugate line segments. LINEAR CONTROL GeoCue Corporation has subsequently developed a system to measure control points and associated conjugate points in orthos and LIDAR Orthos. This system (The Feature Matching and Adjustment System CuePac, or the FMAS CuePac) is intended for assessing the horizontal and vertical accuracy of LIDAR data and the horizontal accuracy of orthos. The system can also compute a 6 parameter linear transform in the horizontal components and a simple shift in Z. The user can optionally apply these corrections to a LIDAR data set using the distributed processing system of the GeoCue software system (GeoCue). The system is primarily designed to use externally provided control data in the form of FMAS Control Features. FMAS supports: Points 2 point line segments Polylines Simple polygons Hierarchical combinations of the above (e.g. traffic stop bar) As we tested the initial version of the system, we realized a particularly valuable workflow of comparing LIDAR data horizontal position by using a known control orthophoto. Basically the flow is as shown in Figure 3 1. During this investigation we recognized the value in developing an algorithm to solve for an affine correction using only parallel line segments (conjugate lines). 1 Of course, one could argue that the wrong elevations have been used in generating the orthos since the LIDAR data are assumed off in X, Y. This is a phenomenon that will be quantified in a future paper.
4 Figure 3. Adjusting LIDAR with metric orthophotos. Conjugate Lines When implementing the flow of Figure 3 we quickly discovered that it is very difficult to accurately digitize point features (such as the end of a paint stripe) in LIDAR intensity images. What looked to be more fruitful was digitizing line segments. However, accurately digitizing the start and end points of the conjugates of the line segments in LIDAR orthos proved unreliable. This is primarily due to the resolution difference between these sources and the fact that often the best control is a segment of a line such as an edge-of-pavement mark. An example of this difficulty is illustrated in Figure 4. On the left is a 1.5 m LIDAR ortho derived from 2 m post spaced LIDAR data. On the right is a USGS 0.33 m natural color ortho. Notice that linear features are clearly visible in the LIDAR ortho. However, precisely delineating the start or end point of a line segment is not possible due to the lower resolution of the LIDAR ortho. Figure 4. LIDAR Ortho (left) compared to a natural color Ortho (right). What did look promising was to digitize a control line segment such that the segments were on the same linear feature but not necessarily conjugate (see Figure 5). Notice that we have digitized a segment along the same linear
5 feature (the white stripe down the west side of a runway) but these segments are not conjugates since they represent different regions of the conjugate line 2. Figure 5. Control segment (vertical orange on right) and Conjugate (green line on left). Note that it is important to realize that the control and conjugate lines are conjugate but not the segments. We then faced the task of developing algorithms to use these semi-conjugate line segments in our point driven system. Thus the statement of the problem is: given a control segment and a conjugate line on which the conjugate segment is located (but whose endpoints are currently unknown), find an algorithm to determine the end points of the conjugate system. Considering the Segments Consider the example of Figure 6. Here we have 4 control segments and 4 target segments. The target segments represent segments of linear features where the linear features are the conjugate features of Control lines. Since these target segments are digitized along conjugate lines but not in conjugate locations, we will refer to these collections as Control Segments (denoted by C ) and Related Segments (denoted by R ). The blue vector V represents the average relative displacement between the control data set and the conjugate data set. C 2 C 1 R 2 R 1 C 3 C 4 R 3 V R 4 Figure 6. Example of Control Segments and Related Segments. 2 Recall that a line is defined by two points and extends to infinity in both directions whereas a line segment terminates at each end by the defining points.
6 We will use the notation defined in Table 1. Table 1. Notation Symbol C i R i c i,k r i,k G i g i,k L i m i Definition The i th control segment The i th related segment The k th vertex (node or point) on the i th control segment The k th vertex (node or point) on the i th related segment The i th conjugate segment The k th vertex (node or point) on the i th conjugate segment The i th line (as opposed to segment) The slope of the i th segment (or line) An approximation or estimate of x It is very important to appreciate that the user, in general, has not attempted to digitize conjugate segments in the target. The only requirement is that the user digitizes a segment on the same linear feature as was used to digitize the control segment. This means that R i is not the conjugate segment of C i. Thus we know that in the case of no error between the control image and the target image,. C i and R i will both lie on the same line but they will not be (in general) coincident. This is why we say that R i is related to C i as opposed to being the conjugate. An illustration of our example, post adjustment, is depicted in Figure 7. Note that even after a fairly good adjustment (as indicated by the diminished length of the displacement vector, V), the related segments are not coincident with the control segments. This is due to the fact that the related segments were digitized on lines conjugate to the lines defined by the control segments but not in conjugate locations along these lines. C 1 R 1 R 2 C 2 C 3 C 4 R 3 V R 4 Figure 7. The appearance of adjusted data without conjugate segments. The Mathematical Approach We know from the specification of the problem that a control segment, C i, defines a line L i and that, with no error, the related segment in the data to be adjusted (R i ) should be coincident with L i. Thus we could use a linear constraint to find adjustment parameters. However, a direct linear constraint will not solve for a scale difference. This makes intuitive sense because scaling something of infinite extent and zero width (a line) yields the same line! Our approach has been to attempt to find the true conjugates of the control segments and then apply standard point adjustment algorithms to these data.
7 Consider a magnified view of the example data. If we can approximate the maximum displacement V m, we know each conjugate point g i,j must be within a distance V m of its control point c i,j (we refer to this as the radius of ambiguity) 3. Furthermore, g i,j must be on the line defined by R i (within digitizing error). This is illustrated in Figure 8. The conjugate points must lie on the orange lines in this figure. With this approximation of the locations of the conjugate points, we will now develop an algorithm to refine he estimates. The approach that we will take is to use orthogonal independence to separately estimate the x and y components of each conjugate point. r = V m C 1,2 g 1, 2 lies somewhere on this segment C 1,1 V m Figure 8. Radius of ambiguity of conjugate points. Consider horizontal and vertical segments as shown in Figure 9. Here we have depicted the radius of ambiguity of each conjugate point and sketched the possible location of each in orange. Now consider the special circumstances that exist for these lines. 3 There is a bit of a subtle point here. If conjugate points are collected by human interaction or an automated technique such as correlation, one cannot make this statement (due to outlier collection errors). However here we are discussing the true and as of yet uncollected conjugate points and thus we can quantify a radius of ambiguity.
8 C 2 R 2 V C 3 Radius of ambiguity R 3 Figure 9. The special case of horizontal and vertical segments. We know, for a horizontal control line, that the y component of the conjugate point is given by the y value of the conjugate line (which is a constant since the line is horizontal). Similarly, the x component of the conjugate of a vertical line is given by the x value of the conjugate vertical line. Thus you will note that by selecting horizontal lines we can remove (to within the error of digitizing) the ambiguity in the y coordinate and by selecting vertical lines we remove the ambiguity in the x coordinate. This leads us to use the slope of the line to motivate the following arguments: (1) (2) where, are the orthogonal components of the standard deviation 4. Note that we are assuming that the slope of the control line is approximately the same as the slope of the conjugate line. In other words horizontal lines provide strength to the y component of the fit and vertical lines to the horizontal fit. Conversely, a horizontal line provides no information about the x component of a conjugate point ) and a vertical line lends no information to the y component of a conjugate point ). An inspection of Figure 9 lends intuitive credence to this. First Pass Conjugate Approximation The previous observations provide the basis for an algorithm that will allow us to approximate the conjugate points (the g i,k ). Our first approximations for the conjugate points are the points defined by the intersection of a line perpendicular to the control point C i,k and intersecting R i. This is illustrated in Figure 10. We denote this first set of approximated conjugate points as where i is the i th related line and j is 1 or 2 representing the end point of the approximated conjugate line segment. Note that these approximations can be as far as V m /2 from their correction locations and thus will require a refinement. 4 One could argue that this should be more properly the sample standard deviation but another subtly occurs here. This will be addressed in a future paper. We do use sample standard deviation in the least squares adjustments.
9 c 1,1 C 1 R 1 c 3,1 c 1,2 C 3 R 3 c 3,2 Figure 10. First approximation to conjugate points. Second Pass Refinement We next need to quantify, for each of these approximated conjugate points. We do this as follows: Let be the unsigned angle between the ith line defined by Ri and the horizontal. Thus. We now want to compute weights for the x and y components of the conjugate points such that (3) since we are going to solve for the initial refinement in a weighted, separable least squares adjustment. Our base model 5 uses: (4) (5) We chose this because the sine and cosine are obviously the slope components with the advantage of not blowing up at the extremes of horizontal and vertical. We needed a much lower weight on diagonal lines and the reciprocal relationship between W(x) and W(y). We used the 4 th power (rather than 2 nd ) because we wanted a rapid off-axis dive in the weight. The weights as a function of are plotted in Figure Weight equations are a topic for further research.
10 1 W x ( m) W y ( m) m Figure 11. Independent coordinate least square weights. We next computed a 6 parameter affine transformation: using a weighted least squares adjustment. The least squares algorithm selected must support independent x and y weights since this is the trick of the first pass adjustment. Note that true conjugate points can be mixed in with this phase of the adjustment (if available) and, indeed, they add strength to the solution. These normal conjugate points should have weights inversely proportionate to their sample standard deviations (if known), normalized to the same weight range as the pseudo conjugate points ( ). If the deviations are unknown, use W = 1.0. The result of the least square solution will be a transform that represents the true adjustment between the control system and the target system. Thus this adjustment can now be used to construct conjugate points for the related segments. If is the six parameter transform that projects the c 1,1 C 2 conjugate space (G) onto the control space (C), then (since we are using a six C 1 R 2 parameter transform you will need to augment the T matrix) projects C onto G. R 1 C 3 We now make a second pass adjustment of the approximated conjugate C 4 R 3 points by the following procedure. Project all of the control points through yielding a result similar to R 4 Figure 12 where we denote the back-projected control V points as green diamonds. Figure 12. Control points (block) projected into conjugate space (green diamonds). (6) (7)
11 Our projected conjugate points ( ) will vary somewhat from the digitized Related lines due to the general measurement error, the fact that the error is not really exactly modeled by a six parameter transform and the distribution of the original pseudo conjugate points. However, our second pass conjugate point g i,k must be on R i (again, to within our digitizing error) since this is the original assumption of collecting conjugate lines. Thus we perform a second pass adjustment by the following procedure. Construct a line through the points c i,j,. Find where this line intersects R i. Make this the new (and final for the purposes of this paper) conjugate point g i,j. This is illustrated, greatly magnified, in Figure 13. c 1,1 g 1,1 C 1 R 1 c 1,2 g 1,2 Figure 13. Constructing refined conjugate points. These new conjugate points are now considered the true conjugates. Final Adjustment The newly computed g i,j are now treated as the true conjugate points to the original control segment end points c i,j. The weights of the deviations of these computed conjugate points are now set such that. These control, conjugate pairs are now used in a second pass, equally weighted least squares adjustment with all other ordinary measurement pairs added to the data set. This yields the final transform T. Any major changes in T between the intermediate adjustment and the final adjustment are indicative of a weak original geometry. In the example of adjusting a LiDAR data set to known linear control, the final process would be to apply the resulting transform T to the each and every LiDAR data point. We do this using a distributed processing function within the GeoCue Enterprise system, allowing the LiDAR files to be processed in parallel. CONCLUSIONS The Ohio Department of Transportation has begun using the horizontal accuracy determination and correction functionality contained in the GeoCue LiDAR FMAS CuePac module. To date, five projects have been analyzed with the software. The current process involves running the algorithms after the strip alignment task has been performed. The control points are loaded and the software automatically drives the view to the next control point in the list for assessment. The features within the LiDAR data are manually digitized in order to facilitate comparison to the control points. While manually digitizing the LiDAR features does create some error, careful selection of features with adequate intensity delineation minimizes the magnitude of these errors. Additionally, as mentioned earlier, avoiding the endpoints of features such as stop bars where LiDAR data may not have sufficiently defined the feature, is recommended. The existing FMAS CuePac is an intuitive and user friendly tool capable of comparing control points and LiDAR intensity features. The software provides the necessary statistics to analyze and adjust the data to remove
12 any global bias. The horizontal accuracy can now be easily quantified. Specifying the horizontal as well as the vertical accuracy of LiDAR data should quickly become the standard. REFERENCES De Berg,M., Van Krevald, M., Overmars, M., Schwarzkopf, O., Computational Geometry, Springer-Verlag, Berlin. Bevington, P., Data Reduction and Error Analysis for the Physical Sciences, McGraw-Hill, New York. Press, W., et al., Numerical Recipes in C, Cambridge University Press, New York.
Summary of Research and Development Efforts Necessary for Assuring Geometric Quality of Lidar Data
American Society for Photogrammetry and Remote Sensing (ASPRS) Summary of Research and Development Efforts Necessary for Assuring Geometric Quality of Lidar Data 1 Summary of Research and Development Efforts
More informationTools, Tips, and Workflows Breaklines, Part 4 Applying Breaklines to Enforce Constant Elevation
Breaklines, Part 4 Applying Breaklines to l Lewis Graham Revision 1.0 In the last edition of LP360 News, we discussed the creation of 3D breaklines. Recall that, for our purposes, a 3D breakline is a vector
More informationU.S. Geological Survey (USGS) - National Geospatial Program (NGP) and the American Society for Photogrammetry and Remote Sensing (ASPRS)
U.S. Geological Survey (USGS) - National Geospatial Program (NGP) and the American Society for Photogrammetry and Remote Sensing (ASPRS) Summary of Research and Development Efforts Necessary for Assuring
More informationLIDAR Quality Levels Part 1
Lewis Graham, CTO GeoCue Group The LIDAR community has been gradually adopting the terminology Quality Level to categorize airborne LIDAR data. The idea is to use a simple scheme to characterize the more
More informationCCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12
Tool 1: Standards for Mathematical ent: Interpreting Functions CCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12 Name of Reviewer School/District Date Name of Curriculum Materials:
More informationVisualization and Analysis of Inverse Kinematics Algorithms Using Performance Metric Maps
Visualization and Analysis of Inverse Kinematics Algorithms Using Performance Metric Maps Oliver Cardwell, Ramakrishnan Mukundan Department of Computer Science and Software Engineering University of Canterbury
More informationTerrestrial Laser Scanning: Applications in Civil Engineering Pauline Miller
Terrestrial Laser Scanning: Applications in Civil Engineering Pauline Miller School of Civil Engineering & Geosciences Newcastle University Overview Laser scanning overview Research applications geometric
More informationEXTREME POINTS AND AFFINE EQUIVALENCE
EXTREME POINTS AND AFFINE EQUIVALENCE The purpose of this note is to use the notions of extreme points and affine transformations which are studied in the file affine-convex.pdf to prove that certain standard
More informationAmarillo ISD Math Curriculum
Amarillo Independent School District follows the Texas Essential Knowledge and Skills (TEKS). All of AISD curriculum and documents and resources are aligned to the TEKS. The State of Texas State Board
More informationIntegrated Mathematics I Performance Level Descriptors
Limited A student performing at the Limited Level demonstrates a minimal command of Ohio s Learning Standards for Integrated Mathematics I. A student at this level has an emerging ability to demonstrate
More informationMathematics High School Geometry An understanding of the attributes and relationships of geometric objects can be applied in diverse contexts
Mathematics High School Geometry An understanding of the attributes and relationships of geometric objects can be applied in diverse contexts interpreting a schematic drawing, estimating the amount of
More informationGeometry. Cluster: Experiment with transformations in the plane. G.CO.1 G.CO.2. Common Core Institute
Geometry Cluster: Experiment with transformations in the plane. G.CO.1: Know precise definitions of angle, circle, perpendicular line, parallel line, and line segment, based on the undefined notions of
More informationCS231A Course Notes 4: Stereo Systems and Structure from Motion
CS231A Course Notes 4: Stereo Systems and Structure from Motion Kenji Hata and Silvio Savarese 1 Introduction In the previous notes, we covered how adding additional viewpoints of a scene can greatly enhance
More information1
Zeros&asymptotes Example 1 In an early version of this activity I began with a sequence of simple examples (parabolas and cubics) working gradually up to the main idea. But now I think the best strategy
More informationTools, Tips and Workflows LP360 Under the Hood - The Conflation Point Cloud Task LP360
Tools, Tips and Workflows LP360 Lewis Graham Revision 1.0 l In the January issue of GeoCue News, we discussed the new Point (PCT) overhaul that we released in LP360 2013.2. The focus of the last article
More information6.1 Leaping Lizards! A Develop Understanding Task
6.1 Leaping Lizards! A Develop Understanding Task Animated films and cartoons are now usually produced using computer technology, rather than the hand drawn images of the past. Computer animation requires
More informationAnswers to practice questions for Midterm 1
Answers to practice questions for Midterm Paul Hacking /5/9 (a The RREF (reduced row echelon form of the augmented matrix is So the system of linear equations has exactly one solution given by x =, y =,
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More informationSmallest Intersecting Circle for a Set of Polygons
Smallest Intersecting Circle for a Set of Polygons Peter Otfried Joachim Christian Marc Esther René Michiel Antoine Alexander 31st August 2005 1 Introduction Motivated by automated label placement of groups
More informationGlossary of dictionary terms in the AP geometry units
Glossary of dictionary terms in the AP geometry units affine linear equation: an equation in which both sides are sums of terms that are either a number times y or a number times x or just a number [SlL2-D5]
More informationMathematics High School Geometry
Mathematics High School Geometry An understanding of the attributes and relationships of geometric objects can be applied in diverse contexts interpreting a schematic drawing, estimating the amount of
More informationTools, Tips and Workflows Thinning Feature Vertices in LP360 LP360,
LP360, 2017.1 l L. Graham 7 January 2017 We have lot of LP360 customers with a variety of disparate feature collection and editing needs. To begin to address these needs, we introduced a completely new
More informationThis blog addresses the question: how do we determine the intersection of two circles in the Cartesian plane?
Intersecting Circles This blog addresses the question: how do we determine the intersection of two circles in the Cartesian plane? This is a problem that a programmer might have to solve, for example,
More informationMathematics GCSE 9 1 Higher Syllabus. Yes. Does the subject set according to ability? Skills Covered. Unit
Mathematics GCSE 9 1 Higher Syllabus Does the subject set according to ability? Unit Unit 1 Unit 2 Unit 3 Unit 4 Yes Skills Covered understand and apply place value correctly. estimate values including
More informationGraphics and Interaction Transformation geometry and homogeneous coordinates
433-324 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationPITSCO Math Individualized Prescriptive Lessons (IPLs)
Orientation Integers 10-10 Orientation I 20-10 Speaking Math Define common math vocabulary. Explore the four basic operations and their solutions. Form equations and expressions. 20-20 Place Value Define
More informationCOMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates
COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates Department of Computer Science and Software Engineering The Lecture outline Introduction Vectors and matrices Translation
More informationAUTOMATIC PHOTO ORIENTATION VIA MATCHING WITH CONTROL PATCHES
AUTOMATIC PHOTO ORIENTATION VIA MATCHING WITH CONTROL PATCHES J. J. Jaw a *, Y. S. Wu b Dept. of Civil Engineering, National Taiwan University, Taipei,10617, Taiwan, ROC a jejaw@ce.ntu.edu.tw b r90521128@ms90.ntu.edu.tw
More informationShould Contours Be Generated from Lidar Data, and Are Breaklines Required? Lidar data provides the most
Should Contours Be Generated from Lidar Data, and Are Breaklines Required? Lidar data provides the most accurate and reliable representation of the topography of the earth. As lidar technology advances
More informationREGISTRATION OF AIRBORNE LASER DATA TO SURFACES GENERATED BY PHOTOGRAMMETRIC MEANS. Y. Postolov, A. Krupnik, K. McIntosh
REGISTRATION OF AIRBORNE LASER DATA TO SURFACES GENERATED BY PHOTOGRAMMETRIC MEANS Y. Postolov, A. Krupnik, K. McIntosh Department of Civil Engineering, Technion Israel Institute of Technology, Haifa,
More informationUNIT 1 GEOMETRY TEMPLATE CREATED BY REGION 1 ESA UNIT 1
UNIT 1 GEOMETRY TEMPLATE CREATED BY REGION 1 ESA UNIT 1 Traditional Pathway: Geometry The fundamental purpose of the course in Geometry is to formalize and extend students geometric experiences from the
More informationGeometric Rectification of Remote Sensing Images
Geometric Rectification of Remote Sensing Images Airborne TerrestriaL Applications Sensor (ATLAS) Nine flight paths were recorded over the city of Providence. 1 True color ATLAS image (bands 4, 2, 1 in
More informationGeometric Constructions
HISTORY OF MATHEMATICS Spring 2005 Geometric Constructions Notes, activities, assignment; #3 in a series. Note: I m not giving a specific due date for this somewhat vague assignment. The idea is that it
More informationTHE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS. Alekseí Yu. Chekunov. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69, 1 (2017), 12 22 March 2017 research paper originalni nauqni rad THE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS Alekseí Yu. Chekunov Abstract. In this
More informationGeometry I Can Statements I can describe the undefined terms: point, line, and distance along a line in a plane I can describe the undefined terms:
Geometry I Can Statements I can describe the undefined terms: point, line, and distance along a line in a plane I can describe the undefined terms: point, line, and distance along a line in a plane I can
More informationIntensity Augmented ICP for Registration of Laser Scanner Point Clouds
Intensity Augmented ICP for Registration of Laser Scanner Point Clouds Bharat Lohani* and Sandeep Sashidharan *Department of Civil Engineering, IIT Kanpur Email: blohani@iitk.ac.in. Abstract While using
More informationChapters 1 9: Overview
Chapters 1 9: Overview Chapter 1: Introduction Chapters 2 4: Data acquisition Chapters 5 9: Data manipulation Chapter 5: Vertical imagery Chapter 6: Image coordinate measurements and refinements Chapters
More informationBasics of Computational Geometry
Basics of Computational Geometry Nadeem Mohsin October 12, 2013 1 Contents This handout covers the basic concepts of computational geometry. Rather than exhaustively covering all the algorithms, it deals
More informationLight: Geometric Optics
Light: Geometric Optics Regular and Diffuse Reflection Sections 23-1 to 23-2. How We See Weseebecauselightreachesoureyes. There are two ways, therefore, in which we see: (1) light from a luminous object
More informationStructural and Syntactic Pattern Recognition
Structural and Syntactic Pattern Recognition Selim Aksoy Department of Computer Engineering Bilkent University saksoy@cs.bilkent.edu.tr CS 551, Fall 2017 CS 551, Fall 2017 c 2017, Selim Aksoy (Bilkent
More informationifp Universität Stuttgart Performance of IGI AEROcontrol-IId GPS/Inertial System Final Report
Universität Stuttgart Performance of IGI AEROcontrol-IId GPS/Inertial System Final Report Institute for Photogrammetry (ifp) University of Stuttgart ifp Geschwister-Scholl-Str. 24 D M. Cramer: Final report
More information1. Introduction. A CASE STUDY Dense Image Matching Using Oblique Imagery Towards All-in- One Photogrammetry
Submitted to GIM International FEATURE A CASE STUDY Dense Image Matching Using Oblique Imagery Towards All-in- One Photogrammetry Dieter Fritsch 1, Jens Kremer 2, Albrecht Grimm 2, Mathias Rothermel 1
More informationThe Geometry of Carpentry and Joinery
The Geometry of Carpentry and Joinery Pat Morin and Jason Morrison School of Computer Science, Carleton University, 115 Colonel By Drive Ottawa, Ontario, CANADA K1S 5B6 Abstract In this paper we propose
More informationCOLOR FIDELITY OF CHROMATIC DISTRIBUTIONS BY TRIAD ILLUMINANT COMPARISON. Marcel P. Lucassen, Theo Gevers, Arjan Gijsenij
COLOR FIDELITY OF CHROMATIC DISTRIBUTIONS BY TRIAD ILLUMINANT COMPARISON Marcel P. Lucassen, Theo Gevers, Arjan Gijsenij Intelligent Systems Lab Amsterdam, University of Amsterdam ABSTRACT Performance
More informationA New Way to Control Mobile LiDAR Data
A New Way to Control Mobile LiDAR Data Survey control has always been a critically important issue when conducting mobile LiDAR surveys. While the accuracies currently being achieved with the most capable
More informationStereo Vision. MAN-522 Computer Vision
Stereo Vision MAN-522 Computer Vision What is the goal of stereo vision? The recovery of the 3D structure of a scene using two or more images of the 3D scene, each acquired from a different viewpoint in
More informationMathematics Curriculum
6 G R A D E Mathematics Curriculum GRADE 6 5 Table of Contents 1... 1 Topic A: Area of Triangles, Quadrilaterals, and Polygons (6.G.A.1)... 11 Lesson 1: The Area of Parallelograms Through Rectangle Facts...
More informationGeoCue Group, Inc Madison Blvd. Suite 202 Madison, AL DEM CuePac User Guide Version 2017.
GeoCue Group, Inc. 9668 Madison Blvd. Suite 202 Madison, AL 35758 1-256-461-8289 www.geocue.com DEM CuePac User Guide Version 2017.1 5 May 2017 NOTICES The material in GeoCue Group, Inc. documents is protected
More informationCommon Core Cluster. Experiment with transformations in the plane. Unpacking What does this standard mean that a student will know and be able to do?
Congruence G.CO Experiment with transformations in the plane. G.CO.1 Know precise definitions of angle, circle, perpendicular line, parallel line, and line segment, based on the undefined notions of point,
More informationDigital Image Processing ERRATA. Wilhelm Burger Mark J. Burge. An algorithmic introduction using Java. Second Edition. Springer
Wilhelm Burger Mark J. Burge Digital Image Processing An algorithmic introduction using Java Second Edition ERRATA Springer Berlin Heidelberg NewYork Hong Kong London Milano Paris Tokyo 12.1 RGB Color
More informationSubstituting a 2 b 2 for c 2 and using a little algebra, we can then derive the standard equation for an ellipse centred at the origin,
Conics onic sections are the curves which result from the intersection of a plane with a cone. These curves were studied and revered by the ancient Greeks, and were written about extensively by both Euclid
More informationChicago is Rural: The Inconsistencies and Absurdities of Street Connectivity Indices Paul Knight Georgia Tech MCRP Option Paper March 1, 2011
1 Chicago is Rural: The Inconsistencies and Absurdities of Street Connectivity Indices Paul Knight Georgia Tech MCRP Option Paper March 1, 2011 Introduction Many contemporary urban design metrics have
More information2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into
2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into the viewport of the current application window. A pixel
More informationFigure 1: Methods for counting nodes and links for use in the calculation of a connectivity index.
Chicago is Rural: The Inconsistencies and Absurdities of Street Connectivity Indices Paul Knight Georgia Institute of Technology CNU 19 Academic Paper Submission December 6, 2010 Abstract Various connectivity
More information5. GENERALIZED INVERSE SOLUTIONS
5. GENERALIZED INVERSE SOLUTIONS The Geometry of Generalized Inverse Solutions The generalized inverse solution to the control allocation problem involves constructing a matrix which satisfies the equation
More informationAalborg Universitet. Published in: Accuracy Publication date: Document Version Early version, also known as pre-print
Aalborg Universitet A method for checking the planimetric accuracy of Digital Elevation Models derived by Airborne Laser Scanning Høhle, Joachim; Øster Pedersen, Christian Published in: Accuracy 2010 Publication
More informationUSING THE ORTHOGONAL PROJECTION FOR PARAMETER INITIALIZATION IN THE 3D RECONSTRUCTION OF DISTANT OBJECTS INTRODUCTION
USING THE ORTHOGONAL PROJECTION FOR PARAMETER INITIALIZATION IN THE 3D RECONSTRUCTION OF DISTANT OBJECTS Keith F. Blonquist, Junior Project Engineer Lidar Pacific Corporation Logan, UT 84341 kfblonquist@lidarpacific.com
More informationTools, Tips, and Workflows LP360 Under the Hood: Cross-sections LP360
LP360 Lewis Graham Revision 1.0 l As we outlined in a previous LP360 Under the Hood, we have added volumetric analysis to LP360. This feature was first released with 2013.2 at the end of last year with
More informationConic Sections. College Algebra
Conic Sections College Algebra Conic Sections A conic section, or conic, is a shape resulting from intersecting a right circular cone with a plane. The angle at which the plane intersects the cone determines
More informationNAEP Released Items Aligned to the Iowa Core: Geometry
NAEP Released Items Aligned to the Iowa Core: Geometry Congruence G-CO Experiment with transformations in the plane 1. Know precise definitions of angle, circle, perpendicular line, parallel line, and
More informationExam 2 Review. 2. What the difference is between an equation and an expression?
Exam 2 Review Chapter 1 Section1 Do You Know: 1. What does it mean to solve an equation? 2. What the difference is between an equation and an expression? 3. How to tell if an equation is linear? 4. How
More informationMath Lab- Geometry Pacing Guide Quarter 3. Unit 1: Rational and Irrational Numbers, Exponents and Roots
1 Jan. 3-6 (4 days) 2 Jan. 9-13 Unit 1: Rational and Irrational Numbers, Exponents and Roots ISTEP+ ISTEP Framework Focus: Unit 1 Number Sense, Expressions, and Computation 8.NS.1: Give examples of rational
More information6 Mathematics Curriculum
New York State Common Core 6 Mathematics Curriculum GRADE GRADE 6 MODULE 5 Table of Contents 1 Area, Surface Area, and Volume Problems... 3 Topic A: Area of Triangles, Quadrilaterals, and Polygons (6.G.A.1)...
More informationMastery. PRECALCULUS Student Learning Targets
PRECALCULUS Student Learning Targets Big Idea: Sequences and Series 1. I can describe a sequence as a function where the domain is the set of natural numbers. Connections (Pictures, Vocabulary, Definitions,
More informationCOMPUTER SIMULATION TECHNIQUES FOR ACOUSTICAL DESIGN OF ROOMS - HOW TO TREAT REFLECTIONS IN SOUND FIELD SIMULATION
J.H. Rindel, Computer simulation techniques for the acoustical design of rooms - how to treat reflections in sound field simulation. ASVA 97, Tokyo, 2-4 April 1997. Proceedings p. 201-208. COMPUTER SIMULATION
More informationT:\Staff Resources\Maths\1 Dept Admin\School website\key Stage 4\GCSE Higher Topic Detail doc 1
Module (time) Contents Notes Number (8) Recognise triangular, square and cube numbers Recall integer squares and corresponding square roots to 15 x 15 Recall the cubes of 2,3,4,5 and 10 Find multiples,
More informationCOMPUTING APPARENT CURVATURE PROGRAM euler_curvature. Euler curvature computation flow chart. Computing Euler curvature
COMPUTING APPARENT CURVATURE PROGRAM euler_curvature Euler curvature computation flow chart Inline dip Crossline dip curvature3d k 1 Strike of k 1 k 2 Strike of k 1 euler_curvature (φ=-90 0 ) (φ=0 0 )
More informationAn Underdetermined Linear System for GPS
An Underdetermined Linear System for GPS Dan Kalman Dan Kalman (kalman@american.edu) joined the mathematics faculty at American University in 1993, following an eight year stint in the aerospace industry
More informationIMAGINE Objective. The Future of Feature Extraction, Update & Change Mapping
IMAGINE ive The Future of Feature Extraction, Update & Change Mapping IMAGINE ive provides object based multi-scale image classification and feature extraction capabilities to reliably build and maintain
More information2.3 Circular Functions of Real Numbers
www.ck12.org Chapter 2. Graphing Trigonometric Functions 2.3 Circular Functions of Real Numbers Learning Objectives Graph the six trigonometric ratios as functions on the Cartesian plane. Identify the
More information8.B. The result of Regiomontanus on tetrahedra
8.B. The result of Regiomontanus on tetrahedra We have already mentioned that Plato s theory that the five regular polyhedra represent the fundamental elements of nature, and in supplement (3.D) to the
More information16.6. Parametric Surfaces. Parametric Surfaces. Parametric Surfaces. Vector Calculus. Parametric Surfaces and Their Areas
16 Vector Calculus 16.6 and Their Areas Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. and Their Areas Here we use vector functions to describe more general
More informationCW High School. Advanced Math A. 1.1 I can make connections between the algebraic equation or description for a function, its name, and its graph.
1. Functions and Math Models (10.00%) 1.1 I can make connections between the algebraic equation or description for a function, its name, and its graph. 4 Pro cient I can make connections between the algebraic
More informationROAD SURFACE STRUCTURE MONITORING AND ANALYSIS USING HIGH PRECISION GPS MOBILE MEASUREMENT SYSTEMS (MMS)
ROAD SURFACE STRUCTURE MONITORING AND ANALYSIS USING HIGH PRECISION GPS MOBILE MEASUREMENT SYSTEMS (MMS) Bonifacio R. Prieto PASCO Philippines Corporation, Pasig City, 1605, Philippines Email: bonifacio_prieto@pascoph.com
More informationCS6100: Topics in Design and Analysis of Algorithms
CS6100: Topics in Design and Analysis of Algorithms Guarding and Triangulating Polygons John Augustine CS6100 (Even 2012): Guarding and Triangulating Polygons The Art Gallery Problem A simple polygon is
More informationMADISON ACADEMY GEOMETRY PACING GUIDE
MADISON ACADEMY GEOMETRY PACING GUIDE 2018-2019 Standards (ACT included) ALCOS#1 Know the precise definitions of angle, circle, perpendicular line, parallel line, and line segment based on the undefined
More informationSection 1: Introduction to Geometry Points, Lines, and Planes
Section 1: Introduction to Geometry Points, Lines, and Planes Topic 1: Basics of Geometry - Part 1... 3 Topic 2: Basics of Geometry Part 2... 5 Topic 3: Midpoint and Distance in the Coordinate Plane Part
More informationINDEPENDENT SCHOOL DISTRICT 196 Rosemount, Minnesota Educating our students to reach their full potential
INDEPENDENT SCHOOL DISTRICT 196 Rosemount, Minnesota Educating our students to reach their full potential MINNESOTA MATHEMATICS STANDARDS Grades 9, 10, 11 I. MATHEMATICAL REASONING Apply skills of mathematical
More informationA METHOD TO PREDICT ACCURACY OF LEAST SQUARES SURFACE MATCHING FOR AIRBORNE LASER SCANNING DATA SETS
A METHOD TO PREDICT ACCURACY OF LEAST SQUARES SURFACE MATCHING FOR AIRBORNE LASER SCANNING DATA SETS Robert Pâquet School of Engineering, University of Newcastle Callaghan, NSW 238, Australia (rpaquet@mail.newcastle.edu.au)
More informationAlgorithm of correction of error caused by perspective distortions of measuring mark images
Mechanics & Industry 7, 73 (206) c AFM, EDP Sciences 206 DOI: 0.05/meca/206077 www.mechanics-industry.org Mechanics & Industry Algorithm of correction of error caused by perspective distortions of measuring
More informationPre AP Geometry. Mathematics Standards of Learning Curriculum Framework 2009: Pre AP Geometry
Pre AP Geometry Mathematics Standards of Learning Curriculum Framework 2009: Pre AP Geometry 1 The content of the mathematics standards is intended to support the following five goals for students: becoming
More informationGEOMETRY. Changes to the original 2010 COS is in red. If it is red and crossed out, it has been moved to another course.
The Geometry course builds on Algebra I concepts and increases students knowledge of shapes and their properties through geometry-based applications, many of which are observable in aspects of everyday
More informationGeometry GEOMETRY. Congruence
Geometry Geometry builds on Algebra I concepts and increases students knowledge of shapes and their properties through geometry-based applications, many of which are observable in aspects of everyday life.
More informationChapter 15 Introduction to Linear Programming
Chapter 15 Introduction to Linear Programming An Introduction to Optimization Spring, 2015 Wei-Ta Chu 1 Brief History of Linear Programming The goal of linear programming is to determine the values of
More informationChapter 5. Transforming Shapes
Chapter 5 Transforming Shapes It is difficult to walk through daily life without being able to see geometric transformations in your surroundings. Notice how the leaves of plants, for example, are almost
More informationLecture 3: Some Strange Properties of Fractal Curves
Lecture 3: Some Strange Properties of Fractal Curves I have been a stranger in a strange land. Exodus 2:22 1. Fractal Strangeness Fractals have a look and feel that is very different from ordinary curves.
More informationIGTF 2016 Fort Worth, TX, April 11-15, 2016 Submission 149
IGTF 26 Fort Worth, TX, April -5, 26 2 3 4 5 6 7 8 9 2 3 4 5 6 7 8 9 2 2 Light weighted and Portable LiDAR, VLP-6 Registration Yushin Ahn (yahn@mtu.edu), Kyung In Huh (khuh@cpp.edu), Sudhagar Nagarajan
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationTRAINING MATERIAL HOW TO OPTIMIZE ACCURACY WITH CORRELATOR3D
TRAINING MATERIAL WITH CORRELATOR3D Page2 Contents 1. UNDERSTANDING INPUT DATA REQUIREMENTS... 4 1.1 What is Aerial Triangulation?... 4 1.2 Recommended Flight Configuration... 4 1.3 Data Requirements for
More informationv Overview SMS Tutorials Prerequisites Requirements Time Objectives
v. 12.2 SMS 12.2 Tutorial Overview Objectives This tutorial describes the major components of the SMS interface and gives a brief introduction to the different SMS modules. Ideally, this tutorial should
More informationNearest approaches to multiple lines in n-dimensional space
Nearest approaches Nearest approaches to multiple lines in n-dimensional space Lejia Han and John C. Bancroft ABSTRACT An efficient matrix approach is described to obtain the nearest approach to nonintersecting
More informationCamera model and multiple view geometry
Chapter Camera model and multiple view geometry Before discussing how D information can be obtained from images it is important to know how images are formed First the camera model is introduced and then
More informationUnit 1. Word Definition Picture. The number s distance from 0 on the number line. The symbol that means a number is greater than the second number.
Unit 1 Word Definition Picture Absolute Value The number s distance from 0 on the number line. -3 =3 Greater Than The symbol that means a number is greater than the second number. > Greatest to Least To
More informationMathematics Standards for High School Geometry
Mathematics Standards for High School Geometry Geometry is a course required for graduation and course is aligned with the College and Career Ready Standards for Mathematics in High School. Throughout
More information8 th Grade Pre Algebra Pacing Guide 1 st Nine Weeks
8 th Grade Pre Algebra Pacing Guide 1 st Nine Weeks MS Objective CCSS Standard I Can Statements Included in MS Framework + Included in Phase 1 infusion Included in Phase 2 infusion 1a. Define, classify,
More information, Geometry, Quarter 1
2017.18, Geometry, Quarter 1 The following Practice Standards and Literacy Skills will be used throughout the course: Standards for Mathematical Practice Literacy Skills for Mathematical Proficiency 1.
More informationcorrelated to the Utah 2007 Secondary Math Core Curriculum Geometry
correlated to the Utah 2007 Secondary Math Core Curriculum Geometry McDougal Littell Geometry: Concepts and Skills 2005 correlated to the Utah 2007 Secondary Math Core Curriculum Geometry The main goal
More informationChapters 1 7: Overview
Chapters 1 7: Overview Chapter 1: Introduction Chapters 2 4: Data acquisition Chapters 5 7: Data manipulation Chapter 5: Vertical imagery Chapter 6: Image coordinate measurements and refinements Chapter
More informationNEW MONITORING TECHNIQUES ON THE DETERMINATION OF STRUCTURE DEFORMATIONS
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 003. NEW MONITORING TECHNIQUES ON THE DETERMINATION OF STRUCTURE DEFORMATIONS D.Stathas, O.Arabatzi, S.Dogouris, G.Piniotis,
More informationCourse of study- Algebra Introduction: Algebra 1-2 is a course offered in the Mathematics Department. The course will be primarily taken by
Course of study- Algebra 1-2 1. Introduction: Algebra 1-2 is a course offered in the Mathematics Department. The course will be primarily taken by students in Grades 9 and 10, but since all students must
More information