Point Cloud Processing

Size: px
Start display at page:

Download "Point Cloud Processing"

Transcription

1 Point Cloud Processing Has anyone seen the toothpaste? Given a point cloud: how do you detect and localize objects? how do you map terrain?

2 What is a point cloud? Point cloud: a set of points in 3-D space just a set of 3-d points Mesh: each point is a vertex of a triangulated face a set of vertices AND connectivity information Point cloud Mesh Images from Course INF 555 slides, Ecole Polytechnique, Paris

3 What is a point cloud? Point cloud: a set of points in 3-D space just a set of 3-d points Mesh: each point is a vertex of a triangulated face a set of vertices AND connectivity information Many depth sensors produce point clouds natively Point cloud But, a mesh contains a lot more information Mesh Images from Course INF 555 slides, Ecole Polytechnique, Paris

4 Time of flight sensors Hokuyo UTM-30LX-EW Scanning Laser Range Finder

5 Time of flight sensors Slide from Course INF 555 slides, Ecole Polytechnique, Paris

6 Time of flight sensors

7 Structured light sensors

8 Slide: John MacCormick, Dickinson University

9 Slide: John MacCormick, Dickinson University

10 Slide: John MacCormick, Dickinson University

11 Slide: John MacCormick, Dickinson University

12 Slide: John MacCormick, Dickinson University

13 Slide: John MacCormick, Dickinson University

14 Calculating surface normals Point cloud Point cloud w/ surface normals (normals are subsampled)

15 Calculating surface normals Question: How do we calculate the surface normal given only points? Answer: 1. Calculate the sample covariance matrix of the points within a local neighborhood of the surface normal 2. Take Eigenvalues/Eigenvectors of that covariance matrix

16 Calculating surface normals Let C denote the set of points in the point cloud Suppose we want to calculate the surface normal at Let denote the r-ball about x is the set of points in the cloud within r of x

17 Calculating surface normals Calculate the sample covariance matrix of the points in

18 Calculating surface normals Length = Eigenvalues of Length = Principle axes of ellipse point in directions of corresponding eigenvectors

19 Calculating surface normals So: surface normal is in the direction of the Eigenvector corresponding to the smallest Eigenvalue of

20 Calculating surface normals: Summary 1. calculate points within r-ball about x: 2. calculate covariance matrix: 3. calculate Eigenvectors: and Eigenvalues: (\lambda_3 is smallest) 4. v_3 is parallel or antiparallel to surface normal

21 Calculating surface normals Important note: the points alone do not tell us the sign of the surface normal

22 Calculating surface normals Important note: the points alone do not tell us the sign of the surface normal

23 Calculating surface normals How large a point neighborhood to use when calculating? Because points can be uneven, don't use k-nearest neighbor. it's important to select a radius r and stick w/ it. which what value of r to use?

24 Calculating surface normals Images from Course INF 555 slides, Ecole Polytechnique, Paris

25 Calculating surface normals Images from Course INF 555 slides, Ecole Polytechnique, Paris

26 Outlier removal Similar approach as in normal estimation: 1. calculate local covariance matrix 2. estimate Eigenvectors/Eigenvalues 3. use that information somehow... Images from Course INF 555 slides, Ecole Polytechnique, Paris

27 Outlier removal If points lie on a line, then is small If points are uniformly random, then Outlier removal: delete all points for which is close to 1 is above a threshold Images from Course INF 555 slides, Ecole Polytechnique, Paris

28 Point cloud registration: ICP Find an affine transformation that aligns two partially overlapping point clouds Images from Course INF 555 slides, Ecole Polytechnique, Paris

29 ICP Problem Statement This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

30 ICP: key idea This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

31 Step 1: center the two point clouds This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

32 Step 2: use SVD to get min t and R This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

33 Step 2: use SVD to get min t and R This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

34 ICP data association problem This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

35 ICP Algorithm Input: two point sets, X and P Output: translation t and rotation R that best aligns pt sets 1. Start with a good alignment 2. Repeat until t and R are small: 3. for every point in X, find its closest neighbor in P 4. find min t and R for that correspondence assignment 5. translate and rotate P by t and R 6. Figure out net translation and rotation, t and R Converges if the point sets are initially well aligned Besl and McKay, 1992

36 ICP example This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

37 ICP Variants This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

38 Selecting points to align This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

39 Normal-space sampling This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

40 Comparison: normal space sampling vs random This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

41 Feature based sampling This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

42 ICP: data association This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

43 ICP: data association This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

44 Closest point matching This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

45 Normal shooting This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

46 Point-to-plane distances This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

47 Closest compatible point This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

48 ICP: summary This slide from: Burgard, Stachniss, Bennewitz, Arras, U. Freiburg

49 Another approach to alignment: RANSAC This slide from: Kavita Bala, Cornell U.

50 RANSAC This slide from: Kavita Bala, Cornell U.

51 How does regression work here?

52 Image alignment problem This slide from: Kavita Bala, Cornell U.

53 Outliers This slide from: Kavita Bala, Cornell U.

54 RANSAC This slide from: Kavita Bala, Cornell U.

55 RANSAC key idea This slide from: Kavita Bala, Cornell U.

56 Counting inliers This slide from: Kavita Bala, Cornell U.

57 Counting inliers This slide from: Kavita Bala, Cornell U.

58 How do we find the best line? This slide from: Kavita Bala, Cornell U.

59 RANSAC This slide from: Kavita Bala, Cornell U.

60 This slide from: Kavita Bala, Cornell U.

61 This slide from: Kavita Bala, Cornell U.

