基于二次函数拟合的太阳影子定位技术研究..docxVIP

基于二次函数拟合的太阳影子定位技术研究马玉雪，令狐林玉，杨东海，孙国伟，买阿丽( 运城学院应用数学系，山西运城 044000)摘 要: 利用影子长度和太阳高度角的关系研究太阳影子定位技术，得到太阳高度角为 45 度时影子和物体一样长; 并通过对太阳影子长度变化的数据进行分析，得出影子长度与时间的二次函数拟合方程，从而根据影长最短的时间点确定拍摄地点的经度; 然后根据太阳方位角与时角，经纬度，高度角之间的关系，求出高度角，从而得到影子和纬度的关系; 最终确定出视频的拍摄地点和拍摄日期。关键词: 太阳高度角; 太阳方位角; 二次函数拟合; 影子定位中图分类号: TP391． 7 文献标志码: A 文章编号: 1008-8008( 2015) 06-0024-05DOI:10.15967/14-1316/g4.2015.06.0070． 引言太阳影子定位技术就是通过分析视频中物体的太阳影子变化，通过测量影子方位及长度与杆的实际长度、经纬度、时间等参数，确定视频的拍摄地点和拍摄日期的一种方法。而物体影子的变化主要源于地球自转和绕太阳的公转，直杆位置及长度的测量是利用 CCD 相机线性成像的交比不变性原理和灭点产生原理，所以通过分析影响影子变化的参数，就可确定视频拍摄的地点和日期。对于太阳影子定位技术考虑以下问题:第一，建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律。第二，当日期已知，直杆高度未知时，根据水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。第三，当日期与直杆高度都未知时，根据水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。第四，根据一根 2 米高的直杆在太阳下的影子变化的视频，建立确定视频拍摄地点的数学模型，同时给出可能的拍摄地点。当拍摄日期未知，根据数学模型确定可能的拍摄地点与日期。收稿日期: 2015-08-161． 直杆影子长度变化曲线1． 1 问题分析影响影子长度的主要因素是杆子高度和高度角。首先确立影响太阳高度角［1］变化的因素并确定高度角的计算公式，然后利用影子长度和太阳高度角的三角函数关系，得到影子长度变化规律。以天安门广场北纬 39 度 54 分 26 秒，东经 116 度 23 分 29 秒位置( 文中记为 P 点) ，3 米直杆在 2015 年 10 月 22 日北京时间 9∶ 00 － 15∶ 00 之间的太阳影子长度的变化曲线为例，说明这些量之间的变换关系。 1． 2 确定影子长度与太阳高度角的关系首先，我们确定太阳高度角［2］θ，收集资料可知太阳高度角θ的变化与太阳赤纬［3］δ、纬度和太阳时角［4］w 三个因素有关:( 1) 太阳赤纬δ与日期 T 有关，其关系式如下: sinδ = 0． 39795cos［0． 98563( T － 173) ］ ( 1) 其中 T = 1 表示1 月1 日，T = 365 表示12 月31日。( 2) 纬度，其数值在 0 至 90 度之间。( 3) 太阳时角 w 是指日面中心的时角，可以由北京时间得到，设北京时间为 t，以 24 小时计，经度为φ，可以得到关系式［3］:基金项目: 山西省高等学校科技项目( ; 运城学院科研基金项目( XK －2014030，XK －2014035) 作者简介: 马玉雪( 1994-) ，女，山西潞城人，运城学院应用数学系 1301 班学生。w = 15° × ［t + ( 120° － φ) ÷ 15° － 12］ ( 2)根据图 1，确定出了这三个因素关系式为［2］:sinθ = sinφsinδ + cosφcosδcosw ( 3) 最后建立函数模型，确立影子长度与太阳高度角θ的函数关系，如图2 所示，设直杆高度为 h 米，影子长度为 l 米，可以得到函数模型为:htanθ =( 4)l图 1 太阳高度角计算示意图图 2 影子长度与太阳高度角的关系图1． 3 画出变化曲线图应用建立的模型，确立了 P 点 3 米高的直杆在2015 年10 月22 日北京时间9: 00 － 15: 00 之间的太阳影子长度随时间变化的二次函数关系，其中由计算可得φ = 116． 39， = 39． 9072，h = 3，T = 265，从而用Matlab 画出了北京时间9: 00 － 15: 00 之间该地的太阳影子长度的变化曲线，如图 3 所示。分析可知，如果为东经 120 度时，北京时间 12 点时影子长度最短，又因为每晚 4 分钟，经度偏西 1 度，而天安门广场经度为东经 116 度 23 分 29 秒，比东经 120 度少 3． 61 度，所以该地的时间要比标准北京时间要晚 14． 44 分钟，所以天安门广场影子最短的时间为 12 点 14 分 26 秒。