Analisi ed approssimazione di alcuni problemi appartenenti alla classe Shape-from-X
|
|
- Timothy Nelson
- 5 years ago
- Views:
Transcription
1 Analisi ed approssimazione di alcuni problemi appartenenti alla classe Shape-from-X Silvia Tozza Assegnista INdAM Unità di ricerca INdAM presso Dip. di Matematica, Sapienza Università di Roma Convegno GNCS 2018, Hotel Belvedere, Montecatini Terme Febbraio 2018
2 Introduction - Shape-from-X Problems Goal: Reconstruction of the shape of an object starting from some kind of data. Shape-from-Shading Data: The grey-level measured in an image of the object Problem Surface Photo Findsurface(s)that givethesameimage(s) Shape-from-Polarization Data: A polarization image of the object.
3 Linear Differen*al Constraints for Photopolarimetric Height Es*ma*on Silvia Tozza William A.P. Smith Dizhong Zhu Ravi Ramamoorthi Edwin R. Hancock
4 Overview Polarimetric images of glossy object in uncontrolled, outdoor illumina#on Degree of polarisa#on Phase angle Unpolarised intensity Polarisa#on image Es#mated ligh#ng Es#mated depth map Texture mapped surface, novel pose A unified PDE system for height estimation Only need to solve large, sparse linear system A polarisation image from multichannel data Arbitrary uncalibrated illuminations (estimated) Remark Uncalibrated two source photometric stereo solvable with two polarisation images.
5 Overview Polarimetric images of glossy object in uncontrolled, outdoor illumina#on Degree of polarisa#on Phase angle Unpolarised intensity Polarisa#on image Es#mated ligh#ng Es#mated depth map Texture mapped surface, novel pose A unified PDE system for height estimation Only need to solve large, sparse linear system A polarisation image from multichannel data Arbitrary uncalibrated illuminations (estimated) Remark Uncalibrated two source photometric stereo solvable with two polarisation images.
6 Outline Introduction a. Polarimetric image capturing b. Polarisation image Multichannel polarisation image estimation Photo-polarimetric height constraints A Unified PDE formulation Two source lighting estimation Numerical Experiments
7 Introduction - Shape from polarisation (SfP) Problem Assumptions: Orthographic projection Refractive index of the surface known Diffuse polarisation model and diffuse reflectance model Dielectric (i.e. non-metallic) material Notations: x = (x, y) is an image point z(x) surface height Light source direction s and viewer direction v, with s v Unit surface normal n(x) formulated via surface gradient n(x) = [ p(x) q(x) 1]T 1 + z(x) 2, where p(x) = x z(x) and q(x) = y z(x).
8 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
9 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
10 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
11 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
12 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
13 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
14 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
15 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
16 Introduction - Polarimetric image capturing Figure: Rotate linear polarising filter in front of camera Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle Figure: Intensity varies sinusoidally
17 Introduction - Polarimetric image capturing The measured intensity at a pixel varies sinusoidally with the polariser angle ϑ j, j {1,..., P}, with P 3: i ϑj (x) = i un (x) ( 1 + ρ(x) cos(2ϑ j 2φ(x)) ) Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle The polarisation image is obtained by decomposing the sinusoid at every pixel into three quantities [Wolff, 1997]. ρ(x) = I max(x) I min (x) I max (x) + I min (x) and i un (x) = I max(x) + I min (x) 2
18 Introduction - Polarisation image 9:% ;9:% <=9:% 8% 0.35 /&)*$+,*-&.%47*#"% !"#$""%&'%(&)*$+,*-&.% /0*,"%1.#)"% 2.(&)*$+,"3%4.5".,+56% Remark Using beam splitters or custom CCDs it is possible to make the required measurements in a single shot.
19 Multichannel polarisation image estimation Colour images (3 channels), polarisation images with two different light source directions (2 channels) or both (6 channels). ρ and φ constant over the channels (they depend only on surface geometry). i un will vary between channels. Multichannel observations in channel c with polariser angle ϑ j i c ϑ j (x) = i c un(x)(1 + ρ(x) cos(2ϑ j 2φ(x))). The system of equations is linear in the unpolarised intensities and, by a change of variables, can be made linear in ρ and φ [Huynh et a., 2010].
20 Multichannel polarisation image estimation We wish to solve a bilinear system and do so in a least squares sense using interleaved alternating minimisation: 1 Fixing ρ and φ and then solve linearly for iun c in each channel 2 Fix the unpolarised intensities and solve linearly for ρ and φ using all channels simultaneously. Concretely, for a single pixel, we solve min C [ I i 1 un (x),..., iun(x) ] C T 2 di, where i 1 un (x),...,ic un (x) C I = (1 + ρ(x) cos(2ϑ 1 2φ(x)))I C., (1 + ρ(x) cos(2ϑ P 2φ(x)))I C d I = [ i 1 ϑ 1 (x),..., i C ϑ 1 (x), i 1 ϑ 2 (x),..., i C ϑ P (x) ] T. and I C denoting the C C identity matrix.
21 Multichannel polarisation image estimation Then, with the unpolarised intensities fixed, we solve for ρ and φ using the following linearisation: min a,b C ρφ [ ] a 2 d b ρφ, where [a b] T = [ρ(x) cos(2φ(x)), ρ(x) sin(2φ(x))] T, iun(x) 1 cos(2ϑ 1) iun(x) 1 sin(2ϑ 1) iϑ 1 1 (x) iun(x) 1... i 1 C ρφ = un(x) cos(2ϑ P ) iun(x) 1 sin(2ϑ P ) i iun(x) 2 cos(2ϑ 1) iun(x) 2 ϑ 1, d sin(2ϑ ρφ = P (x) i 1 un(x) 1) iϑ 2 1 (x) i 2. un(x)... iun(x) C cos(2ϑ P ) iun(x) C sin(2ϑ P ) iϑ C P (x) iun(x) C We estimate ρ and φ from the linear parameters using φ(x) = 1 2 atan2(b, a) and ρ(x) = a 2 + b 2.
22 Multichannel polarisation image estimation The multichannel result is visibly less noisy than the single channel performance. Input Single channel estimation Input Multichannel estimation Figure: Multichannel polarisation image estimation. Left to right: an image from the input sequence; phase angle (φ) and degree of polarisation (ρ) estimated from a single channel; phase angle (φ) and degree of polarisation (ρ) estimated from three colour channels and two light source directions.
23 Photo-polarimetric height constraints - I Diffuse polarisation model ρ d (x) = 4 cos(θ(x)) where η is the refractive index. ) 2 sin(θ(x)) ( 2 η 1 η η 2 sin(θ(x)) 2 sin(θ(x)) ( 2 η + 1 η ) η2 + 2, Typical values for dielectrics η [1.4, 1.6]. We assume η = 1.5 for the rest of the talk. Degree of polarisation constraint We rearrange the previous equation arriving to cos(θ(x)) = n(x) v = f (ρ d (x), η) = 2 ρ + 2 η 2 ρ 2 η 2 + η 4 + ρ η 2 ρ 2 η 4 ρ 2 4 η 3 ρ (ρ 1) (ρ + 1) + 1 η 4 ρ η 4 ρ + η η 2 ρ η 2 ρ 2 η 2 + ρ ρ + 1 where we drop the dependency of ρ d on x for brevity.
24 Photo-polarimetric height constraints - II Lambertian reflectance model where γ(x) is the albedo. Shading constraint i un (x) = γ(x) cos(θ i ) = γ(x)n(x) s, (1) Writing n(x) in terms of z, (1) can be rewritten as follows: i un (x) = γ(x) z(x) s + s z(x), with s = (s 1, s 2 ). (2) 2 Remark If s v (a configuration physically impossible to achieve precisely) then this equation provides no more information than the degree of polarisation. Hence, we assume s v.
25 Photo-polarimetric height constraints - III The second degree of freedom of the surface normal direction is the azimuth angle and this is related to the phase angle of the sinusoid Intensity /4 /2 3 /4 5pi/4 3 /2 7 /4 2 Polariser angle The azimuth angle is either equal to the phase angle or shifted by pi.
26 Photo-polarimetric height constraints - III The polarisation cue restricts n(x) to two possible directions!"#$#%&'&(&)*"+! *,-&./)0#& Collinearity condition: n(x) [cos(φ(x)) sin(φ(x)) 0] T = 0. (3) Since the nonlinear normalisation term is always 0, we can write (3) in terms of the surface gradient arriving to Phase angle constraint p(x) cos(φ(x)) + q(x) sin(φ(x)) = 0
27 Linearisation of height constraints Combining the degree of polarisation and the shading constraints, we can arrive at a linear equation DOP ratio constraint z(x) ṽ + v 3 f (ρ d (x), η) = γ(x) z(x) s + s 3. (4) i un (x) We can rewrite (4) arriving to the following PDE: b(x) z(x) = h(x), (5) where b(x) :=b (f,iun) = i un (x)ṽ γ(x)f (ρ d (x), η) s, (6) h(x) :=h (f,iun) = i un (x)v 3 γ(x)f (ρ d (x), η) s 3. (7) with ṽ = (v 1, v 2 ) and s = (s 1, s 2 ).
28 Linearisation of height constraints Let us consider two unpolarised intensities, i un,1, i un,2, taken from two different light source directions, s, t. By applying the shading constraint twice, once for each light source, we get Intensity ratio constraint i un,2 ( z(x) s + s 3 ) = i un,1 ( z(x) t + t 3 ). (8) we can rewrite (8) as a PDE in the form of (5) with b(x) := b (i un,1,i un,2 ) = i un,2 (x) s i un,1 (x) t, (9) where t = (t 1, t 2 ), and h(x) := h (i un,1,i un,2 ) = i un,2 (x)s 3 i un,1 (x) t 3. (10)
29 A Unified PDE formulation B(x) z(x) = h(x), where B : Ω R J 2, h : Ω R J 1, Ω is the reconstruction domain, J = 2, 3 or 4 depending on the cases. Single light and polarisation formulation [Smith et al., 2016] [ ] (f,i b un) (f,iun) B = 1 b 2, h = [h (f,iun), 0] T, cos φ sin φ with b (f,iun) and h (f,iun) defined by (6) and (7). Here, we need a single polarisation image uniform γ(x)
30 A Unified PDE formulation B(x) z(x) = h(x), where B : Ω R J 2, h : Ω R J 1, Ω is the reconstruction domain, J = 2, 3 or 4 depending on the cases. Proposed 1: Albedo invariant formulation B(x) = [ (i b un,1,i un,2 ) 1 b (i ] un,1,i un,2 ) 2, h(x) = cos φ sin φ [ ] h (i un,1,i un,2 ), 0 where b (i un,1,i un,2 ) and h (i un,1,i un,2 ) defined as in (9) and (10). Here, we need two unpolarised images taken from two different light source direction, s and t.
31 A Unified PDE formulation B(x) z(x) = h(x), where B : Ω R J 2, h : Ω R J 1, Ω is the reconstruction domain, J = 2, 3 or 4 depending on the cases. Proposed 2: Phase invariant formulation B = [b (f,i un,1), b (f,i un,2), b (i un,1,i un,2 ) ] T, and h = [h (f,i un,1), h (f,i un,2), h (i un,1,i un,2 ) ] T, Here, we need two unpolarised images taken from s and t. Knowledge of the albedo map s, t, v non-coplanar to have B not singular. Note that f = f (ρ d (x), η) is the same for the two required images.
32 A Unified PDE formulation B(x) z(x) = h(x), where B : Ω R J 2, h : Ω R J 1, Ω is the reconstruction domain, J = 2, 3 or 4 depending on the cases. Proposed 3: Most constrained formulation b (f,i un,1) 1 b (f,i un,1) 2 B = cos φ b (f,i un,2) 1 b (f,i un,2) 2 b (i un,1,i un,2 ) 1 b (i un,1,i un,2 ) 2 sin φ, h = Here, we combine all the previous constraints. h (f,i un,1) h (f,i un,2) h (i un,1,i un,2 ). We require known albedo. Nevertheless, it is possible to first apply proposed method 1, estimate the albedo and then re-estimate surface height. 0
33 Height estimation via linear least squares We discretize the gradient via finite differences, arriving to Discrete linear system in z Az = h, A = BG, with A R JM M, M the number of pixels, J = 2, 3, or 4 depending on the cases. G R 2M M the matrix of finite difference gradients B R JM 2M is the discrete per-pixel version of B(x), h R JM 1 is the discrete per-pixel version of h(x), z R M 1 the vector of the unknown height values. The system is large but sparse. A is a full-rank matrix for each choice of B that comes from the proposed formulations. A related to [Smith et al., 2016] is full-rank except in one case: s 1, s 2 0, s 1 = s 2 and φ = π/4 at least in one pixel.
34 Two source lighting estimation It is possible to estimate both light source directions simultaneously, in an albedo invariant manner. From the intensity ratio, we have 1 equation per pixel and 6 unknowns. Hp: The intensity of the light source remains constant in each colour channel across the two images the length of the light source vectors is arbitrary. So, we constrain them to unit length. In spherical coordinates: (θ s, α s ) and (θ t, α t ), such that [s 1, s 2, s 3 ] = [cos α s sin θ s, sin α s sin θ s, cos θ s ] and [t 1, t 2, t 3 ] = [cos α t sin θ t, sin α t sin θ t, cos θ t ]. This reduces the number of unknowns to four. We have two possible surface normal directions at each pixel and therefore two possible gradients: p(x) ± cos φ(x) tan θ(x), q(x) ± sin φ(x) tan θ(x).
35 Two source lighting estimation The residuals at pixel x j in channel c are given by either: r j,c(θ s, α s, θ t, α t) =i c un,1(x j)( p(x j)t 1 q(x j)t 2 + t 3) i c un,2(x j)( p(x j)s 1 q(x j)s 2 + s 3) or q j,c(θ s, α s, θ t, α t) =i c un,1(x j)(p(x j)t 1 + q(x j)t 2 + t 3) i c un,2(x j)(p(x j)s 1 + q(x j)s 2 + s 3). Minimisation problem for light source direction estimation min min[r 2 θ s,α s,θ t,α t j,c(θ s, α s, θ t, α t ), qj,c(θ 2 s, α s, θ t, α t )]. j,c This optimisation is non-convex, but even with a random initialisation, it almost always converges to the global minimum. Convex/concave ambiguity (s, t) is a solution (Ts, Tt) is also a solution (with T = diag([ 1, 1, 1])).
36 Numerical Tests - Synthetic case Input Input (uniform albedo) (varying albedo) Prop. 1 Prop. 2 Prop. 3 Prop. 1+3 Varying albedo Uniform albedo [Smith et al., 2016] Ground truth height Figure: Qualitative results on synthetic data. S. Tozza - INdAM Analisi ed approssimazione di problemi appartenenti alla classe Shape-from-X
37 Numerical Tests - Synthetic case Table: HeightS. and Tozzasurface - INdAMnormal Analisi errors ed approssimazione on synthetic di problemi data. appartenenti alla classe Shape-from-X Setting Uniform albedo, known lighting Uniform albedo, estimated lighting Unknown albedo, known lighting Unknown albedo, estimated lighting Method σ = 0% σ = 0.5% σ = 2% Height Normal Height Normal Height Normal (pix) (deg) (pix) (deg) (pix) (deg) [Smith et al., 2016] Prop Prop Prop Prop [Smith et al., 2016] Prop Prop Prop Prop [Smith et al., 2016] Prop Prop Prop Prop [Smith et al., 2016] Prop Prop Prop Prop
38 Numerical Tests - Real objects Input Estimated Normals Estimated Surface Estimated Albedo [Smith et al., 2016] Figure: Qualitative results on real objects with varying albedo obtained by using Prop. 1+3 and comparison to [Smith et al., 2016]. S. Tozza - INdAM Analisi ed approssimazione di problemi appartenenti alla classe Shape-from-X
39 Numerical Tests - Real objects Input Estimated Normals Estimated Albedo Figure: Qualitative results on real objects with varying albedo obtained by using Prop. 1+3 and comparison to [Smith et al., 2016]. S. Tozza - INdAM Analisi ed approssimazione di problemi appartenenti alla classe Shape-from-X
40 Numerical Tests - Real objects Estimated Surface [Smith et al., 2016] Figure: Qualitative results on real objects with varying albedo obtained by using Prop. 1+3 and comparison to [Smith et al., 2016].
41 Conclusions We proposed a unified PDE formulation for recovering height from photo-polarimetric data We proposed a variety of methods that use different combinations of linear constraints The used equations are linear, so depth estimation is simply a linear least squares problem. We proposed a more robust way to estimate a polarisation image from multichannel data We showed how to estimate lighting from two source photo-polarimetric images Together, our methods provide uncalibrated, albedo invariant shape estimation with only two light sources.
42 Work in progress/future Perspectives 1 Moving to a perspective projection 2 Considering more complex reflectance models 3 Exploiting better the information available in specular reflection and polarisation 4 Allowing mixtures of the two polarisation models (diffuse and specular) 5 Adding multiview polarisation images From a single polarisation image [Smith et al., 2016], estimating the albedo [Smith et al., submitted].
43 References Wolff, L.B., Polarization vision: a new sensory approach to image understanding, Image Vision Comput., 15(2): 81 93, Huynh, C.P., Robles-Kelly, A., Hancock, E., Shape and refractive index recovery from single-view polarisation images, In: Proc. CVPR, pp , W. A.P. Smith, R. Ramamoorthi, S. Tozza, Linear depth estimation from an uncalibrated, monocular polarisation image, Lecture Notes in Computer Science 9912, pp , Springer S. Tozza, W. A.P. Smith, D. Zhu, R. Ramamoorthi, E. R. Hancock, Linear Differential Constraints for Photo-polarimetric Height Estimation, 2017 IEEE International Conference on Computer Vision (ICCV), pp , W. A.P. Smith, R. Ramamoorthi, S. Tozza, Height-from-Polarisation with Unknown Lighting or Albedo, IEEE Transactions on Pattern Analysis and Machine Intelligence, submitted.
44 Shape-from-Shading: Works in progress Resolution of the SfS problem via Filtered schemes (in collaboration with M. Falcone, G. Paolucci (Sapienza)) Multi-view/SfS/PS via minimization techniques (in collaboration with Y. Quèau (TUM Munich) and Jean-Denis Durou (IRIT, Toulouse)) Resolution of the SfS problem via Learning? (in collaboration with E. Rodolà (Sapienza))
Linear Differential Constraints for Photo-polarimetric Height Estimation
Linear Differential Constraints for Photo-polarimetric Height Estimation Silvia Tozza Sapienza - Università di Roma tozza@matuniroma1it Ravi Ramamoorthi UC San Diego ravir@csucsdedu William A P Smith University
More informationWHEN unpolarised light is reflected by a surface it
IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, VOL. X, NO. X, AUGUST 2018 1 Height-from-Polarisation with Unknown Lighting or Albedo William A. P. Smith, Member, IEEE, Ravi Ramamoorthi,
More informationGeneral Principles of 3D Image Analysis
General Principles of 3D Image Analysis high-level interpretations objects scene elements Extraction of 3D information from an image (sequence) is important for - vision in general (= scene reconstruction)
More informationImage Processing 1 (IP1) Bildverarbeitung 1
MIN-Fakultät Fachbereich Informatik Arbeitsbereich SAV/BV (KOGS) Image Processing 1 (IP1) Bildverarbeitung 1 Lecture 20: Shape from Shading Winter Semester 2015/16 Slides: Prof. Bernd Neumann Slightly
More informationShape from Diffuse Polarisation
Shape from Diffuse Polarisation Gary Atkinson and Edwin Hancock Department of Computer Science University of York, York, YO1 5DD, UK. atkinson@cs.york.ac.uk Abstract When unpolarised light is reflected
More informationPhotometric Stereo. Lighting and Photometric Stereo. Computer Vision I. Last lecture in a nutshell BRDF. CSE252A Lecture 7
Lighting and Photometric Stereo Photometric Stereo HW will be on web later today CSE5A Lecture 7 Radiometry of thin lenses δa Last lecture in a nutshell δa δa'cosα δacos β δω = = ( z' / cosα ) ( z / cosα
More informationLambertian model of reflectance I: shape from shading and photometric stereo. Ronen Basri Weizmann Institute of Science
Lambertian model of reflectance I: shape from shading and photometric stereo Ronen Basri Weizmann Institute of Science Variations due to lighting (and pose) Relief Dumitru Verdianu Flying Pregnant Woman
More informationAnd if that 120MP Camera was cool
Reflectance, Lights and on to photometric stereo CSE 252A Lecture 7 And if that 120MP Camera was cool Large Synoptic Survey Telescope 3.2Gigapixel camera 189 CCD s, each with 16 megapixels Pixels are 10µm
More informationLambertian model of reflectance II: harmonic analysis. Ronen Basri Weizmann Institute of Science
Lambertian model of reflectance II: harmonic analysis Ronen Basri Weizmann Institute of Science Illumination cone What is the set of images of an object under different lighting, with any number of sources?
More informationPhotometric Stereo. Photometric Stereo. Shading reveals 3-D surface geometry BRDF. HW3 is assigned. An example of photometric stereo
Photometric Stereo Photometric Stereo HW3 is assigned Introduction to Computer Vision CSE5 Lecture 6 Shading reveals 3-D surface geometry Shape-from-shading: Use just one image to recover shape. Requires
More informationShape and Refractive Index from Single-View Spectro-Polarimetric Images
Shape and Refractive Index from Single-View Spectro-Polarimetric Images Cong Phuoc Huynh 1, Antonio Robles-Kelly 1,2,3 and Edwin R. Hancock 4 1 National ICT Australia (NICTA), Locked Bag 8001, Canberra
More informationSome Illumination Models for Industrial Applications of Photometric Stereo
Some Illumination Models for Industrial Applications of Photometric Stereo Yvain Quéau and Jean-Denis Durou Université de Toulouse, IRIT, UMR CNRS 5505, Toulouse, France ABSTRACT Among the possible sources
More informationX.media.publishing. 3D Computer Vision. Efficient Methods and Applications. von Christian Wöhler. 1. Auflage
X.media.publishing 3D Computer Vision Efficient Methods and Applications von Christian Wöhler 1. Auflage 3D Computer Vision Wöhler schnell und portofrei erhältlich bei beck-shop.de DIE FACHBUCHHANDLUNG
More informationAnnouncements. Photometric Stereo. Shading reveals 3-D surface geometry. Photometric Stereo Rigs: One viewpoint, changing lighting
Announcements Today Photometric Stereo, next lecture return to stereo Photometric Stereo Introduction to Computer Vision CSE152 Lecture 16 Shading reveals 3-D surface geometry Two shape-from-x methods
More informationRe-rendering from a Dense/Sparse Set of Images
Re-rendering from a Dense/Sparse Set of Images Ko Nishino Institute of Industrial Science The Univ. of Tokyo (Japan Science and Technology) kon@cvl.iis.u-tokyo.ac.jp Virtual/Augmented/Mixed Reality Three
More informationShape and Refractive Index Recovery from Single-View Polarisation Images
Shape and Refractive Index Recovery from Single-View Polarisation Images Cong Phuoc Huynh 1 Antonio Robles-Kelly 1, Edwin Hancock 3 1 School of Engineering, Australian National University, Canberra ACT
More informationCS4670/5760: Computer Vision Kavita Bala Scott Wehrwein. Lecture 23: Photometric Stereo
CS4670/5760: Computer Vision Kavita Bala Scott Wehrwein Lecture 23: Photometric Stereo Announcements PA3 Artifact due tonight PA3 Demos Thursday Signups close at 4:30 today No lecture on Friday Last Time:
More informationStereo 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 informationInverse Rendering with a Morphable Model: A Multilinear Approach
ALDRIAN, SMITH: INVERSE RENDERING OF FACES WITH A MORPHABLE MODEL 1 Inverse Rendering with a Morphable Model: A Multilinear Approach Oswald Aldrian oswald@cs.york.ac.uk William A. P. Smith wsmith@cs.york.ac.uk
More informationRadiance. Pixels measure radiance. This pixel Measures radiance along this ray
Photometric stereo Radiance Pixels measure radiance This pixel Measures radiance along this ray Where do the rays come from? Rays from the light source reflect off a surface and reach camera Reflection:
More informationMultiple View Geometry
Multiple View Geometry CS 6320, Spring 2013 Guest Lecture Marcel Prastawa adapted from Pollefeys, Shah, and Zisserman Single view computer vision Projective actions of cameras Camera callibration Photometric
More informationDirect Shape Recovery from Photometric Stereo with Shadows
Direct Shape Recovery from Photometric Stereo with Shadows Roberto Mecca, Aaron Wetzler, Ron Kimmel and Alfred Marcel Bruckstein Department of Computer Science Technion - Israel Institute of Technology,
More informationA Factorization Method for Structure from Planar Motion
A Factorization Method for Structure from Planar Motion Jian Li and Rama Chellappa Center for Automation Research (CfAR) and Department of Electrical and Computer Engineering University of Maryland, College
More informationEpipolar geometry contd.
Epipolar geometry contd. Estimating F 8-point algorithm The fundamental matrix F is defined by x' T Fx = 0 for any pair of matches x and x in two images. Let x=(u,v,1) T and x =(u,v,1) T, each match gives
More informationHigh Quality Shape from a Single RGB-D Image under Uncalibrated Natural Illumination
High Quality Shape from a Single RGB-D Image under Uncalibrated Natural Illumination Yudeog Han Joon-Young Lee In So Kweon Robotics and Computer Vision Lab., KAIST ydhan@rcv.kaist.ac.kr jylee@rcv.kaist.ac.kr
More informationOther approaches to obtaining 3D structure
Other approaches to obtaining 3D structure Active stereo with structured light Project structured light patterns onto the object simplifies the correspondence problem Allows us to use only one camera camera
More information3D Shape Recovery of Smooth Surfaces: Dropping the Fixed Viewpoint Assumption
IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE, VOL., NO., 1 3D Shape Recovery of Smooth Surfaces: Dropping the Fixed Viewpoint Assumption Yael Moses Member, IEEE and Ilan Shimshoni Member,
More informationDIFFUSE-SPECULAR SEPARATION OF MULTI-VIEW IMAGES UNDER VARYING ILLUMINATION. Department of Artificial Intelligence Kyushu Institute of Technology
DIFFUSE-SPECULAR SEPARATION OF MULTI-VIEW IMAGES UNDER VARYING ILLUMINATION Kouki Takechi Takahiro Okabe Department of Artificial Intelligence Kyushu Institute of Technology ABSTRACT Separating diffuse
More informationCapturing light. Source: A. Efros
Capturing light Source: A. Efros Review Pinhole projection models What are vanishing points and vanishing lines? What is orthographic projection? How can we approximate orthographic projection? Lenses
More informationA Novel Illumination-Invariant Loss for Monocular 3D Pose Estimation
A Novel Illumination-Invariant Loss for Monocular 3D Pose Estimation Srimal Jayawardena Marcus Hutter Nathan Brewer Australian National University srimal(dot)jayawardena(at)anu(dot)edu(dot)au http://users.cecs.anu.edu.au/~srimalj
More informationBIL Computer Vision Apr 16, 2014
BIL 719 - Computer Vision Apr 16, 2014 Binocular Stereo (cont d.), Structure from Motion Aykut Erdem Dept. of Computer Engineering Hacettepe University Slide credit: S. Lazebnik Basic stereo matching algorithm
More informationStructure from Motion. Prof. Marco Marcon
Structure from Motion Prof. Marco Marcon Summing-up 2 Stereo is the most powerful clue for determining the structure of a scene Another important clue is the relative motion between the scene and (mono)
More informationLigh%ng and Reflectance
Ligh%ng and Reflectance 2 3 4 Ligh%ng Ligh%ng can have a big effect on how an object looks. Modeling the effect of ligh%ng can be used for: Recogni%on par%cularly face recogni%on Shape reconstruc%on Mo%on
More informationPhotometric stereo. Recovering the surface f(x,y) Three Source Photometric stereo: Step1. Reflectance Map of Lambertian Surface
Photometric stereo Illumination Cones and Uncalibrated Photometric Stereo Single viewpoint, multiple images under different lighting. 1. Arbitrary known BRDF, known lighting 2. Lambertian BRDF, known lighting
More information3-D Stereo Using Photometric Ratios
3-D Stereo Using Photometric Ratios Lawrence B. Wolff and Elli Angelopoulou Computer Vision Laboratory Department of Computer Science, The Johns Hopkins University Baltimore, MD 21218 Abstract. We present
More informationComputer Vision I - Algorithms and Applications: Multi-View 3D reconstruction
Computer Vision I - Algorithms and Applications: Multi-View 3D reconstruction Carsten Rother 09/12/2013 Computer Vision I: Multi-View 3D reconstruction Roadmap this lecture Computer Vision I: Multi-View
More informationDirect Surface Reconstruction using Perspective Shape from Shading via Photometric Stereo
Direct Surface Reconstruction using Perspective Shape from Shading via Roberto Mecca joint work with Ariel Tankus and Alfred M. Bruckstein Technion - Israel Institute of Technology Department of Computer
More informationSegmentation and Tracking of Partial Planar Templates
Segmentation and Tracking of Partial Planar Templates Abdelsalam Masoud William Hoff Colorado School of Mines Colorado School of Mines Golden, CO 800 Golden, CO 800 amasoud@mines.edu whoff@mines.edu Abstract
More informationPhotometric Stereo with Auto-Radiometric Calibration
Photometric Stereo with Auto-Radiometric Calibration Wiennat Mongkulmann Takahiro Okabe Yoichi Sato Institute of Industrial Science, The University of Tokyo {wiennat,takahiro,ysato} @iis.u-tokyo.ac.jp
More informationToday. Stereo (two view) reconstruction. Multiview geometry. Today. Multiview geometry. Computational Photography
Computational Photography Matthias Zwicker University of Bern Fall 2009 Today From 2D to 3D using multiple views Introduction Geometry of two views Stereo matching Other applications Multiview geometry
More informationOptic Flow and Basics Towards Horn-Schunck 1
Optic Flow and Basics Towards Horn-Schunck 1 Lecture 7 See Section 4.1 and Beginning of 4.2 in Reinhard Klette: Concise Computer Vision Springer-Verlag, London, 2014 1 See last slide for copyright information.
More informationA Survey of Light Source Detection Methods
A Survey of Light Source Detection Methods Nathan Funk University of Alberta Mini-Project for CMPUT 603 November 30, 2003 Abstract This paper provides an overview of the most prominent techniques for light
More information3D Motion Analysis Based on 2D Point Displacements
3D Motion Analysis Based on 2D Point Displacements 2D displacements of points observed on an unknown moving rigid body may provide information about - the 3D structure of the points - the 3D motion parameters
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationAnnouncement. Lighting and Photometric Stereo. Computer Vision I. Surface Reflectance Models. Lambertian (Diffuse) Surface.
Lighting and Photometric Stereo CSE252A Lecture 7 Announcement Read Chapter 2 of Forsyth & Ponce Might find section 12.1.3 of Forsyth & Ponce useful. HW Problem Emitted radiance in direction f r for incident
More informationPolarimetric Multi-View Stereo
Polarimetric Multi-View Stereo Zhaopeng Cui, Jinwei Gu 1, Boxin Shi 3, Ping Tan and Jan Kautz 1 1 Nvidia Research, USA Simon Fraser University, Canada 3 National Institute of Advanced Industrial Science
More informationRecovering illumination and texture using ratio images
Recovering illumination and texture using ratio images Alejandro Troccoli atroccol@cscolumbiaedu Peter K Allen allen@cscolumbiaedu Department of Computer Science Columbia University, New York, NY Abstract
More informationA Shape from Shading Approach for the Reconstruction of Polyhedral Objects using Genetic Algorithm
A Shape from Shading Approach for the Reconstruction of Polyhedral Objects using Genetic Algorithm MANOJ KUMAR RAMA BHARGAVA R. BALASUBRAMANIAN Indian Institute of Technology Roorkee, Roorkee P.O. Box
More informationSpecularities Reduce Ambiguity of Uncalibrated Photometric Stereo
CENTER FOR MACHINE PERCEPTION CZECH TECHNICAL UNIVERSITY Specularities Reduce Ambiguity of Uncalibrated Photometric Stereo Ondřej Drbohlav Radim Šára {drbohlav,sara}@cmp.felk.cvut.cz O. Drbohlav, and R.
More informationcalibrated coordinates Linear transformation pixel coordinates
1 calibrated coordinates Linear transformation pixel coordinates 2 Calibration with a rig Uncalibrated epipolar geometry Ambiguities in image formation Stratified reconstruction Autocalibration with partial
More informationAssignment #2. (Due date: 11/6/2012)
Computer Vision I CSE 252a, Fall 2012 David Kriegman Assignment #2 (Due date: 11/6/2012) Name: Student ID: Email: Problem 1 [1 pts] Calculate the number of steradians contained in a spherical wedge with
More informationFace Re-Lighting from a Single Image under Harsh Lighting Conditions
Face Re-Lighting from a Single Image under Harsh Lighting Conditions Yang Wang 1, Zicheng Liu 2, Gang Hua 3, Zhen Wen 4, Zhengyou Zhang 2, Dimitris Samaras 5 1 The Robotics Institute, Carnegie Mellon University,
More informationScene Segmentation Using Polarisation Information
Scene Segmentation Using Polarisation Information Nitya Subramaniam Submitted for the degree of Master of Science by Research University of York Computer Science April 2010 Abstract This thesis draws its
More informationStructure from Motion
Structure from Motion Outline Bundle Adjustment Ambguities in Reconstruction Affine Factorization Extensions Structure from motion Recover both 3D scene geoemetry and camera positions SLAM: Simultaneous
More informationA Machine learning approach for Shape From Shading
2nd International Conference on Signal, Image, Vision and their Applications (SIVA13) November 18-20, 2013 - Guelma, Algeria. A Machine learning approach for Shape From Shading Lyes ABADA Laboratory of
More informationGlobal Illumination and the Rendering Equation
CS294-13: Special Topics Lecture #3 Advanced Computer Graphics University of California, Berkeley Handout Date??? Global Illumination and the Rendering Equation Lecture #3: Wednesday, 9 September 2009
More informationCHAPTER 9. Classification Scheme Using Modified Photometric. Stereo and 2D Spectra Comparison
CHAPTER 9 Classification Scheme Using Modified Photometric Stereo and 2D Spectra Comparison 9.1. Introduction In Chapter 8, even we combine more feature spaces and more feature generators, we note that
More informationStereo CSE 576. Ali Farhadi. Several slides from Larry Zitnick and Steve Seitz
Stereo CSE 576 Ali Farhadi Several slides from Larry Zitnick and Steve Seitz Why do we perceive depth? What do humans use as depth cues? Motion Convergence When watching an object close to us, our eyes
More informationGeneralised Perspective Shape from Shading with Oren-Nayar Reflectance
JU ET AL.: GENERALISED PERSPECTIVE SFS WITH OREN-NAYAR REFLECTANCE 1 Generalised Perspective Shape from Shading with Oren-Nayar Reflectance Yong Chul Ju 1 ju@tu-cottbus.de Silvia Tozza 2 tozza@mat.uniroma1.it
More informationLecture 22: Basic Image Formation CAP 5415
Lecture 22: Basic Image Formation CAP 5415 Today We've talked about the geometry of scenes and how that affects the image We haven't talked about light yet Today, we will talk about image formation and
More informationOn Optimal Light Configurations in Photometric Stereo
O. Drbohlav and M. Chantler: On Optimal Light Configurations in Photometric Stereo. In ICCV 25: Proceedings of the 1th IEEE International Conference on Computer Vision, vol. II, pp. 177 1712, Beijing,
More informationComputer Vision Lecture 17
Computer Vision Lecture 17 Epipolar Geometry & Stereo Basics 13.01.2015 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Announcements Seminar in the summer semester
More informationCOMP 558 lecture 16 Nov. 8, 2010
Shading The term shading typically refers to variations in irradiance along a smooth Lambertian surface. Recall that if a surface point is illuminated by parallel light source from direction l, then the
More informationComputer Vision Lecture 17
Announcements Computer Vision Lecture 17 Epipolar Geometry & Stereo Basics Seminar in the summer semester Current Topics in Computer Vision and Machine Learning Block seminar, presentations in 1 st week
More informationStructure from Small Baseline Motion with Central Panoramic Cameras
Structure from Small Baseline Motion with Central Panoramic Cameras Omid Shakernia René Vidal Shankar Sastry Department of Electrical Engineering & Computer Sciences, UC Berkeley {omids,rvidal,sastry}@eecs.berkeley.edu
More informationPerception and Action using Multilinear Forms
Perception and Action using Multilinear Forms Anders Heyden, Gunnar Sparr, Kalle Åström Dept of Mathematics, Lund University Box 118, S-221 00 Lund, Sweden email: {heyden,gunnar,kalle}@maths.lth.se Abstract
More informationSpecular Reflection Separation using Dark Channel Prior
2013 IEEE Conference on Computer Vision and Pattern Recognition Specular Reflection Separation using Dark Channel Prior Hyeongwoo Kim KAIST hyeongwoo.kim@kaist.ac.kr Hailin Jin Adobe Research hljin@adobe.com
More informationApplications of Light Polarization in Vision
Applications of Light Polarization in Vision Lecture #18 Thanks to Yoav Schechner et al, Nayar et al, Larry Wolff, Ikeuchi et al Separating Reflected and Transmitted Scenes Michael Oprescu, www.photo.net
More informationConstructing a 3D Object Model from Multiple Visual Features
Constructing a 3D Object Model from Multiple Visual Features Jiang Yu Zheng Faculty of Computer Science and Systems Engineering Kyushu Institute of Technology Iizuka, Fukuoka 820, Japan Abstract This work
More informationToward a Stratification of Helmholtz Stereopsis
Toward a Stratification of Helmholtz Stereopsis Todd E Zickler Electrical Engineering Yale University New Haven CT 652 zickler@yaleedu Peter N Belhumeur Computer Science Columbia University New York NY
More informationPolarimetric Multi-View Stereo
Polarimetric Multi-View Stereo Zhaopeng Cui 1 Jinwei Gu Boxin Shi 3 Ping Tan 1 Jan Kautz 1 Simon Fraser University NVIDIA Research 3 Artificial Intelligence Research Center, National Institute of AIST
More informationDD2429 Computational Photography :00-19:00
. Examination: DD2429 Computational Photography 202-0-8 4:00-9:00 Each problem gives max 5 points. In order to pass you need about 0-5 points. You are allowed to use the lecture notes and standard list
More informationCSC418 / CSCD18 / CSC2504
5 5.1 Surface Representations As with 2D objects, we can represent 3D objects in parametric and implicit forms. (There are also explicit forms for 3D surfaces sometimes called height fields but we will
More informationPartial Calibration and Mirror Shape Recovery for Non-Central Catadioptric Systems
Partial Calibration and Mirror Shape Recovery for Non-Central Catadioptric Systems Nuno Gonçalves and Helder Araújo Institute of Systems and Robotics - Coimbra University of Coimbra Polo II - Pinhal de
More informationEfficient Photometric Stereo on Glossy Surfaces with Wide Specular Lobes
Efficient Photometric Stereo on Glossy Surfaces with Wide Specular Lobes Hin-Shun Chung Jiaya Jia Department of Computer Science and Engineering The Chinese University of Hong Kong {hschung,leojia}@cse.cuhk.edu.hk
More information3D and Appearance Modeling from Images
3D and Appearance Modeling from Images Peter Sturm 1,Amaël Delaunoy 1, Pau Gargallo 2, Emmanuel Prados 1, and Kuk-Jin Yoon 3 1 INRIA and Laboratoire Jean Kuntzmann, Grenoble, France 2 Barcelona Media,
More informationTwo-View Geometry (Course 23, Lecture D)
Two-View Geometry (Course 23, Lecture D) Jana Kosecka Department of Computer Science George Mason University http://www.cs.gmu.edu/~kosecka General Formulation Given two views of the scene recover the
More informationCamera Parameters Estimation from Hand-labelled Sun Sositions in Image Sequences
Camera Parameters Estimation from Hand-labelled Sun Sositions in Image Sequences Jean-François Lalonde, Srinivasa G. Narasimhan and Alexei A. Efros {jlalonde,srinivas,efros}@cs.cmu.edu CMU-RI-TR-8-32 July
More informationUnderstanding Variability
Understanding Variability Why so different? Light and Optics Pinhole camera model Perspective projection Thin lens model Fundamental equation Distortion: spherical & chromatic aberration, radial distortion
More informationPassive 3D Photography
SIGGRAPH 99 Course on 3D Photography Passive 3D Photography Steve Seitz Carnegie Mellon University http:// ://www.cs.cmu.edu/~seitz Talk Outline. Visual Cues 2. Classical Vision Algorithms 3. State of
More informationOn Shape and Material Recovery from Motion
On Shape and Material Recovery from Motion Manmohan Chandraker NEC Labs America Abstract. We present a framework for the joint recovery of the shape and reflectance of an object with dichromatic BRDF,
More informationStructure from motion
Structure from motion Structure from motion Given a set of corresponding points in two or more images, compute the camera parameters and the 3D point coordinates?? R 1,t 1 R 2,t 2 R 3,t 3 Camera 1 Camera
More informationToward a Stratification of Helmholtz Stereopsis
Appears in Proc. CVPR 2003 Toward a Stratification of Helmholtz Stereopsis Todd E. Zickler Electrical Engineering Yale University New Haven, CT, 06520 zickler@yale.edu Peter N. Belhumeur Computer Science
More informationComputer Vision Systems. Viewing Systems Projections Illuminations Rendering Culling and Clipping Implementations
Computer Vision Systems Viewing Systems Projections Illuminations Rendering Culling and Clipping Implementations Viewing Systems Viewing Transformation Projective Transformation 2D Computer Graphics Devices
More informationStructure from Motion. Introduction to Computer Vision CSE 152 Lecture 10
Structure from Motion CSE 152 Lecture 10 Announcements Homework 3 is due May 9, 11:59 PM Reading: Chapter 8: Structure from Motion Optional: Multiple View Geometry in Computer Vision, 2nd edition, Hartley
More informationA Combinatorial Transparent Surface Modeling from Polarization Images
A Combinatorial Transparent Surface Modeling from Polarization Images Mohamad Ivan Fanany 1, Kiichi Kobayashi 1, and Itsuo Kumazawa 2 1 NHK Engineering Service Inc., 1-10-11 Kinuta Setagaya-ku Tokyo, Japan,
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 informationReflectivity Function based Illumination and Sensor Planning for Industrial Inspection
Reflectivity Function based Illumination and Sensor Planning for Industrial Inspection Marc M. Ellenrieder, Christian Wöhler and Pablo d Angelo a a DaimlerChrysler AG, Research & Technology, REI/AI, P.O.
More informationPhotometric Stereo, Shape from Shading SfS Chapter Szelisky
Photometric Stereo, Shape from Shading SfS Chapter 12.1.1. Szelisky Guido Gerig CS 6320, Spring 2012 Credits: M. Pollefey UNC CS256, Ohad Ben-Shahar CS BGU, Wolff JUN (http://www.cs.jhu.edu/~wolff/course600.461/week9.3/index.htm)
More informationLecture 15: Shading-I. CITS3003 Graphics & Animation
Lecture 15: Shading-I CITS3003 Graphics & Animation E. Angel and D. Shreiner: Interactive Computer Graphics 6E Addison-Wesley 2012 Objectives Learn that with appropriate shading so objects appear as threedimensional
More informationNonrigid Surface Modelling. and Fast Recovery. Department of Computer Science and Engineering. Committee: Prof. Leo J. Jia and Prof. K. H.
Nonrigid Surface Modelling and Fast Recovery Zhu Jianke Supervisor: Prof. Michael R. Lyu Committee: Prof. Leo J. Jia and Prof. K. H. Wong Department of Computer Science and Engineering May 11, 2007 1 2
More informationCS223b Midterm Exam, Computer Vision. Monday February 25th, Winter 2008, Prof. Jana Kosecka
CS223b Midterm Exam, Computer Vision Monday February 25th, Winter 2008, Prof. Jana Kosecka Your name email This exam is 8 pages long including cover page. Make sure your exam is not missing any pages.
More informationShadow Graphs and Surface Reconstruction
Shadow Graphs and Surface Reconstruction Yizhou Yu and Johnny T. Chang Department of Computer Science University of Illinois at Urbana-Champaign Urbana, IL 61801, USA {yyz,jtchang}@uiuc.edu Abstract. We
More informationThe main problem of photogrammetry
Structured Light Structured Light The main problem of photogrammetry to recover shape from multiple views of a scene, we need to find correspondences between the images the matching/correspondence problem
More informationComputer Vision I Name : CSE 252A, Fall 2012 Student ID : David Kriegman Assignment #1. (Due date: 10/23/2012) x P. = z
Computer Vision I Name : CSE 252A, Fall 202 Student ID : David Kriegman E-Mail : Assignment (Due date: 0/23/202). Perspective Projection [2pts] Consider a perspective projection where a point = z y x P
More informationA Canonical Framework for Sequences of Images
A Canonical Framework for Sequences of Images Anders Heyden, Kalle Åström Dept of Mathematics, Lund University Box 118, S-221 00 Lund, Sweden email: andersp@maths.lth.se kalle@maths.lth.se Abstract This
More informationComputer Vision I. Dense Stereo Correspondences. Anita Sellent 1/15/16
Computer Vision I Dense Stereo Correspondences Anita Sellent Stereo Two Cameras Overlapping field of view Known transformation between cameras From disparity compute depth [ Bradski, Kaehler: Learning
More informationShape, Albedo, and Illumination from a Single Image of an Unknown Object Supplementary Material
Shape, Albedo, and from a Single Image of an Unknown Object Supplementary Material Jonathan T. Barron and Jitendra Malik UC Berkeley {barron, malik}@eecs.berkeley.edu 1. Gradient Norm We will detail how
More informationAnalysis of Planar Light Fields from Homogeneous Convex Curved Surfaces Under Distant Illumination
Analysis of Planar Light Fields from Homogeneous Convex Curved Surfaces Under Distant Illumination Ravi Ramamoorthi and Pat Hanrahan {ravir,hanrahan}@graphics.stanford.edu http://graphics.stanford.edu/papers/planarlf/
More informationIdentifying Car Model from Photographs
Identifying Car Model from Photographs Fine grained Classification using 3D Reconstruction and 3D Shape Registration Xinheng Li davidxli@stanford.edu Abstract Fine grained classification from photographs
More information