62 This slide from: Kavita Bala, Cornell U.

63 Using RANSAC to Fit a Sphere

64 Using RANSAC to Fit a Sphere

65 Using RANSAC to Fit a Sphere Radius? Center?

66 Using RANSAC to Fit a Sphere How generate candidate spheres? How score spheres?

67 Using RANSAC to Fit a Sphere How generate candidate spheres? 1. sample a point How score spheres?

68 Using RANSAC to Fit a Sphere How generate candidate spheres? 1. sample a point 2. estimate surface normal How score spheres?

69 Using RANSAC to Fit a Sphere radius How generate candidate spheres? 1. sample a point 2. estimate surface normal 3. sample radius How score spheres?

70 Using RANSAC to Fit a Sphere radius How generate candidate spheres? 1. sample a point 2. estimate surface normal 3. sample radius 4. estimate center to be radius distance from sampled point along surface normal How score spheres?

71 Using RANSAC to Fit a Sphere radius How generate candidate spheres? 1. sample a point 2. estimate surface normal 3. sample radius 4. estimate center to be radius distance from sampled point along surface normal How score spheres? 1. count num pts within epsilon of candidate sphere surface

72 Using RANSAC to Fit a Cylinder How generate candidate cylinders?

73 Using RANSAC to Fit a Cylinder How generate candidate cylinders? 1. sample two pts

74 Using RANSAC to Fit a Cylinder How generate candidate cylinders? 1. sample two pts 2. estimate surface normal at both pts

75 Using RANSAC to Fit a Cylinder How generate candidate cylinders? 1. sample two pts 2. estimate surface normal at both pts 3. get sample axis by taking cross product between two normals

76 Using RANSAC to Fit a Cylinder How generate candidate cylinders? 1. sample two pts 2. estimate surface normal at both pts 3. get sample axis by taking cross product between two normals 4. project points onto plane orthogonal to axis 5. fit a circle using a method similar to what we did for the sphere.

77 Using RANSAC to Fit a Cylinder 3xn matrix of pts projected onto plane How generate candidate cylinders? 1. sample two pts 2. estimate surface normal at both pts 3. get sample axis by taking cross product between two normals 4. project points onto plane orthogonal to axis 5. fit a circle using a method similar to what we did for the sphere. 3x1 unit vector in direction of axis 3xn matrix of pts in 3d space

78 RANSAC: the parameters This slide from: Kavita Bala, Cornell U.

79 RANSAC: how many rounds? This slide from: Kavita Bala, Cornell U.

80 RANSAC: how many parameters to sample? This slide from: Kavita Bala, Cornell U.

81 RANSAC Summary This slide from: Kavita Bala, Cornell U.

82 Hough transform This slide from: Kavita Bala, Cornell U.

83 Hough transform This slide from: Kavita Bala, Cornell U.

84 Hough transform This slide from: Kavita Bala, Cornell U.

85 Hough transform This slide from: Kavita Bala, Cornell U.

86 Hough transorm This slide from: Kavita Bala, Cornell U.

Introduction to Mobile Robotics Iterative Closest Point Algorithm. Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras

Introduction to Mobile Robotics Iterative Closest Point Algorithm. Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras Introduction to Mobile Robotics Iterative Closest Point Algorithm Wolfram Burgard, Cyrill Stachniss, Maren Bennewitz, Kai Arras 1 Motivation 2 The Problem Given: two corresponding point sets: Wanted: translation

More information

INFO0948 Fitting and Shape Matching

INFO0948 Fitting and Shape Matching INFO0948 Fitting and Shape Matching Renaud Detry University of Liège, Belgium Updated March 31, 2015 1 / 33 These slides are based on the following book: D. Forsyth and J. Ponce. Computer vision: a modern

More information

Surface Registration. Gianpaolo Palma

Surface Registration. Gianpaolo Palma Surface Registration Gianpaolo Palma The problem 3D scanning generates multiple range images Each contain 3D points for different parts of the model in the local coordinates of the scanner Find a rigid

More information

3D Perception. CS 4495 Computer Vision K. Hawkins. CS 4495 Computer Vision. 3D Perception. Kelsey Hawkins Robotics

3D Perception. CS 4495 Computer Vision K. Hawkins. CS 4495 Computer Vision. 3D Perception. Kelsey Hawkins Robotics CS 4495 Computer Vision Kelsey Hawkins Robotics Motivation What do animals, people, and robots want to do with vision? Detect and recognize objects/landmarks Find location of objects with respect to themselves

More information

Scan Matching. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics

Scan Matching. Pieter Abbeel UC Berkeley EECS. Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Scan Matching Pieter Abbeel UC Berkeley EECS Many slides adapted from Thrun, Burgard and Fox, Probabilistic Robotics Scan Matching Overview Problem statement: Given a scan and a map, or a scan and a scan,

More information

5.2 Surface Registration

5.2 Surface Registration Spring 2018 CSCI 621: Digital Geometry Processing 5.2 Surface Registration Hao Li http://cs621.hao-li.com 1 Acknowledgement Images and Slides are courtesy of Prof. Szymon Rusinkiewicz, Princeton University

More information

Algorithm research of 3D point cloud registration based on iterative closest point 1

Algorithm research of 3D point cloud registration based on iterative closest point 1 Acta Technica 62, No. 3B/2017, 189 196 c 2017 Institute of Thermomechanics CAS, v.v.i. Algorithm research of 3D point cloud registration based on iterative closest point 1 Qian Gao 2, Yujian Wang 2,3,

More information

