一种残膜回收机防缠绕挑膜装置的制 一种秧草收获机用电力驱动行走机构

一种基于时间同步的车载侧后视图像生成方法与流程

2023-02-04 13:36:17 来源:中国专利 TAG:


1.本发明涉及汽车领域,尤其涉及一种基于时间同步的车载侧后视图像生成方法。


背景技术:

2.现有侧后视应用如雨后春笋般迅速发展,普遍采用车载摄像头直接替代后视镜,将各个位置摄像头采集到的鱼眼图像显示在屏幕的对应区域。但是在直观感受上,各个摄像头的视场图像是独立不连贯的,驾驶员观看屏幕的主观感受较差。此外由于摄像头是电子产品,具备一定的失效性,外置摄像头仍然存在被灰尘和泥浆污染的情况,使得画面某些区域模糊或者遮挡,造成驾驶感受的降低甚至出现意外。但现阶段,摄像头被污染仍然是依赖于人眼进行判断,尚未有摄像头被污染自动判断的方法,这导致摄像头出现污染影响成像时需要人判断校正,如果驾驶员不及时发现处理,依赖被污染摄像头的判断可能会发生安全事故。因此,现有技术仍然存在较多的缺陷,需要进行改进。


技术实现要素:

3.基于现有技术中存在的缺陷,本发明提供一种基于时间同步的车载侧后视图像生成方法,至少包括:
4.步骤s1,对多个车载以太网摄像头进行图像采集的时间进行时间同步,并对采集图像进行畸变矫正;
5.步骤s2,获取冗余车载以太网摄像头拍摄的图像,在进行图像融合后获得融合图像,根据融合图像对应原始畸变图像判断选择未被污染的车载以太网摄像头作为工作摄像头;
6.步骤s3,获取不同安装位置的工作摄像头拍摄图像进行畸变矫正中的差异矫正单应性矩阵h
da
posk、矫正融合单应性矩阵h
af
pos,并对工作摄像头拍摄的画面进行图像映射,得到安装位置工作摄像头对应的融合平面的图像;
7.步骤s4,融合图像上的不同安装位置工作摄像头对应图像互相覆盖区域的各个像素值根据信息权重进行覆盖区域的图像融合。
8.一种基于时间同步的车载侧后视图像生成方法,进一步地,车载以太网摄像头通过以太网eavb协议传输图像内容和控制信息,内置rtc时钟,智能主机通过车载以太网总线每隔一段时间修正每个车载以太网摄像头的rtc时钟。
9.一种基于时间同步的车载侧后视图像生成方法,进一步地,车载以太网摄像头进行图像采集的时间进行时间同步具体包括:
10.车载以太网摄像头计算各自与主处理器的链路传播延迟p
delay

11.主处理器不定期对车载以太网摄像头发送主处理器本地时间报文时,车载以太网摄像头解析出报文中记录的主处理器本地时间戳tm,再加上传播链路延迟p
delay
,获取收到报文后此时主处理器的本地rtc时间戳tc;
12.同时,车载太网摄像头将各自本地的rtc时钟时间戳修改为tc。
13.一种基于时间同步的车载侧后视图像生成方法,进一步地,畸变矫正至少包括第一次畸变矫正,其中,第一次畸变矫正至少包括:通过冗余车载以太网摄像头获取垂直于地面的预设靶图成像后的畸变图像,利用光学畸变参数映射表对畸变图像进行第一畸变矫正后获取第一畸变矫正图像。
14.一种基于时间同步的车载侧后视图像生成方法,进一步地,畸变矫正至少包括第二次畸变矫正,其中,第二次畸变矫正至少包括:建立差异矫正平面坐标系,分别获取预设靶图的特征点在第一畸变矫正图像的对应的第一坐标点以及在差异矫正平面对应的第二坐标点,利用透视变换,求解单应性矩阵h
da
pos,利用单应性矩阵h
da
pos进行反变化对第一畸变的图像进行第二次畸变矫正;
15.获取单应性矩阵h
da
pos并利用反透视变换将建立差异矫正平面对应的所有第一坐标点在第一畸变矫正图像对应第二坐标点,获取第二坐标点在第一畸变矫正图像对应的像素值,并将像素值赋值给第一坐标点。
16.一种基于时间同步的车载侧后视图像生成方法,进一步地,在步骤s2中,针对安装在同一位置的冗余车载以太网摄像头,将同步采集后的原始摄像头的畸变图像分为若干区域,每个区域独立进行清晰度得分计算,一个摄像头得到一组各个区域的清晰度得分df(f)
vh
,每个冗余摄像头k各自对应一组清晰度得分dfk(f)
vh

17.一种基于时间同步的车载侧后视图像生成方法,进一步地,在步骤s2中,针对安装在同一位置的冗余车载以太网摄像头,将矫正后的车载以太网摄像头画面映射到实际显示应用的融合画面,得到各个车载以太网摄像头安装位置差异矫正后对应融合画面的矫正融合单应性矩阵h
af
pos,并设定融合画面的区域信息权重函数vld
(x,y)
,并根据该安装位置各个冗余车载以太网摄像头k清晰度系数dfk(f)
vh
及其对应的差异矫正单应性矩阵h
da
posk得到信息权重infk(f),冗余摄像头中信息权重最高的摄像头为工作摄像头,其余为待机摄像头。
18.一种基于时间同步的车载侧后视图像生成方法,进一步地,覆盖区域包括:右侧和后侧工作摄像头覆盖区域融合图像q
rr(x,y)
,左侧和后侧工作摄像头覆盖区域融合图像q
rl(x,y)

19.摄像头覆盖区域融合图像包括相应区域的工作摄像头在最终融合平面上的图像和信息权重的加权平均值,融合计算方法包括:
[0020][0021][0022]
其中,左侧、后侧、右侧工作摄像头的信息权重leftinf(f)
(x,y)
,rearinf(f)
(x,y)
,rightinf(f)
(x,y)

