视觉检测系统报告

视觉检测系统报告

年春季学期研究生课程考核（阅读报告、研究报告）考核科目:视觉测量系统学所在院（系）:电气工程及自动化学院学生所在学科:仪器科学与技术学生姓名:***学

号:10S001***学生类别:工学硕士考核结果: 阅卷人: 视觉测量系统课程报告第一部分视觉测量系统发展现状综述机器视觉自起步发展到现在，已有15年的发展历史。应该说机器视觉作为一种应用系统，其功能特点是随着工业自动化的发展而逐渐完善和发展的。

目前全球整个视觉市场总量大概在60～70亿美元，是按照每年

8、8%的增长速度增长的。而在中国，这个数字目前看来似乎有些庞大，但是随着加工制造业的发展，中国对于机器视觉的需求将承上升趋势。

一、机器视觉的定义及特点简言之，机器视觉就是用机器代替人眼来做测量和判断。机器视觉系统是指通过机器视觉产品（即图像摄取装置，分CMOS和CCD两种）将被摄取目标转换成图像信号，传送给专用的图像处理系统，根据像素分布和亮度、颜色等信息，转变成数字化信号；图像系统对这些信号进行各种运算来抽取目标的特征，进而根据判别的结果来控制现场的设备动作。

机器视觉系统的特点是提高生产的柔性和自动化程度。在一些不适合于人工作业的危险工作环境或人工视觉难以满足要求的场合，常用机器视觉来替代人工视觉；同时在大批量工业生产过程中，用人工视觉检查产品质量效率低且精度不高，用机器视觉检测方法可以大大提高生产效率和生产的自动化程度。而且机器视觉易于实现信息集成，是实现计算机集成制造的基础技术。

正是由于机器视觉系统可以快速获取大量信息，而且易于自动处理，也易于同设计信息以及加工控制信息集成，因此，在现代自动化生产过程中，人们将机器视觉系统广泛地用于工况监视、成品检验和质量控制等领域。在中国，这种应用也在逐渐被认知，且带来最直接的反应就是国内对于机器视觉的需求将越来越多。

二、机器视觉在国内外的应用现状在国外，机器视觉的应用普及主要体现在半导体及电子行业，其中大概40%～50%都集中在半导体行业。具体如PCB印刷电路：各类生产印刷电路板组装技术、设备；单、双面、多层线路板，覆铜板及所需的材料及辅料；辅助设施以及耗材、油墨、药水药剂、配件；电子封装技术与设备；丝网印刷设备及丝网周边材料等。SMT表面贴装：SMT工艺与设备、焊接设备、测试仪器、返修设备及各种辅助工具及配件、SMT材料、贴片剂、胶粘剂、焊剂、焊料及防氧化油、焊膏、清洗剂等；再流焊机、波峰焊机及自动化生产线设备。电子生产加工设备：电子元件制造设备、半导体及集成电路制造设备、元

器件成型设备、电子工模具。机器视觉系统还在质量检测的各个方面已经得到了广泛的应用，并且其产品在应用中占据着举足轻重的地位。除此之外，机器视觉还用于其他各个领域。

而在中国，以上行业本身就属于新兴的领域，再加之机器视觉产品技术的普及不够，导致以上各行业的应用几乎空白，即便是有，也只是低端方面的应用。目前在我国随着配套基础建设的完善，技术、资金的积累，各行各业对采用图像和机器视觉技术的工业自动化、智能化需求开始广泛出现，国内有关大专院校、研究所和企业近两年在图像和机器视觉技术领域进行了积极思索和大胆的尝试，逐步开始了工业现场的应用。其主要应用于制药、印刷、矿泉水瓶盖检测等领域。这些应用大多集中在如药品检测分装、印刷色彩检测等。真正高端的应用还很少，因此，以上相关行业的应用空间还比较大。当然、其他领域如指纹检测等等领域也有着很好的发展空间。

