太阳影子定位

1、太阳影子定位摘要：太阳影子定位技术就是通过分析物体影子变化，确定物体所在经纬度的技术。本文根据几何关系，结合穷举算法和相似度算法对给定大致经纬度范围内的物体进行精确经纬度定位的研究。我们首先通过太阳高度角与各个参数之间的关系，确立了影子长度与经纬度、日期、时间以及杆高之间的函数关系。接着，我们探讨建立了日期已知的经纬度定位模型，优化了之前的影长变化模型，将其结合了遍历算法。先通过影子顶点变化曲线特性分析缩小遍历范围,再遍历了不同杆长所有可能的经纬度。利用欧式距离算法衡量影长与遍历计算由的影长之间的相似度，欧氏距离最小的地点为最可能的直杆所处地点。最后，我们探讨了日期未知下的经纬度定位模型。对之

2、前的模型进行了进一步的优化：引入影子方位角变化量作为描述影子方位随时间变化特性的参数，与影长结合共同描述影子变化特性。以影长相似度最高作为目标进行一次遍历，接着以影子方位角变化量相似度最高作为目标，对一次遍历结果进行二次遍历，以求得与其相似度最高的地点的经纬度和对应的日期，并用欧式距离检验。关键词：经纬度、遍历算法、欧氏距离DOI:10.16640/ki.37-1222/t.2016.09.2181影长变化模型的建立我们设直杆影子长度为1,直杆长度为Zo太阳高度角是太阳射到地面的平行光线与地面的夹角，设其为h,范围为(0°,90°)。因为太阳光射到地面为平行光，又假设地面平

3、整，固定直杆垂直于地面，则有影长与直杆长度z和不同时刻的太阳高度角h的三角关系。其中，h是太阳高度角，为当地纬度，规定北纬度数为正，南纬度数为负，所以范围为(-90,90°),6为太阳赤纬角，即太阳直射点所在的纬度，范围为(-23°26',23°26'),Q为太阳时角。由图1可以看由，太阳赤纬角就是太阳直射点所在的纬度。太阳在地球上的直射点的纬度变化以年为周期，直射点在北纬23°26'到南纬23°26'的范围内移动，太阳赤纬由日期确定，夏至时太阳照射在北回归线，即北纬2326',冬至时太阳照射在南回归线，

4、即南纬2326'。一天之中赤纬角的变化可忽略不计。因此固定日期的赤纬角6可以当做常数。赤纬角的计算公式可以查阅资料(资料名称见参考文献2)得到，如下：8=0.006918-0.399912Xcos(B)+0.0070257Xsin(B)-0.006758Xcos(2B)+0.000907Xsin(20)-0.002697Xcos(30)+0.00148xsin(30)(1-3)其中，B为日角，n为日数，从1月1日开始起n=1,比如1月3日n=3o因此只要已知日期(年月日)，就可以算由具体的赤纬角值。由图2可看由，太阳时角Q是指从地球上的观测点天球子午圈沿天赤道量至太阳所在时圈的角距离。

5、子午圈和天赤道量属于地理学中赤道坐标系的概念，在本文中不过多展开。解释一下太阳时角就是任一时刻太阳与直杆的连线在赤道面上的投影与观测点正南方向的夹角。设该地方的经度为a,西经为负值，东经为正值，a范围为(-180°,180°)。设该地方的北京时间为t(当前时刻，单位为小时，例：12点36分即12.6小时)，对于观测点来说，不同时刻太阳相对于自己的方位也不一样。这导致太阳时角的具体值也不同。所以太阳时角Q与北京时间t和当地经度a有关。查阅文献2可知太阳时角Q的计算公式：其中S表示当地的真太阳时，即当地太阳位于正南向的瞬时，即当地的正午。综上公式1-1、1-2、1-3,如果已知

6、奥地经度a和北京时间t,就可以求由此时刻的太阳时角Q;如果已知该地的日期，就可以求生该地的太阳赤位角8,如果再知道该地纬度,就可以求生太阳高度角h,结合杆高z可以求生影长lo通过公式1-1、1-2、1-3可得建立影长变化模型：2有日期下经纬度定位模型的建立由于只要给由任意杆长就可以遍历所有经纬度得到最有可能的地点，为了计算便捷，我们只给由杆长为2米和3米两种情况下，最有可能的经纬度。2.1 影长的计算若已知条件中直接给由了影长则可以直接进行step2,若影长给由的是坐标(x,y),则可以算由影长，然后将经纬度范围内的每个经纬点这个时刻的影长也算由来与之对比即可。2.2 遍历范围的缩小由于遍历算