[0023]
right(q
x
、qy)、rear(q
x
、qy)、left(q
x
、qy)分别为右侧、后侧、左侧工作摄像头在最终融合平面上的图像。
[0024]
一种基于时间同步的车载侧后视图像生成方法,进一步地,左侧、后侧、右侧工作摄像头的信息权重leftinf(f)
(x,y)
,rearinf(f)
(x,y)
,rightinf(f)
(x,y)
的计算公式如下:
[0025]
leftinf(f)
(x,y)
=leftdfk(f)
(x,y)
×
vld
(x,y)
[0026]
rearinf(f)
(x,y)
=reardfk(f)
(x,y)
×
vld
(x,y)
[0027]
rightinf(f)
(x,y)
=rightdfk(f)
(x,y)
×
vld
(x,y)
[0028]
其中reardfk(f)
(x,y)
为后侧工作摄像头的清晰度得分;leftdfk(f)
(x,y)
为左侧工作摄像头的信息权重,rightdfk(f)
(x,y)
为右侧工作摄像头的清晰度得分,vld
(x,y)
为信息权重函数。
[0029]
一种基于时间同步的车载侧后视图像生成方法,进一步地,清晰度得分的获取具体包括:将获取的原始畸变图像按边长为b1划分为vnum
×
hnum个小块imb(f)
vh
,统计每个小块内部的清晰度分值var(cf)
vh
,连续统计n帧图像的每一个小块内部的清晰度分数进行累加,获得清晰度得分df(f)
vh

[0030][0031]
其中,f为当前图像帧序号,resh、resv分别为图像水平分辨率和垂直分辨率,b1为小块的边长,v∈[1,vnum],h∈[1,hnum],其中小块的边长,v∈[1,vnum],h∈[1,hnum],其中
[0032]
一种基于时间同步的车载侧后视图像生成方法,进一步地,清晰度分值var(cf)
vh
计算方法包括:
[0033][0034]
其中,imb(f)
vh
表示垂直序号为v,水平序号为h的小块,f为当前图像的帧序号,bl为小块的边长,i和j分别表示小块图像内的像素坐标,i∈[1,b1],j∈[1,bl];
[0035]
设定清晰度分值阈值thre,thre根据bl大小和统计帧数n设定,thre
min
范围为8~12;
[0036]
清晰度分值阈值thre:
[0037]
thre=(thre
min
×
bl)2×n[0038]
将各个小块统计的清晰度分值与清晰度分值阈值比较,大于等于清晰度分值阈值,认为该区域图像清晰度满足要求。
[0039]
一种基于时间同步的车载侧后视图像生成方法,进一步地,融合图像的获取至少包括:
[0040]
步骤s100,获取预设靶图矩形roi区域4个顶点在融合矫正图像中的像素坐标,按照顺时针方向分别为(resq
x1
,resq
y1
)、(resq
xk
,resq
yk
)、(resq
x(t-k 1)
,resq
y(t-k 1)
)、(resq
xt
,resq
yt
);
[0041]
步骤s101,设观察点在三维坐标系下的o点,o点在水平面g的投影为g,g点为三维坐标系的原点,g点位于靶图平面的后摄像头成像的垂直中心线上,以g点为原点建立靶图坐标系;
[0042]
步骤s102,根据角点顺序按照顺时针方向分别为c1、c2、c3、c4,求解靶图坐标系下右冗余摄像头、左冗余摄像头、后冗余摄像头拍摄靶图成像后的roi区域以观察点o观测的融合图像的像素坐标namefc(i)x,namefc(i)y,name=rgt,rer,lft;i=1,2,3,4;
[0043]
步骤s103,设左右摄像头安装位置距离其对应靶图的像素距离为disz,将右冗余摄像头、左冗余摄像头、后冗余摄像头成像后的roi区域的以观察点o观测的融合图像的像素坐标转化为roi区域的顶点坐标namec(i)x,namec(i)y,其中,name=right,rear,left;i=1,2,3,4。
[0044]
一种基于时间同步的车载侧后视图像生成方法,进一步地,根据融合图像获取对应原始畸变图像至少包括:
[0045]
步骤s201,根据观测的融合图像的像素坐标和靶图矩形roi区域4个顶点的像素坐标namec(i)x,namec(i)y,其中,name=right,rear,left;i=1,2,3,4,利用透视变换进行求解h
af
pos,获取多摄像头融合图像的单应性矩阵h
af
pos和多摄像头的差异矫正图像,其中,融合图像的像素坐标记为namefc(i)x,namefc(i)y,其中,name=rgt,rer,lft;i=1,2,3,4.;
[0046]
步骤s202,重复步骤s201,依次获得不同位置的右、后、左冗余摄像头对应的单应性矩阵h
af
pos和多摄像头的差异矫正图像;
[0047]
步骤s203,将摄像头最终融合平面图像(q
x
、qy)、h
af
pos利用透视变换,获取各个摄像头的最终融合图像各像素点坐标对应的差异矫正图像坐标(r
x
,ry);
[0048]
步骤s204,差异矫正图像坐标(r
x
,ry)和对应摄像头的单应性矩阵h
da
pos利用透视变换,获取各个摄像头的差异矫正图像坐标对应的畸变矫正图像坐标(resp
x
,respy);
[0049]
步骤s205,根据鱼眼相机的畸变矫正原理,获取畸变矫正图像坐标(resp
x
,respy)相对应的原始畸变图像的坐标(resp

x
,resp
′y);
[0050]
步骤s206,重复步骤s203至s205,获取不同摄像头对应的原始畸变图像坐标(resp

xk
,resp

yk
),其中,k为摄像头的编号。
[0051]
一种基于时间同步的车载侧后视图像生成方法,进一步地,信息权重infk(f)计算方法包括:
[0052][0053]
其中,k表示冗余摄像头的编号,resqh、resqv分别为融合后图像水平分辨率和垂直分辨率,x,y表示融合图像像素坐标,f表示帧的序列号;
[0054]
dfk(f)
(x,y)
表示冗余摄像头k的清晰度得分,vld
(x,y)
表示融合画面的区域信息权重函数。
[0055]
一种基于时间同步的车载侧后视图像生成方法,进一步地,区域信息权重函数vld
(x,y)
获取包括:
[0056]
vld
(x,y)
=vld
(x)
×
vld
(y)
[0057]
其中,垂直信息权重函数vld
(y)
,水平信息权重函数vld
(x)