Using Augmented Measurements to Improve the Convergence of ICP. Jacopo Serafin and Giorgio Grisetti

Using Augmented Measurements to Improve the Convergence of ICP. Jacopo Serafin and Giorgio Grisetti Jacopo Serafin and Giorgio Grisetti Point Cloud Registration We want to find the rotation and the translation that maximize the overlap between two point clouds Page 2 Point Cloud Registration We want

More information

Image Stitching. Slides from Rick Szeliski, Steve Seitz, Derek Hoiem, Ira Kemelmacher, Ali Farhadi

Image Stitching. Slides from Rick Szeliski, Steve Seitz, Derek Hoiem, Ira Kemelmacher, Ali Farhadi Image Stitching Slides from Rick Szeliski, Steve Seitz, Derek Hoiem, Ira Kemelmacher, Ali Farhadi Combine two or more overlapping images to make one larger image Add example Slide credit: Vaibhav Vaish

More information

Chapter 4. Clustering Core Atoms by Location

Chapter 4. Clustering Core Atoms by Location Chapter 4. Clustering Core Atoms by Location In this chapter, a process for sampling core atoms in space is developed, so that the analytic techniques in section 3C can be applied to local collections

More information

Structured Light II. Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov

Structured Light II. Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov Structured Light II Johannes Köhler Johannes.koehler@dfki.de Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov Introduction Previous lecture: Structured Light I Active Scanning Camera/emitter

More information

Structured Light II. Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov

Structured Light II. Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov Structured Light II Johannes Köhler Johannes.koehler@dfki.de Thanks to Ronen Gvili, Szymon Rusinkiewicz and Maks Ovsjanikov Introduction Previous lecture: Structured Light I Active Scanning Camera/emitter

More information

Photo by Carl Warner

Photo by Carl Warner Photo b Carl Warner Photo b Carl Warner Photo b Carl Warner Fitting and Alignment Szeliski 6. Computer Vision CS 43, Brown James Has Acknowledgment: Man slides from Derek Hoiem and Grauman&Leibe 2008 AAAI

More information

Instance-level recognition

Instance-level recognition Instance-level recognition 1) Local invariant features 2) Matching and recognition with local features 3) Efficient visual search 4) Very large scale indexing Matching of descriptors Matching and 3D reconstruction

More information

Field-of-view dependent registration of point clouds and incremental segmentation of table-tops using time-offlight

Field-of-view dependent registration of point clouds and incremental segmentation of table-tops using time-offlight Field-of-view dependent registration of point clouds and incremental segmentation of table-tops using time-offlight cameras Dipl.-Ing. Georg Arbeiter Fraunhofer Institute for Manufacturing Engineering

More information

3D Computer Vision. Structured Light II. Prof. Didier Stricker. Kaiserlautern University.

3D Computer Vision. Structured Light II. Prof. Didier Stricker. Kaiserlautern University. 3D Computer Vision Structured Light II Prof. Didier Stricker Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz http://av.dfki.de 1 Introduction

More information

3D Modeling from Range Images

3D Modeling from Range Images 1 3D Modeling from Range Images A Comprehensive System for 3D Modeling from Range Images Acquired from a 3D ToF Sensor Dipl.-Inf. March 22th, 2007 Sensor and Motivation 2 3D sensor PMD 1k-S time-of-flight

More information

3D Environment Reconstruction

3D Environment Reconstruction 3D Environment Reconstruction Using Modified Color ICP Algorithm by Fusion of a Camera and a 3D Laser Range Finder The 2009 IEEE/RSJ International Conference on Intelligent Robots and Systems October 11-15,

More information

LOAM: LiDAR Odometry and Mapping in Real Time

LOAM: LiDAR Odometry and Mapping in Real Time LOAM: LiDAR Odometry and Mapping in Real Time Aayush Dwivedi (14006), Akshay Sharma (14062), Mandeep Singh (14363) Indian Institute of Technology Kanpur 1 Abstract This project deals with online simultaneous

More information

Intensity Augmented ICP for Registration of Laser Scanner Point Clouds

Intensity 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 information

Gaussian and Mean Curvature Planar points: Zero Gaussian curvature and zero mean curvature Tangent plane intersects surface at infinity points Gauss C

Gaussian and Mean Curvature Planar points: Zero Gaussian curvature and zero mean curvature Tangent plane intersects surface at infinity points Gauss C Outline Shape Analysis Basics COS 526, Fall 21 Curvature PCA Distance Features Some slides from Rusinkiewicz Curvature Curvature Curvature κof a curve is reciprocal of radius of circle that best approximates

More information

Instance-level recognition

Instance-level recognition Instance-level recognition 1) Local invariant features 2) Matching and recognition with local features 3) Efficient visual search 4) Very large scale indexing Matching of descriptors Matching and 3D reconstruction

More information

Automatic Image Alignment (feature-based)

Automatic Image Alignment (feature-based) Automatic Image Alignment (feature-based) Mike Nese with a lot of slides stolen from Steve Seitz and Rick Szeliski 15-463: Computational Photography Alexei Efros, CMU, Fall 2006 Today s lecture Feature

More information

#$ % $ $& "$%% " $ '$ " '

#$ % $ $& $%%  $ '$  ' ! " This section of the course covers techniques for pairwise (i.e., scanto-scan) and global (i.e., involving more than 2 scans) alignment, given that the algorithms are constrained to obtain a rigid-body

More information

Iterative Closest Point Algorithm in the Presence of Anisotropic Noise