三、中国机器视觉未来发展趋势在机器视觉赖以普及发展的诸多因素中，有技术层面的，也有商业层面的，但制造业的需求是决定性的。制造业的发展，带来了对机器视觉需求的提升；也决定了机器视觉将由过去单纯的采集、分析、传递数据，判断动作，逐渐朝着开放性的方向发展，这一趋势也预示着机器视觉将与自动化更进一步的融合。

需求决定产品，只有满足需求的产品才有生存的空间，这是不变的规律。机器视觉也是如此。

未来，中国机器视觉发展主要表现为以下一些特性：

1、随着产业化的发展对机器视觉的需求将呈上升趋势机器视觉发展空间较大的部分在半导体和电子行业，而据我国相关数据显示，全球集成电路产业复苏迹象明显；与此同时，全球经济衰退使我国集成电路产业获取了市场优势、成本优势、人才回流等优势；国家加大对集成电路产业这一战略领域的规划力度，“信息化带动工业化”，走“新兴工业化道路”为集成电路产业带来了巨大的发展机遇，特别是高端产品和创新产品市场空间巨大，设计环节、国家战略领域、3C应用领域、传统产业类应用领域成为集成电路产业未来几年的重点投资领域。

此外，中国已成为全球集成电路的一个重要需求市场。据相关数据显示，2002年我国集成电路市场需求规模为11

35、5亿元人民币，占世界市场规模的

9、76%。2002年中国集成电路市场总销量为2

83、2亿块，总销售额为11

35、5亿元，同比增长

26、2%。中国已成为近年来世界半导体投资的热点。在全国许多地区，特别是长江三角洲地区，都有新的IC制造线和封装测试线投资兴建，IC设计公司的数量每年成倍增长。在产业政策的引导下，上海、北京、天津和深圳等地出现投资IC的好势头：天津Motorola投资15亿美元，月投

2、5万片的8英寸芯片生产线和上海中芯国际投资14亿美元，月投8英寸芯片硅片

4、2万片的项目已经投入运行。另外，中国半导体行业协会最新调研数据表明，2000年6月到2002年8月两年间，中国IC 产业的投资总额约300亿元，相当于过去40年的投资总和。全国IC设计单位数量两年之间翻两番，已激增到389家，收入过亿元的达7～8家；专业测试公司已有10家左右，我国的IC测试业初具雏形。

就以上数据显示，中国的半导体和电子市场已初具规模，而如此强大的半导体产业将需要高质量的技术做后盾。同时他对于产品的高质量、高集成度的要求将越来越高。恰巧，机器视觉将能帮助他们解决以上的问题，因此该行业将是机器视觉最好的用武之地。同时，对于机器视觉的需求将蒸蒸日上。

2、统一开放的标准是机器视觉发展的原动力目前国内有近数家机器视觉产品厂商，与国外机器视觉产品相比，国内产品最大的差距并不单纯是在技术上，而且还包括品牌和知识产权上。另一现状是目前国内的机器视觉产品主要以代理国外品牌为主，以此来逐渐朝着自主研发产品的路线靠近，起步较晚。未来，机器视觉产品的好坏不能够通过单一因素来衡量，应该逐渐按照国际化的统一标准判定，随着中国自动化的逐渐开放，将带领与其相关的产品技术也逐渐开放。因此，依靠封闭的技术难以促进整个行业的发展，只有形成统一而开放的标准才能让更多的厂商在

相同的平台上开发产品，这也是促进中国机器视觉朝国际化水平发展的原动力。

3、基于嵌入式的产品将取代板卡式产品从产品本身看，机器视觉会越来越趋于依靠PC技术，并且与数据采集等其他控制和测量的集成会更紧密。且基于嵌入式的产品将逐渐取代板卡式产品，这是一个不断增长的趋势。主要原因是随着计算机技术和微电子技术的迅速发展，嵌入式系统应用领域越来越广泛，尤其是其具备低功耗技术的特点得到人们的重视。另外，嵌入式操作系统绝大部分是以C语言为基础的，因此使用C高级语言进行嵌入式系统开发是一项带有基础性的工作，使用高级语言的优点是可以提高工作效率，缩短开发周期，更主要的是开发出的产品可靠性高、可维护性好、便于不断完善和升级换代等。因此，嵌入式产品将会取代板卡式产品。