[0058][0059][0060]
其中,resqh为最终融合图像的像素水平分辨率;
[0061]
其中,ofst_h1和ofst_h2,low_val均为人为设定阈值,low_val为主观认为的最低权重,ofst_h1是主观认定的最低权重区域,ythre为根据经验设定的阈值;
[0062]
ofst_vt表示最上边像素到右侧相机拍摄到靶图right camera roi的点(resp
x1
,resp
y1
)的垂直像素距离。
[0063]
一种基于时间同步的车载侧后视图像生成方法,优选的,其中,第一次畸变矫正包括:
[0064]
步骤s1,获取车载以太网摄像头的光学畸变参数,光学畸变参数包括摄像头镜头指定cmos芯片条件下的实验模组进行实际实验环境图像采集的离散视场角度θ与其对应像高γ的映射表;
[0065]
步骤s2,根据小孔成像原理,建立实际拍摄平面u各个像素点对应到成像面i上鱼眼图像的各个点的空间坐标关系,建立离散视场角度θ与物高λ的第一映射函数;
[0066]
λ=tan(θ)
[0067]
步骤s3,根据光学畸变参数映射表,构建像高γ与离散视场角度θ第二映射函数θ=g(γ):
[0068]
步骤s4,根据成像几何关系,求解实际拍摄平面u与对应的矫正成像面、畸变成像面之间的矫正缩放因子su和畸变缩放因子si;
[0069]
步骤s5:利用第一映射函数、第二映射函数以及矫正缩放因子su和畸变缩放因子si,建立矫正后的无畸变图像的像素级坐标resp
x
,respy与畸变图像的像素级坐标resp

x
,resp
′y的对应第三映射函数,并将resp

