数字图像处理-光度立体视觉.ppt

信息视觉处理,光度立体视觉 Photometric Stereo,Materials from:,DC&CV Lab.CSIE NTU,双向反射分布函数（BRDF）,The bidirectional reflectance distribution function is the fraction of incident light emitted in one direction when the surface is illuminated from another direction.ratio of the scene radiance to the scene irradiance,Bi-directional Reflectance Function,:Polar angle,:Azimuth angle,双向反射分布函数（BRDF）,:polar angle between surface normal and lens center:azimuth angle of the sensor:emitting from:incident to:irradiance of the incident light at the illuminated surface:radiance of the reflected light:ratio of the scene radiance to the scene irradiance,Differential reflectance model:,双向反射分布函数（BRDF）,定义：输出方向的辐射度与输入方向的辐照度的比率；BRDF可以从零（在该方向没有反射光）变化至无穷大（输出方向的单位辐射度来自于输入方向任意小的辐射度）；BRDF在输入和输出方向是对称的。（著名的Helmholtz互易原理）；尽管对于某些输入与输出的角度，BDRF可以较大，但对大多数来说它不能大。事实上平均值必须相当小。,source,surfacereflection,surface,incidentdirection,bodyreflection,Body Reflection:Diffuse ReflectionMatte AppearanceNon-Homogeneous MediumClay,paper,etc,Surface Reflection:Specular ReflectionGlossy AppearanceHighlightsDominant for Metals,Image Intensity=Body Reflection+Surface Reflection,Mechanisms of Reflection,Body Reflection:Diffuse ReflectionMatte AppearanceNon-Homogeneous MediumClay,paper,etc,Surface Reflection:Specular ReflectionGlossy AppearanceHighlightsDominant for Metals,Many materials exhibitboth Reflections:,Example Surfaces,viewingdirection,surfaceelement,normal,incidentdirection,Lambertian BRDF is simply a constant:,albedo,Surface appears equally bright from ALL directions!(independent of),Surface Radiance:,Commonly used in Vision and Graphics!,source intensity,source intensity I,Diffuse Reflection and Lambertian BRDF,Diffuse Reflection and Lambertian BRDF,CANT perceive the shape of the snow covered terrain!,CAN perceive shape in regions lit by the street lamp!WHY?,White-out:Snow and Overcast Skies,Assume Lambertian Surface with Albedo=1(no absorption)Assume Sky radiance is constant Substituting in above Equation:,Radiance of any patch is the same as Sky radiance!(white-out condition),Diffuse Reflection from Uniform Sky,source intensity I,surfaceelement,normal,incidentdirection,viewingdirection,specular/mirror direction,Mirror BRDF is simply a double-delta function:,Valid for very smooth surfaces.All incident light energy reflected in a SINGLE direction(only when=).,Surface Radiance:,specular albedo,Specular Reflection and Mirror BRDF,Observed Image Color=a x Body Color+b x Specular Reflection Color,R,G,B,Klinker-Shafer-Kanade 1988,Color of Source(Specular reflection),Color of Surface(Diffuse/Body Reflection),Does not specify any specific model forDiffuse/specular reflection,Combing Specular and Diffuse:Dichromatic Reflection,Diffuse and Specular Reflection,diffuse,specular,diffuse+specular,Photometric Stereo,Image Intensity and 3D Geometry,Shading as a cue for shape reconstructionWhat is the relation between intensity and shape?Reflectance Map,Surface Normal,surface normal,Surface Normal,Gradient Space,Reflectance Map,Relates image irradiance I(x,y)to surface orientation(p,q)for given source direction and surface reflectanceLambertian case:,Image irradiance:,Lambertian case,Reflectance Map,Lambertian case,iso-brightnesscontour,Note:is maximum when,Reflectance Map,Glossy surfaces(Torrance-Sparrow reflectance model),diffuse term,specular term,Reflectance Map,Shape from a Single Image?,Given a single image of an object with known surface reflectance taken under a known light source,can we recover the shape of the object?Given R(p,q)(pS,qS)and surface reflectance)can we determine(p,q)uniquely for each image point?,NO,Solution,Take more imagesPhotometric stereoAdd more constraintsShape-from-shading,Photometric Stereo,We can write this in matrix form:,Image irradiance:,Lambertian case:,Photometric Stereo,Solving the Equations,More than Three Light Sources,Get better results by using more lights,Least squares solution:,Solve for as before,Color Images,The case of RGB imagesget three sets of equations,one per color channel:Simple solution:first solve for n using one channelThen substitute known n into above equations to getOr combine three channels and solve for n,Computing light source directions,Trick:place a chrome sphere in the scenethe location of the highlight tells you the source direction,For a perfect mirror,light is reflected about N,Specular Reflection-Recap,We see a highlight when Then s is given as follows:,Computing the Light Source Direction,Can compute N by studying this figureHints:use this equation:can measure c,h,and r in the image,N,rN,C,H,c,h,Chrome sphere that has a highlight at position h in the image,image plane,sphere in 3D,Depth from Normals,Get a similar equation for V2Each normal gives us two linear constraints on zcompute z values by solving a matrix equation,积分方法,积分的方法可以分为局部方法(local approach)和全局方法(global approach)两类。局部方法比较容易操作，计算效率也比较高，可是它会传递误差，因此对光度立体得到的法向数据要求很精确。,局部和全局积分方法,局部方法：Coleman和 Jain提出了一种扫描算法，是从点P相邻两点的法向矢量来计算P点的深度。Healey 和 Jain是从P 点的相邻8个点的法向来恢复P点的深度的。Wu 和 Li 用格林理论和多路径积分去恢复相对高度。最简单的局部积分方法就是利用梯度从不同路径积分，然后去平均全局方法：把表面积分当作一个优化问题。因为他把梯度数据当成全局数据，所以这种算法比较抗噪声。Ikeuchi用最小二乘法从针图来估计表面深度。Horn、Frankot and Chellappa、Simchony等也提出了一些相应的全局算法。,积分方法的比较,真实形状,不同方法积分结果,Possion方程优化方法对尺度的保持比较好,Frankot Chellapp频域优化的方法对形态保持的比较好,不同方法积分结果,仿射变换方法,加权能量方程,Limitations,Big problemsDoesnt work for shiny things,semi-translucent thingsShadows,inter-reflectionsSmaller problemsCamera and lights have to be distantCalibration requirementsmeasure light source directions,intensitiescamera response function,Trick for Handling Shadows,Weight each equation by the pixel brightness:,Gives weighted least-squares matrix equation:,Solve for as before,由多个视图获取表面法线与发射率,一个球面组成的5幅图像，是使用正交投影从一个视角得到的。,由多个视图获取表面法线与发射率,表面的反射系数,由多个视图获取表面法线与发射率,法线向量场,由法线获取表面形状,通过积分求形状,实际例子,乒乓球图像：,(a)(b)(c)(d),(e)(f)(g)(h),实际例子,恢复结果：,Original Images,Results-Shape,Shallow reconstruction(effect of interreflections),Accurate reconstruction(after removing interreflections),Results-Albedo,No Shading Information,Original Images,Results-Shape,Results-Albedo,Results,Estimate light source directionsCompute surface normalsCompute albedo valuesEstimate depth from surface normalsRelight the object(with original texture and uniform albedo),成像系统,传统的光度立体视觉视觉需要光源从多个方向依次对物体照明，并拍摄多幅图像，才能完成光度立体的计算需要。近些年来，有学者提出了彩色光度立体的概念，即分别用红、绿、蓝三色光从三个方向对物体进行光照，在一幅彩色图像中分别进行三色光的光度计算，进而恢复立体信息。,最新进展,最新进展,三色光同时记录的光度信息是彼此分离的，因为照相机对三色光的光度记录是分开的，对于Bayer相机来说，虽然没有分别记录，但是CCD上每一个像素点的只对一种色光感应，因此也能够进行光度彩色光度立体。,进一步学习材料,David A.Forsyth等著,林学訚 等译.计算机视觉一种现代方法M.北京：电子工业出版社，2004.6.Dense Photometric stereo,http:/www.cse.ust.hk/cktang/sample_pub/dense_ps_pami06.pdf,