4、标准化、一体化解决方案也将是机器视觉的必经之路另外，由于机器视觉是自动化的一部分，没有自动化就不会有机器视觉，机器视觉软硬件产品正逐渐成为协作生产制造过程中不同阶段的核心系统，无论是用户还是硬件供应商都将机器视觉产品作为生产线上信息收集的工具，这就要求机器视觉产品大量采用“标准化技术”，直观的说就是要随着自动化的开放而逐渐开放，可以根据用户的需求进行二次开发。当今，自动化企业正在倡导软硬一体化解决方案，机器视觉的厂商在未来5-6年内也应

该不单纯是只提供产品的供应商，而是逐渐向一体化解决方案的系统集成商迈进。

在未来的几年内，随着中国加工制造业的发展，对于机器视觉的需求也逐渐增多；随着机器视觉产品的增多，技术的提高，国内机器视觉的应用状况将由初期的低端转向高端。由于机器视觉的介入，自动化将朝着更智能、更快速的方向发展。另外，由于用户的需求是多样化的，且要求程度也不相同。那么，个性化方案和服务在竞争中将日益重要，即用特殊定制的产品来代替标准化的产品也是机器视觉未来发展的一个取向。机器视觉的应用也将进一步促进自动化技术向智能化发展。第二部分图像传感器(相机)及其产品的发展现状综述目前主流的影像传感器（图像传感器）主要有CCD和CMOS两种。是数码相机、数码摄像头等产品的核心部件。CCD是电荷藕合器件图像传感器的简称，CMOS是互补性氧化金属半导体的简称，它们都采用高感光度的半导体材料制成。能把光线转变成电荷，然后转为信号。两者的生产工艺和所使用的设备和计算机芯片差不多，所以目前很多主流CMOS/CCD的厂商同时也是芯片制造商。

随着近几年半导体制造工艺的成熟，CCD/CMOS产品成本逐渐降低，主要体现在制造工艺上从

5、2微米逐步向45nm演进，使得单位面积的CCD/CMOS能够承载更多的像素单位，目前市场上800万像素的数码相机已经很

常见。在成本上跟电脑内存差不多，容量18个月翻一番，价格却在持续下降。

最近几年数码影像产品的价格已经跌倒普通消费者能够接受的水平，向一般家庭普及。传统相机由于技术的局限，已经无法取得突破性的进步，在成本和性能都被采用影像传感器的新兴影像设备所取代，目前传统影像设备已经逐渐从市场上消失，仅在高端市场上坚守最后一块阵地。

目前影像传感技术已经普及到国民经济的各个部门，比如医疗影像、军事侦查、卫星拍摄、气象预报、安防产业、光学检测、工业自动化控制、指纹检查等。和现代电子计算机技术、软件技术、光学技术的有机结合将促进影像传感器的进一步普及，更多的新应用将呈现在人们的眼前。

CMOS技术发展迅速有可能成为未来主流

当前CCD在成像质量上优于CMOS，所以在数码相机、医疗影像、卫星拍摄等对分辨率要求较高的领域CCD是主流，而另一方面由于CMOS采用标准的半导体生产工艺，生产成本低，耗电少，普遍用在手机和电脑摄像头。由于技术差异，CCD较多得用于静态影像，如拍照；CMOS则擅长于动态影像，如视频监控；另外CMOS 还可以把一些周边电路集成到芯片中，在便携式设备中大量使用。

对比两种技术的未来发展趋势，笔者认为未来CMOS将是主流。我们可以从以下几个方面来进行思考：

一、CMOS影像传感器技术演进速度远远快于CCD技术。CCD 技术的像素从

5、2微米演进至

1、7微米耗时12年之久，而CMOS技术完成这一变化则只用了3年时间。

二、从市场的扩张速度来看，CMOS也高于CCD，xx年预计图像传感器市场交货总额为

7、12亿美元，总体市场份额中CCD出货量从40%减少到30%左右。

