环境因素对空三摄影影响实验装置和测评方法与流程

1.本发明涉及空中三角测量领域，特别涉及环境因素对空中三摄影纹理数据的影响与评测方法。


背景技术：

2.目前，对于同一目标三维成像，有很多测量方法，但是，各种测量方法都没有给出测量精度的评估方法。对于空中三角测量来说，拍摄图像质量评估是三维数据重建的关键因素。由于用于空中三角测量的飞行器会受到各种干扰，干扰因素，对飞行器拍摄数据质量会造成影响，目前，没有方法和手段评估环境因素对空中三角测量拍摄数据的影响。
3.激光点云与空三点云广泛应用于地形数据获取。申请号202010736839.7《一种三维激光点云和全景图像的配准方法》将激光三维点云与全景图像二维数据进行配置。申请号：20210090259.5《综合机载和车载三维激光点云以及街景影像的三维建筑物精细几何重建方法》综合了机载和车载三维激光点云以及街景影像的三维建筑物精细几何重建。申请号201510228288.7《空中三角测量的方法与装置》给出了一种空中三角测量的方法与装置。上述三个技术方案均有所在领域的技术突破，但与本发明的技术不同，本发明用于模拟环境因素并评估环境因素对空三纹理影响。


技术实现要素：

4.本发明的目的是提供一种实验装置，根据实验装置采集的数据，正确评估环境因素对目标拍摄的影响的方法，还原真实场景的重要信息。
5.本发明的目的是这样达到的：一种环境因素对空三摄影影响实验装置。实验装置由评测模具组合、倾斜摄影无人机、光源、风源组成；评测模具组合由十五个评测模具构成。含球形评测模具、正多面体评测模具。球形模具包括：球形分辨率评测模具，球形色度评测模具，球形灰度评测模具。正多面体评测模具包括：正四面体分辨率评测模具，正四面体色度评测模具，正六面体分辨率评测模具，正六面体色度评测模具，正十二面体分辨率评测模具，正十二面体色度评测模具，正二十面体分辨率评测模具，正二十面体色度评测模具；正四面体灰度评测模具；正六面体灰度评测模具；正十二面体灰度评测模具；正二十面体灰度评测模具。十五个评测模具采用任意方式摆放在一个平面。
6.正多面体评测模具各面有内接圆。
7.球形评测模具为半径为r的球体，球形评测模具球体的外表面称为球面，球体按照上、下、左、右、前、后划分为上端评测球面、下端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面；将每一个评测球面划分为64个区域，称为评测区。
8.球形评测模具球体的评测区与正多面体评测模具各面内接圆上分别涂上分辨率评测图样、色度评测图样和灰度评测图样。
9.光源由三脚架、反光面罩、可调灯座、试验灯构成；三脚架用于固定反光面罩；可调灯座安装在反光面罩上，用于将实验灯放置在合适位置，并为试验灯提供电源以及电源连
接座；反光面罩用于将实验灯的光反射，照到被拍摄的模具上；试验灯用于发出需要照度和色温的光，所发的光被反光面罩照射到被拍摄摸具上。
10.风源由风机三脚架、轴流风机、风机调速器、风速仪构成。
11.风机调速器连接电源，轴流风机电源连接线连接到风机调速器，由风机调速器控制轴流风机的风速；轴流风机安装在风机三脚架上；风速仪用于测量轴流风机出口的风速。
12.对于球形分辨率评测模具，在各评测球区的64个色区，涂分辨率评测图样，分辨率评测图样为黑白相间的颜色。
13.对于球形色度评测模具，在评测区涂色的64个色区，涂色度评测图样，色度评测图样为1-64号不同颜色。
14.对于球形灰度评测模具，在各评测球区的64个色区涂灰度评测图样，灰度评测图样为灰色。
15.正多面体分辨率评测模具中，正四面体分辨率评测模具为正四面体、正六面体分辨率评测模具为正六面体、正十二面体分辨率评测模具为正十二面体、正二十面体分辨率评测模具为正二十面体，各个正多面体各面的内接圆涂有分辨率评测图样；分辨率评测图样为黑白相间的颜色。
16.正多面体色度评测模具中，正四面体色度评测模具为正四面体、正六面体色度评测模具为正六面体、正十二面体色度评测模具为正十二面体、正二十面体色度评测模具为正二十面体，各个正多面体各面的内接圆涂有色度评测图样,色度评测图样为1-64号不同颜色。
17.正多面体灰度评测模具中，正四面体灰度评测模具为正四面体，正六面体色度评测模具为正六面体，正十二面体色度评测模具为正十二面体，正二十面体色度评测模具为正二十面体，各个正多面各面的内接圆涂有灰度评测图样，灰度评测图样为灰色。
18.光源的反光面罩形状为一个旋转抛物面，含罩体、固定柄、反光层、遮光罩固定沿和遮光板。反光层为反光材料，贴在反光面罩内侧；设反光面罩所在的旋转抛物面顶点为a，直线ao为反光面罩的旋转抛物面轴线；过轴线的平面与旋转抛物面相交线为抛物线，设该抛物线的焦点为o；过o点与直线ao垂直的平面与反光面罩的旋转抛物面相交线为圆，该圆称为焦点平面圆，在焦点平面圆上，等距离分布三个调节杆安装槽，相邻两个调节杆安装槽之间弧度为120度：每个调节杆安装槽两侧，各有四个调节杆定位孔；遮光罩固定沿为圆环形，连接在反光面罩底端，其所在平面与直线ao垂直；在遮光罩固定沿上有两个遮光罩固定孔，两个遮光罩固定孔的连线通过遮光罩固定沿的圆心。
19.可调灯座由试验灯安装板、调节杆组成；试验灯安装板中间位置，钻有若干个试验灯安装螺孔，用于安装试验灯；在试验灯安装板圆盘外侧，焊接有三个调节杆安装口；调节杆安装口为u形，由左侧安装支架、右侧安装支架、底部固定板构成；左侧安装支架、右侧安装支架、底部固定板为一个整体；左侧安装支架、右侧安装支架形状相同，为矩形和半圆形构成，半圆形一侧有一个定位固定孔88。
20.调节杆由调节杆柱状端和调节杆扁平端构成，调节杆柱状端长度远远比调节杆扁平端长度长；调节杆柱状端与调节杆扁平端构成一个整体；在调节杆扁平状外侧中部，有一个扁平端固定孔。扁平端固定孔直径与定位固定孔直径相同；调节杆装入调节杆定位管定位安装。
21.调节杆定位管分为定位管本体、定位管左侧限定座、定位管右侧限定座，在定位管右侧限定座上端设置有调节杆紧固螺孔，调节杆紧固螺钉将调节杆紧固。调节杆紧固螺钉由紧固螺钉外丝与紧固螺钉手柄组合而成，紧固螺钉外丝为外丝螺钉，与调节杆紧固螺孔匹配；紧固螺钉手柄用于旋转调节杆紧固螺钉。
22.定位管左侧限定座、定位管右侧限定座安装在调节定位杆限定板内，安装在调节定位杆限定板内，调节定位杆限定板用于将调节杆定位管限制在反光面罩；两端为限定板端侧；中间为限定板中间部分。限定板中间部分为拱形，用于限制定位管右侧限定座、定位管左侧限定座的活动范围。
23.调节杆与可调灯座之间安装时，调节杆扁平端插入到左侧安装支架和右侧安装支架中间，调节杆安装螺栓依次穿过右侧安装支架定位固定孔、扁平端固定孔、左侧安装支架定位固定孔后，用调节杆安装螺帽和调节杆安装螺栓拧紧，使得调节杆固定在调节杆安装口，且可灵活转动。
24.调节杆与反光面罩之间的安装时，调节杆安装在定位管本体内。
25.调节杆定位管的安装：调节定位杆限定板将调节杆定位管的定位管右侧限定座和定位管左侧限定座固定在反光面罩上；固定时，定位管本体嵌入调节杆安装槽，调节杆定位管的定位管右侧限定座和定位管左侧限定座嵌入调节定位杆限定板的限定板中间部分，每个调节定位杆限定板两端的限定板固定孔与同一侧的调节杆定位孔对齐；每个调节定位杆限定板均用两对螺栓螺帽固定；限定板安装螺栓穿过限定板固定孔、调节杆定位孔后，用螺帽拧紧，将调节定位杆限定板固定；调节杆定位管的定位管右侧限定座和定位管左侧限定座被调节定位杆限定板限制，能够绕着定位管右侧限定座和定位管左侧限定座的轴线转动。
26.可调灯座与调节杆的安装：三个调节杆的调节杆扁平端通过调节杆安装口固定在可调灯座上，三个调节杆的调节杆柱状端分别穿过三个定位管本体，被调节杆紧固螺钉固定在调节杆定位管上；调节杆紧固螺钉拧松的时候，调节杆的调节杆柱状端在调节杆定位管内滑动；调节杆紧固螺钉拧紧的时候，调节杆的调节杆柱状端与调节杆定位管固定。
27.遮光版由遮光板和遮光板折叠支架构成；遮光板边框由一个边框侧边、两个边框支撑边、摄影用柔光布构成；边框侧边和边框支撑板均为长条形，长度相同，长度大于遮光罩固定沿的外径，两个边框支撑边在遮光板边框两端，与遮光板边框垂直连接；在两个边框支撑边的中间均有支撑边固定孔。
28.摄影用柔光布裁剪为正方形，边长与边框侧边相同，其中三个边粘贴在边框侧边和两个边框支撑边，构成遮光板。
29.遮光板折叠支架由若干跟折叠杆体连接而成；相邻折叠杆体采用螺栓与螺帽紧固，相邻折叠杆体之间的折叠杆体连接孔称为折叠杆体邻间连接孔；螺栓穿过相邻折叠杆的折叠杆体邻间连接孔后，用螺帽紧固；最外的孔为折叠支架固定孔，两个折叠支架固定孔一端对应支撑边固定孔一端对应遮光板固定孔；遮光板折叠支架一端，螺栓穿过遮光板固定孔与折叠支架固定孔，用螺帽紧固，使得遮光板折叠支架固定在遮光罩固定沿上；遮光板折叠支架另一端，螺栓穿过支撑边固定孔与折叠支架固定孔，用螺帽紧固，使得遮光板折叠支架固定在边框支撑边上采用两个遮光板折叠支架，将两个边框支撑边连接到遮光罩固定沿。
30.球形评测模具为半径为r的球体；在球体上设置球体z线、球体x线、球体y线；z线31、球体x线32、球体y线33穿过球心，彼此之间相互垂直；z线为垂直方向，球体z线与球面有两个交点：与上边球面的交点称为球体顶端；与下边球面的交点称为球体底端。
31.设球体x线为左右方向；设球体y线为前后方向。
32.球体x线与球面有两个交点：与左边球面的交点称为球体左端；与右边球面的交点称为球体右端。
33.球体y线与球面有两个交点：与前边球面的交点称为球体前端；与后边球面的交点称为球体后端；在球体顶端下方，距离球体顶端0.28r，且与球体z线垂直的平面将球面分割为上下两部分，分割平面上边部分称为上端评测球面。
34.在球体底端上方，距离球体底端0.28r，且与球体z线垂直的平面将球面分割为上下两部分，分割平面下边部分称为下端评测球面。
35.在球体左端右方，距离球体左端0.28r，且与球体x线垂直的平面将球面分割为左右两部分，分割平面左边部分称为左端评测球面。
36.在球体右端左方，距离球体右端0.28r，且与球体x线垂直的平面将球面分割为左右两部分，分割平面右边部分称为右端评测球面。
37.在球体前端后方，距离球体前端0.28r，且与球体y线垂直的平面将球面分割为前后两部分，分割平面前边部分称为前端评测球面。
38.在球体后端前方，距离球体后端0.28r，且与球体y线垂直的平面将球面分割为前后两部分，分割平面后边部分称为后端评测球面。
39.上端评测球面、下端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面面积、形状、涂色均相同，只是朝向不同。
40.球形评测模具的上端评测球面、下端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面均有1-64个色区，64个色区为测评区。
41.在上端评测球面将球体z线垂直的底部小圆称为上端底部小圆，将上端底部小圆等分成64段等长的弧，共计64个分段点，称为等弧分段点；在上端分测球面上，连接球体顶端与64个等弧分段点的圆弧，将上端分辨率评测球面分割成64个扇面，将64个扇面称为64个区域，按顺时针方向将64个区域依次称为1号色区，2号色区
……
,64号色区，64个色区为测评区。
42.对于球形分辨率评测模具，在各评测球区的64个色区，分别涂色黑白相间的颜色；黑色颜色分量分别为：红色，0；绿色，0；蓝色，0；白色颜色分量分别为：红色，255；绿色，255；蓝色，255。对于球形色度评测模具，在评测区涂色的64个色区，分别涂色红色、绿色、蓝色分量不同的1-64号色。
43.对于球形灰度评测模具，在各评测球区的64个色区均涂色为灰色，第n号色区，所涂颜色为：红色分量＝蓝色分量＝黄色分量＝4(n-1)。
44.对于正四面体分辨率评测模具为正四面体、正六面体分辨率评测模具为正六面体、正十二面体分辨率评测模具为正十二面体、正二十面体分辨率评测模具为正二十面体，各个正多面体各面的内接圆涂有分辨率评测图样，评测图样为分别涂色黑白相间的颜色；黑色颜色分量分别为：红色，0；绿色，0；蓝色，0；白色颜色分量分别为：红色，255；绿色，255；蓝色，255。
45.对于正四面体色度评测模具为正四面体、正六面体色度评测模具为正六面体、正十二面体色度评测模具为正十二面体、正二十面体色度评测模具为正二十面体，各个正多面体各面的内接圆涂有色度评测图样，评测图样为红色、绿色、蓝色分量不同的1-64号色。
46.对于正四面体灰度评测模具为正四面体，正六面体色度评测模具为正六面体，正十二面体色度评测模具为正十二面体，正二十面体色度评测模具为正二十面体，正二十面体各面的内接圆上，各个正多面各面的内接圆涂有灰度评测图样，第n号色区，所涂颜色为：红色分量＝蓝色分量＝黄色分量＝4(n-1)。
47.无人机操控包括正摄影、倾斜摄影，正摄时：
48.一机位：评测模具组合外接圆圆心正上方；
49.二无人机高度：拍摄图片刚好包含所有模具的高度；
50.三无人机操控：人工操控；
51.四：无人机离风源距离：保证无人机安全的最近距离；
52.倾斜摄影时：
53.一倾斜角度：45度
54.二机位：评测模具组合外接圆圆心与水平面45度方向；
55.三无人机高度：拍摄图片刚好包含所有模具的高度；
56.四无人机操控：人工操控；
57.五：无人机离风源距离：保证无人机安全的最近距离；
58.风源控制包括风速选择和轴流风机朝向；
59.风速选择：选择任意风速。
60.轴流风机朝向：轴流风机轴向对准无人机；轴流风机轴向与水平面的角度设置为任意角度。
61.光源控制：光照角度、照度、光阴；
62.光照角度：反光面罩轴向对准模具外接圆圆心，反光面罩轴向与水平面的角度为任意角度；
63.反光面罩距离调节为光线覆盖评测模具组合外接圆区域；
64.照度采用任意照度；
65.光阴：光阴效果由遮光板产生；通过遮光板遮罩面积，产生任意光阴效果；
66.所述无人机操控、风源控制、光源控制在测评是采用任意组合。
67.测评控制组合为；
68.一光源组合：
69.光强组合：四种：最大照度；1/2照度；1/4照度；1/8照度；
70.光阴组合：四种：遮光板采用如下方式：不遮罩；遮罩1/4；遮罩1/2；遮罩3/4；
71.光照角度组合：二种：反光面罩轴向与水平面的角度为：45度；90度；
72.光源采用上述光强、光阴、光照角度组合，总共组合形态：32种；
73.二风源组合：
74.风速组合：七种：风速1：0米/秒；风速2：1.5米/秒；风速3：3.3米/秒；风速4：5.4米/秒；风速5：7.9米/秒；风速6：10.7米/秒；风速7：13.8米/秒；
75.风向组合：三种：轴流风机轴向与水平面的角度分别为0度，45度，-45度；
76.风源组合方式：计21种。
77.三摄影方式组合
78.两种：正射，倾斜
79.四合计组合：32
×
21
×
2＝1344种。
80.在各种评价组合下拍摄数据，针对每种评价组合进行评测，评测参数计算过程包括数据预处理、像素序列提取、黑白分辨率、灰度评价参数计算和色彩评价参数计算，其中，色彩评价参数计算含色相h、亮度l、饱和度s计算和色相偏差，亮度偏差，饱和度偏差计算。
81.数据预处理：
82.一评测区域
83.对于球形评测模具：针对上端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面能够拍出完整评测图像的区域作为评测区域；
84.对于正四面体、正六面体、正十二面体、正二十面体评测模具：针对能够拍摄完整评测图像的面作为评测区域；
85.二数据预处理
86.以平行于照片底边为其中一条边，做矩形，提取评测图像，矩形的四条边为评测图像的切线；设宽度方向包含的像素为n个，高度方向包含的像素为m个,且n大于等于m；用i表行，j表示列；i为整数，取值为1、2、3、
……
、m；j为整数，取值为1、2、3、
……
、n；
87.图片左上角为第1行第1列；右下角为第m行第n列第i行，第j列的像素颜色表达为：红色分量：ri(i,j)；绿色分量：gi(i,j)；蓝色分量：bi(i,j)；
88.如果m等于n：对评测图像不需进行转换；
89.如果m小于n，按如下方法转换图像，转换后的图像，宽度方向为n个像素，高度方向为n个像素；k为整数，为1、2、3、
……
、n；
90.转换后的图片左上角为第1行第1列；右下角为第n行第n列；
91.转换后的第k行，第j列的像素颜色表达为：红色分量：ro(k,j)；绿色分量：go(k,j)；蓝色分量：bo(k,j)；
92.对所提取的评测图像进行如下转换：
93.(1)第一步：计算m除以(n-m)，设x为m除以(n-m)的整数部分；
94.(2)第二步：q＝0；p＝0；r＝0；
95.(3)第三步：q＝q 1；p＝p 1；
96.ro(q,j)＝ri(q-r,j)；go(q-r,j)＝gi(q,j)；bo(q,j)＝bi(q-r,j)；
97.j取值为1、2、3、
……
、n；
98.(4)第四步：判断q是否等于n；等于n，进入第八步；不等于n，则进入第五步；
99.(5)第五步：判断p是否等于x；p等于x，进入第六步；p不等于x，转入第三步；(6)第六步：p＝0；q＝q 1；r＝r 1；
100.ro(q,j)＝(ri(q-1,j) ri(q,j))/2；
101.go(q,j)＝(gi(q-1,j) gi(q,j))/2；
102.bo(q,j)＝(bi(q-1,j) bi(q,j))/2；
103.j取值为1、2、3、
……
、n；
104.(7)第七步：判断q是否等于n；等于n，进入第八步；不等于n，则进入第三步；
105.(8)第八步：转换结束；
106.对转换后的彩色像素进行灰度计算，计算得到的图像用igr表示：
107.gr(k，j)＝ro(k，j)
×
0.299 g o(k，j)
×
0.587 b o(k，j)
×
0.114；
108.计算得到的值取整；
109.像素序列提取方法：设n除以2，取整为n50；n除以4，取整为n25；提取半径rn的像素序列,rn小于n50，行号用row表示，列号用col表示：像素序列用a表示，a＝1,2,3,4,5，
……
110.(1)第一步a＝1；temp＝0；xa＝0；
111.(2)第二步：col(a)＝n50
–
rn xa；
112.(3)第三步：计算对vx取整数
113.(4)第四步：判断vx-temp是否大于1；小于等于1，row(a)＝n50-vx，进入第七步；大于1，
114.则令t＝vx-temp，进入第五步；
115.(5)第五步：row(a)＝n50-temp；a＝a 1；t＝t-1；
116.(6)第六步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp 1；进入第五步；否则，进入第七步；
117.(7)第七步：temp＝vx,a＝a 1；xa＝xa 1；
118.(8)第八步：如果xa大于rn,则转到第九步；否则转到第二步；
119.(9)第九步：col(a)＝n50
–
rn xa计算对vx取整数；
120.(10)第十步：判断temp-vx是否大于1；小于等于1，row(a)＝n50-vx，进入第十三步；大于1，则令t＝temp-vx，进入第十一步；
121.(11)第十一步：row(a)＝n50-temp；a＝a 1；t＝t-1；
122.(12)第十二步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp-1；进入第十一步；否则，进入第十三步；
123.(13)第十三步：temp＝vx,a＝a 1；xa＝xa 1；
124.(14)第十四步：如果xa大于2rn,则转到第十五步；否则转到第九步；
125.(15)第十五步：xa＝xa-1；
126.(16)第十六步：col(a)＝n50
–
rn xa计算对vx取整数；
127.(17)第十七步：判断vx-temp是否大于1；小于等于1，row(a)＝n50 vx，进入第二十步；大于1，则令t＝vx-temp，进入第十八步；
128.(18)第十八步：row(a)＝n50 temp；a＝a 1；t＝t-1；
129.(19)第十九步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp 1；进入第十八步；否则，进入第二十步；
130.(20)第二十步：temp＝vx,a＝a 1；xa＝xa-1
131.(21)第二十一步：如果xa等于rn,则转到第二十二步；否则转到第十六步；
132.(22)第二十二步：col(a)＝n50
–
rn xa，计算对vx取整数；
133.(23)第二十三步：判断temp-vx是否大于1；小于等于1，row(a)＝n50- vx，进入第二十六步；大于1，则令t＝temp-vx，进入第二十四步；
134.(24)第二十四步：row(a)＝n50 temp；a＝a 1；t＝t-1；
135.(25)第二十五步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp-1；进入第二十四步；否则,进入第二十六步；
136.(26)第二十六步：temp＝vx,a＝a 1；xa＝xa-1；
137.(27)第二十七步：如果xa大于零,则转到第二十二步；否则转到第二十八步；
138.(28)第二十八步：cnta＝a；程序结束。
139.黑白分辨率：
140.(1)第一步：rn＝80；
141.(2)第二步：根据像素序列提取方法，计算像素序列的行列序列：col(a)，row(a)，计算cnta；cnta为序列的长度；
142.(3)第三步：从灰度图像igr中，提取像素序列，像素序列用si(b)表示，b＝0、1、2、3、
……
、cnta；si(b)＝igr(row(b),col(b))；
143.(4)第四步：令his(c)＝0,c＝0、1、2、3、
……
、255；d＝1；
144.(5)第五步：his(si(d))＝his(si(d)) 1；
145.(6)第六步：d＝d 1；
146.(7)第七步：判断d是否小于等于cnta？是：进入第五步；否：进入第八步；
147.(8)第八步：suml＝0，sumr＝0；cntl＝0；cntr＝255；
148.(9)第九步：suml＝suml his(cntl)；sumr＝sumr his(cntr)；
149.(10)第十步：判断suml是否大于sumr，是则cntr＝cntr-1，sumr＝sumr his(cntr)；否：则cntl＝cntl 1，suml＝suml his(cntl)；
150.(11)第十一步：判断cntl是否大于cntr：是，则进入第十二步；否，进入第十步；
151.(12)第十二步：d＝1,maxlv＝his(0)；maxl＝0；
152.(13)第十三步：判断his(d)是否大于maxlv；是：maxlv＝his(d),maxl＝d；进入第十四步；否：进入第十四步；
153.(14)第十四步：d＝d 1；
154.(15)第十五步：判断d是否大于cntl，如果是，则进入第十六步；如果否，则进入第十三步；
155.(16)第十六步：d＝254,maxrv＝his(255)；maxr＝255；
156.(17)第十七步：判断his(d)是否大于maxrv；是：maxrv＝his(d),maxr＝d；进入第十八步；否：进入第十八步；
157.(18)第十八步：d＝d-1，进入第十九步；
158.(19)第十九步：判断d是否大于cntl，否，则进入第二十步；是，则进入第十七步；
159.(20)第二十步：d＝1；stepl＝0.5(maxl cntl)；stepr＝0.5(maxr cntl)；f＝0；cntinv＝0；(21)第二十一步：si(d)是否小于stepl？是：f＝-1进入第二十六步；否，进入第二十二步；
160.(22)第二十二步：si(d)是否大于stepr？是：f＝1进入第二十三步；否，进入第二十一步；
161.(23)第二十三步：d＝d 1
162.(24)第二十四步：判断d是否大于cnta？是，进入第二十九步；否，进入第二十五步；
163.(25)第二十五步：si(d)是否小于stepl？是：f＝-1，cntinv＝cntinv 1；进入第二十六步；否，进入第二十三步；
164.(26)第二十六步：d＝d 1；
165.(27)第二十七步：判断d是否大于cnta？是，进入第二十九步；否，进入第二十八步；
166.(28)第二十八步：si(d)是否大于stepr？如果是：f＝1，cntinv＝cntinv 1；进入第二十三步；如果否，进入第二十六步；
167.(29)第二十九步：判断cntinv是否大于62；如果是进入第三十步；如果否，则rn＝rn 1，进入第二步；
168.(30)第三十步：黑白分辨率为rn对应的像素值
169.(31)第三十一步：程序结束。
170.灰度评价参数计算:
171.对于灰度评测摸具对应的图片，进行灰度评价参数计算；灰度评价参数指标有：灰度完整性，用wzgr表示；各色区颜色分量偏差，红色分量用rpcgr(a)表示，绿色分量用gpcgr(a)表示，蓝色分量用bpcgr(a)表示；色区灰度偏差，用grpc(a)表示；
172.(1)第一步：采用人工方法判别64个色区；
173.(2)第二步：人工在每个色区选择一个矩形区域；提取矩形区域每个像素的红色分量，绿色分量，蓝色分量；
174.(3)第三步：分别计算每个色区的的红色分量，绿色分量，蓝色分量平均值，用rgr(a),ggr(a),bgr(a)表示，a＝1、2、3、
……
、64，
175.(4)第四步：计算：grgr(a)＝rgr(a)
×
0.299 ggr(a)
×
0.587 bgr(a)
×
0.114，a＝1、2、3、
……
、64，
176.(5)第五步：计算灰度完整性，用wzgr表示,wzgr＝(grgr(64)-grgr(1))/252；
177.(6)第六步：计算各色区颜色分量校正值，红色分量用rjzgr(a)表示，绿色分量用gjzgr(a)表示，蓝色分量用bjzgr(a)表示，a＝1、2、3、
……
、64，
178.rjzgr(a)＝252
×
(rgr(a)
–
rgr(1))/(rgr(64)
–
rgr(1))
179.gjzgr(a)＝252
×
(ggr(a)
–
ggr(1))/(ggr(64)
–
ggr(1))
180.bjzgr(a)＝252
×
(bgr(a)
–
bgr(1))/(bgr(64)
–
bgr(1))
181.(7)第七步：计算各色区颜色分量偏差，红色分量用rpcgr(a)表示，绿色分量用gpcgr(a)表示，蓝色分量用bpcgr(a)表示，a＝1、2、3、
……
、64，
182.rpcgr(a)＝|rjzgr(a)-4(a-1)|；符号“||”表示取绝对值
183.(8)第八步：计算各色区灰度校正值，用grjz(a)表示，a＝1、2、3、
……
、64，
184.grjzgr(a)＝(grgr(a)
–
grgr r(1))/wzgr；
185.(9)第九步：计算各色区灰度偏差，用grpc(a)表示，a＝1、2、3、
……
、64，
186.grpc(a)＝|grjzgr(a)-4(a-1)|；符号“||”表示取绝对值。
187.色彩评价参数计算：
188.对于色度评测模具对应的图片，进行色彩评价参数计算；色彩评价参数指标有：色相偏差，亮度偏差，饱和度偏差。
189.色相h、亮度l、饱和度s计算：找出评测像素红色分量r、绿色分量g、蓝色分量b中间的最大值与最小值，设其为max、min，则rgb颜色空间转hsl颜色空间的计算方法如下:
190.(1)色相h的计算:
191.当max＝min时，h＝0；当max＝r，且g大于等于b时，h＝60
×
(g-b)/(max-min)；
192.当max＝r，且g小于b时，h＝360 60
×
(g-b)/(max-min)；
193.当max＝g时：h＝120 60
×
(b-r))/(max-min)；
194.当max＝b时：h＝240 60
×
(r-g))/(max-min)；
195.(2)亮度l计算：l＝(max min)/2
196.(3)饱和度s计算：
197.当l＝0或者max＝min时：s＝0；
198.当l大于零且小于等于0.5时：s＝0.5
×
(max-min)/l
199.当l大于0.5时：s＝0.5
×
(max-min)/(1-l)
200.色相偏差，亮度偏差，饱和度偏差计算方法：
201.(1)第一步：采用人工方法判别64个色区；
202.(2)第二步：人工在每个色区选择一个矩形区域；提取矩形区域每个像素的红色分量，绿色分量，蓝色分量；
203.(3)第三步：分别计算每个色区的的红色分量，绿色分量，蓝色分量平均值，用rgr(a),ggr(a),bgr(a)表示，a＝1、2、3、
……
、64；
204.(4)第四步：根据表1“评测区涂色”，计算1-64号色区；的标准色相、标准饱和度、标准亮度l
205.(5)根据rgr(a),ggr(a),bgr(a)，计算测量色相h、测量亮度l、测量饱和度s；(6)根据每个色区的测量色相与标准色相，测量饱和度与标准饱和度，测量亮度l-标准亮度l计算色相偏差、饱和度偏差、亮度偏差，色相偏差＝测量色相h-标准色相，饱和度偏差＝测量饱和度s-标准饱和度，亮度偏差＝测量亮度l-标准亮度l。
206.本发明的积极效果是：
207.1、提出一种用于评估空中三角测量纹理数据生成质量评估的设备，该设备可以模拟飞行器的飞行环境和照相机拍摄环境，以采集各种环境条件下照相机拍摄数据。
208.2、通过本发明方法，给出无人机倾斜摄影风速、关照、阴影等不同环境下纹理数据的成像质量。
附图说明
209.图1是本发明的环境因素对空三摄影影响实验装置结构示意图。
210.图2是实验装置的评测模具组合示意图。
211.图3是球形模的球体大小与各部位命名规则示意图。
212.图4是球形分辨率评测模具一个评测区分辨率评测图样示意图。
213.图5实验装置的光源结构示意图。
214.图6是反光面罩结构示意图。
215.图7反光面罩俯视图。
216.图8是可调灯座试验灯安装板俯视图。
217.图9是可调灯座调节杆安装口正视图。
218.图10是调节杆安装口侧视图。
219.图11是调节杆结构示意图。
220.图12是调节杆定位管结构示意图。
221.图13是调节杆紧固螺钉示意图。
222.图14是调节杆限定板俯视图。
223.图15是调节杆限定板正视图。
224.图16是调节杆与可调灯座之间的安装示意图。
225.图17是调节杆与反光面罩之间的安装示意图。
226.图18是调节杆与反光面罩之间的安装正视图。
227.图19是遮光板边框示意图。
228.图20是遮光板折叠支架示意图。
229.图21遮光板折叠支架结构示意图。
230.图22是实验装置的风源结构示意图。
231.图中，1倾斜摄影无人机，2光源，3风源，4评测模具组合，11球形分辨率评测模具，12球形色度评测模具，13正四面体分辨率评测模具，14正四面体色度评测模具，15正六面体分辨率评测模具，16正六面体色度评测模具，17正十二面体分辨率评测模具，18正十二面体色度评测模具，19正二十面体分辨率评测模具，20正二十面体色度评测模具，21球形灰度评测模具，22正四面体灰度评测模具，23正六面体灰度评测模具，24正十二面体灰度评测模具，25正二十面体灰度评测模具，31球体z线，32球体x线，33球体y线，34球体顶端，35球体底端，36球体右端，37球体左端，38球体后端，39球体前端，50三脚架，51反光面罩，52可调灯座，53试验灯，60罩体，61固定柄，62-1、62-2灯座调节杆，63-1、63-2遮光罩固定沿，64反光层，70焦点平面圆，71遮光罩固定沿内侧，72遮光罩固定沿外侧，73-1、73-2、73-3调节杆安装槽，74-1a、74-1b、74-1c、74-1d、74-2a、74-2b、74-2c、74-2d、74-31a、74-3b、74-3c、74-3d调节杆定位孔，75-1、75-2遮光板固定孔，80试验灯安装板，81-1、81-2、81-3调节杆安装口，82-1、82-2、82-3、82-4试验灯安装螺孔，85左侧安装支架，86右侧安装支架，87底部固定板，80、88-1、88-2定位固定孔，90调节杆，91调节杆柱状端，92调节杆扁平端，93扁平端固定孔，95调节杆定位管，96定位管本体，97定位管右侧限定座，98定位管左侧限定座，99调节杆紧固螺孔，100调节杆紧固螺钉，101紧固螺钉外丝，102紧固螺钉手柄，105-1右侧限定板固定孔、105-2左侧限定板固定孔，106-1限定板右端侧，106-2限定板左端侧，107限定板中间部分，110调节杆安装螺栓，111调节杆安装螺栓头部，112调节杆安装螺帽，73调节杆安装槽，106-3、106-4限定板端侧，116-1、116-2、116-3、116-4限定板安装螺栓头部，115-1、115-2限定板安装螺栓，117-1、117-2、限定板安装螺帽，120边框侧边，121-1、121-2边框支撑边，122-1、122-2支撑边固定孔，130折叠杆体，131-1、131-2折叠杆体连接孔；133-1～133-5折叠杆体邻间连接孔，150风机三脚架，151轴流风机，152风机调速器，153风速仪。
具体实施方式
232.参见附图。实验装置由评测模具组合4、倾斜摄影无人机1、光源2、风源3组成。
233.评测模具组合由十五个评测模具构成，含球形评测模具、正多面体评测模具；球形模具包括：球形分辨率评测模具11，球形色度评测模具12，球形灰度评测模具21；正多面体评测模具包括：正四面体分辨率评测模具13，正四面体色度评测模具14，正六面体分辨率评
测模具15，正六面体色度评测模具16，正十二面体分辨率评测模具17，正十二面体色度评测模具18，正二十面体分辨率评测模具19，正二十面体色度评测模具20；正四面体灰度评测模具22；正六面体灰度评测模具23；正十二面体灰度评测模具24；正二十面体灰度评测模具25。十五个评测模具采用任意方式摆放在一个平面。实施例采用图2分布方式。
234.图3给出了球体大小与各部位命名规则。球体半径用r表示；球体的外表面称为球面。
235.球体z线31：指穿过球形的垂直线；球体z线与球面有两个交点：与上边球面的交点称为球体顶端34；与下边球面的交点称为球体底端35。
236.设球体x线为左右方向；设球体y线为前后方向,z线为垂直方向。
237.球体x线与球面有两个交点：与左边球面的交点称为球体左端37；与右边球面的交点称为球体右端36；球体y线与球面有两个交点：与前边球面的交点称为球体前端39；与后边球面的交点称为球体后端38。
238.在球体顶端34下方，距离球体顶端0.28r，且与球体z线垂直的平面将球面分割为上下两部分，分割平面上边部分称为上端评测球面；
239.在球体底端35上方，距离球体底端0.28r，且与球体z线垂直的平面将球面分割为上下两部分，分割平面下边部分称为下端评测球面；
240.在球体左端37右方，距离球体左端0.28r，且与球体x线垂直的平面将球面分割为左右两部分，分割平面左边部分称为左端评测球面；
241.在球体右端36左方，距离球体右端0.28r，且与球体x线垂直的平面将球面分割为左右两部分，分割平面右边部分称为右端评测球面；
242.在球体前端39后方，距离球体前端0.28r，且与球体y线垂直的平面将球面分割为前后两部分，分割平面前边部分称为前端评测球面；
243.在球体后端38前方，距离球体后端0.28r，且与球体y线垂直的平面将球面分割为前后两部分，分割平面后边部分称为后端评测球面；
244.上端评测球面、下端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面面积、形状、涂色均相同，只是朝向不同。以上端评测球面进行阐述；
245.在上端评测球面，将球体z线垂直的底部小圆称为上端底部小圆。将上端底部小圆等分成64段等长的弧，共计64个分段点，称为等弧分段点。在上端测评球面上，连接球体顶端与64个等弧分段点的圆弧，将上端分辨率评测球面分割成64个扇面，将64个扇面称为64个区域，按顺时针方向将64个区域依次称为1号色区，2号色区
……
,64号色区。
246.图4是球形分辨率评测模具一个评测区分辨率评测图样示意图。
247.在64个区域上，分别涂色黑白相间的颜色。黑色颜色分量分别为：红色，0；绿色，0；蓝色，0。白色颜色分量分别为：红色，255；绿色，255；蓝色，255。
248.对于球形灰度评测模具，在各评测球区的64个色区涂灰度评测图样，灰度评测图样为灰色。球形色度评测模具在评测区涂色的64个色区，涂色度评测图样，色度评测图样为1-64号不同颜色，具体色区如表1。
249.表1
[0250][0251]
图4亦为正多面体分辨率评测模具一个内接圆评测区分辨率评测图样示意图。
[0252]
正多面体分辨率评测模具中，正四面体分辨率评测模具为正四面体、正六面体分辨率评测模具为正六面体、正十二面体分辨率评测模具为正十二面体、正二十面体分辨率评测模具为正二十面体，各个正多面体各面的内接圆涂有分辨率评测图样；将各个正多面体各面的内接圆等分为64个扇形，64个扇形分别涂色黑白相间的颜色；黑色颜色分量分别为：红色，0；绿色，0；蓝色，0；白色颜色分量分别为：红色，255；绿色，255；蓝色，255。
[0253]
正多面体色度评测模具中，正四面体色度评测模具为正四面体、正六面体色度评测模具为正六面体、正十二面体色度评测模具为正十二面体、正二十面体色度评测模具为正二十面体，各个正多面体各面的内接圆等分为64个扇形，64个扇形涂有色度评测图样,色
度评测图样为如表11-64号不同颜色。
[0254]
对于正四面体灰度评测模具为正四面体，正六面体色度评测模具为正六面体，正十二面体色度评测模具为正十二面体，正二十面体色度评测模具为正二十面体，正二十面体各面的内接圆等分为64个扇形，64个扇形涂有灰度评测图样，第n号色区，所涂颜色为：红色分量＝蓝色分量＝黄色分量＝4(n-1)。
[0255]
参见附图5～图15。光源2由三脚架50、反光面罩51、可调灯座52、试验灯53构成。三脚架固定反光面罩。灯座安装在反光面罩上，用于将实验灯放置在合适位置，并为实验灯提供电源以及电源连接座；反光面罩用于将实验灯的光反射，照到被拍摄的模具上；试验灯用于发出需要照度和色温的光，所发的光照射到被拍摄摸具上。
[0256]
本实施例三脚架采用现成产品：意大利曼富图摄影包公司生产，型号：mkbfrc4gtxp-bh试验灯采用现成产品：深圳市神牛摄影器材有限公司，型号：sl150ii bi升级版。
[0257]
反光面罩51形状为一个旋转抛物面，由受力性能好，不易变形的材料制作，实施例材料选择铝合金面板，内侧贴有摄影反光布。反光面罩51由罩体60、固定柄61、反光层(64)、遮光罩固定沿63-1、63-2和遮光板.反光层64为反光材料，贴在反光面罩(51)内侧；设反光面罩所在的旋转抛物面顶点为a，直线ao为反光面罩的旋转抛物面轴线；过轴线的平面与旋转抛物面相交线为抛物线，设该抛物线的焦点为o；过o点与直线ao垂直的平面与反光面罩的旋转抛物面相交线为圆，该圆称为焦点平面圆70，在焦点平面圆上，等距离分布三个调节杆安装槽73-1、73-2、73-3，相邻两个调节杆安装槽之间弧度为120度：每个调节杆安装槽两侧，有四个调节杆定位孔74-1a、74-1b、74-1c、74-1d、74-2a、74-2b、74-2c、74-2d、74-31a、74-3b、74-3c、74-3d；遮光罩固定沿63-1、63-2为圆环形，连接在反光面罩51底端，其所在平面与直线ao垂直；在遮光罩固定沿63-1、63-2上有两个遮光罩固定孔75-1、75-2，两个遮光罩固定孔的连线通过遮光罩固定沿的圆心。
[0258]
参见附图8。可调灯座52由试验灯安装板80，调节杆组成。
[0259]
试验灯安装板80安装试验灯，由力学性能好，易成形材料构成，实施例采用铝合金材料；做出圆盘状；在试验灯安装板中间位置，钻有若干个试验灯安装螺孔，用于安装试验灯；实施例钻有四个试验灯安装螺孔82-1、82-2、82-3、82-4；在试验灯安装板圆盘外侧，焊接三个调节杆安装口81-1、81-2、81-3。
[0260]
调节杆安装口如图9、10。调节杆安装口为u形，由左侧安装支架85、右侧安装支架86、底部固定板87构成；左侧安装支架85、右侧安装支架86、底部固定板为一个整体，实施例采用铝合金材料；左侧安装支架、右侧安装支架形状相同，为矩形和半圆形构成，半圆形一侧有定位固定孔88-1。
[0261]
调节杆结构如图11.调节杆由调节杆柱状端91和调节杆扁平端92构成，调节杆柱状端长度远远比调节杆扁平端长度长；调节杆柱状端与调节杆扁平端构成一个整体；在调节杆扁平状外侧中部，有一个扁平端固定孔。扁平端固定孔直径与定位固定孔直径相同；调节杆90装入调节杆定位管95定位安装。
[0262]
调节杆定位管95结构如图12。调节杆定位管95分为定位管本体96、定位管左侧限定座98、定位管右侧限定座97，由力学性能良好材料构成，本实施例采用铝合金。在定位管右侧限定座上端设置有调节杆紧固螺孔，调节杆紧固螺钉100将调节杆91紧固。调节杆紧固
螺钉100结构如图13。调节杆紧固螺钉由紧固螺钉外丝101与紧固螺钉手柄102组合而成，实施例采用铝合金制作。紧固螺钉外丝为外丝螺钉，与调节杆紧固螺孔匹配；紧固螺钉手柄用于旋转调节杆紧固螺钉。
[0263]
参见附图14、15。定位管左侧限定座98、定位管右侧限定座97安装在调节定位杆限定板内。调节定位杆限定板用于将调节杆定位管限制在反光面罩；两端为限定板端侧106-1、106-2；中间为限定板中间部分107；限定板中间部分为拱形，用于限制定位管右侧限定座、定位管左侧限定座的活动范围。调节定位杆限定板由力学性能良好的材料制作，本实施例采用不锈钢。
[0264]
见附图16.调节杆与可调灯座之间的安装：调节杆扁平端92插入到左侧安装支架85和右侧安装支架86中间，调节杆安装螺栓110依次穿过右侧安装支架定位固定孔105-1、扁平端固定孔93、左侧安装支架定位固定孔105-2后，用调节杆安装螺帽112和调节杆安装螺栓110拧紧，使得调节杆固定在调节杆安装口，且可以灵活转动。
[0265]
参见附图17、18.调节杆与反光面罩之间的安装：调节杆与反光面罩之间的安装时，调节杆91安装在定位管本体96内。
[0266]
调节杆定位管95的安装：调节定位杆限定板将调节杆定位管95的定位管右侧限定座97和定位管左侧限定座98固定在反光面罩上；固定时，定位管本体96嵌入调节杆安装槽，调节杆定位管的定位管右侧限定座97和定位管左侧限定座98嵌入调节定位杆限定板的限定板中间部分107，每个调节定位杆限定板两端的限定板固定孔105-1、105-2与同一侧的调节杆定位孔对齐；每个调节定位杆限定板均用两对螺栓螺帽固定；限定板安装螺栓穿过限定板固定孔、调节杆定位孔后，用螺帽拧紧，将调节定位杆限定板固定；调节杆定位管的定位管右侧限定座97和定位管左侧限定座98被调节定位杆限定板限制，能够绕着定位管右侧限定座97和定位管左侧限定座98的轴线转动。
[0267]
可调灯座与调节杆的安装：三个调节杆的调节杆扁平端92通过调节杆安装口固定在可调灯座上，三个调节杆的调节杆柱状端91分别穿过三个定位管本体96，被调节杆紧固螺钉固定在调节杆定位管上；调节杆紧固螺钉拧松的时候，调节杆的调节杆柱状端91在调节杆定位管内滑动；调节杆紧固螺钉拧紧的时候，调节杆的调节杆柱状端91与调节杆定位管固定。
[0268]
参见附图19-21.遮光板由遮光板边框和遮光板折叠支架构成；遮光板边框由一个边框侧边120、两个边框支撑边121-1、121-2、摄影用柔光布构成；边框侧边和边框支撑板均为长条形，长度相同，长度大于遮光罩固定沿的外径，两个边框支撑边在遮光板边框两端，与遮光板边框垂直连接；由力学性能良好的材料制作，实施例采用铝合金材料，在两个边框支撑边的中间均有支撑边固定孔122-1、122-2。摄影用柔光布裁剪为正方形，边长与边框侧边相同，其中三个边粘贴在边框侧边和两个边框支撑边，构成遮光板。折叠杆体为条状，采用力学性能良好的材料，实施例采用铝合金；折叠杆体两端有折叠杆体连接孔。
[0269]
遮光板折叠支架由若干跟折叠杆体130-1～130-6连接而成；相邻折叠杆体采用螺栓与螺帽紧固，133-1～133-5为相邻折叠杆体之间的折叠杆体邻间连接孔；螺栓穿过相邻折叠杆的折叠杆体邻间连接孔后，用螺帽紧固；最外的孔为折叠支架固定孔132-1、132-2，两个折叠支架固定孔一端对应支撑边固定孔122-1、122-2，一端对应遮光板固定孔75-1、75-2；遮光板折叠支架一端，螺栓穿过遮光板固定孔与折叠支架固定孔，用螺帽紧固，使得
遮光板折叠支架固定在遮光罩固定沿上；遮光板折叠支架另一端，螺栓穿过支撑边固定孔与折叠支架固定孔132-1、132-2，用螺帽紧固，使得遮光板折叠支架固定在边框支撑边120上；采用两个遮光板折叠支架，将两个边框支撑边连接到遮光罩固定沿。
[0270]
通过控制相邻折叠之间的夹角，控制折叠杆的长度；通过控制折叠杆长度，控制遮光板.遮盖光线的面接以及区域。
[0271]
参见附图22。风源由风机三脚架150、轴流风机151、风机调速器152、风速仪153构成。风机调速器连接电源，轴流风机电源连接线连接到风机调速器，风机调速器控制轴流风机的风速；轴流风机安装在风机三脚架上；风速仪测量轴流风机出口的风速。本实施例中，风机三脚架：意大利曼富图摄影包公司生产，型号：mkbfrc4gtxp-bh；轴流风机：艾美特电器(深圳)有限公司，型号：dpt10-20a；风机调速器：上海德力西开关有限公司：产品型号：dty-4000w；风速仪：深圳市句茂源科技有限公司，产品型号：gm8902 。
[0272]
本发明的评测模具中，正多面体评测模具设有内接圆，分辨率评测模具内接圆上的评测图样与球形分辨率评测模具评测图样一致，为黑白相间评测图样，其红色，0；绿色，0；蓝色，0；白色颜色分量分别为：红色，255；绿色，255；蓝色，255。色度评测图样与球形色度评测模具评测图样一致，为表1中所列64中评测图样。灰色评测模具内接圆上的评测图样与球形灰色评测模具评测图样一致，为灰色。
[0273]
本发明在实际评测时，无人机操控正摄影或倾斜摄影。正摄时：
[0274]
一机位：评测模具组合，外接圆圆心正上方；
[0275]
二无人机高度：拍摄图片刚好包含所有模具的高度；
[0276]
三无人机操控：人工操控；
[0277]
四：无人机离风源距离：保证无人机安全的最近距离；
[0278]
倾斜摄影时：
[0279]
一倾斜角度：45度
[0280]
二机位：评测模具组合各外接圆圆心与水平面45度方向；
[0281]
三无人机高度：拍摄图片刚好包含所有模具的高度；
[0282]
四无人机操控：人工操控；
[0283]
五：无人机离风源距离：保证无人机安全的最近距离；
[0284]
风源控制包括风速选择和轴流风机朝向；
[0285]
风速选择：选择任意风速。本实施例选择如下几种风速测试：
[0286]
风速1：0米/秒；风速2：1.5米/秒；风速3：3.3米/秒；风速4：5.4米/秒；
[0287]
风速5：7.9米/秒；风速6：10.7米/秒；风速7：13.8米/秒。
[0288]
轴流风机朝向：轴流风机轴向对准无人机；轴流风机轴向与水平面的角度设置为任意角度。实施例中，轴流风机轴向与水平面的角度分别为0度，45度，-45度。
[0289]
光源控制：光照角度、照度、光阴
[0290]
光照角度：反光面罩轴向对准模具外接圆圆心，反光面罩轴向与水平面的角度为任意角度；反光面罩距离调节为光线覆盖评测模具各组合外接圆区域；
[0291]
照度采用任意照度；
[0292]
光阴：光阴效果由遮光板产生；通过遮光板遮罩面积，产生任意光阴效果；
[0293]
所述无人机操控、风源控制、光源控制在测评是采用任意组合。
[0294]
实施例组合类型如下：
[0295]
一光源组合：
[0296]
光强组合：四种：最大照度；1/2照度；1/4照度；1/8照度；
[0297]
光阴组合：四种：遮光板采用：不遮罩；遮罩1/4；遮罩1/2；遮罩3/4；
[0298]
光照角度组合：二种：反光面罩轴向与水平面的角度为：45度；90度；
[0299]
光源采用上述光强、光阴、光照角度组合，总共组合形态：32种。
[0300]
二风源组合
[0301]
风速组合：七种：风速1：0米/秒；风速2：1.5米/秒；风速3 3.3米/秒；风速4：5.4米/秒；风速5：7.9米/秒；风速6：10.7米/秒；风速7：13.8米/秒。
[0302]
风向组合：三种：轴流风机轴向与水平面的角度分别为0度，45度，-45度。
[0303]
风源组合方式：21种。
[0304]
三摄影方式组合：两种：正射，倾斜
[0305]
四合计组合方式：32
×
21
×
2＝1344种
[0306]
在各种评价组合下拍摄数据，针对每种评价组合进行评测，评测参数计算过程包括数据预处理以、像素序列提取、黑白分辨率、灰度评价参数计算和色彩评价参数计算，其中，色彩评价参数计算含色相h、亮度l、饱和度s计算和色相偏差，亮度偏差，饱和度偏差计算；
[0307]
数据预处理：
[0308]
一评测区域
[0309]
对于球形评测模具：针对上端评测球面、左端评测球面、右端评测球面、前端评测球面、后端评测球面能够拍出完整评测图像的区域作为评测区域；
[0310]
对于正四面体、正六面体、正十二面体、正二十面体评测模具：针对能够拍摄完整评测图像的面作为评测区域。
[0311]
二数据预处理以平行于照片底边为其中一条边，做矩形，提取评测图像，矩形的四条边为评测图像的切线；设宽度方向包含的像素为n个，高度方向包含的像素为m个,且n大于等于m；用i表行，j表示列；i为整数，取值为1、2、3、
……
、m；j为整数，取值为1、2、3、
……
、n；
[0312]
图片左上角为第1行第1列；右下角为第m行第n列第i行，第j列的像素颜色表达为：红色分量：ri(i,j)；绿色分量：gi(i,j)；蓝色分量：bi(i,j)；
[0313]
如果m等于n：对评测图像不需进行转换；
[0314]
如果m小于n，按如下方法转换图像，转换后的图像，宽度方向为n个像素，高度方向为n个像素；k为整数，为1、2、3、
……
、n；
[0315]
转换后的图片左上角为第1行第1列；右下角为第n行第n列；
[0316]
转换后的第k行，第j列的像素颜色表达为：红色分量：ro(k,j)；绿色分量：go(k,j)；蓝色分量：bo(k,j)；
[0317]
对所提取的评测图像进行如下转换：
[0318]
(1)第一步：计算m除以(n-m)，设x为m除以(n-m)的整数部分；
[0319]
(2)第二步：q＝0；p＝0；r＝0；
[0320]
(3)第三步：q＝q 1；p＝p 1；
[0321]
ro(q,j)＝ri(q-r,j)；go(q-r,j)＝gi(q,j)；bo(q,j)＝bi(q-r,j)；
[0322]
j取值为1、2、3、
……
、n；
[0323]
(4)第四步：判断q是否等于n；等于n，进入第八步；不等于n，则进入第五步；
[0324]
(5)第五步：判断p是否等于x；p等于x，进入第六步；p不等于x，转入第三步；(6)第六步：p＝0；q＝q 1；r＝r 1；
[0325]
ro(q,j)＝(ri(q-1,j) ri(q,j))/2；
[0326]
go(q,j)＝(gi(q-1,j) gi(q,j))/2；
[0327]
bo(q,j)＝(bi(q-1,j) bi(q,j))/2；j取值为1、2、3、
……
、n；
[0328]
(7)第七步：判断q是否等于n；等于n，进入第八步；不等于n，则进入第三步；
[0329]
(8)第八步：转换结束。
[0330]
对转换后的彩色像素进行灰度计算，计算得到的图像用igr表示：
[0331]
gr(k，j)＝ro(k，j)
×
0.299 go(k，j)
×
0.587 bo(k，j)
×
0.114；
[0332]
计算得到的值取整；
[0333]
像素序列提取方法：设n除以2，取整为n50；n除以4，取整为n25；提取半径rn的像素序列,rn小于n50，行号用row表示，列号用col表示：像素序列用a表示，a＝1,2,3,4,5，
……
[0334]
(1)第一步a＝1；temp＝0；xa＝0；
[0335]
(2)第二步：col(a)＝n50
–
rn xa；
[0336]
(3)第三步：计算对vx取整数；
[0337]
(4)第四步：判断vx-temp是否大于1；小于等于1，row(a)＝n50-vx，进入第七步；大于1，则令t＝vx-temp，进入第五步；
[0338]
(5)第五步：row(a)＝n50-temp；a＝a 1；t＝t-1；
[0339]
(6)第六步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp 1；进入第五步；否则，进入第七步；
[0340]
(7)第七步：temp＝vx,a＝a 1；xa＝xa 1；
[0341]
(8)第八步：如果xa大于rn,则转到第九步；否则转到第二步；
[0342]
(9)第九步：col(a)＝n50
–
rn xa，计算对vx取整数；
[0343]
(10)第十步：判断temp-vx是否大于1；小于等于1，row(a)＝n50-vx，进入第十三步；大于1，则令t＝temp-vx，进入第十一步；
[0344]
(11)第十一步：row(a)＝n50-temp；a＝a 1；t＝t-1；
[0345]
(12)第十二步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp-1；进入第十一步；否则，进入第十三步；
[0346]
(13)第十三步：temp＝vx,a＝a 1；xa＝xa 1；
[0347]
(14)第十四步：如果xa大于2rn,则转到第十五步；否则转到第九步；
[0348]
(15)第十五步：xa＝xa-1
[0349]
(16)第十六步：col(a)＝n50
–
rn xa，计算对vx取整数；
[0350]
(17)第十七步：判断vx-temp是否大于1；小于等于1，row(a)＝n50 vx，进入第二十步；大于1，则令t＝vx-temp，进入第十八步；
[0351]
(18)第十八步：row(a)＝n50 temp；a＝a 1；t＝t-1；
[0352]
(19)第十九步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp 1；进入第十八步；否则，进入第二十步；
[0353]
(20)第二十步：temp＝vx,a＝a 1；xa＝xa-1
[0354]
(21)第二十一步：如果xa等于rn,则转到第二十二步；否则转到第十六步；
[0355]
(22)第二十二步：col(a)＝n50
–
rn xa，计算对vx取整数；
[0356]
(23)第二十三步：判断temp-vx是否大于1；小于等于1，row(a)＝n50- vx，进入第二十六步；大于1，则令t＝temp-vx，进入第二十四步；
[0357]
(24)第二十四步：row(a)＝n50 temp；a＝a 1；t＝t-1；
[0358]
(25)第二十五步：判断t是否大于1，大于1，col(a)＝n50
–
rn xa，temp＝temp-1；进入第二十四步；否则,进入第二十六步；
[0359]
(26)第二十六步：temp＝vx,a＝a 1；xa＝xa-1；
[0360]
(27)第二十七步：如果xa大于零,则转到第二十二步；否则转到第二十八步；
[0361]
(28)第二十八步：cnta＝a；程序结束；
[0362]
黑白分辨率：
[0363]
(1)第一步：rn＝80；
[0364]
(2)第二步：根据像素序列提取方法，计算像素序列的行列序列：col(a)，row(a)，计算cnta；cnta为序列的长度；
[0365]
(3)第三步：从灰度图像igr中，提取像素序列，像素序列用si(b)表示，b＝0、1、2、3、
……
、cnta；si(b)＝igr(row(b),col(b))
[0366]
(4)第四步：令his(c)＝0,c＝0、1、2、3、
……
、255；d＝1；
[0367]
(5)第五步：his(si(d))＝his(si(d)) 1；
[0368]
(6)第六步：d＝d 1
[0369]
(7)第七步：判断d是否小于等于cnta？是：进入第五步；否：进入第八步；
[0370]
(8)第八步：suml＝0，sumr＝0；cntl＝0；cntr＝255；
[0371]
(9)第九步：suml＝suml his(cntl)；sumr＝sumr his(cntr)；
[0372]
(10)第十步：判断suml是否大于sumr，是则cntr＝cntr-1，sumr＝sumr his(cntr)；否：
[0373]
则cntl＝cntl 1，suml＝suml his(cntl)；
[0374]
(11)第十一步：判断cntl是否大于cntr：是，则进入第十二步；否，进入第十步；
[0375]
(12)第十二步：d＝1,maxlv＝his(0)；maxl＝0；
[0376]
(13)第十三步：判断his(d)是否大于maxlv；是：maxlv＝his(d),maxl＝d；进入第十四步；否：进入第十四步；
[0377]
(14)第十四步：d＝d 1；
[0378]
(15)第十五步：判断d是否大于cntl，如果是，则进入第十六步；如果否，则进入第十三步；
[0379]
(16)第十六步：d＝254,maxrv＝his(255)；maxr＝255；
[0380]
(17)第十七步：判断his(d)是否大于maxrv；是：maxrv＝his(d),maxr＝d；进入第
十八步；否：进入第十八步；
[0381]
(18)第十八步：d＝d-1，进入第十九步；
[0382]
(19)第十九步：判断d是否大于cntl，否，则进入第二十步；是，则进入第十七步；
[0383]
(20)第二十步：d＝1；stepl＝0.5(maxl cntl)；stepr＝0.5(maxr cntl)；f＝0；cntinv＝0；(21)第二十一步：si(d)是否小于stepl？是：f＝-1进入第二十六步；否，进入第二十二步；
[0384]
(22)第二十二步：si(d)是否大于stepr？是：f＝1进入第二十三步；否，进入第二十一步；
[0385]
(23)第二十三步：d＝d 1
[0386]
(24)第二十四步：判断d是否大于cnta？是，进入第二十九步；否，进入第二十五步；
[0387]
(25)第二十五步：si(d)是否小于stepl？是：f＝-1，cntinv＝cntinv 1；进入第二十六步；否，进入第二十三步；
[0388]
(26)第二十六步：d＝d 1；
[0389]
(27)第二十七步：判断d是否大于cnta？是，进入第二十九步；否，进入第二十八步；
[0390]
(28)第二十八步：si(d)是否大于stepr？如果是：f＝1，cntinv＝cntinv 1；进入第二十三步；如果否，进入第二十六步；
[0391]
(29)第二十九步：判断cntinv是否大于62；如果是进入第三十步；如果否，则rn＝rn 1，进入第二步；
[0392]
(30)第三十步：黑白分辨率为rn对应的像素值；
[0393]
(31)第三十一步：程序结束。
[0394]
灰度评价参数计算:
[0395]
对于灰度评测摸具对应的图片，进行灰度评价参数计算；灰度评价参数指标有：灰度完整性，用wzgr表示；各色区颜色分量偏差，红色分量用rpcgr(a)表示，绿色分量用gpcgr(a)表示，蓝色分量用bpcgr(a)表示；色区灰度偏差，用grpc(a)表示；
[0396]
(1)第一步：采用人工方法判别64个色区；
[0397]
(2)第二步：人工在每个色区选择一个矩形区域；提取矩形区域每个像素的红色分量，绿色分量，蓝色分量；
[0398]
(3)第三步：分别计算每个色区的的红色分量，绿色分量，蓝色分量平均值，用rgr(a),ggr(a),bgr(a)表示，a＝1、2、3、
……
、64，
[0399]
(4)第四步：计算grgr(a)＝rgr(a)
×
0.299 ggr(a)
×
0.587 bgr(a)
×
0.114，a＝1、2、3、
…
64，
[0400]
(5)第五步：计算灰度完整性，用wzgr表示,wzgr＝(grgr(64)-grgr(1))/252；
[0401]
(6)第六步：计算各色区颜色分量校正值，红色分量用rjzgr(a)表示，绿色分量用gjzgr(a)表示，蓝色分量用bjzgr(a)表示，a＝1、2、3、
……
、64，
[0402]
rjzgr(a)＝252
×
(rgr(a)
–
rgr(1))/(rgr(64)
–
rgr(1))；
[0403]
gjzgr(a)＝252
×
(ggr(a)
–
ggr(1))/(ggr(64)
–
ggr(1))；
[0404]
bjzgr(a)＝252
×
(bgr(a)
–
bgr(1))/(bgr(64)
–
bgr(1))；
[0405]
(7)第七步：计算各色区颜色分量偏差，红色分量用rpcgr(a)表示，绿色分量用gpcgr(a)表示，蓝色分量用bpcgr(a)表示，a＝1、2、3、
……
、64，
[0406]
rpcgr(a)＝|rjzgr(a)-4(a-1)|；符号“||”表示取绝对值
[0407]
(8)第八步：计算各色区灰度校正值，用grjz(a)表示，a＝1、2、3、
……
、64，
[0408]
grjzgr(a)＝(grgr(a)
–
grgr r(1))/wzgr；
[0409]
(9)第九步：计算各色区灰度偏差，用grpc(a)表示，a＝1、2、3、
……
、64，
[0410]
grpc(a)＝|grjzgr(a)-4(a-1)|；符号“||”表示取绝对值。
[0411]
色彩评价参数计算：
[0412]
对于色度评测模具对应的图片，进行色彩评价参数计算；色彩评价参数指标有：色相偏差，亮度偏差，饱和度偏差；
[0413]
色相h、亮度l、饱和度s计算：找出评测像素红色分量r、绿色分量g、蓝色分量b中间的最大值与最小值，设其为max、min，则rgb颜色空间转hsl颜色空间的计算方法如下:
[0414]
(1)色相h的计算:
[0415]
当max＝min时，h＝0；当max＝r，且g大于等于b时，h＝60
×
(g-b)/(max-min)；
[0416]
当max＝r，且g小于b时，h＝360 60
×
(g-b)/(max-min)；
[0417]
当max＝g时：h＝120 60
×
(b-r))/(max-min)；
[0418]
当max＝b时：h＝240 60
×
(r-g))/(max-min)；
[0419]
(2)亮度l计算：l＝(max min)/2
[0420]
(3)饱和度s计算
[0421]
当l＝0或者max＝min时：s＝0；
[0422]
当l大于零且小于等于0.5时：s＝0.5
×
(max-min)/l
[0423]
当l大于0.5时：s＝0.5
×
(max-min)/(1-l)
[0424]
色相偏差，亮度偏差，饱和度偏差计算方法：
[0425]
(1)第一步：采用人工方法判别64个色区；
[0426]
(2)第二步：人工在每个色区选择一个矩形区域；提取矩形区域每个像素的红色分量，绿色分量，蓝色分量；
[0427]
(3)第三步：分别计算每个色区的的红色分量，绿色分量，蓝色分量平均值，用rgr(a),ggr(a),bgr(a)表示，a＝1、2、3、
……
、64；
[0428]
(4)第四步：根据表1“评测区涂色”，计算1-64号色区；的标准色相、标准饱和度、标准亮度l
[0429]
(5)根据rgr(a),ggr(a),bgr(a)，计算测量色相h、测量亮度l、测量饱和度s；(6)根据每个色区的测量色相与标准色相，测量饱和度与标准饱和度，测量亮度l-标准亮度l计算色相偏差、饱和度偏差、亮度偏差，
[0430]
色相偏差＝测量色相h-标准色相，
[0431]
饱和度偏差＝测量饱和度s-标准饱和度，
[0432]
亮度偏差＝测量亮度l-标准亮度l。