x
,resp
′y坐标位置的像素值赋值给resp
x
,respy的表示的坐标位置,从而获取畸变矫正后的图像。
[0070]
有益效果:
[0071]
1.本发明的技术方案通过一种基于时间同步的侧后视图像生成系统,通过采用以太网摄像头,组合了以太网同步技术和cmos的被动触发采集,实现了各个摄像头在同一时刻的采集,使得各个摄像头采集的图像真正的具有了对三维世界的同时摄影,从根本上避免了行车过程中由于各个摄像头采集不同时而拍摄的实际物体的位置发生改变导致的图像变化,从而解决了对多摄像头图像融合带来的融合边界图像不可控误差的根本原因。此外通过冗余摄像头同步工作,通过算法实时选择画面信息量更大的摄像头作为工作摄像头,从而使得摄像头能尽可能多的提供驾驶所需要的图像信息,提升了驾驶的安全性,基于平面贴图的多图像融合技术更迎合了人眼的主观观察感受。
[0072]
2.本发明的技术方案通过车载以太网摄像头,再加上eavb协议以及对应的时间同
步修正,能够保持冗余摄像头采集的图像在同一时刻,有助于融合图像的拼接,提高成像效果。
[0073]
3.本发明提供的技术方案中,在判断摄像头是否被污染的情况,通过引入信息权重,并且创造性的定义了信息权重等于清晰度得分和融合画面区域信息的权重函数。在计算信息权重的过程中,并未将摄像头所有的区域都进行考虑,通过预设靶图,建立靶图坐标系,通过图像融合后确定最终的有效区域,计算有效区域的清晰度得分和区域信息权重,使其计算结果更准确,避免不参加融合摄像头形成图像的部分被考虑进来。即对形成安全不产生影响的图像不考虑进来。
[0074]
4.本发明提供的技术方案中,创造性提出的清晰度得分的计算方法,并图像区域进行分块并计算每小块的像素差异,然后连续计算多帧图像的差异,设定阈值进行比较,使得计算结果更科学,准确。
附图说明
[0075]
以下附图仅对本发明做示意性说明和解释,并不限定本发明的范围。
[0076]
图1为本发明一实施例中车载以太网摄像头与主处理的时间链路示意图;
[0077]
图2为本发明一实施例中lens配型指定cmos的情况下的实验模组进行实际实验环境图像采集得出的离散视场角度与其对应像高的映射表。
[0078]
图3为本发明一实施例中鱼眼摄像头的成像模型原理示意图。
[0079]
图4为本发明一实施例中鱼眼摄像头拍摄的未畸变矫正的棋盘格图片。
[0080]
图5为本发明一实施例中经过第一次畸变矫正的棋盘格图片。
[0081]
图6为本发明一实施例中经过第二次畸变矫正的图片。
[0082]
图7为本发明一实施例中连续统计n帧图像的每一个小块内部的水平垂直的清晰度分值示意图。
[0083]
图8为本发明一实施例中摄像头表面被严重污染以及对污染区域进行染色的示意图,图8a表示表面被严重污染原图像,图8b表示原图经过计算处理并染色的图像。
[0084]
图9为本发明一实施例中摄像头表面轻微污染以及对污染区域进行染色的示意图,图9a表示表面轻微污染原图像,图9b表示原图经过计算处理并染色的图像。
[0085]
图10为本发明一实施例中靶图t矩形roi区域4个顶点在畸变矫正图像中的像素坐标示意图。
[0086]
图11为本发明一实施例中左、后、右摄像头形成的靶图摆放在同一平面并垂于水平面的示意图。
[0087]
图12为本发明一实施例中以观察点为视角对靶图进行观察的示意图。
[0088]
图13为本发明一实施例中左、后、右摄像头拍摄画面进行融合后覆盖区域显示的示意图。
[0089]
图14为本发明一实施例中工作摄像头安装在预设安装位置并朝向指定视场后拍摄的图像经过在融合平面形成的图像示意图。
[0090]
图15为本发明一实施例中将图像中信息权重显示1染色为纯白,显示为0染色为纯黑形成的示意图。
[0091]
图16为本发明一实施例中各个区域冗余摄像头(未污染和污染)的图像。
[0092]
图17为本发明一实施例中各个区域冗余摄像头最终融合图像加上贴图的最终图像。
具体实施方式
[0093]
为了对本文的技术特征、目的和效果有更加清楚的理解,现对照附图说明本发明的具体实施方式,在各图中相同的标号表示相同的部分。为使图面简洁,各图中的示意性地表示出了与本发明相关部分,而并不代表其作为产品的实际结构。另外,为使图面简洁便于理解,在有些图中具有相同结构或功能的部件,仅示意性地绘示了其中的一个,或仅标出了其中的一个。
[0094]
关于控制系统,功能模块、应用程序(app)本领域技术人员熟知的是,其可以采用任何适当的形式,既可以是硬件也可以是软件,既可以是离散设置的多个功能模块,也可以是集成到一个硬件上的多个功能单元。作为最简单的形式,所述控制系统可以是控制器,例如组合逻辑控制器、微程序控制器等,只要能够实现本技术描述的操作即可。当然,控制系统也可以作为不同的模块集成到一个物理设备上,这些都不偏离本发明的基本原理和保护范围。
[0095]
本发明中“连接”,即可包括直接连接、也可以包括间接连接、通信连接、电连接,特别说明除外。
[0096]
本文中所使用的术语仅为了描述特定实施方案的目的并且不旨在限制本公开。如本文中所使用地,单数形式“一个”、“一种”、以及“该”旨在也包括复数形式,除非上下文明确地另作规定。还将理解的是,当在说明书中使用时,术语“包括”和/或“包含”是指存在有所陈述的特征、数值、步骤、操作、元件和/或组分,但是并不排除存在有或额外增加一个或多个其它的特征、数值、步骤、操作、元件、组分和/或其组成的群组。作为在本文中所使用的,术语“和/或”包括列举的相关项的一个或多个的任何和全部的组合
[0097]
应当理解,此处所使用的术语“车辆”或“车辆的”或其它类似术语一般包括机动车辆,例如包括运动型多用途车辆(suv)、公共汽车、卡车、各种商用车辆的乘用汽车,包括各种舟艇、船舶的船只,航空器等等,并且包括混合动力车辆、电动车辆。
[0098]
具体地,本发明提供一种基于时间同步的车载侧后视图像生成方法,应用于具备多摄像头的汽车,多摄像头通过车载以太网总线连接智能主机,智能主机控制多摄像头拍摄图像的融合拼接形成最终的侧后视图。
[0099]
现有技术的解决方案中,如车载环视系统,在汽车前后左右安装四个设摄像头,实现360环视系统,但是同一区域仅安装一个摄像头,当该区域的摄像头被污染后,其无法实现环视系统。
[0100]
因此,在本实施例中,同一区域安装有多个摄像头,其中一个摄像头被污染了,另外的摄像头可以启动来工作。但多个摄像头并没有主次之分,在都没有被污染的情况下,任意其中一个工作即可,若被污染,选择未被污染的一种摄像头进行工作即可。
[0101]
优先的,冗余摄像头的个数为2个;
[0102]
优先的,为了实现多个设想摄像头采集的图像能够保持同步,以便后续的拼接,本实施例采用车载以太网摄像头,类型镜头类型为鱼眼相机。
[0103]
一种基于时间同步的车载侧后视图像生成方法,具体包括:
[0104]
步骤s1,对多个车载以太网摄像头进行图像采集的时间进行时间同步,并对采集图像进行畸变矫正;
[0105]
步骤s2,获取冗余车载以太网摄像头拍摄的图像,在进行图像融合后获得融合图像,根据融合图像对应原始畸变图像判断选择未被污染的车载以太网摄像头作为工作摄像头;
[0106]
步骤s3,获取不同安装位置的工作摄像头拍摄图像进行畸变矫正中的差异矫正单应性矩阵h
da
posk、矫正融合单应性矩阵h
af
pos,并对工作摄像头拍摄的画面进行图像映射,得到安装位置工作摄像头对应的融合平面的图像;
[0107]
步骤s4,融合图像上的不同安装位置工作摄像头对应图像互相覆盖区域的各个像素值根据信息权重进行覆盖区域的图像融合。
[0108]
具体地,车载以太网设摄像头的时间同步过程如下:
[0109]
冗余车载以太网摄像头同步采集模块负责多路摄像头同步采集图像数据并传递图像数据到鱼眼畸变矫正模块。
[0110]
具体的,冗余车载以太网摄像头同步采集模块包括至少2个广角高清以太网摄像头,每个摄像头通过以太网线与其他模块连接。
[0111]
车载以太网摄像头以eavb协议对采集到的图像数据进行以太网封包并通过以太网连接线传输并进行时间同步,不同的车载以太网摄像头根据同步时间修正自身的rtc时钟,并在统一时间戳进行一帧图像的采集。
[0112]
具体的,采集的图像数据可以采用ieee 1722协议进行以太网封包和传输,其中avtp payload之前的数据统称包头,包头包含可定制控制信息,例如通过stream id来区分摄像头,同一摄像头的以太网封包内的avtp payload内部为该车载以太网摄像头顺序发送的图像数据段。
[0113]
冗余摄像头间的时间同步,具体的,可以通过ieee 802.1as协议进行时间同步如图1所示。
[0114]
例如单个车载以太网摄像头通过以太网向主处理器发送pdelay_req报文,请求测量以太网信号链路传播延迟,当pdelay_req报文离开车载以太网摄像头mac层时,车载以太网摄像头通过其本地rtc时钟记录此时的时间戳t1,同时当主处理器mac层接收到pdelay_req报文时,通过主处理器的本地rtc时钟记录此时的时间戳t2;主处理器随后通过车载以太网向以太网摄像头发送pdelay_resp报文回复主处理器接收到的pdelay_req报文的本地时间戳t2,并且通过主处理器本地rtc时钟记录pdelay_resp报文从主处理器mac层发出的时间t3,随后通过pdelay_reso_follow_up报文将t3发送给车载以太网摄像头;车载以太网摄像头记录下其mac层接收到pdelay_resp报文时候的本地rtc时钟的时间戳t4,并通过解析pdelay_resp_follow_up报文搭载的pdelay_resp报文从主处理器mac层发送的时间戳t3。以太网mac层以下链路的传输延迟基本满足对称关系,即报文从以太网摄像头到主处理器的传输时间和主处理器到以太网摄像头的传输时间是相同的,则链路传播延迟p
delay
可通过以下公式计算:
[0115][0116]
不同以太网摄像头通过此方式得到各自与主处理器的链路传播延迟p
delav1