三、CMOS已经在成像的通透性、对实物的色彩还原能力等方面迎头赶上，总体性能的超越似乎只是时间问题。

四、目前CMOS已经主导手机拍照市场。CCD在过去三年已经是成熟的技术，由于其技术的独特性，改进困难，在成本上已经很难与CMOS竞争，由于其必须采用高于10伏的电压，功率较高，很难集成到一些便捷式产品中。

五、CCD技术主要掌握在日本人的手中，阻碍了技术的传播，而几乎一般的半导体制造商都可以制造CMOS，所以成本下降很快，也有利于技术进步。

因此，CMOS取代CCD是市场的发展趋势，06年的时候，美光影像集团新市场开发总监Paul Gallagher更是断言，CMOS技术将最终取代CCD技术，在未来4～5年内，CCD市场将不断下滑，预计在10年后CCD市场份额将仅剩10%左右。

全球市场竞争格局和近几年发展趋势

目前CCD市场的主要厂商为索尼、夏普、松下等日本厂商，后来进入图像传感器领域的厂商基本上都采用CMOS技术，因为它采用标准的半导体工艺，技术门槛相对校低，并且技术进步快，被大多数公司看好。索尼牢牢占据着CCD第一大供应商的位置。

跟CCD浓重的日本色彩不同，CMOS厂商则散布于美欧日等发达国家，以xx年为例，CMOS供应商按产值的排名如下：美国OmniVision、美国美光半导体、欧洲意法半导体、日本东芝、日本索尼。根据表二所示，去年OmniVision以

7、25亿美元的销售收入跃居第一位，增长了34%。美光、意法半导体、东芝均出现了不同程度的下滑。索尼在积极补课CMOS，索尼的熊本新厂于07年开始完工并量产，把索尼拉进了CMOS前五。

全文结束》》年9月日本的调查公司TechnoSystemResearch 统计指出，xx年CMOS传感器的全球市场规模达4亿3440万个，为20xx年的

2、2倍；金额则为13亿2300万美元（约1480亿日圆），为20xx年的

1、7倍。其市场规模在数量上已超越CCD。根据市场调研公司In-Stat 07年提供的数据，虽然数码相机增长强劲，但手机仍是图像传感器的主要市场，占总出货量的75%以上。据市场调研公司IC Insights数据，xx年可拍照手机的出货量将从xx年的

3、65亿部增长到

4、75亿部，而且可拍照手机在总体手机市场中所占的比例也在不断上升。据IC Insights的报告，xx年可拍照手机占总体手机出货量的比例将达54%左右，而xx年是45%。

自2002至xx年全球图像传感器市场大约以35%的速度在发展，主要得益于拍照手机的全面普及，还有数码相机需求的飙涨。于是乎，在这种形势下，各大厂商纷纷增加产量，而xx年却迎来了市场的寒冬。年，CCD和CMOS图像传感器合并销售下降7%，主要由于上半年摄像手机库存清理，市场需求增长放缓，以及工业和办公设备应用中的销售疲软，比如制造系统中图像机器、复印机和扫描仪等。CMOS销售下滑12%，CCD增长为零。

ICInsights公司今年四月份发布的分析报告指出，xx年图像传感器市场将出现反弹。预计xx年图像传感器市场规模将达76亿美元，比xx年增长10%。ICInsights分析员RobLineback表示，预计xx年CMOS图像传感器销售将上升19%，达到44亿美元；同时，预计xx年CCD市场下降1%，为32亿美元。不过该报告在撰写的时候可能没考虑到下半年欧美会出现严重金融危机，全球需求萎缩的情况，实际情况要等明年的分析报告。根据信产部提供的数据今年整体手机销量严重下滑，多家国产手机厂商生存艰难，国外品牌也被迫降价，深圳山寨机库存严重，一系列消息，都让人对今年的CMOS市场感到担忧。

国内现状与思考