7、法计算量巨大，所以我们考虑结合题目给由的北京时间数据和影子顶点变化曲线的特性来缩小遍历的经纬度范围。(1)估算大致经度范围。根据日夜状态与一定的背景信息先行缩小经度范围。(2)估算大致纬度范围。查阅地理学文献3可以知道对于处在不同纬度的直杆，一天当中的影长轨迹可能有两种情况。以北半球的夏季为例，北半球的夏季太阳直射点在北半球(赤道与北回归线之间)，那么位于直射点以南的地区(包括南半球)，影子轨迹就如左图；位于直射点以北的地区，一天的影子轨迹就如右图。2.3 将经纬度范围代入遍历程序我们考虑先在较大经纬度范围内进行粗遍历，先找由最可能点所在的大致范围，再缩小经纬度范围，缩小步长，进行精确定位。(

8、1)缩小可能的经纬度范围：由于初始遍历的经纬范围比较大，为防止程序计算量过大导致计算时间过长，我们可以先选取较大的经纬度变化步长，。(2)精确确定经纬度：在得到粗略的经纬度之后可以缩小步长为0.01°,求生较为精确的地理坐标。3无日期下经纬度定位模型的建立上一个模型已经分析，只有已知奥地的日期、纬度、经度a、北京时间t和杆高Z可以求生影长1,在无日期下，如果沿用之前的模型，则需要遍历确定日期和杆长两个量。3.1 沿用之前模型计算不同日期下的经纬度坐标我们在之前的遍历模型外加上一层遍历365天的循环,然后重新遍历。在进过试验后可以看由来，沿用之前的模型仍能求生一年365天每天的最有可能

9、的经纬度位置。但是无法确定附件给生的具体日期。因此我们加入地理学中太阳方位角的概念来确定日期。一天中物体影子的变化由影长变化和影子方位(和正南方向的夹角)变化共同确定，通过上述遍历法求得的结果很多且相似度很高是由于只考虑了影长随时间的变化的相似度，没有考虑影子方位随时间变化的相似度。上述模型中缺少描述影子方位变化参数。由此，我们引入极坐标系，来讨论影子方位角的变化值，以此确定影子方位变化，对模型进行优化。3.2 日期的确定(1)极坐标系的建立。我们考虑将给由影长坐标的平面直角坐标系转化为极坐标系。以直杆底端为原点，平面直角坐标系的x轴正方向为极轴正方向，取逆时针方向为角度的正方向，建立极坐标系

10、：与原平面直角坐标系转化关系如下：(3-2)其中p=1(影长)，e为影子方向与极轴正方向所成的夹角。我们考虑利用上述极坐标系求生附件2中各个时刻的e(2)太阳方位角A的引入。固定直杆底端为极坐标系的原点，变化极轴的方向，特别地，若以正南方向为极轴正方向，建立极坐标系。此时，有e=a,a为太阳方位角(直杆影子与正南方向的夹角)，太阳方位角A与太阳高度角h间关系如下：(3-2)上式中，6为太阳赤纬，Q为太阳时角，计算方法之前已经给由所以，通过求解得到的每一天相似度最高地点的经纬度、杆高、日期和观测时间，即可算由每一天相似度最高的地点，任意观测时刻的太阳方位角。但是，由于建立极坐标系时极轴方向的不同

11、，通过计算得到太阳方位角A与利用附件原始平面直角坐标系转化得到的极坐标系中的8是不同的（当且仅当极轴方向为正南方向时e=a）。因此，我们考虑引入影子与极轴正方向夹角e关于时间的变化量ae作为描述影子方位随时间变化特性的参数，定义ae为影子方位角变化量。我们定义ae为影长坐标给由的观测后一时刻的e与前一时刻e的差值，由上述分析，太阳方位角a即为极轴方向为直杆正南方向时的e,所以应当有ae=AA。实际计算时，由于ae与aa不可能严格相等，我们考虑沿用欧式距离算法作为判定ae与aa相似度的算法，对求由来的一年365天每天的最有可能的经纬度进行二次穷举，算由这365个经纬地点方位角的变化量AA与附件所给的影子与极轴正方向夹角A9的欧氏距离，筛选由最小的欧氏距离，得到经纬度。此时的经纬度即最有可能的经纬度地点，这个经纬度对应的日数即可推算由当天日期。第二次欧氏距离计算公式如下：(3-2)遍历365个日期对应的经纬度，求取A9与AA相