p
delay2
、......、p
delayn
,则当主处理器不定期对各个以太网摄像头发送主处理器本地时间报文时,各个以太网摄像头解析出报文中记录的主处理器本地时间戳tm,再加上各自的传播链路延迟p
delay
,就得到了收到报文后此时主处理器的本地rtc时间戳tc,同时各个以太网摄像头将各自本地的rtc时钟时间戳修改为tc,则实现了所有以太网摄像头与主处理器的时间同步。
[0117]
各个摄像头根据时间同步后的本地rtc时钟,每间隔预设好的固定时间ts进行一帧图像的采集,则实现了冗余以太网摄像头的图像同步采集,彻底避免了不同摄像头采集时间差异导致的图像空间误差。
[0118]
在步骤s1中,畸变矫正至少包括第一次畸变矫正,其中,第一次畸变矫正至少包括:通过冗余车载以太网摄像头获取垂直于地面的预设靶图成像后的畸变图像,利用光学畸变参数映射表对畸变图像进行第一畸变矫正后获取获得第一畸变矫正图像;
[0119]
第一次畸变矫正具体包括:
[0120]
车载以太网摄像头基本选择鱼眼镜头,镜头lens供应商对该款镜头会提供光学畸变参数表,光学畸变参数表为lens配型指定cmos的情况下的实验模组进行实际实验环境图像采集得出的离散视场角度与其对应像高的映射表如图2。
[0121]
获取车载以太网摄像头的光学畸变参数,光学畸变参数包括摄像头镜头指定cmos芯片条件下的实验模组进行实际实验环境图像采集的离散视场角度θ与其对应像高γ的映射表;
[0122]
根据小孔成像原理,建立实际拍摄平面u各个像素点对应到成像面i上鱼眼图像的各个点的空间坐标关系,建立离散视场角度θ与物高λ的第一映射函数;
[0123]
具体地,参见图3,
[0124]
物高为λ,图1中光学畸变参数表的角度θ单位为度,像高γ单位为mm。
[0125]
设lens光心到实际拍摄平面距离为1,满足关系,则:
[0126]
λ=tan(θ),(γ-》θ)
ꢀꢀꢀ
公式(1)
[0127]
θ=arctan(λ),(θ-》γ)
ꢀꢀꢀ
公式(2)
[0128]
第一映射函数:λ=tan(θ)
[0129]
根据光学畸变参数映射表,构建像高γ与离散视场角度θ第二映射函数θ=g(γ),
[0130]
具体地,在本实施例中,假设:
[0131]
对于任意角度θ,满足θ(n)≤θ≤θ(n 1);在畸变参数表中,查找θ(n)、θ(n 1)相对应的γ(n)与γ(n 1),则与θ对应的γ应满足满足γ(n)≤γ≤γ(n 1);
[0132]
其中n为像素点的序号,θ(n)、θ(n 1)表示相邻离散角度,γ(n),γ(n 1)相邻像高,θ(n)、θ(n 1)、γ(n)、γ(n 1)根据已知θ或γ从光学畸变参数表中查找获得;
[0133]
当γ为已知时,
[0134][0135]
当θ为已知时,
[0136][0137]
根据成像几何关系,求解实际拍摄平面u与对应的矫正成像面、畸变成像面之间的
矫正缩放因子su和畸变缩放因子si;
[0138]
具体包括:
[0139]
设此时op连线与水平方向夹角为ψ,ψ的取值范围为0~180度,则p点在实际拍摄平面u上面的水平坐标p
x
和垂直坐标py和p’点在成像面i上面的水平坐标p

x
和垂直坐标p
′y满足关系,:
[0140]
p
x
=λ
×
cos(ψ),py=λ
×
sin(ψ)
ꢀꢀꢀ
公式(5)
[0141]
p

x
=γ
×
cos(ψ),p
′y=γ
×
sin(ψ)
ꢀꢀꢀ
公式(6)
[0142]
选择实际以太网摄像头的cmos成像面水平方向最大宽度2γ
max
和实际鱼眼畸变图像最大水平分辨率reshi、最大垂直分辨率resvi和畸变矫正图像最大水平分辨率reshu、最大垂直分辨率resvu,带入式1,式2,即可得到无畸变实际拍摄图像的λ
max
和畸变缩放因子si和矫正缩放因子su,满足关系:
[0143][0144][0145]
利用第一映射函数、第二映射函数以及矫正缩放因子su和畸变缩放因子si,建立矫正后的无畸变图像的像素级坐标resp
x
,respy与畸变图像的像素级坐标resp

x
,resp
′y的对应第三映射函数,并将resp

x
,resp
′y坐标位置的像素值赋值给resp
x
,respy的表示的坐标位置,从而获取畸变矫正后的图像。
[0146]
具体包括:
[0147]
成像面鱼眼畸变图像的像素级坐标resp

x
,resp
′y和实际拍摄面无畸变图像的像素级坐标resp
x
,respy满足关系,
[0148][0149][0150][0151][0152]
将公式(9)与公式(11)的等式两边分别相除,将公式(1),公式(5)、公式(6)、公式(7)和公式(8)代入,获取resp
x
与resp

x
的对应第三映射函数:
[0153][0154]
同理:将公式(10)与公式(12)的等式两边分别相除,将公式(1),公式(5)、公式(6)、公式(7)和公式(8)代入,获取respy与resp
′y的对应第三映射函数:
[0155]
[0156]
其中,γ为p点的对应的像高,γ
max
摄像头的cmos成像面水平方向最大宽度的一半;为点p与拍摄平面u的中心o点的连线与水平方向的夹角,reshi、resvi表示实际鱼眼畸变图像最大水平分辨率和最大垂直分辨率,reshu、resvu表示畸变矫正图像最大水平分辨率最大水平分辨率和最大垂直分辨率。
[0157]
以某款cmos水平长度4mm为例,成像面刚刚覆盖完cmos水平方向时,成像的鱼眼畸变图片水平分辨率reshi为1280个像素,此时op在水平方向上对应的γ
max
为2;
[0158]
设光学畸变参数表最小精度满足γ(n)≤γ≤γ(n 1)时的γ(n)为1.95;
[0159]
γ(n 1)为2.03,θ(n)为55.5,θ(n 1)为56,
[0160]
则带入式1,式2可得此时实际拍摄平面u上面op在水平方向上对应点p的λ
max
为1.4721,
[0161]
设实际拍摄无畸变图像水平分辨率reshu为1280像素,则根据式4可得si为320,su为434.75。
[0162]
畸变矫正的图像即实际拍摄面u图像,其上任意像素点坐标(resp
x
,respy)带入式公式(14)和公式(15),则可以求取到鱼眼畸变图像即成像面i图像对应的像素点坐标(resp

x
,resp
′y),由于图像像素坐标为整数,计算得到的对应像素点坐标为浮点数(resp
x
,respy),因此实际p点的像素值可由p

点的周围像素值进行估值获取;
[0163]
具体的估算方法如下:
[0164]
(resp

x
,resp
′y)为浮点时,其周围距离最近的4个像素点p0、p1、p2、p3坐标分别为为
[0165]
获取p0、p1、p2、p3对应的像素值分别为val
p0
、val
p1
、val
p2
、val
p3
;计算p0、p1、p2、p3相应的权重系数:
[0166][0167][0168][0169][0170]
根据权重系数计算出像素值val
p