中国在图像传感器领域目前还处于起步阶段，仍有很长的路要走，由于图像传感器属于半导体产业的一部分，因此，影像传感器的发展取决于国内芯片制造技术的进步。目前具有台湾背景的中芯国际、和舰科技、台积电等厂商在大陆的芯片厂已经颇具规模，欧美的因特尔、意法半导体也相继在中国大连和深圳设厂。根据最新的消息今明两年台湾有可能开放12英寸晶圆及0、13微米制程甚至90纳米制程技术到大陆投资，势必会进一步促进中国芯片产业的发展。在技术选择方面，无疑CMOS必选的方向，CCD技术一直掌握日本人手中，欧美发达国家厂商多年来都无法在该领域分一杯羹，我们应该吸取教训，选对方向。年IBM向中芯国际授权45纳米大批量CMOS技术，xx年10月中芯国际宣布成功开发0、11微米CMOS 图像传感器工艺技术，在此工艺下生产的器件，其分辨率、暗光噪声和相对照度都将得到增强。无疑，中国本土厂商进军影像传感器产业的号角已经吹响！第三部分摄像机模型的Matlab实现在机器视觉中，刚体变换经常用于计算一个经过旋转和平移后的新坐标。比较常用的旋转矩阵的表示形式有欧拉角表示法和四元素表示法。

一、欧拉角表示法用旋转矩阵表示刚体的旋转变换简化了许多运算，但它需要９个元素来完全描述这种旋转变换。被称为欧拉角的三个角度能很好描述刚体的旋转变换。相应的Matlab程序如下：x=[0123 0123 0123 0123];y=[0 0 0 0111122223333];z=[1231231231231231];%定义任意16个点坐标

p=[x;y;z];wx=pi/6;wy=pi/6;wz=pi/6; %分别设定偏转、俯仰、

侧倾角度值

r11=cos(wz)*cos(wy);r12=cos(wz)*sin(wy)*sin(wx)-

sin(wz);r13=cos(wz)*sin(wy)*cos(wx)-

sin(wz)*sin(wx);r21=sin(wz)*cos(wy);r22=sin(wz)*sin(wy)*s in(wx)+cos(wz)*cos(wx);r23=sin(wz)*sin(wy)*cos(wx)+cos(wz )*sin(wx);r31=-

sin(wy);r32=cos(wy)*sin(wx);r33=cos(wy)*cos(wx);R=[r11,r1 2,r13;r21,r22,r23;r31,r32,r33];%求解欧拉角表示法的旋转矩

阵t=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; %设定平移

量

pz=R*p+t;xz=[pz(1),pz(4),pz(7),pz(10),pz(13),pz(16),pz(19 ),pz(22),pz(25),pz(28),pz(31),pz(34),pz(37),pz(40),pz(43) ,pz(46)];yz=[pz(2),pz(5),pz(8),pz(11),pz(14),pz(17),pz(20 ),pz(23),pz(26),pz(29),pz(32),pz(35),pz(38),pz(41),pz(44) ,pz(47)];zz=[pz(3),pz(6),pz(9),pz(12),pz(15),pz(18),pz(21 ),pz(24),pz(27),pz(30),pz(33),pz(36),pz(39),pz(42),pz(45) ,pz(48)];%经坐标变换后的坐标

figure;subplot(1,2,1);grid,stem3(x,y,z);axis([-

1414]);subplot(1,2,2);grid,stem3(xz,yz,zz);axis([-1414]);首先定义16个点的三维坐标，分别为[001]，[102]，[203]，

[301]，[012]，[113]，[211]，[312]，[023]，[121]，[222]，[323]，[031]，[132]，[233]，[333]。而后设定偏转、俯仰、侧倾角度值，本程序将三个角度值都设为30度，平移量设为0。程序仿真后的图形如下：

左侧图为原坐标，右侧图为经变换后的坐标。得出变换后的16个点坐标分别为［0、1250，0、6495，0、7500］，［

1、0000，

1、7321，

1、0000］，［

1、8750，

2、8146，

1、2500］，［

2、3750，

1、9486，-0、7500］，［-0、0335，

2、1740，

1、9330］，［0、8415，

3、2566，

2、1830］，［

1、3415，

2、3905，0、1830］，［

2、2465，

