A Preliminary Study of
Image-base Indoor
Navigation with
Panoramic Camera
Motivation
Department of Geomatics, National Cheng Kung University
2
Motivation
Introduction
Department of Geomatics, National Cheng Kung University
4
Image-based method perform well in indoor environments
Input image
Output result
www.ctimes.com.tw www.apk3.com
Flow chart
Establish the control images database Find the corresponding images Measure the coordinate of corresponding points Calculate the position and orientation Control points Back-resection Bundle adjustment SURF RANSAC ExperimentDepartment of Geomatics, National Cheng Kung University
6
Limitation of FOV Few conjugate points
1.Maybe get the bad intersection geometry
2.Need lots of images to establish the control image database
Frame image
!?
What will happened if we replace the frame image by panoramic image?
Panoramic image
• Ladybug 5 • 8000 x 4000
• Horizontal field of view: 360 degrees • Vertical field of view: about 150 degrees • Contain lot of image feature information
Department of Geomatics, National Cheng Kung University
• Calculate the exterior orientation – Single image resection
– SPI bundle adjustment
• Collection of a set control images
Images
Position information
Single image resection
Department of Geomatics, National Cheng Kung University
Bundle adjustment
SPI E(m) N(m) h(m) ω(˚) φ(˚) κ(˚) S1 109.06 239.63 -3.77 4.09 1.59 -97.25 S5 106.18 258.76 -2.29 5.93 3.52 90.31 S6 103.45 236.77 -4.31 5.06 0.54 73.98 SPI E(m) N(m) h(m) ω(˚) φ(˚) κ(˚) S2 109.68 267.22 -3.75 -1.57 -0.05 -99.34 S3 110.37 282.46 -3.66 -0.22 -3.26 -92.13 S4 107.57 290.49 -3.77 2.72 2.28 87.17 Check point Control image Control pointImage matching methods
• Manual matching
– Search conjugate points
– Compute the vectors of conjugate points • Auto-matching with panoramic image
– Detect feature points and matching (SURF)
– Eliminate error matching conjugate points (RANSAC)
Department of Geomatics, National Cheng Kung University
Experiment introduction
• Case I
– Movement test
– Distance between 2 station is 0.95 meters – A1, A2, A3, A4, A5
A1
A5
A2 A3
Trajectory Trajectory
Movement result
Department of Geomatics, National Cheng Kung University
14
Experiment introduction
• Case II
– Rotational test
– Rotation angle is 90 degrees clockwise – B1, B2, B3, B4 B1
B2 B4
Rotation result
Rotational direction
Department of Geomatics, National Cheng Kung University
16
Experiment introduction
• Case III
– Oblique test
– Tilt angle is about 5 degrees – C1, C2, C3
Oblique result
downward upward
Before: red image After: blue image
Department of Geomatics, National Cheng Kung University
Epipolar geometry
Department of Geomatics, National Cheng Kung University
20 epipolar plane
Epipolar plane: plane defined by the camera centers and word point
Epipolar line: set of points that project to the same point in left image, when projected to right image forms a line
• The Essential matrix can be computed directly from the camera coordinates with singular value decomposition
• 𝑥𝑥𝑐𝑐 𝑦𝑦𝑐𝑐 𝑧𝑧𝑐𝑐 𝑒𝑒𝑒𝑒1121 𝑒𝑒𝑒𝑒1222 𝑒𝑒𝑒𝑒1323 𝑒𝑒31 𝑒𝑒32 𝑒𝑒33 𝑥𝑥𝑞𝑞 𝑦𝑦𝑞𝑞 𝑧𝑧𝑞𝑞 = 0 • 𝐴𝐴 = 𝑥𝑥𝑠𝑠𝑥𝑥𝑞𝑞 𝑥𝑥𝑠𝑠𝑦𝑦𝑞𝑞 𝑥𝑥𝑠𝑠𝑧𝑧𝑞𝑞 𝑦𝑦𝑠𝑠𝑥𝑥𝑞𝑞 𝑦𝑦𝑠𝑠𝑦𝑦𝑞𝑞 𝑦𝑦𝑠𝑠𝑧𝑧𝑞𝑞 𝑧𝑧𝑠𝑠𝑥𝑥𝑞𝑞 𝑧𝑧𝑠𝑠𝑦𝑦𝑞𝑞 𝑧𝑧𝑠𝑠𝑧𝑧𝑞𝑞 • 𝑥𝑥 = 𝑒𝑒𝑒𝑒1121 𝑒𝑒𝑒𝑒1222 𝑒𝑒𝑒𝑒1323 𝑒𝑒31 𝑒𝑒32 𝑒𝑒33
• One conjugate point can provide one equation
Movement result
Department of Geomatics, National Cheng Kung University
Oblique result
Department of Geomatics, National Cheng Kung University
Movement Rotation Oblique
Affine(Avg.) 31.8% 38.3% 35%
Essential(Avg.) 77.8% 96% 93%
Total umber of pairs(Avg.) 445 1871 1252
Search ability of RANSAC
• The Essential matrix can be decomposed into R and t with SVD
• Assume projection matrix of control image
– P = [I | 0]
• There are 4 possible solutions of query image
– PA=[Ra |t] – PB=[Ra |-t] – PC=[Rb |t] – PD=[Rb |-t]
Department of Geomatics, National Cheng Kung University
26
θ P p’ 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑃𝑃𝐵𝐵 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑝𝑝𝑝𝐵𝐵 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑃𝑃𝐵𝐵 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑝𝑝𝑝𝐵𝐵 𝜃𝜃 = cos−1( 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑃𝑃𝐵𝐵 ∙ 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑝𝑝𝑝𝐵𝐵 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑃𝑃𝐵𝐵 𝑅𝑅𝐵𝐵𝑀𝑀𝑟𝑟𝑝𝑝𝑝𝐵𝐵 ) θ ≈ 0˚
Possible solution
𝑧𝑧𝐵𝐵 𝑦𝑦𝐵𝐵 𝑥𝑥𝐵𝐵 𝑌𝑌𝑀𝑀 𝑋𝑋𝑀𝑀 𝑍𝑍𝑀𝑀Department of Geomatics, National Cheng Kung University
28
Conclusion
• RANSAC is an algorithm which we can change model for the computation case by case
• Coplanarity constraint is useful for spherical panoramic image • The vectors in the panoramic image is more complicated than
frame images
– Movement: Radial transpire and gather – Rotation: straight
– Oblique: swirl
• Matching ability of SURF:
rotation(~1800)>oblique(~1200)>movement (~450)
Department of Geomatics, National Cheng Kung University
Conclusion
• Essential matrix is not only suitable for the RANSAC model but is good at calculating the relative orientation between two panorama images.
• The way which is used for judging the ambiguity in frame images is not suitable for panoramic image.
• With the angle between vectors, we can get the correct orientation in 4 possible solutions.
• Try another feature detection and matching method such as
SIFT.(Maybe it can provide more robust and reliable conjugate points.)
• Develop the method to find the corresponding control images with query image. (A possible solution is do the image
matching and record the number of conjugate points)
• Use other panoramic cameras to test the robustness of our method.
• Calculate the scale of control and query image.( If two
corresponding control images are found, this mission can be solved.)
Department of Geomatics, National Cheng Kung University
32
• Berthold K.P. Horn, Recovering Baseline and Oreintation from ‘Essential’ Matrix
• Cansm Yildiz, An Implementation on Recognizing Panoramas
• David Nister, An Efficient Solution to the Five-Point Relative Pose Problem
• Daesik Kim, 3D Reconstruction
• Henrik Stewenius, Christopher Engels, David Nister, Recent Developments on Direct Relative Orientation
• Herbert Bay, Andreas Ess, Tinne Tuytelaars, Luc Van Gool, Speed-Up Robust Features
• Richard Hartley and Andrew Zisserman, Multiple View Geometry in computer vision
• Richard I. Hartley, In Defence of the 8-point Algorithm
• Van Vinh Nguyen, Jin Guk Kim, and Jong Weon Lee, Panoramic
Image-References
Thank you
Department of Geomatics, National Cheng Kung University