1、1关于太阳影子定位的模型研究摘要太阳影子定位技术是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点与日期的一种方法。如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面。而本文主要研究的就是通过太阳影子的变化来确定物体地点和日期的问题。针对问题一，要确定直杆太阳影子长度的变化曲线，本文用太阳高度角、赤纬角、时角、经度和纬度之间的关系建立函数模型，利用 MATLAB 拟合出影子长度随时间的变化曲线。针对问题二，要求根据直杆影子顶点坐标求出直杆所在地的位置，根据经纬度等变量之间的关系，建立函数模型： 2)3015cos()51.0cos()51.0sin()( hlTfxg s通过 MAT
2、LAB 中的 fsolve 方法求解出直杆所在地的经纬度 )3150.9,5.108(可知直杆所在地为海南省昌化港。并利用 MATLAB 拟合得到了直杆影长与时间之间的函数，进而验证了该模型的合理性。针对问题三，要求根据直杆影子顶点坐标求出直杆所在地的位置与日期，我们对问题二的模型进行反推得到附件 2 直杆所在经纬度 ，为 6 月 1 日的新疆阿克苏；附件三直杆所)7398.40,96.8(2在经纬度 ，为 4 月 24 日的蒙古赛音山达附近。)65.43,019.(针对问题四，要求确定视频的拍摄位置与日期，我们用 CAD 制图软件处理截取图像中相对杆长与相对影长，结合问题二中模型得到视频的拍
3、摄地点为山西吕梁或内蒙古包头。对于拍摄时间未知的问题，利用类似方法得到直杆所在地为 6 月 30 日的内蒙古鄂尔多斯或 6 月 29 日的内蒙古包头。关键词：太阳影子定位； 最小二乘法；太阳高度角；赤纬角 1、问题提出如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面。太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。根据太阳影子定位技术建立模型，解决下列问题：1.建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应用建立的模型画出 2015 年 10 月 22 日北京时间9:00-15:00 之间天安门广场（北纬 东经 ）3 米高
4、“265439“291直杆的太阳影子长度的变化曲线。2.根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点。将模型应用于附件 1 的影子顶点坐标数据中，给出若干个可能的地点。3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直杆所处的地点和日期。将模型分别应用于附件 2和附件 3 的影子顶点坐标数据，给出若干个可能的地点与日期。4附件 4 为一根高度大约 2 米的直杆在太阳下的影子变化的视频。建立确定视频拍摄地点的数学模型，并用模型给出若干个可能3的拍摄地点。如果拍摄日期未知，能否根据视频确定出拍摄地点与日期。二、问题分析针对问题一，要求建立
5、影子长度变化的模型，分析影子长度关于各个参数的变化规律，通过应用太阳高度角、赤纬角、时角、经度和纬度之间的关系建立模型，利用 MATLAB 拟合出影子长度随时间变化的曲线。针对问题二，要求根据直杆影子顶点坐标求出直杆所在地的位置。对于所求问题，我们建立关于高度角、赤纬角、经纬度等变量之间的函数模型，采用 MATLAB 中 fsolve 的方法求解出直杆可能所在的位置。接着我们用 MATLAB 最小二乘法拟合曲线验证了该模型的合理性。针对问题三，要求根据直杆影子顶点坐标求出直杆所在地的位置与日期。我们在问题二模型的基础上，重新设置变量建立了新的函数模型，用 MATLAB 求解出直杆可能所在位置与
6、日期。针对问题四，要求根据一根已知高度的直杆在太阳下影子变化的视频确定视频的拍摄位置与日期。我们截取了 14 个时刻的视频图片，用 CAD 制图软件测出直杆的相对影长与杆长(见附件 10)，并建立函数模型，得到视频拍摄的位置与日期。三、模型假设1.忽略地球大气层折射效应对太阳光的影响。42.忽略视觉对杆长与影长数据收集的影响。3.忽略视频拍摄角度对相对影长与相对杆长的影响。4、定义与符号说明符号 含义H太阳高度角纬度太阳赤纬角时角n时间序号（从 1 月 1 开始计时)T所求地方时l直杆长h影长sT已知地方时f所求地点经度t时间5、模型的建立与求解5.1 问题一模型的建立与求解5图 1-1 赤纬
7、角示意图 图 1-2 高度角示意图对于赤纬角 2365)284(sin45.23n（1）由题意知，直杆影长变化的时间是 2015 年北京时间 10 月 22 日，可取 ，由 (1) 得295n 107.2对于地方时 3)12(5T（2）地球上经度每相隔 ，地方时相差 4 分钟，故有地方时 与已1 T知地方时 的关系sT4分 钟经 度 差 sT6（3）直杆所在地点是北纬 39 度 54 分 26 秒,东经 116 度 23 分 29 秒，所以可以由（2）得出从 9:00 到 15:00 各个时刻的 。对于太阳高度角 2coscossinsinH（4）直杆影长 、杆长 、太阳高度角 几何关系 :h
8、lHlhcot（5） 然后，我们利用 MATLAB 软件得到了直杆影长在 中随时间变化的曲线，如图 1-3 所示（源文件见附表0:15:9（1）。7图 1-3 影长随时间的变化曲线5.2 问题二模型的建立与求解5.21 模型的建立与求解由题意得测量日期是 2015 年 4 月 18 日，可算出日期序号，由（1）得：08n 51.0设直杆长 ，影长 ，而 与影子顶点坐标 有：lhyx,2yx（6）结合 (5)得：822tan1sinhlH（7）由(2)、(3)可得：3015sTf（8）结合公式（4）得到：)3015cos(cossin2 sTfhl （9）下面我们假设函数： 2)3015cos(
9、)51.0cos()51.0sin()( hlTfxg s（10） 而后在附件 1 中选取六组数据分别代入（10），利用 MATLAB 中最小二乘法 来求出杆长 、经度 与纬度 ，从而确定直杆所6 lf处的地点。所得情况如图 2-1(源程序见附表 (2)）。序号 北京时间 X 坐标（米） Y 坐标（米） 影长（米） 经/纬度（度）14:42 1.0365 0.4973 1.14962582615:00 1.2448 0.5311 1.353364049115:30 1.6438 0.5892 1.74620591108.5990/19.315014:42 1.0365 0.4973 1.149
10、62582615:09 1.3568 0.5483 1.463399853215:21 1.516 0.5715 1.620144515108.7485/19.176914:45 1.0699 0.5029 1.18219897615:15 1.4349 0.5598 1.540231817315:39 1.7801 0.6074 1.880875001108.6433/19.27054 14:48 1.1038 0.5085 1.215296955 108.5369/19.3582915:00 1.2448 0.5311 1.35336404915:24 1.5577 0.5774 1.66
11、127061314:51 1.1383 0.5142 1.24905105215:12 1.3955 0.5541 1.501481622515:42 1.8277 0.6135 1.927918447108.7285/19.219214:54 1.1732 0.5198 1.2831953415:18 1.4751 0.5657 1.579853316615:33 1.6882 0.5952 1.790050915108.3289/19.5196图 2-1每组数据得到的经纬度数值基本一致。所以附件 1 直杆的经纬度坐标可取 ，使)3150.9,5.108(用谷歌地球得知，此地处于海南省昌化港(图 2-2)。图 2-25.22 模型的验证根据附件 1 中某直杆在水平地面上的太阳影子的顶点坐标，可以通过 Excel 软件求出直杆在在个时间段的影子长度(图 2-3)。10图 2-3通过使用最小二乘法 ，由 MATLAB 拟合出从 14:42 到 15:42 里1影长随时间的变化曲线，如图 2-4（源程序见附表(3)）。
