CN101702056A - 基于立体图像对的立体图像显示方法 - Google Patents

基于立体图像对的立体图像显示方法 Download PDF

Info

Publication number
CN101702056A
CN101702056A CN200910185622A CN200910185622A CN101702056A CN 101702056 A CN101702056 A CN 101702056A CN 200910185622 A CN200910185622 A CN 200910185622A CN 200910185622 A CN200910185622 A CN 200910185622A CN 101702056 A CN101702056 A CN 101702056A
Authority
CN
China
Prior art keywords
image
stereo
point
picture
width
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN200910185622A
Other languages
English (en)
Other versions
CN101702056B (zh
Inventor
吴华夏
沈威
武维生
张涛
董戴
李园
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AVIC Huadong Photoelectric Co Ltd
Original Assignee
Anhui East China Institute of Optoelectronic Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Anhui East China Institute of Optoelectronic Technology filed Critical Anhui East China Institute of Optoelectronic Technology
Priority to CN2009101856225A priority Critical patent/CN101702056B/zh
Publication of CN101702056A publication Critical patent/CN101702056A/zh
Application granted granted Critical
Publication of CN101702056B publication Critical patent/CN101702056B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Processing (AREA)

Abstract

本发明属于光栅自由立体成像技术领域,涉及一种基于立体图像对的立体图像显示方法,本方法的具体步骤如下:步骤一,首先通过立体摄像机采集设备获取立体图像对;步骤二,对上述采集到的图像进行矫正和重采样,获得满足纯平移关系的两幅新图像;步骤三,对重采样的两幅图像进行密集匹配,建立点与点之间的对应关系;步骤四,由重采样的两幅图像插补出给定视点下的新图像;步骤五,对所得新图像进行反矫正,获得新视点下的立体图像;最后,将该图像在自由立体显示器上显示。本方法可以广泛应用于光栅LCD类型的自由立体显示设备中,通过本发明应用到具体实施里,可以得到很好的立体效果和立体观感。

Description

基于立体图像对的立体图像显示方法
【技术领域】
本发明涉及一种基于立体图像对的立体图像显示处理方法,让使用者不佩戴立体眼镜,裸眼就能在立体显示器上看到指定视点下的立体图像。
【背景技术】
人看物体之所以是三维的,是因为人有两个眼睛,并且两眼有一定的间距,物体在两眼视网膜上产生两幅有细微差别的图象,经人脑处理后合成为一幅三维图象。立体成像技术是利用人眼立体成像原理,人眼看物体时是从不同角度看到两个稍有差别的图像,大脑将这两个具有视差的图像合成后形成立体视觉。现有的立体成像方法主要有1、偏振镜法,2、滤色镜法,3、视差立体法,4、体视镜法。偏振镜法,也是最常用的方法,原理就是利用光的不同偏振角度,让两个镜片透过不同的光,将两幅有细微差别的图象分别投射到人眼中,给人以三维立体感。但是,此种佩戴立体眼镜的观赏方式常常让观赏者有不舒适感,特别是佩戴眼镜的观众,观赏立体影像较为困难。滤色镜法仅限于黑白立体图像。视差立体法,不需要戴立体眼镜,可以直接看到立体图像,但只适合彩色印刷,不适合在电子计算机上使用。体视镜法,需要通过两套相同的光学系统观察立体图像,一般应用于显微镜上。因此需要一种不必佩戴立体眼镜就能观察立体影像的方式,以满足上述观众的需求。
【发明内容】
本发明所要解决的技术问题在于采用基于立体图像对的立体图像显示方法,让观赏者不需要佩戴立体眼镜,裸眼即可观赏到高效率,高保真的立体图像画面。
本发明是通过以下技术方案解决上述技术问题的:
一种基于立体图像对的立体图像显示方法,包括如下步骤:
步骤一:首先通过立体摄像机对准同一对象以不同的角度所拍摄的两幅图像,即立体图像对;
步骤二:对所述采集到的立体图像对进行矫正和重采样,通过对所述立体图像对进行极线校正,获得满足纯平移关系的两幅新图像,以实现双目图像的致密匹配;
步骤三:对所述经过重采样的立体图像对进行密集匹配,建立点与点之间的对应关系;
步骤四:由所述经过重采样的立体图像对插补出给定视点下的新图像;
步骤五,对所述新图像进行反矫正,获得新视点下的体式图像;最后,将该图像在自由立体显示器上显示。
现有的立体摄像机具有两个间距约为人眼间距的镜头,在拍摄时分别记录下由于间距产生的两幅有细微差别的图象,模拟人眼,经过计算机合成后也就成为三维图象。所述立体图像对可以利用两部摄像获取设备,如数码相机或者计算机模拟的相机,在同一水平面,按照直线或者弧线路径,所述获取装置镜头要对准同一对象以不同的角度所拍摄的两幅图像。
本发明的优点在于:通过本发明的应用,可以得到很好的立体效果和立体观感,该方法可以广泛应用于光栅LCD类型的自由立体显示设备中,对于大力推动立体显像设备的产业化,具有重要的实用意义。
【附图说明】
下面参照附图结合实施例对本发明作进一步的描述。
图1是本发明的立体图像对体式显示处理方法流程图,
图2是本发明的立体图像对生成图像的图像插补示意图,
图3是本发明的密集匹配中卷积运算示意图,
图4是本发明的Is在m产生空洞,
图5是本发明的Is在m发生重叠。
【具体实施方式】
现将有关发明的详细内容及技术说明,配合附图说明如下:
参见图1、图2,图2是本发明的立体图像对生成体式图像的坐标系示意图。在图中我们对视点做了如下定义:设摄像机坐标系关于世界坐标系的旋转矩阵为R,平移变量为T,则称摄像机位于视点(R|T)。设11,12分别为某场景在两个不同视点(I|O)和(R|T)下的图像,其中I为单位矩阵(以第一个摄像机所在的坐标系为参考坐标系,R和T分别为第二个摄像机相对于第一个摄像机的旋转和平移)。如图1所示,由原图像I1,I2,求摄像机在新视点cs(Rs,Ts)下的图像Is(λ,β)的过程就称为图像插补。
1>对原图像对进行矫正和重采样,使得校正后的图像对满足纯平移关系:它们的视差仅在平行于X轴方向上,这样进行匹配时仅需要沿着一条水平线找对应点就可以了。进行矫正和重采样的过程如下:
1.1>首先将图像坐标原点移到图像的中心:设图像中心坐标为x0,其齐次坐标表示为(u,v,1)T,可构造变换
Figure G2009101856225D0000031
图像I2的对极点e2经此变换后为e′2=Te2
1.2>接着构造旋转矩阵R,使对极点e′2被旋转到x轴上的点(f,0,1)T。令e′2=(e′2x,e′2y,1)T
Figure G2009101856225D0000032
则R可选择为
Figure G2009101856225D0000033
对极点e′2经此变换后为
Figure G2009101856225D0000034
最后,构造变换矩阵G,映射对极点e″2到无穷远点,G可选择为
G = 1 0 0 0 1 0 - 1 / f 0 1
且将对极点e″2变换到无穷远点e′″2=Ge″2=(f,0,0)T
综合以上各步可知,合成变换H2=GRT把原图像的对极点e2变换到无穷远点(f,0,0)T
1.3>接下来就是把一个映射施加到第一幅图像上,使得两幅图像的对极线相互匹配。
设原图像为I1,I2,需要对两幅图像进行重采样,分别施加两个射影变换H1,H2作用于I1,I2。令I1,I2是具有基本矩阵F=[e′]xM的两幅图像(e′是F的左零矢量,M=[e′]xF+e′vT,其中v为任意三维列矢量),其中H1=(1+H2e2aT)H2M。
利用上述方法可以求得变换H1,分别用H1,H2矫正I1,I2,使得重采样后的图像对应的对极线平行于x轴,而对应点的y坐标相等。
利用上述的射影变换H1和H2矫正两幅原图像,并分别进行重采样就可获得两幅新图像。
2>密集匹配,对两幅图像进行了矫正和重采样,下面就是对这两幅图像进行密集匹配步骤如下:
2.1>对第一幅图像扫描线上的点进行一阶差分滤波,检测出角点,这里采用如下一阶高斯差分滤波器对图像进行离散卷积运算:
h ( 1 ) ( n ) = - 2 cn e - n 2 2 σ 2
其中c是归一化常数,它可由公式
Figure G2009101856225D0000042
求得,σ通常取常数1,在实际使用中,h(n)可用一个长度为2N+1的窗口内的值表示,N取2的情况见图3所示。
2.2>由第一幅图像上某一扫描线上检测到的角点,利用灰度相关匹配算法[6]在第二幅图像的同一扫描线上寻找对应的角点。
2.3>用仿射变换建立较为准确的点与点的对应关系。设(u1,v1),(u′1,v′1)与(u2,v2),(u′2,v′2)为两对对应的角点,其中(u1,v1),(u2,v2)是I1上的点,(u′1,v′1),(u′2,v′2)是I2上的点,它们在同一扫描线上,v1=v1′=v2=v2′。设(u,v)为(u1,v1),(u2,v2)之间的任一点,u1≤u≤u2,于是有
Figure G2009101856225D0000043
其中u′=u+δ,v′=v,点(u′,v′)在(u′1,v′1),(u′2,v′2)之间,它与(u,v)是一对对应点。这样就完成了点与点的对应。
3>图像插补,该方法的步骤如下:
3.1>确定插补图像像素位置
为确定图像像素位置,首先构造一个如下形式的Warp函数:
W(U1,s)=sU1+(1-s)f(U1)
其中s∈[0,1],U1表示图像I1上点,f表示图像I1,I2之间的对应关系,W表示插补图像上的点。为了使图像是新视点下的像,进一步对Warp函数进行改进,即令U1=(u1,v1)T∈I1和U2=(u2,v2)T∈I2为一对对应点,Us=(us,vs)T表示对应的插补图像点,那么此时Warp函数的形式为:
Us=sU1+(1-s)U2
将对应点坐标代入上式,得
u s = u 2 + ( u 1 - u s ) s v s = v 2 + ( v 1 - v 2 ) s
3.2>确定插补图像像素灰度
假定两个视点下所对应的原图像I1,I2光照条件不变,从而两幅原图像I1,I2中的对应点的灰度值差别不大,因此插补图像的灰度可以用下面的算法近似估计。令
Figure G2009101856225D0000052
设对应点U1∈I1,U′1=f(U1)∈I2的灰度值分别为r(U1),r(U′1),
Figure G2009101856225D0000053
上的对应点W(U1,s)的灰度值为r(W(U1,s))=sr(U1)+(1-s)r(U′1)综合上述算式得出原图像I1,I2的插补图像为
Figure G2009101856225D0000054
4>反矫正得到实际的体式图像
通过上述方法得到了插补图像,但这样的图像并不是新视点下的真实图像,因此需要通过变换矩阵Hs=sH1+(1-s)H2对该图像进行反矫正处理,最终得到新视点下的体式图像。
用上述方法我们可以获得新视点下的体式的插补图像,但是在插补过程中无法避免会出现空洞现像(参见图4)和重叠现像(参见图5)。为了得到较好的插补结果,空洞中的像素灰度值可以利用邻近插值法获得,重叠部分的像素灰度可取该像素上多个灰度的平均值。
目前,通过立体图像配合光栅技术显示方式是比较好的技术之一,随着科技的进步,薄、轻的平板显示器已经代替CRT显示器,使得我们能将立体显示技术应用于液晶显示器成为可能。通过液晶自由立体显示器,我们不仅能看到一般立体静画影像,还能看到交互式的实时立体影像。通过结合立体图像对的体式显示技术,可以将大量实时采集的图像转化成任意视点的体式图像,给立体显示器更大的使用空间。

Claims (6)

1.一种基于立体图像对的立体图像显示方法,其特征是:所述方法包括如下步骤:
步骤一:首先通过立体摄像机对准同一对象以不同的角度所拍摄的两幅图像,即立体图像对;
步骤二:对所述采集到的立体图像对进行矫正和重采样,通过对所述立体图像对进行极线校正,获得满足纯平移关系的两幅新图像;
步骤三:对所述经过重采样的立体图像对进行密集匹配,建立点与点之间的对应关系;
步骤四:由所述经过重采样的立体图像对插补出给定视点下的新图像;
步骤五,对所述新图像进行反矫正,获得新视点下的体式图像;
最后,将该图像在自由立体显示器上显示。
2.根据权利要求1所述的一种基于立体图像对的立体图像显示方法其特征是:所述方法步骤二中进行矫正和重采样的过程如下:
1>首先将图像坐标原点移到图像的中心:设图像中心坐标为x0,其齐次坐标表示为(u,v,1)T,可构造变换
Figure F2009101856225C0000011
图像I2的对极点e2经此变换后为e′2=Te2
2>接着构造旋转矩阵R,使对极点e′2被旋转到x轴上的点(f,0,1)T。令e′2=(e′2x,e′2y,1)T
Figure F2009101856225C0000012
则R可选择为
Figure F2009101856225C0000013
对极点e′2经此变换后为
Figure F2009101856225C0000014
最后,构造变换矩阵G,映射对极点e″2到无穷远点,G可选择为
G = 1 0 0 0 1 0 - 1 / f 0 1
且将对极点e″2变换到无穷远点e″′2=Ge″2=(f,0,0)T
通过以上步骤,合成变换H2=GRT把原图像的对极点e2变换到无穷远点(f,0,0)T
3.根据权利要求1所述的一种基于立体图像对的立体图像显示方法,其特征是:所述方法步骤二中所述立体图像对的密集匹配步骤如下:
设原图像为I1,I2,需要对两幅图像进行重采样,分别施加两个射影变换H1,H2作用于I1,I2。令I1,I2是具有基本矩阵F=[e′]xM的两幅图像(e′是F的左零矢量,M=[e′]xF+e′vT,其中v为任意三维列矢量),其中
H1=(1+H2e2aT)H2M;
利用上述方法可以求得变换H1,分别用H1,H2矫正I1,I2,使得重采样后的图像对应的对极线平行于x轴,而对应点的y坐标相等;
利用上述的射影变换H1和H2矫正两幅原图像,并分别进行重采样就可获得两幅新图像。
4.根据权利要求1所述的一种基于立体图像对的立体图像显示方法其特征是:所述方法步骤四对所述经过矫正和重采样的立体图像对进行密集匹配步骤如下:
1>对第一幅图像扫描线上的点进行一阶差分滤波,检测出角点,这里采用如下一阶高斯差分滤波器对图像进行离散卷积运算:
h ( 1 ) ( n ) = - 2 cne - n 2 2 σ 2
其中c是归一化常数,它可由公式
Figure F2009101856225C0000022
求得,σ通常取常数1,其中2N-3<n<2N+3,N为自然数。
2>由第一幅图像上某一扫描线上检测到的角点,利用灰度相关匹配算法[6]在第二幅图像的同一扫描线上寻找对应的角点。
3>用仿射变换建立较为准确的点与点的对应关系。设(u1,v1),(u′1,v′1)与(u2,v2),(u′2,v′2)为两对对应的角点,其中(u1,v1),(u2,v2)是I1上的点,(u′1,v′1),(u′2,v′2)是I2上的点,它们在同一扫描线上,v1=v′1=v2=v′2。设(u,v)为(u1,v1),(u2,v2)之间的任一点,u1≤u≤u2,于是有
δ = ( u 1 ′ - u 1 ) + ( u 2 ′ - u 2 ) - ( u 1 ′ - u 1 ) u 2 - u 1 ( u - u 1 )
其中u′=u+δ,v′=v,点(u′,v′)在(u′1,v′1),(u′2,v′2)之间,它与(u,v)是一对对应点。这样就完成了点与点的对应。
5.根据权利要求1所述的一种基于立体图像对的立体图像显示方法,其特征是:所述方法步骤四中所述立体图像对的插补步骤如下:
1>确定插补图像像素位置
为确定图像像素位置,首先构造一个如下形式的Warp函数:
W(U1,s)=sU1+(1-s)f(U1)
其中s∈[0,1],U1表示图像I1上点,f表示图像I1,I2之间的对应关系,W表示插补图像上的点;
为了使图像是新视点下的像,进一步对Warp函数进行改进,即令U1=(u1,v1)T∈I1和U2=(u2,v2)T∈I2为一对对应点,Us=(us,vs)T表示对应的插补图像点,那么此时Warp函数的形式为:
Us=sU1+(1-s)U2
将对应点坐标代入上式,得
u s = u 2 + ( u 1 - u s ) s v s = v 2 + ( v 1 - v 2 ) s
2>确定插补图像像素灰度:
假定两个视点下所对应的原图像I1,I2光照条件不变,从而两幅原图像I1,I2中的对应点的灰度值差别不大,因此插补图像的灰度可以用下面的算法近似估计。令
Figure F2009101856225C0000041
设对应点U1∈I1,U′1=f(U1)∈I2的灰度值分别为r(U1),r(U′1),
Figure F2009101856225C0000042
上的对应点W(U1,s)的灰度值为r(W(U1,s))=sr(U1)+(1-s)r(U′1)
综合上述算式得出原图像I1,I2的插补图像为
Figure F2009101856225C0000043
6.根据权利要求5所述的一种基于立体图像对的立体图像显示方法,其特征是:所述像素的灰度值用邻近插值法获得,重叠部分的像素灰度取该像素上多个灰度的平均值。
CN2009101856225A 2009-11-25 2009-11-25 基于立体图像对的立体图像显示方法 Active CN101702056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101856225A CN101702056B (zh) 2009-11-25 2009-11-25 基于立体图像对的立体图像显示方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101856225A CN101702056B (zh) 2009-11-25 2009-11-25 基于立体图像对的立体图像显示方法

Publications (2)

Publication Number Publication Date
CN101702056A true CN101702056A (zh) 2010-05-05
CN101702056B CN101702056B (zh) 2011-07-20

Family

ID=42156976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101856225A Active CN101702056B (zh) 2009-11-25 2009-11-25 基于立体图像对的立体图像显示方法

Country Status (1)

Country Link
CN (1) CN101702056B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608765A (zh) * 2011-01-20 2012-07-25 承景科技股份有限公司 立体图像显示方法
CN102801991A (zh) * 2011-05-27 2012-11-28 翔德电子科技(深圳)有限公司 获取3d影像的方法
CN102842118A (zh) * 2012-07-17 2012-12-26 刘怡光 一种鲁棒的立体像对矫正新方法
CN104519330A (zh) * 2013-09-26 2015-04-15 深圳市云立方信息科技有限公司 一种对3d视频识别的方法及装置
CN106023230A (zh) * 2016-06-02 2016-10-12 辽宁工程技术大学 一种适合变形图像的稠密匹配方法
TWI618394B (zh) * 2012-09-26 2018-03-11 Lg伊諾特股份有限公司 利用影像特性之立體匹配裝置
CN109313809A (zh) * 2017-12-26 2019-02-05 深圳配天智能技术研究院有限公司 一种图像匹配方法、装置及存储介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6331295A (ja) * 1986-07-25 1988-02-09 Canon Inc 立体画像信号処理装置
CN1231456A (zh) * 1999-01-23 1999-10-13 中国科学院固体物理研究所 用计算机图对合成方法制作立体图片
WO2008080281A1 (fr) * 2006-12-28 2008-07-10 Nuctech Company Limited Procédé et système d'imagerie radiologique pour balayage à double observation

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608765A (zh) * 2011-01-20 2012-07-25 承景科技股份有限公司 立体图像显示方法
CN102608765B (zh) * 2011-01-20 2015-06-03 承景科技股份有限公司 立体图像显示方法
CN102801991A (zh) * 2011-05-27 2012-11-28 翔德电子科技(深圳)有限公司 获取3d影像的方法
CN102842118A (zh) * 2012-07-17 2012-12-26 刘怡光 一种鲁棒的立体像对矫正新方法
TWI618394B (zh) * 2012-09-26 2018-03-11 Lg伊諾特股份有限公司 利用影像特性之立體匹配裝置
CN104519330A (zh) * 2013-09-26 2015-04-15 深圳市云立方信息科技有限公司 一种对3d视频识别的方法及装置
CN104519330B (zh) * 2013-09-26 2018-04-06 深圳市云立方信息科技有限公司 一种对3d视频识别的方法及装置
CN106023230A (zh) * 2016-06-02 2016-10-12 辽宁工程技术大学 一种适合变形图像的稠密匹配方法
CN106023230B (zh) * 2016-06-02 2018-07-24 辽宁工程技术大学 一种适合变形图像的稠密匹配方法
CN109313809A (zh) * 2017-12-26 2019-02-05 深圳配天智能技术研究院有限公司 一种图像匹配方法、装置及存储介质
CN109313809B (zh) * 2017-12-26 2022-05-31 深圳配天智能技术研究院有限公司 一种图像匹配方法、装置及存储介质

Also Published As

Publication number Publication date
CN101702056B (zh) 2011-07-20

Similar Documents

Publication Publication Date Title
CN101702056B (zh) 基于立体图像对的立体图像显示方法
CN107105213B (zh) 立体图像显示装置
CN103974055B (zh) 3d照片生成系统及方法
Stern et al. Perceivable light fields: Matching the requirements between the human visual system and autostereoscopic 3-D displays
CN108174182A (zh) 三维跟踪式裸眼立体显示视区调整方法及显示系统
US20160295194A1 (en) Stereoscopic vision system generatng stereoscopic images with a monoscopic endoscope and an external adapter lens and method using the same to generate stereoscopic images
CN107277495A (zh) 一种基于视频透视的智能眼镜系统及其透视方法
CN105611278A (zh) 防裸眼3d观看眩晕感的图像处理方法及系统和显示设备
CN104247411B (zh) 校正由立体显示的调节效应引起的失真误差的方法和装置
CN104618706A (zh) 一种分时实现多人多角度全息立体显示的移动终端及方法
CN105791800A (zh) 立体显示系统及立体显示方法
CN102547350B (zh) 一种基于梯度光流算法合成虚拟视点方法及立体显示装置
CN105809654A (zh) 目标对象跟踪方法、装置和立体显示设备及方法
CN112929636B (zh) 3d显示设备、3d图像显示方法
US20140347444A1 (en) Method and device for stereo base extension of stereoscopic images and image sequences
US20220060680A1 (en) Head mounted display apparatus
Hwang et al. Instability of the perceived world while watching 3D stereoscopic imagery: a likely source of motion sickness symptoms
CN109428987A (zh) 一种头戴式全景360度立体摄像装置及摄像处理方法
CN105657401A (zh) 一种裸眼3d显示方法、系统及裸眼3d显示装置
CN204377059U (zh) 基于软镜的立体显示系统
CN106991715A (zh) 基于光场采集的光栅棱柱三维显示渲染方法
WO2007017834A3 (en) Disparity value generator
CN108259888A (zh) 立体显示效果的测试方法及系统
CN110602478A (zh) 一种三维显示装置及系统
JP2017098596A (ja) 画像生成方法及び画像生成装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: AVIC HUADONG PHOTOELECTRIC CO., LTD.

Free format text: FORMER OWNER: HUADONG PHOTOELECTRIC TECHNIQUE INSTITUTE OF ANHUI PROVINCE

Effective date: 20110317

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20110317

Address after: 241002 Yijiang District, Wuhu, China Hi tech Development Zone, science and Technology Park

Applicant after: Avic Huadong Photoelectric Co., Ltd.

Address before: 241002 Yijiang District, Wuhu, China Hi tech Development Zone, science and Technology Park

Applicant before: Huadong Photoelectric Technique Institute of Anhui Province

C14 Grant of patent or utility model
GR01 Patent grant