3、4731，0、4330］，［-0、1920，

3、6986，

3、1160］，［0、3080，

2、8325，

1、1160］，［

1、1830，

3、9151，

1、3360］，［

2、0580，

4、9976，

1、6160］，［-0、7255，

3、2745，

2、0490］，［0、1495，

4、3571，

2、2990］，［

1、0245，

5、4396，

2、5490］,[

1、5245

4、5736 0、5490]。

二、四元素表示法四元素是一个四元矢量，可用来描述坐标旋转，实践证明，对于旋转与定位求解问题，用四元素法可以很好的数值解。由以上16个点的前后变换坐标，可得出四元素法的

变换矩阵。程序如下：t1=[t(1),t(4),t(7); t(2),t(5),t(8);

t(3),t(6),t(9)];p1=[p(1),p(4),p(7); p(2),p(5),p(8);

p(3),p(6),p(9)];pz1=[pz(1),pz(4),pz(7); pz(2),pz(5),pz(8); pz(3),pz(6),pz(9)];Rq=(pz1-t1)*pinv(p1);%求逆矩阵%加入噪

声干扰变换后的坐标

pz1pz2=[pz1(1),pz1(2),pz1(3),pz1(4),pz1(5),pz1(6),pz1(7),

pz1(8),pz1(9)];pz3=pz2+rand(1,9);pz4=[pz3(1),pz3(4),pz3(7 ); pz3(2),pz3(5),pz3(8); pz3(3),pz3(6),pz3(9)];Rq1=(pz4-

t1)*pinv(p1);欧拉角法的变换矩阵为四元素法求得的变换矩阵

同样为对空间坐标加入rand产生的随机噪声后，将程序运行3次，求得的变换矩阵分别为由数据可分析出，噪声误差对变换矩

阵的影响还是很大的，因此在用变换坐标求相应的变换矩阵时，

应尽量提高测量精度，减少误差带来的影响。

三、摄像机线性模型摄像机通过成像透镜将三维场景投影到

摄像机二维像平面上，这个投影可用成像变换描述，即摄像机成

像模型。摄像机模型可将空间点的世界坐标变换为图像坐标及像

素坐标。线性模型程序如下：U=1316;V=1035;dX=0、

0000068;dY=0、0000068;u0=U/2;v0=V/2; %设定像素UV，像素尺

寸dXdY及u0v0%给出空间点的世界坐标Xw=5;Yw=5;Zw=5;%给出外

部参数wx=0;wy=0;wz=0; %分别设定偏转、俯仰、侧倾角度值

r11=cos(wz)*cos(wy);r12=cos(wz)*sin(wy)*sin(wx)-

sin(wz);r13=cos(wz)*sin(wy)*cos(wx)-

sin(wz)*sin(wx);r21=sin(wz)*cos(wy);r22=sin(wz)*sin(wy)*s in(wx)+cos(wz)*cos(wx);r23=sin(wz)*sin(wy)*cos(wx)+cos(wz )*sin(wx);r31=-

sin(wy);r32=cos(wy)*sin(wx);r33=cos(wy)*cos(wx);R=[r11,r1 2,r13;r21,r22,r23;r31,r32,r33];t=[0;0;0];M2=[R,t;0 0 01];%求摄像机坐标系下的齐次坐标

camora_xyz=M2*[Xw;Yw;Zw;1];%设定焦距，求图像坐标f=0、012;P=[f 0 0 0;0 f 0 0;0 01 0];image_xyz=P*camora_xyz;%求空间点的图像坐标

X=image_xyz(1)/image_xyz(3);Y=image_xyz(2)/image_xyz(3);%设定内参,求像素坐标M1=[1/dX 0 u0;01/dY v0;0

01];pixel_xyz=M1*[X;Y;1];u=pixel_xyz(1);v=pixel_xyz(2);给出空间点世界坐标，设定外部参数（欧拉角表示法三角均为0度，位移为0），便可求得点的图像坐标及像素坐标。给定9个点的空间世界坐标分别为[1，2，3]，[2，2，3]，[-1，-3，3]，[1，-2，4]，[-5，5，5]，[0、5，0、6，0、7]，[30，31，32]，[0，0，40]，[100，120，130]。运行程序后，得相应的图像坐标分别为[0、0040，0、0080]，[0、0080，0、0080]，[-0、0024，-0、0072]，[0、0030，-0、0060]，[-0、0120，0、0120]，[0、0086，0、0103]，[0、0113，0、0116]，[0，0]，[0、0092，0、0111]。像素坐标分别为[12

46、2，16

94、0]，[18

34、5，16

94、0]，

