第二课时 纬线和纬度 经线和经度 利用经纬网定位
学习目标：1、学会在地球仪上识别经线、纬线、两极、地轴、本初子午线和赤道。
2、通过学生观察地球仪，使学生知道经线和纬线，经度和纬度的特点。
3. 能熟练地利用经纬网确定某一地点的地理位置
学习重难点：经纬线的特征和经纬度的分布特征。利用经纬网确定某一地点的地理位置。
学习过程：
一、自主学习
学习任务一：阅读教材《经线和纬线》图文内容并结合地球仪，掌握经纬线的特点及经纬度的标度规律。
纬线和纬度
1． 什么是纬线？纬线是什么形状的？
2． 各纬线长度的变化有什么规律？
3． 找出最长的一条纬线，这条纬线把地球分为哪两个半球？
4． 纬度数值的排列有什么规律？(学生在图中标出主要纬线的度数)
5． 纬度最大值的纬线有什么特征？该最大值出现在哪个地区？
6． 请同学描画出0°、30°、60°三条纬线，找一位学生在图上标出低纬中纬高纬三个区域。
7． 在地球仪上找出全部位于北半球的大洲、大部分位于南半球的大洲和地跨赤道两侧的大洲。
8． 在地球仪上找出分别找出位于低纬、中纬、高纬三个区域的国家各一个。
经线和经度
1、 经线在形状、提示方向、长度的分布三个方面与纬线有什么不同？
2、 经度以哪条经线为起点的？
3、 经度的标度范围是多少度？
4、 对比纬度的变化规律，说明地球仪上经度的变化规律有什么特点。
5、 东西半球是如何划分的？这样划分的意义是什么？
6、 在地球仪上找出主要位于东半球和西半球的大洲，并说出它们的名称。
用途、如何利用经纬网确定位置。）
基于经纬度的物体太阳光线投影长度的研究
基于经纬度的物体太阳光线投影长度的研究
摘要
本文研究的问题是如何通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期。
建模整体思路是，先建立一系列分析用到的物理量，设定一些假设和约束条件，使得问题求解有可行性，之后对这些物理量进行演绎。
建模使用的软件平台主要是matlab，分析用到的主要参量是太阳赤纬、时角、高度角、方位角、纬度，分析过程当中用到的方法有，建立物理概念，明确物理意义，比如引用天球坐标系的概念，在天球坐标系的基础上进行物理分析，通过对建立的参变量进行物理关系的推导，形成公式体系进行求解，对题目所给予的影子坐标数据进行适当变换处理，使用matlab进行合理的拟合，对于用公式法和方程法没法顺利解决的问题使用穷举法作为解题的补充，对于视频中坐标的取法用到了坐标转换的思想。
cossin其中主要公式有 1.sinA cosh
H2.tanh L
3. cosAsinhsinsin coshcos
4. sinh=cosΩcosφcosδ+sinφsinδ
第一问，通过物理量变换，先求出高度角，进而得到影子长度与时间变化关系。 第二问，拟合点求经度，取点套公式求纬度。
第三问，方程思想，过程复杂，采用穷举法近似实现求解。
第四问，难点在于通过视频分析，得到影子端点的变化坐标，进而将问题转化成第二问，已知日期（太阳赤纬），时间（时角），求解经度纬度。
关键词：天球坐标系 物理量演绎分析 matlab数据拟合分析 二元方程组近似穷举法 坐标转换思想
1.问题重述与分析
如何确定视频的拍摄地点和拍摄日期是视频数据分析的重要方面，太阳影子定位技术就是通过分析视频中物体的太阳影子变化，确定视频拍摄的地点和日期的一种方法。
1. 建立影子长度变化的数学模型，分析影子长度关于各个参数的变化规律，并应
用你们建立的模型画出2015年10月22日北京时间9:00-15:00之间天安门广
场（北纬39度54分26秒,东经116度23分29秒）3米高的直杆的太阳影子长
度的变化曲线。
分析：模型的参数有经度（地方时），纬度，日期（太阳赤纬）
如果能够根据这三个变量建立相关模型，则地球上任意地点任意时刻的物体影子的形状和方位都能够确定
2. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直
杆所处的地点。将你们的模型应用于附件1的影子顶点坐标数据，给出若干个
可能的地点。
分析：这属于一个模型的逆过程，根据已经得到的影子的轨迹形状、日期来推断地点
3. 根据某固定直杆在水平地面上的太阳影子顶点坐标数据，建立数学模型确定直
杆所处的地点和日期。将你们的模型分别应用于附件2和附件3的影子顶点坐
标数据，给出若干个可能的地点与日期。
分析：第三问与第二问的不同在于第二问有具体的日期，而第三问中并没有具体的日期这就为求解带来了一定的不确定性和难度
4．
（1）附件4为一根直杆在太阳下的影子变化的视频，并且已通过某种方式估计出直杆的高度为2米。请建立确定视频拍摄地点的数学模型，并应用你们的模型给出若干个可能的拍摄地点。
（2）如果拍摄日期未知，你能否根据视频确定出拍摄地点与日期？
分析：根据视频提取某一时刻的影子的长度，视角之间的转换关系，方向的确定都是值得分析的地方
2.模型约定与假设
本文采用如下假定：
1．太阳光线视为平行光
2．研究地面上的杆子的时候地面视为平的
3．一年365天 一天24h 南北回归线纬度为23°26′
4.本文采用天球坐标系
5.宏观上设地球为光滑标准球体 不考虑大气层折射影响
6．数据中的时间在处理的时候都应
该处理成当地的地方时
7.地球的自转和公转都是匀速运
动，其中，公转为圆形轨道
3.总体模型及相关概念建立 1）地平坐标系（天球坐标系之一）
过观测者O (天球中心)的铅垂
线﹐延伸后与天球交于两点﹐朝上
的一点Z 称为天顶﹐朝下的一点
Z ’称为天底（右下图）
过天顶Z 和天体作一垂直圈﹐
它与地平圈交于垂足D 点﹐则天体 在地平坐标系中的第一坐标就是大圆弧D 或极距 Z 。D =h 称为地平纬度﹐又称地平高度﹐简称高度﹔而Z= 称为天顶距。地平高度也可以用平面角OD 来量度﹐而天顶距也可以用平面角OZ 来量度。天球上与地平圈相平行的小圆称为地平纬圈﹐也称平行圈。同一地平纬圈上任意点的地平高度都是相同的﹐因此可以称为等高圈。南点S 与垂足D 之间的大圆弧SD =a ﹐是地平坐标系中的第二坐
标﹐称为地平经度或天文方位角﹐简称方位角。
方
位角也可以用平面角SOD 来量度﹐天文学中习惯
从南点起按顺时针方向量度。以地平圈为基圈﹑子
午圈为主圈﹑南点为主点的坐标系称为地平坐标
系。由于周日视运动﹐天体的地平坐标不断发生变
化。另一方面﹐对不同的观测者﹐由于铅垂线方向
的不同﹐就有不同的地平坐标系﹐同一天体也就
有不同的地平坐标。 【1】
2）首先为了定义阐释诸多物理量，我们可以建立
一个以观测者为原点的空间直角坐标系，将地球视
为一个完美的球体，在观察者的位置我们可以做一
个切面，记为H，y轴经过此面，y轴亦即东西方向，
而在此平面内的过观察者的与y轴垂直的直线方向
即正南正北，正北方向可以记为p，XOY平面为与赤
道面平行的平面，太阳光的方向用l来表示，l与
xoy平面的交角记为δ，而O-Z轴的方向，正是地球
极轴指向。由地理概念易知，Ω（ω）是时角。以
上都是为了求解太阳的高度角和方位角设置的参变
量。 【2】
3）δ（太阳赤纬）的计算
选取12月22日为基准点，从这一天起到以后某一天，地球在其公转轨道上上走过一个转角，这个转角的大小记为α。其中由地理学知识容易知道，12月22日这一天太阳赤纬为南纬23°26′
将360° 365等分，易得α的计算公式α=0.9863（d2-d1）其中d1是12月22日的日期数，d2-d1的整体含义是待求日期与冬至日的日期差。 【3】
4）查资料得太阳赤纬δ
sinδ=0.39795cos[0.98563(n-173)]（n为此时与1月1日距离的天数） 【4】
5）时角
由地理学知识，Ω=15*T+k（其中T为格林尼治时间，二十四小时制，k为观测点的时间，东经为正，西经为负） 【5】
6）太阳高度角
由图片易知，太阳高度角即为向量n和向量l的夹角的余角，地理意义上即为太阳光线与地平面的夹角。
如图平面过y轴，平面与oz轴的夹角为φ，由解析几何的知识可以知道平面H在空间
中的单位法向量的表达式为n=（cosφ,0,sinφ）,根据太阳的直射纬度（太阳赤纬）
以及视角Ω可以得到单位向量l=（cosΩcosδ,sinΩsinδ,sinδ），设太阳高度角为
n,lh由余弦公式可以知道sinh=cos,综合可以知道
sinh=cosΩcosφcosδ+sinφsinδ
【6】
7)方位角的定义：太阳光的单位向量在地平面上的投影线与正南方向的射线，按照从正南方向顺时针到投影线的顺序所构成的角成为方位角。逆时针为负数,顺时针为正. 太阳的方位角A公式经过推导可知是
cossinsinA cosh
cosAsinhsinsin 【7】 coshcos
若规定东西方向为x方向则，南北方向为y方向，则tanA=x/y
从而,太阳下，竹竿影子的端点关于一系列参数的方程就得到了，为接下来的一些工作奠定了基础。
4．具体模型建立与求解
4.1第一问建模与求解
分析：第一问当中，确定太阳长度的变化曲线即确定太阳影子长度和时间的函数关系，利用高度角、杆长、影子长度的关系，将求解影子长度转化成求高度角，通过太阳赤纬，地理纬度换算，时角换算来得到太阳高度角随时间的变化关系，进而求出影长的变化，同时给出了，影子变化的实际轨迹图，非常直观。
如图，在被投影物体长度一定的情况下，影子长度的变化主要取决于夹角A的值，而夹角A与太阳高度角B相等，于是，决定影子长度的唯一因素就是太阳高度角。而太阳高度角主要由三个因素影响。1.太阳赤纬。2.地理纬度。3.地方时。
三个因素的求解方法：
1.太阳赤纬。太阳赤纬的变化是由于地球公转所引起的，因此，将公转的一周365等分，并得到此时与12月22日的差值天数，可以计算出此时的公转角=0.9863d,在通过三角关系，可得到太阳赤纬与公转角之间的关系为：sin=-sin23.5cos=-0.39795cos,以上是我们的初步想法。而为了方便起见，我
们引用他人的资料得：sin=0.39795cos[0.98563(n-173)]（n为此时与1月1日距离的天数）
2.地理纬度。给定的坐标，无须计算。
3.地方时与北京时间的换算。由于经度的不同，每个地区的地方时有差异，而北京时间较为普遍，因此，我们需要通过给定的北京时间来进行当地地方时的换算。地球自转平均1度需要4分钟，因此，换算成地方时只需要计算与120度地区相差的度数即可。通过地方时，我们可以确定时角，即以12时为0时角，每小时变化360/24=15度。
已知上述的三个因素，我们可以通过空间向量中夹角公式，得出太阳高度角的表达式：
sinh = costcosφcosδ + sinφsinδ（h为太阳高度角，t为时角，δ为太阳赤纬，φ为地理纬度），通过太阳高度角，我们可以进而通过三角函数求解影子的长度。
通过上述分析，我们可以用matlab实现此问题，并得到实现此类问题的程序，以题
目中要求的数据为例，
1. 我们计算太阳赤纬，根据公式需要用到当前日期与元旦的差值
2. 换算纬度
3. 根据经度求与北京的时间差
4. 求解地方时
5. 换算成时角
6. 求解太阳高度角
7. 计算影子长度
8. 做出图线
程序参见程序1
得到的图线如图所示：
注意，此处的时间是北京时间，而并不是地方时间，所以对于地方时来说，时间并不是关于正午对称的，所以得到的图像并不对称，但是显然，图像的最低点对应的时间正是以北京时间来表示的当地正午时间，为12：14左右。
更进一步，可以得到影子的轨迹图，如下：
竹竿影子的轨迹图
4.2第二问建模与求解
分析：第二问的处理思路是将经度、纬度的求解分开。对于经度，采用matlab进行数据多项式拟合的方法，得出二次曲线，进而得到当地正午的北京时间，通过两个正午时间差的换算，得到当地的经度。对于纬度，通过高度角的余弦公式求解。这里如果运用方程组的思想的话，只需要一个点（正午以外）即可解决，因为加上其他物理量之间的制约关系，刚好可以解决。
1． 经度的求算：
首先,
），再根据长度与时间的点，进行拟合，由于函数曲线先减后增，考虑二次函数拟合，用matlab拟合得到多项式为：
y=0.1489*x^2-3.7519*x+24.1275
现在，将此多项式处理，求得最小值时x的值（即正午的北京时间）;
12.5987
于是，在北京时间12.5987时，当地达到正午。
通过时差换算，可算出当地经度：
(12-12.5987)*60/4+120=111.0195
即经度为东经111.02度
详见程序2、3、4
2. 纬度的求算：
首先计算太阳赤纬 a=asin(0.39795*cos(0.98563*(n-173)))
之后从给出的数据中取出一组
x=1.7337 y=0.6013;
得到方位角
A=atan(x/y)
算出此时的时角
t=(15+36/60-12.5987)*15*pi/180
进而带入公式求得高度角的余弦，进而求出能够使公式成立的纬度
φ=30.57
即北纬30.57度
综上，可以得到最终答案 即东经111.0195°北纬30.57°
4.3 第三问的建模与求解
分析：第三问，两组数据较第二问而言，缺少了时间，即太阳赤纬的没法直接求算，这里如果用方程组的思想，从给出的数据当中取两个点（非正午，非关于正午对称）即可解决，而由于得出的方程比较复杂，所以我们采取另外一种简单有效，精度稍低的方式——穷举法，对纬度和日期（太阳赤纬）进行穷举，从当中挑选出最接近正确答案的一组，由此得到答案。
建模的思路如下
先求经度
1. 将数据导入matlab进行拟合（二次拟合）
根据横纵坐标求长度，对时间进行换算
2. 对时间和影子长度两个参数进行拟合求出二次函数的系数
3. 求出影子最短的时候对应的时间和长度
4. 通过时差换算得到当地的经度 (12-15.2141)*60/4+120=71.7885 即E71.7885°
同理得到附件三中经度为E108.9405°
接下来进行纬度和日期的计算，纬度和日期可以根据太阳方位角和高度角公式，计算得到。首先考虑
cossin1. sinA cosh
2. sinhcoscoscossinsin
然而以上计算方法较为复杂且容易出错，此处考虑用穷举法来计算。
由于一年有365天，因此从元旦开始，考虑365个数，进行穷举循环，从而计算出太阳赤炜，如下：
d=asin(0.39795*cos(0.98563*(i-173)*pi/180)) （i为1到365的整数）
由于纬度范围有限（南北90度），考虑北纬从0-pi/2之间9000个值。
取一组数据
T=13:29 x=-0.8215 y=0.4001 来进行穷举计算。
此时的时角为 ((13+29/60)-15.2141)*15*pi/180
方位角可以由坐标值得到，为arctan（-0.8215/0.4001）
因此，可以用如下方法。
对从i=1:365的每一天，从角度为零到九十度开始扫描，寻求能够使目标表达式 sqrt(1-(cos(a)*cos(d)*cos(-0.4531)+sin(a)*sin(d))^2)-cos(d)*sin(-0.4531)/(-0.8990)最小的解。
得到的结果
x =64 282 y =599 599
同理，第三组数据也可有如下方法算得。
以这一组数据为标准计算
T=13:39 X=1.7589 Y=3.2907
x =95 251 y = 3257 3257
即第三问附件二
时间为3月5日或10月9日，坐标为东经71.79度，北纬5.99度
同理，也可能为9月5日或4月9日，坐标为东经71.79度，南纬5.99度
第三问附件三
时间为4月5日或9月8日，坐标为东经108.94度，北纬32.57度
同理也可以为10月5日或3月8日，坐标为东经108.94度，南纬32.57度。
程序参见程序5、6、7
4.4第四问的建模和求解
第四问是前面几问的综合应用，而第四问与前面几问的区别在于，第四问需要从视频中提取坐标信息，进而归结成第三问的问题。为了缩小研究范围，从视频中显示的植被，日期，建筑风格，人种等等可以判断应该是北半球的亚洲国家。至于从视频中提取坐标等信息的方法，限于时间和知识，在此不能完美解答。
总之，第四问在思路上就是通过视角，坐标系的转换，将图中信息转化成前面的坐标等，进而求解。
5.误差分析
误差