[0171]
val
p
=val
p0
×
w0 val
p1
×
w1 val
p2
×
w2 val
p4
×
w3
[0172]
按照畸变矫正图像的图像分辨率遍历p点像素坐标(resp
x
,respy)并计算出对应的鱼眼畸变图像p点的像素坐标(resp

x
,resp
′y),将p点的像素值赋为p点插值计算出的val
p
,则完成了鱼眼基本图像的第一次畸变矫正。鱼眼畸变图像如图4,畸变矫正图像如图5。
[0173]
冗余摄像头安装位置必然存在空间差异,因此其拍摄视场图像必然存在差异,且不同的以太网摄像头在生产过程中由于生产工艺精度也必然导致成像面与cmos平面存在一定的光心偏移和焦平面角度误差,这都导致按照既定畸变矫正数学模型将鱼眼畸变图像还原为畸变矫正图像的误差。
[0174]
冗余摄像头差异矫正将冗余的不同以太网摄像头的畸变矫正图像roi区域差异进
行最大限度的消除。
[0175]
因此,需要在第一次畸变矫正图像的基础上进行第二次畸变矫正;
[0176]
畸变矫正至少包括第二次畸变矫正,其中,第二次畸变矫正至少包括:第二次畸变矫正,其中,第二次畸变矫正至少包括:建立差异矫正平面坐标系,分别获取预设靶图的特征点在第一畸变矫正图像的对应的第一坐标点以及在差异矫正平面对应的第二坐标点,利用透视变换,求解单应性矩阵h
da
pos,利用单应性矩阵h
da
pos进行反变化对第一畸变的图像进行第二次畸变矫正。
[0177]
获取单应性矩阵h
da
pos并利用反透视变换将建立差异矫正平面对应的所有第一坐标点在第一畸变矫正图像对应第二坐标点,获取第二坐标点在第一畸变矫正图像对应的像素值,并将像素值赋值给第一坐标点。
[0178]
具体包括:
[0179]
在差异矫正平面坐标系,预设靶图的特征点(t
xn
,t
yn
)在在差异矫正平面对应的第一坐标(r
xn
,r
yn
)的映射函数为:
[0180]rxn
=t
xn
×
sc ofst
x
[0181]ryn
=t
yn
×
sc ofsty[0182]
其中,sc为平面坐标尺度变换因子,ofst
x
和ofsty分别为水平和垂直像素偏移量;
[0183]
步骤s101,通过冗余摄像头获取垂直于地面的预设靶图,获取预设靶图中多个特征点在靶图平面的实际物理坐标(t
xn
,t
yn
);
[0184]
步骤s102,根据实际物理坐标(t
xn
,t
yn
),求解经过畸变校正后实际拍摄面u中与实际物理坐标(t
xn
,t
yn
)相对应的特征点的像素坐标(resp
xn
,resp
yn
);
[0185]
步骤s103,建立差异矫正平面坐标系,求解与靶图中多个特征点的坐标映射到差异矫正平面坐标系下相对应的特征点坐标(r
xn
,r
yn
);
[0186]
步骤s104,根据透视变换,求解(r
xn
,r
yn
)至(resp
xn
,resp
yn
)的单应性矩阵h
da
pos:
[0187]
透视变换:
[0188][0189]hda
pos:
[0190][0191]
步骤s105,再次利用透视变换变化,将差异矫正图像所有像素坐标(r
xi
,r
yi
)带入式7替换(r
xn
,r
yn
),利用步骤s104中求出的h
da
pos,计算出对应的畸变矫正后图像的对应坐标(resp
xi
,resp
yi
)。
[0192]
然后利用如第一畸变矫正中采用的像素值估算方法,计算差异矫正图像所有像素坐标对应的畸变矫正后的像素坐标,获取畸变矫正图像中的像素值并赋值给对应差异矫正图像中的坐标,形成最终的差异矫正图像如图6所示。
[0193]
具体地,在行车过程中或者雨雪等天气影响下,尽管摄像头表面通过憎水膜的物理手段可以极大的降低污物粘着的概率,但仍然有一定可能沾上污渍,污渍沾染具有突发
性,附着可变化的特征,这会减少以太网摄像头获取的图像信息使得驾驶员对真实环境情况感知受阻。清晰度判决计算出摄像头各个区域的清晰度得分,从而判断出摄像头各个区域的污染程度。从而确定工作摄像头和非工作摄像头。
[0194]
将单个摄像头采集的原始图像按边长为b1划分为若干小块imb,每个小块内部按照传统清晰度计算算法计算其清晰度得分。清晰度是表征图像信息含量的尺度,图像高频区域和中频区域的信息占据了绝大部分的图像信息,清晰度计算方式之一通过计算相邻像素的差异性,能体现出信息含量的尺度,当图像清晰时相邻像素差异性大,图像模糊时相邻像素差异性小。
[0195]
具体的,当前帧序号f的水平分辨率为resh、垂直分辨率为resv的图像imb(f),按照边长为bl可以划分为水平hnum、垂直vnum个小正方形块imb(f)
vh
,(v∈[1,vnum],h∈[1,hnum]),其中小块内每个像素的像素值为imb(f)
vh
(i,j),(i∈[1,b1],j∈[1,b1]),则每个小块垂直和水平方向的清晰度分值var(f)
vh
计算如下:
[0196][0197]
其中,imb(f)
vh
表示垂直序号为v,水平序号为h的小块,f为当前图像的帧序号,bl为小块的边长,i和j分别表示小块图像内的像素坐标,i∈[1,bl],j∈[1,b1];
[0198]
连续统计n帧图像的每一个小块内部的水平垂直的清晰度分值,如图7所示,分别对帧内划分的每一个小块imb(f)
vh
,(v∈[1,vnum],h∈[1,hnum])累加统计其在之前连续n帧图像的方差,得到清晰度得分df(f)
vh
,,各个冗余摄像头对应的清晰度得分记为dfk(f)
vh