[3

50、0588，-5

41、3235]，[10

99、2，-3

64、8529]，[-1106、7，22

82、2]，[19

18、5，20

30、1]，[23

12、4，22

27、1]，[658，5

17、5]，[xx、5，21

46、5]。

四、摄像机非线性模型实际上，由于实际的镜头并不是理想的透视成像，而是带有不同程度的畸变，加入畸变参数的模型即为非线性模型。程序如下：U=1316;V=1035;dX=0、

0000068;dY=0、0000068;u0=U/2;v0=V/2; %设定像素UV，像素尺寸dXdY及u0v0%给出空间点的世界坐标Xw=5;Yw=5;Zw=5;%给出外部参数wx=0;wy=0;wz=0; %分别设定偏转、俯仰、侧倾角度值

r11=cos(wz)*cos(wy);r12=cos(wz)*sin(wy)*sin(wx)-

sin(wz);r13=cos(wz)*sin(wy)*cos(wx)-

sin(wz)*sin(wx);r21=sin(wz)*cos(wy);r22=sin(wz)*sin(wy)*s in(wx)+cos(wz)*cos(wx);r23=sin(wz)*sin(wy)*cos(wx)+cos(wz )*sin(wx);r31=-

sin(wy);r32=cos(wy)*sin(wx);r33=cos(wy)*cos(wx);R=[r11,r1 2,r13;r21,r22,r23;r31,r32,r33];t=[0;0;0];M2=[R,t;0 0 01];%求摄像机坐标系下的齐次坐标

camora_xyz=M2*[Xw;Yw;Zw;1];%设定焦距，求图像坐标f=0、012;P=[f 0 0 0;0 f 0 0;0 01 0];image_xyz=P*camora_xyz;%求空间点的图像坐标

fale_X=image_xyz(1)/image_xyz(3);fale_Y=image_xyz(2)/imag e_xyz(3);%加入非线性畸变影响k1=0、

001;r2=fale_X^2+fale_Y^2;derta_x=fale_X*k1*r2;derta_y=fal e_Y*k1*r2;X=fale_X+derta_x;Y=fale_Y+derta_y;%设定内参,求像素坐标M1=[1/dX 0 u0;01/dY v0;0

01];pixel_xyz=M1*[X;Y;1];u=pixel_xyz(1);v=pixel_xyz(2);设畸变系数k=0、001，上述9个点的图像坐标分别为：[0、0040，0、0080]，[0、0080，0、0080]，[-0、0024，-0、0072]，[0、0030，-0、0060]，[-0、0120，0、0120]，[0、0086，0、0103]，[0、0113，0、0116]，[0，0]，[0、0092，0、0111]。像素坐标分别为[12

46、2，16

94、0]，[18

34、5，16

94、0]，

[3

50、0588，-5

41、3235]，[10

99、2，-3

64、8529]，[-1106、7，22

82、2]，[19

18、5，20

30、1]，[23

12、4，22

27、1]，[658，5

17、5]，[xx、5，21

46、5]。由数据可得，当畸变系数较小时，对坐标几乎没有影响，但畸变系数较大或图像处于透镜边缘时，对坐标还是有影响的。第四部分单目运动系统光点排序问题

一、问题描述

基于立体靶标单目运动系统中，靶标上有LED光点11个，沿光轴方向呈高低分布，其中6个低点、5个高点，且垂直投影到像面为一条直线，从左至右6个低点的编号分别为（1）（2）（3）、（9）（10）（11），5个高点为（4）（5）（6）（7）