Iterative Closest Point Algorithm in the Presence of Anisotropic Noise Iterative Closest Point Algorithm in the Presence of Anisotropic Noise L. Maier-Hein, T. R. dos Santos, A. M. Franz, H.-P. Meinzer German Cancer Research Center, Div. of Medical and Biological Informatics

More information

Efficient SLAM Scheme Based ICP Matching Algorithm Using Image and Laser Scan Information

Efficient SLAM Scheme Based ICP Matching Algorithm Using Image and Laser Scan Information Proceedings of the World Congress on Electrical Engineering and Computer Systems and Science (EECSS 2015) Barcelona, Spain July 13-14, 2015 Paper No. 335 Efficient SLAM Scheme Based ICP Matching Algorithm

More information

Scanning Real World Objects without Worries 3D Reconstruction

Scanning Real World Objects without Worries 3D Reconstruction Scanning Real World Objects without Worries 3D Reconstruction 1. Overview Feng Li 308262 Kuan Tian 308263 This document is written for the 3D reconstruction part in the course Scanning real world objects

More information

CSE 252B: Computer Vision II

CSE 252B: Computer Vision II CSE 252B: Computer Vision II Lecturer: Serge Belongie Scribes: Jeremy Pollock and Neil Alldrin LECTURE 14 Robust Feature Matching 14.1. Introduction Last lecture we learned how to find interest points

More information

3D Point Cloud Processing

3D Point Cloud Processing 3D Point Cloud Processing The image depicts how our robot Irma3D sees itself in a mirror. The laser looking into itself creates distortions as well as changes in intensity that give the robot a single

More information

Geometric Modeling in Graphics

Geometric Modeling in Graphics Geometric Modeling in Graphics Part 10: Surface reconstruction Martin Samuelčík www.sccg.sk/~samuelcik samuelcik@sccg.sk Curve, surface reconstruction Finding compact connected orientable 2-manifold surface

More information

Association-Matrix-Based Sample Consensus Approach for Automated Registration of Terrestrial Laser Scans Using Linear Features

Association-Matrix-Based Sample Consensus Approach for Automated Registration of Terrestrial Laser Scans Using Linear Features Association-Matrix-Based Sample Consensus Approach for Automated Registration of Terrestrial Laser Scans Using Linear Features Kaleel Al-Durgham and Ayman Habib Abstract This paper presents an approach

More information

Three-Dimensional Sensors Lecture 6: Point-Cloud Registration

Three-Dimensional Sensors Lecture 6: Point-Cloud Registration Three-Dimensional Sensors Lecture 6: Point-Cloud Registration Radu Horaud INRIA Grenoble Rhone-Alpes, France Radu.Horaud@inria.fr http://perception.inrialpes.fr/ Point-Cloud Registration Methods Fuse data

More information

The Curse of Dimensionality

The Curse of Dimensionality The Curse of Dimensionality ACAS 2002 p1/66 Curse of Dimensionality The basic idea of the curse of dimensionality is that high dimensional data is difficult to work with for several reasons: Adding more

More information

Elective in A.I. Robot Programming

Elective in A.I. Robot Programming Elective in A.I. Robot Programming 2014/2015 Prof: Daniele Nardi Perception with RGB D sensors Jacopo Serafin Point Cloud Library (PCL), Surface Normals, Generalized ICP Contact Jacopo Serafin Ph.D. Student

More information

Rigid ICP registration with Kinect

Rigid ICP registration with Kinect Rigid ICP registration with Kinect Students: Yoni Choukroun, Elie Semmel Advisor: Yonathan Aflalo 1 Overview.p.3 Development of the project..p.3 Papers p.4 Project algorithm..p.6 Result of the whole body.p.7

More information

3D Terrain Sensing System using Laser Range Finder with Arm-Type Movable Unit

3D Terrain Sensing System using Laser Range Finder with Arm-Type Movable Unit 3D Terrain Sensing System using Laser Range Finder with Arm-Type Movable Unit 9 Toyomi Fujita and Yuya Kondo Tohoku Institute of Technology Japan 1. Introduction A 3D configuration and terrain sensing

More information

Multi-stable Perception. Necker Cube

Multi-stable Perception. Necker Cube Multi-stable Perception Necker Cube Spinning dancer illusion, Nobuuki Kaahara Fitting and Alignment Computer Vision Szeliski 6.1 James Has Acknowledgment: Man slides from Derek Hoiem, Lana Lazebnik, and

More information

10/03/11. Model Fitting. Computer Vision CS 143, Brown. James Hays. Slides from Silvio Savarese, Svetlana Lazebnik, and Derek Hoiem

10/03/11. Model Fitting. Computer Vision CS 143, Brown. James Hays. Slides from Silvio Savarese, Svetlana Lazebnik, and Derek Hoiem 10/03/11 Model Fitting Computer Vision CS 143, Brown James Hays Slides from Silvio Savarese, Svetlana Lazebnik, and Derek Hoiem Fitting: find the parameters of a model that best fit the data Alignment:

More information

SIMULTANEOUS REGISTRATION OF MULTIPLE VIEWS OF A 3D OBJECT Helmut Pottmann a, Stefan Leopoldseder a, Michael Hofer a

SIMULTANEOUS REGISTRATION OF MULTIPLE VIEWS OF A 3D OBJECT Helmut Pottmann a, Stefan Leopoldseder a, Michael Hofer a SIMULTANEOUS REGISTRATION OF MULTIPLE VIEWS OF A 3D OBJECT Helmut Pottmann a, Stefan Leopoldseder a, Michael Hofer a a Institute of Geometry, Vienna University of Technology, Wiedner Hauptstr. 8 10, A