[0199][0200]
设定清晰度分值阈值thre,thre根据b1大小和统计帧数n设定。具体的,认为相邻像素值差异在thre
min
认为是画面开始模糊,则清晰度分值阈值thre满足下式,根据经典的灰阶辨认尺度,thre
min
可选择8~12。
[0201]
thre=(thre
min
×
bl)2×n[0202]
将各个小块统计的累计清晰度分值与清晰度分值阈值比较,大于等于清晰度分值阈值则认为该区域图像清晰度满足要求。
[0203]
为方便直观观察,将小于清晰度分值阈值的小块染色为黑色,大于等于清晰度分值阈值的小块染色为白色,表面污染严重摄像头及其对应小块统计染色如图8所示,图8a表示未被染色的图像,图8b表示经过统计染色的图像。
[0204]
对比摄像头及其对应小块统计染色图9所示,图9表示表示摄像头表面仅轻微污染。图9a表示未被染色的图像,图9b表示经过统计染色的图像。
[0205]
具体地,需要说明的是,并不是所有的摄像头区域都对于行车有最重大的安全意义,比如天空区域相对地面区域而言,这部分的信息对与行车安全而言的更少,摄像头此部
分区域如果被污渍遮挡相对于摄像头地面区域被遮挡的话危害性更小。摄像头作为外部器件,总有一定概率受到污染,因此当地面区域没有被污染的情况下,该摄像头采集的图像对应行车安全还是可以采纳的。同时安装在不同拍摄位置的冗余摄像头最后选取到融合图像中的区域也有所不同,安装在某一拍摄位置的摄像头内一些区域的数据根本不会采纳到最后的融合图像中,因此相对于最终的融合系统而言,这部分图像是冗余的可以舍弃的,也就意味着该区域被污染对整体的融合图像没有影响。有效内容判决模块将摄像头最终显示区域权重配合清晰度得分进行统计,得出各个冗余摄像头的融合信息含量得分,根据各个冗余摄像头融合信息含量得分选择得分最高的摄像头为工作摄像头,其余摄像头为待机摄像头。
[0206]
清晰度得分计算方法得到的清晰度得分小块的块坐标定义为(v,h),v∈[1,vnum],h∈[1,hnum],其水平和垂直像素坐标分别为resph=h
×
bl respv=v
×
bl。
[0207]
设靶图特征点为水平k个,垂直j个,k
×
j=t,根据冗余摄像头差异矫正模块步骤,可知靶图t矩形roi区域4个顶点在融合矫正图像中的像素坐标(resq
x1
,resq
y1
)、(resq
xk
,resq
yk
)、(resq
x(t-k 1)
,resq
y(t-k 1)
)、(resq
xt
,resq
yt
),如图10所示。
[0208]
不同摄像头安装区域分别对应一个靶图,靶图摆放在同一平面垂直于水平面,各靶图相对位置在靶图摆放像素平面q上,3个靶图在水平方向上居中,如图11。其中ofst_right是像素屏幕q最左边像素到右侧相机拍摄到靶图roi的点(resq
x1
,resq
y1
)的水平像素距离,ofst_vt是q最上边像素到右侧相机拍摄到靶图right camera roi的点(resq
x1
,resq
y1
)的垂直像素距离,ofst_vby是q最上边像素到右侧相机拍摄到靶图right camera roi的点(resq
xt
,resq
yt
)的垂直像素距离;
[0209]
ofst_rear是像素屏幕q最左边像素到后侧相机拍摄到靶图rear camera roi的点(resq
x1
,resq
y1
)的水平像素距离;
[0210]
ofst_left是像素屏幕q最左边像素到左侧相机拍摄到靶图left camera roi的点(resq
x1
,resq
y1
)的水平像素距离。
[0211]
设观察点在三维坐标系下的o点,其在水平面g的投影为g,如图12。其中o点离地平面的y轴像素距离为viewy,到靶图的z轴像素距离为viewz,靶图上各点离o点的x轴的像素距离为view
x
。则设g点为三维坐标系的原点,图12中角点顺序从左到右从上到下分别为c1、c2、c3、c4,则各个靶图的roi区域观察到的融合图像q的角点坐标在此坐标系下分别为:
[0212]
右侧摄像头roi区域,
[0213][0214]
rgtfc1y=ofst_vb-ofst_vt
[0215][0216]
rgtfc1y=ofst_vb-ofst_vt
[0217][0218]
rgtfc3y=0
[0219][0220]
rgtfc4y=0
[0221]
后侧摄像头roi区域,
[0222][0223]
rerfc1y=ofst_vb-ofst_vt
[0224][0225]
rerfc2y=ofst_vb-ofst_vt
[0226][0227]
rerfc3y=0
[0228][0229]
rerfc4y=0
[0230]
左侧摄像头roi区域,
[0231][0232]
lftfc1y=ofst_vb-ofst_vt
[0233][0234]
lftfc2y=ofst_vb-ofst-vt
[0235][0236]
lftfc3y=0
[0237][0238]
lftfc4y=0
[0239]
则根据透视原理,设左右摄像头安装位置距离其对应靶图的像素距离为disz,则最终融合图像中各靶图矩形roi区域的4个顶点c(i),i=1、2、3、4,在观察点o看到的融合图像中的像素坐标为的关系为:
[0240]
右侧摄像头roi区域,
[0241][0242]
后侧摄像头roi区域,式17
[0243][0244]
左侧摄像头roi区域,式18
[0245][0246]
设各个安装位置的差异矫正图像与多个摄像头融合图像的单应性矩阵为h
af
pos,pos=left,right,rear。
[0247]
右侧单个摄像头差异矫正图像与多个摄像头融合图像存在单应性矩阵,将(rgtfc(i)
x
,rgtfc(i)y),i=1、2、3、4替代(r
xn
,r
yn
),(rightc
x(i)
,rightc
y(i)
),i=1,2,3,4替代(f
xn
,f
yn
),带入式8,即可求得该右侧摄像头的差异矫正图像和多摄像头融合图像的单应性矩阵h
af
right;
[0248]
后侧单个摄像头差异矫正图像与多个摄像头融合图像存在单应性矩阵,将(rerfc(i)
x
,rerfc(i)y),i=1、2、3、4替代(r
xn
,r
yn
),rearc
x(i)
,rearc
y(i)
),i=1,2,3,4,带入式8,即可求得该右侧摄像头的差异矫正图像和多摄像头融合图像的单应性矩阵h
af
rear;
[0249]
左侧单个摄像头差异矫正图像与多个摄像头融合图像存在单应性矩阵,将(lftfc(i)
x
,lftfc(i)y),i=1、2、3、4替代(r
xn
,r
yn
),(leftc
x(i)
,leftc
y(i)
),,i=1,2,3,4,带入式8,即可求得该左侧摄像头的差异矫正图像和多摄像头融合图像的单应性矩阵h
af
left;
[0250]
设最终融合图像q的像素分辨率为水平resqh、垂直resqv,则点的坐标范围为:
[0251][0252]
qy,y∈[ofst_vb,ofst_vb-resqv]
[0253][0254]
resqv=ofst_vb ofst_va
[0255]
上述给出融合图像qx和qy的取值范围;
[0256]
以右侧摄像头为例进行说明:
[0257]
将右摄像头最终融合平面图像right(q
x
、qy)、h
af
pos利用透视变换,获取右侧各个摄像头的最终融合图像各像素点坐标对应的差异矫正图像坐标(r
x
,ry),公式如下:
[0258][0259]
将(r
x
,ry)和右侧各个摄像头的单应性矩阵h
da
posk,利用透视变换,获取各个摄像头的差异矫正图像坐标对应的畸变矫正图像坐标(resp
x
,respy),公式如下:
[0260][0261]
根据鱼眼相机的畸变矫正原理,获取畸变矫正图像坐标(resp
x
,respy)相对应的原始畸变图像的坐标(resp

x
,resp
′y);
[0262]
重复上述步骤,获取不同摄像头对应的原始畸变图像坐标(resp

xk
,resp

yk
),其中,k为摄像头的编号。
[0263]
根据原始畸变图像(resp

x
,resp
′y)获得对应的图像区域,并从图像区域进行分块,能够获得右侧各个摄像头k最终融合图像各像素点rightk(q
x
、qy)对应的清晰度得分right(dfk(f)
vh
),记为rightdfk(f)
x,y