More information

Lecture 6 Introduction to Numerical Geometry. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2018

Lecture 6 Introduction to Numerical Geometry. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2018 Lecture 6 Introduction to Numerical Geometry Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2018 Outline Introduction Basic concepts in geometry Discrete geometry Metric for discrete

More information

Segmentation of point clouds

Segmentation of point clouds Segmentation of point clouds George Vosselman INTERNATIONAL INSTITUTE FOR GEO-INFORMATION SCIENCE AND EARTH OBSERVATION Extraction of information from point clouds 1 Segmentation algorithms Extraction

More information

Miniature faking. In close-up photo, the depth of field is limited.

Miniature faking. In close-up photo, the depth of field is limited. Miniature faking In close-up photo, the depth of field is limited. http://en.wikipedia.org/wiki/file:jodhpur_tilt_shift.jpg Miniature faking Miniature faking http://en.wikipedia.org/wiki/file:oregon_state_beavers_tilt-shift_miniature_greg_keene.jpg

More information

Processing 3D Surface Data

Processing 3D Surface Data Processing 3D Surface Data Computer Animation and Visualisation Lecture 12 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing

More information

Feature Matching and RANSAC

Feature Matching and RANSAC Feature Matching and RANSAC Recognising Panoramas. [M. Brown and D. Lowe,ICCV 2003] [Brown, Szeliski, Winder, CVPR 2005] with a lot of slides stolen from Steve Seitz, Rick Szeliski, A. Efros Introduction

More information

HSM3D: Feature-Less Global 6DOF Scan-Matching in the Hough/Radon Domain