[0264]
同理,后侧和左侧各个摄像头d最终融合平面图像rear(q
x
、qy)和h
af
rear、left(q
x
、qy)和h
af
left并按照上述步骤即可得到后侧各个摄像头k最终融合图像各像素点rear k(q
x
、qy)对应的清晰度系数rear(dfk(f)
vh
),记为reardfk(f)
x,y
;左侧各个摄像头k最终融合图像各像素点leftk(q
x
、qy)对应的清晰度系数left(dfk(f)
vh
),记为leftdfk(f)
x,y

[0265]
以左侧污染摄像头为例,其最终融合平面图像对应的清晰度系数leftdfk(f)按照阈值thre判断如图13所示;
[0266]
各安装位置摄像头区域信息权重为摄像头实际安装位置拍摄到的融合图像中驾驶员主观观察概率最大区域,即后、侧后方地面及车身延展线以下区域,该区域的获取方式为,将清洁的摄像头安装在指定安装位置并朝向指定视场,得到该摄像头在融合平面上的图像,优选的,选择ofst_vt以下区域,如图14灰色区域.
[0267]
根据透视原理,离融合图像中心区域贴图越近,则信息有效性越大。设定垂直信息权重函数vld
(y)
,水平信息权重函数vld
(x)
,则信息权重函数vld
(x,y)
[0268]
垂直信息权重函数vld
(y)

[0269][0270]
水平信息权重函数vld
(x)

[0271][0272]
ofst_h1和ofst_h2,low_val均为人为设定,low_val为主观认为的最低权重,ofst_h1是主观认定的最低权重区域,ofst1到ofst2之间的区域为权重渐变区域;
[0273]
信息权重函数vld
(x,y)

[0274]
vld
(x,y)
=vld
(x)
×
vld
(y)
[0275]
将信息权重1染色为纯白,0染色为纯黑,则vld
(x,y)
如图15所示
[0276]
单个摄像头k在最终融合图像的当前帧序号f对应的信息权重infk(f)由清晰度系数dfk(f)
(x,y)
和区域信息权重vld
(x,y)
共同决定:
[0277][0278]
其中,k表示冗余摄像头的编号,resqh、resqv分别为融合后图像水平分辨率和垂直分辨率,x,y表示融合图像像素坐标,f表示帧的序列号;
[0279]
dfk(f)
(x,y)
表示冗余摄像头k的清晰度得分,vld
(x,y)
表示融合画面的区域信息权重函数。
[0280]
选择相同安装区域内所有冗余摄像头中信息权重得分最高的一个摄像头作为工作摄像头,其余摄像头作为待机摄像头。融合图像选择各个安装区域中的工作摄像头的最终平面融合图像进行融合。
[0281]
具体地,各个安装区域的摄像头对应的最终融合平面的图像在最终融合平面上存在覆盖区域,覆盖区域图像进行加权融合,像素值加权方式根据信息权重进行计算,其中leftinf(f)
(x,y)
为左侧工作摄像头的信息权重,leftdfk(f)
(x,y)
为左侧工作摄像头的清晰度系数;rearinf(f)
(x,y)
为后侧工作摄像头的信息权重,reardfk(f)
(x,y)
为后侧工作摄像头的清晰度系数;rightinf(f)
(x,y)
为右侧工作摄像头的信息权重,rightdfk(f)
(x,y)
为右侧工作摄像头的清晰度系数,q
rr(x,y)
为右侧和后侧工作摄像头覆盖区域融合图像,q
rl(x,y)
为左侧和后侧工作摄像头覆盖区域融合图像,right(q
x
、qy)、rear(q
x
、qy)、left(q
x
、qy)分别为右侧、后侧、左侧工作摄像头在最终融合平面上的图像。
[0282]
leftinf(f)
(x,y)
=leftdfk(f)
(x,y)
×
vld
(x,y)
[0283]
rearinf(f)
(x,y)
=reardfk(f)
(x,y
)
×
vld
(x,y)
[0284]
rightinf(f)
(x,y)
=rightdfk(f)
(x,y)
×
vld
(x,y)
[0285][0286][0287]
根据实际显示效果,人为定义贴图区域q
c(x,y)
,最终融合图像该区域的像素替换为贴图像素。
[0288]
各个区域冗余摄像头(未污染和污染)如图16所示;
[0289]
各个区域冗余摄像头(未污染)最终融合图像加上贴图效果如图17所示。
[0290]
以上所述的仅是本发明的优选实施方式,本发明不限于以上实施例。本领域的技术人员可以清楚,该实施例中的形式不局限于此,同时可调整方式也不局限于此。可以理解,本领域技术人员在不脱离本发明的基本构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。
再多了解一些

本文用于创业者技术爱好者查询,仅供学习研究,如用于商业用途,请联系技术所有人。

发表评论 共有条评论
用户名: 密码:
验证码: 匿名发表

相关文献