HSM3D: Feature-Less Global 6DOF Scan-Matching in the Hough/Radon Domain HSM3D: Feature-Less Global 6DOF Scan-Matching in the Hough/Radon Domain Andrea Censi Stefano Carpin Caltech 3D data alignment how-to The data is 3D, but the sensor motion is 2D? (e.g., 3D scanner on a

More information

Instance-level recognition part 2

Instance-level recognition part 2 Visual Recognition and Machine Learning Summer School Paris 2011 Instance-level recognition part 2 Josef Sivic http://www.di.ens.fr/~josef INRIA, WILLOW, ENS/INRIA/CNRS UMR 8548 Laboratoire d Informatique,

More information

Mosaics. Today s Readings

Mosaics. Today s Readings Mosaics VR Seattle: http://www.vrseattle.com/ Full screen panoramas (cubic): http://www.panoramas.dk/ Mars: http://www.panoramas.dk/fullscreen3/f2_mars97.html Today s Readings Szeliski and Shum paper (sections

More information

Measuring and Visualizing Geometrical Differences Using a Consumer Grade Range Camera

Measuring and Visualizing Geometrical Differences Using a Consumer Grade Range Camera Measuring and Visualizing Geometrical Differences Using a Consumer Grade Range Camera Gijs Ramaekers Supervisor: Jonny Gustafsson Abstract Measuring geometrically complex parts is currently very labor

More information

A Multi-Resolution ICP with Heuristic Closest Point Search for Fast and Robust 3D Registration of Range Images

A Multi-Resolution ICP with Heuristic Closest Point Search for Fast and Robust 3D Registration of Range Images A Multi-Resolution ICP with Heuristic Closest Point Search for Fast and Robust 3D Registration of Range Images Timothée Jost and Heinz Hügli Institute of Microtechnology, University of Neuchâtel, Breguet

More information

Surface Reconstruction with MLS

Surface Reconstruction with MLS Surface Reconstruction with MLS Tobias Martin CS7960, Spring 2006, Feb 23 Literature An Adaptive MLS Surface for Reconstruction with Guarantees, T. K. Dey and J. Sun A Sampling Theorem for MLS Surfaces,

More information

Overview. Efficient Simplification of Point-sampled Surfaces. Introduction. Introduction. Neighborhood. Local Surface Analysis

Overview. Efficient Simplification of Point-sampled Surfaces. Introduction. Introduction. Neighborhood. Local Surface Analysis Overview Efficient Simplification of Pointsampled Surfaces Introduction Local surface analysis Simplification methods Error measurement Comparison PointBased Computer Graphics Mark Pauly PointBased Computer

More information

Unconstrained 3D-Mesh Generation Applied to Map Building

Unconstrained 3D-Mesh Generation Applied to Map Building Unconstrained 3D-Mesh Generation Applied to Map Building Diego Viejo and Miguel Cazorla Robot Vision Group Departamento de Ciencia de la Computación e Inteligencia Artificial Universidad de Alicante E-03690,

More information

Processing 3D Surface Data

Processing 3D Surface Data Processing 3D Surface Data Computer Animation and Visualisation Lecture 17 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing

More information

DETECTION AND ROBUST ESTIMATION OF CYLINDER FEATURES IN POINT CLOUDS INTRODUCTION

DETECTION AND ROBUST ESTIMATION OF CYLINDER FEATURES IN POINT CLOUDS INTRODUCTION DETECTION AND ROBUST ESTIMATION OF CYLINDER FEATURES IN POINT CLOUDS Yun-Ting Su James Bethel Geomatics Engineering School of Civil Engineering Purdue University 550 Stadium Mall Drive, West Lafayette,

More information

Feature Detectors and Descriptors: Corners, Lines, etc.

Feature Detectors and Descriptors: Corners, Lines, etc. Feature Detectors and Descriptors: Corners, Lines, etc. Edges vs. Corners Edges = maxima in intensity gradient Edges vs. Corners Corners = lots of variation in direction of gradient in a small neighborhood

More information

Geometric Registration for Deformable Shapes 3.3 Advanced Global Matching

Geometric Registration for Deformable Shapes 3.3 Advanced Global Matching Geometric Registration for Deformable Shapes 3.3 Advanced Global Matching Correlated Correspondences [ASP*04] A Complete Registration System [HAW*08] In this session Advanced Global Matching Some practical

More information

Robust Range Image Registration using a Common Plane

Robust Range Image Registration using a Common Plane VRVis Technical Report 1 Robust Range Image Registration using a Common Plane Joachim Bauer bauer@icg.vrvis.at Konrad Karner karner@vrvis.at Andreas Klaus klaus@vrvis.at Roland Perko University of Technology

More information

A 3D Point Cloud Registration Algorithm based on Feature Points

A 3D Point Cloud Registration Algorithm based on Feature Points International Conference on Information Sciences, Machinery, Materials and Energy (ICISMME 2015) A 3D Point Cloud Registration Algorithm based on Feature Points Yi Ren 1, 2, a, Fucai Zhou 1, b 1 School

More information

Uncertainties: Representation and Propagation & Line Extraction from Range data

Uncertainties: Representation and Propagation & Line Extraction from Range data 41 Uncertainties: Representation and Propagation & Line Extraction from Range data 42 Uncertainty Representation Section 4.1.3 of the book Sensing in the real world is always uncertain How can uncertainty

More information

DIRSAC: A Directed Sample and Consensus Algorithm for Localization with Quasi Degenerate Data

DIRSAC: A Directed Sample and Consensus Algorithm for Localization with Quasi Degenerate Data DIRSAC: A Directed Sample and Consensus Algorithm for Localization with Quasi Degenerate Data Chris L Baker a,b, William Hoff c a National Robotics Engineering Center (chris@chimail.net) b Neya Systems

More information

L17. OCCUPANCY MAPS. NA568 Mobile Robotics: Methods & Algorithms

L17. OCCUPANCY MAPS. NA568 Mobile Robotics: Methods & Algorithms L17. OCCUPANCY MAPS NA568 Mobile Robotics: Methods & Algorithms Today s Topic Why Occupancy Maps? Bayes Binary Filters Log-odds Occupancy Maps Inverse sensor model Learning inverse sensor model ML map

More information

Geometric Modeling Summer Semester 2012 Point-Based Modeling

Geometric Modeling Summer Semester 2012 Point-Based Modeling Geometric Modeling Summer Semester 2012 Point-Based Modeling 3D Scanning MLS-Surfaces Reconstruction & Registration Overview... Topics: Subdivision Surfaces Implicit Functions Variational Modeling Point-Based

More information

Object Detection in Point Clouds Using Conformal Geometric Algebra

Object Detection in Point Clouds Using Conformal Geometric Algebra Adv. Appl. Clifford Algebras 27 (2017), 1961 1976 c 2017 The Author(s) This article is an open access publication 0188-7009/031961-16 published online February 16, 2017 DOI 10.1007/s00006-017-0759-1 Advances

More information

GEOMETRIC TOOLS FOR COMPUTER GRAPHICS

GEOMETRIC TOOLS FOR COMPUTER GRAPHICS GEOMETRIC TOOLS FOR COMPUTER GRAPHICS PHILIP J. SCHNEIDER DAVID H. EBERLY MORGAN KAUFMANN PUBLISHERS A N I M P R I N T O F E L S E V I E R S C I E N C E A M S T E R D A M B O S T O N L O N D O N N E W

More information

Prof. Noah Snavely CS Administrivia. A4 due on Friday (please sign up for demo slots)

Prof. Noah Snavely CS Administrivia. A4 due on Friday (please sign up for demo slots) Robust fitting Prof. Noah Snavely CS111 http://www.cs.cornell.edu/courses/cs111 Administrivia A due on Friday (please sign up for demo slots) A5 will be out soon Prelim is coming up, Tuesday, / Roadmap

More information

Two-view geometry Computer Vision Spring 2018, Lecture 10

Two-view geometry Computer Vision Spring 2018, Lecture 10 Two-view geometry http://www.cs.cmu.edu/~16385/ 16-385 Computer Vision Spring 2018, Lecture 10 Course announcements Homework 2 is due on February 23 rd. - Any questions about the homework? - How many of

More information

Midterm Examination CS 534: Computational Photography

Midterm Examination CS 534: Computational Photography Midterm Examination CS 534: Computational Photography November 3, 2016 NAME: Problem Score Max Score 1 6 2 8 3 9 4 12 5 4 6 13 7 7 8 6 9 9 10 6 11 14 12 6 Total 100 1 of 8 1. [6] (a) [3] What camera setting(s)

More information

Homography estimation

Homography estimation RANSAC continued Homography estimation x w? ~x img H~x w Homography estimation? x img ~x w = H 1 ~x img Homography estimation (0,0) (1,0) (6.4,2.8) (8.0,2.9) (5.6, 4.0) (7.8, 4.2) (0,1) (1,1) Ah =0s.tkhk

More information

Intersecting Simple Surfaces. Dr. Scott Schaefer

Intersecting Simple Surfaces. Dr. Scott Schaefer Intersecting Simple Surfaces Dr. Scott Schaefer 1 Types of Surfaces Infinite Planes Polygons Convex Ray Shooting Winding Number Spheres Cylinders 2/66 Infinite Planes Defined by a unit normal n and a point

More information

Homographies and RANSAC

Homographies and RANSAC Homographies and RANSAC Computer vision 6.869 Bill Freeman and Antonio Torralba March 30, 2011 Homographies and RANSAC Homographies RANSAC Building panoramas Phototourism 2 Depth-based ambiguity of position

More information

Correspondence. CS 468 Geometry Processing Algorithms. Maks Ovsjanikov

Correspondence. CS 468 Geometry Processing Algorithms. Maks Ovsjanikov Shape Matching & Correspondence CS 468 Geometry Processing Algorithms Maks Ovsjanikov Wednesday, October 27 th 2010 Overall Goal Given two shapes, find correspondences between them. Overall Goal Given

More information

Surface Reconstruction. Gianpaolo Palma

Surface Reconstruction. Gianpaolo Palma Surface Reconstruction Gianpaolo Palma Surface reconstruction Input Point cloud With or without normals Examples: multi-view stereo, union of range scan vertices Range scans Each scan is a triangular mesh

More information

Bounding Volume Hierarchies

Bounding Volume Hierarchies Tutorial: Real-Time Collision Detection for Dynamic Virtual Environments Bounding Volume Hierarchies Stefan Kimmerle WSI/GRIS University of Tübingen Outline Introduction Bounding Volume Types Hierarchy

More information

CITS 4402 Computer Vision

CITS 4402 Computer Vision CITS 4402 Computer Vision Prof Ajmal Mian Lecture 12 3D Shape Analysis & Matching Overview of this lecture Revision of 3D shape acquisition techniques Representation of 3D data Applying 2D image techniques

More information

Motion Tracking and Event Understanding in Video Sequences

Motion Tracking and Event Understanding in Video Sequences Motion Tracking and Event Understanding in Video Sequences Isaac Cohen Elaine Kang, Jinman Kang Institute for Robotics and Intelligent Systems University of Southern California Los Angeles, CA Objectives!

More information

Computer Graphics Ray Casting. Matthias Teschner

Computer Graphics Ray Casting. Matthias Teschner Computer Graphics Ray Casting Matthias Teschner Outline Context Implicit surfaces Parametric surfaces Combined objects Triangles Axis-aligned boxes Iso-surfaces in grids Summary University of Freiburg

More information

Dynamic Geometry Processing

Dynamic Geometry Processing Dynamic Geometry Processing EG 2012 Tutorial Will Chang, Hao Li, Niloy Mitra, Mark Pauly, Michael Wand Tutorial: Dynamic Geometry Processing 1 Articulated Global Registration Introduction and Overview

More information

Mobile Point Fusion. Real-time 3d surface reconstruction out of depth images on a mobile platform

Mobile Point Fusion. Real-time 3d surface reconstruction out of depth images on a mobile platform Mobile Point Fusion Real-time 3d surface reconstruction out of depth images on a mobile platform Aaron Wetzler Presenting: Daniel Ben-Hoda Supervisors: Prof. Ron Kimmel Gal Kamar Yaron Honen Supported

More information

ECE 470: Homework 5. Due Tuesday, October 27 in Seth Hutchinson. Luke A. Wendt

ECE 470: Homework 5. Due Tuesday, October 27 in Seth Hutchinson. Luke A. Wendt ECE 47: Homework 5 Due Tuesday, October 7 in class @:3pm Seth Hutchinson Luke A Wendt ECE 47 : Homework 5 Consider a camera with focal length λ = Suppose the optical axis of the camera is aligned with

More information

3D Photography: Stereo

3D Photography: Stereo 3D Photography: Stereo Marc Pollefeys, Torsten Sattler Spring 2016 http://www.cvg.ethz.ch/teaching/3dvision/ 3D Modeling with Depth Sensors Today s class Obtaining depth maps / range images unstructured

More information

Model Fitting. Introduction to Computer Vision CSE 152 Lecture 11

Model Fitting. Introduction to Computer Vision CSE 152 Lecture 11 Model Fitting CSE 152 Lecture 11 Announcements Homework 3 is due May 9, 11:59 PM Reading: Chapter 10: Grouping and Model Fitting What to do with edges? Segment linked edge chains into curve features (e.g.,

More information

Accurate 3D Face and Body Modeling from a Single Fixed Kinect

Accurate 3D Face and Body Modeling from a Single Fixed Kinect Accurate 3D Face and Body Modeling from a Single Fixed Kinect Ruizhe Wang*, Matthias Hernandez*, Jongmoo Choi, Gérard Medioni Computer Vision Lab, IRIS University of Southern California Abstract In this

More information

Fitting. Fitting. Slides S. Lazebnik Harris Corners Pkwy, Charlotte, NC

Fitting. Fitting. Slides S. Lazebnik Harris Corners Pkwy, Charlotte, NC Fitting We ve learned how to detect edges, corners, blobs. Now what? We would like to form a higher-level, more compact representation of the features in the image by grouping multiple features according

More information

3D Photography: Active Ranging, Structured Light, ICP

3D Photography: Active Ranging, Structured Light, ICP 3D Photography: Active Ranging, Structured Light, ICP Kalin Kolev, Marc Pollefeys Spring 2013 http://cvg.ethz.ch/teaching/2013spring/3dphoto/ Schedule (tentative) Feb 18 Feb 25 Mar 4 Mar 11 Mar 18 Mar

More information

Range Image Registration Based on Circular Features

Range Image Registration Based on Circular Features Range Image Registration Based on Circular Features Cecilia Chao Chen Ioannis Stamos Dept. of Computer Science Dept. of Computer Science Graduate Center / CUNY Hunter College / CUNY New York, NY 10016

More information

Model Fitting, RANSAC. Jana Kosecka

Model Fitting, RANSAC. Jana Kosecka Model Fitting, RANSAC Jana Kosecka Fitting: Overview If we know which points belong to the line, how do we find the optimal line parameters? Least squares What if there are outliers? Robust fitting, RANSAC

More information

Personal Navigation and Indoor Mapping: Performance Characterization of Kinect Sensor-based Trajectory Recovery

Personal Navigation and Indoor Mapping: Performance Characterization of Kinect Sensor-based Trajectory Recovery Personal Navigation and Indoor Mapping: Performance Characterization of Kinect Sensor-based Trajectory Recovery 1 Charles TOTH, 1 Dorota BRZEZINSKA, USA 2 Allison KEALY, Australia, 3 Guenther RETSCHER,

More information

FARO Scanning Plugin

FARO Scanning Plugin FARO Scanning Plugin For compatibility with Geomagic products, see Release Notes for. Document version E. Copyright 2005, Raindrop Geomagic, Inc. The FARO scanner is a seven-axis measurement device from

More information

Network Traffic Measurements and Analysis

Network Traffic Measurements and Analysis DEIB - Politecnico di Milano Fall, 2017 Introduction Often, we have only a set of features x = x 1, x 2,, x n, but no associated response y. Therefore we are not interested in prediction nor classification,

More information

Fundamentals of 3D. Lecture 4: Debriefing: ICP (kd-trees) Homography Graphics pipeline. Frank Nielsen 5 octobre 2011

Fundamentals of 3D. Lecture 4: Debriefing: ICP (kd-trees) Homography Graphics pipeline. Frank Nielsen 5 octobre 2011 INF555 Fundamentals of 3D Lecture 4: Debriefing: ICP (kd-trees) Homography Graphics pipeline Frank Nielsen 5 octobre 2011 nielsen@lix.polytechnique.fr ICP: Iterative Closest Point Algorithm at a glance

More information

Feature Based Registration - Image Alignment

Feature Based Registration - Image Alignment Feature Based Registration - Image Alignment Image Registration Image registration is the process of estimating an optimal transformation between two or more images. Many slides from Alexei Efros http://graphics.cs.cmu.edu/courses/15-463/2007_fall/463.html

More information

Recognition: Face Recognition. Linda Shapiro EE/CSE 576

Recognition: Face Recognition. Linda Shapiro EE/CSE 576 Recognition: Face Recognition Linda Shapiro EE/CSE 576 1 Face recognition: once you ve detected and cropped a face, try to recognize it Detection Recognition Sally 2 Face recognition: overview Typical

More information

Registration of Low Cost Maps within Large Scale MMS Maps

Registration of Low Cost Maps within Large Scale MMS Maps 2016 IEEE 18th International Conference on High Performance Computing and Communications; IEEE 14th International Conference on Smart City; IEEE 2nd International Conference on Data Science and Systems

More information

Robot Mapping. TORO Gradient Descent for SLAM. Cyrill Stachniss

Robot Mapping. TORO Gradient Descent for SLAM. Cyrill Stachniss Robot Mapping TORO Gradient Descent for SLAM Cyrill Stachniss 1 Stochastic Gradient Descent Minimize the error individually for each constraint (decomposition of the problem into sub-problems) Solve one

More information

Instance-level recognition II.

Instance-level recognition II. Reconnaissance d objets et vision artificielle 2010 Instance-level recognition II. Josef Sivic http://www.di.ens.fr/~josef INRIA, WILLOW, ENS/INRIA/CNRS UMR 8548 Laboratoire d Informatique, Ecole Normale

More information

3D Models and Matching

3D Models and Matching 3D Models and Matching representations for 3D object models particular matching techniques alignment-based systems appearance-based systems GC model of a screwdriver 1 3D Models Many different representations

More information

Final Exam Assigned: 11/21/02 Due: 12/05/02 at 2:30pm

Final Exam Assigned: 11/21/02 Due: 12/05/02 at 2:30pm 6.801/6.866 Machine Vision Final Exam Assigned: 11/21/02 Due: 12/05/02 at 2:30pm Problem 1 Line Fitting through Segmentation (Matlab) a) Write a Matlab function to generate noisy line segment data with

More information

Agenda. Rotations. Camera models. Camera calibration. Homographies

Agenda. Rotations. Camera models. Camera calibration. Homographies Agenda Rotations Camera models Camera calibration Homographies D Rotations R Y = Z r r r r r r r r r Y Z Think of as change of basis where ri = r(i,:) are orthonormal basis vectors r rotated coordinate

More information

Stereo and Epipolar geometry

Stereo and Epipolar geometry Previously Image Primitives (feature points, lines, contours) Today: Stereo and Epipolar geometry How to match primitives between two (multiple) views) Goals: 3D reconstruction, recognition Jana Kosecka

More information

SOLUTION FREQUENCY-BASED PROCEDURE FOR AUTOMATED REGISTRATION OF TERRESTRIAL LASER SCANS USING LINEAR FEATURES INTRODUCTION

SOLUTION FREQUENCY-BASED PROCEDURE FOR AUTOMATED REGISTRATION OF TERRESTRIAL LASER SCANS USING LINEAR FEATURES INTRODUCTION SOLUTION FREQUENCY-BASED PROCEDURE FOR AUTOMATED REGISTRATION OF TERRESTRIAL LASER SCANS USING LINEAR FEATURES ABSTRACT Kaleel Al-Durgham, Ayman Habib, Mehdi Mazaheri Department of Geomatics Engineering,

More information