欢迎访问喜蛋文章网
你的位置:首页 > 经典文章 > 文章正文

高分一号影像的像元的太阳天顶角如何算

时间: 2021-07-16 23:57:19 | 来源: 喜蛋文章网 | 编辑: admin | 阅读: 99次

高分一号影像的像元的太阳天顶角如何算

太阳天顶角是什么

180

天顶角的定义是什么?谢谢

地球绕太阳公转,由于地轴的倾斜,地轴与轨道平面始终保持着大概66`34'的夹角,这样,才引起太阳直射点在南北纬23`26’之间往返移动,并决定了太阳可能直射的范围:春,秋分日,太阳直射赤道---即直射点的纬度为0`;冬至日,太阳直射南回归线--即直射点的纬度为23`26’S;夏至日,太阳直射北回归线--即直射点的纬度为23`26’N。如果某地的纬度已经知道,依据下面的公式就可以计算出此地的正午太阳高度角的大小:
H=90`-纬差(纬差是指某地的地理纬度与当日直射点所在纬度之间的差值)
例如:A地的纬度为40`N,求A地夏至日的正午太阳高度。
夏至日太阳直射的纬度为23`26’N,与A地的纬度差=40`-23`26’=16`34’,那么H=90`-16`34’=73`26’。
如果求A地冬至日的正午太阳高度。冬至日太阳直射的纬度为23`26’S,与A地的纬度差=40`+23`26’=63`26’,那么H=90`-63`26’=26`34’。

有根据某地方的经纬度和时间计算月球天顶角方位角的具体公式吗?

感谢楼主的信任:
计算的基础请详细阅读下面的文章,特别是无论如何也必须记下由赤道坐标向地平坐标的换算(后面的计算过程还会用到黄道坐标向赤道坐标的运算)。
http://wenku.baidu.com/link?url=PxcaPTDzTLsesQuECFQjpUrFAEbD7Z5JLm1WyFz6xkb0pz65eGdDkkbyR3qpOCtDAFtsiyFZ4YYlMXOtNIh_8fqeHQFJpJEmNZkmYNumTbG

好,现在我们知道了换算的公式了,需要的三个参量是:1.观测点的纬度Φ 2.天体的时角t 3.天体的赤纬δ。
第一个数字很容易得到,观测点的纬度是固定的,所以可以看作一个常数。但接下来两个,麻烦大了,先说相对简单一些的赤纬:
根据一个简化的模型,月球在黄道上匀速运动,那我们可以根据月亮的黄经推出月亮的赤经赤纬。月亮的黄经可以根据太阳的黄经和月亮和太阳的相对位置来确定:
举一个例子,比如春分上弦,太阳的黄经是0,月球的黄经是90°,所以可以推算月球的黄经是90°,对应6h,赤纬23.5°。
从例子中我们可以看到,我们需要的是公历日期来推算太阳的位置,然后再根据月球的相对太阳的位置(也就是需要知道上一次朔的时间)来推算月球相对于太阳的位置,再通过黄道和赤道的坐标的换算来推算月球的赤经赤纬。
但是,如果只是这样处理,会有比较大的误差(可以达到10°),第一个是月球在白道上运动,黄白交角5°9′;第二个是月球的轨道并不是正圆,存在一定的偏心率,所以在黄道上不是匀速运动。前者可以通过查阅月球过黄白交点的时间来校正,后者可以通过考察月球过近地点的时间来校正。我想楼主到这里肯定是听得云里雾里了吧?就不再展开了。

好的,最后一个参数,时角。其实在前面的计算中,得到的赤纬的数值之外,同样得到了赤经的数值。
在文章中已经提到,时角等于=恒星时-天体的赤经。而恒星时由日期决定,换而言之我们可以通过查阅日期对应的恒星时,扣除赤经来得到天体的时角。

现在三个参数都有了,代入一开始提到的公式就能计算出赤经赤纬了。我在这里介绍的,大约是东周时期对月球位置计算的精度(如果考虑了提到的两点误差差不多到了东汉的水准),我想对于绝大多数的朋友早就是每个字都认识的不明觉厉系列了吧?
现代天文学对于月球的运动考察了数以百计的参数,如果能说出其中十个以上,那就是功力了得的爱好者。趁这个机会,向天文学家表示由衷的敬意。

太阳天顶角的计算公式

有谁知道如何在已知经度,纬度,时间的情况下,计算太阳天顶角。

太阳高度角和太阳天顶角互为余角。

1日地距离

  地球绕太阳公转的轨道是椭圆形的,太阳位于椭圆两焦点中的一个。发自太阳到达地球表面的辐射能量与日地间距离的平方成反比,因此,一个准确的日地距离值R就变得十分重要了。日地平均距离R0,又称天文单位,

1天文单位=1.496×108km

或者,更准确地讲等于149597890±500km。日地距离的最小值(或称近日点)为0.983天文单位,其日期大约在1月3日;而其最大值(或称远日点)为1.017天文单位,日期大约在7月4日。地球处于日地平均距离的日期为4月4日和10月5日。

  由于日地距离对于任何一年的任何一天都是精确已知的,所以这个距离可用一个数学表达式表述。为了避免日地距离用具体长度计量单位表示过于冗长,一般均以其与日地平均距离比值的平方表示,即ER=(r/r0)2,也有的表达式用的是其倒数,即r0/r,这并无实质区别,只是在使用时,需要注意,不可混淆。

  我们得到的数学表达式为

  ER=1.000423+0.032359sinθ+0.000086sin2θ-0.008349cosθ+0.000115cos2θ(1)

式中θ称日角,即

θ=2πt/365.2422

(2)

这里t又由两部分组成,即

t=N-N0

(3)

式中N为积日,所谓积日,就是日期在年内的顺序号,例如,1月1日其积日为1,平年12月31日的积日为365,闰年则为366,等等。

  N0=79.6764+0.2422×(年份-1985)

-INT〔(年份-1985)/4〕

(4)

2太阳赤纬角

  地球绕太阳公转的轨道平面称黄道面,而地球的自转轴称极轴。极轴与黄道面不是垂直相交,而是呈66.5°角,并且这个角度在公转中始终维持不变。正是由于这一原因形成了每日中午时刻太阳高度的不同,以及随之而来的四季的变迁。太阳高度的变化可以从图1中形象地看到。图中日地中心的连线与赤道面间的夹角每天(实际上是每一瞬间)均处在变化之中,这个角度称为太阳赤纬角。它在春分和秋分时刻等于零,而在夏至和冬至时刻有极值,分别为正负23.442°。

 由于太阳赤纬角在周年运动中任何时刻的具体值都是严格已知的,所以它(ED)也可以用与式(1)相类似的表达式表述,即:

  ED=0.3723+23.2567sinθ+0.1149sin2θ

-0.1712sin3θ-0.758cosθ+0.3656cos2θ

+0.0201cos3θ(5)

式中θ的含义与式(1)中的相同。

3时差

  真正的太阳在黄道上的运动不是匀速的,而是时快时慢,因此,真太阳日的长短也就各不相同。但人们的实际生活需要一种均匀不变的时间单位,这就需要寻找一个假想的太阳,它以均匀的速度在运行。这个假想的太阳就称为平太阳,其周日的持续时间称平太阳日,由此而来的小时称为平太阳时。
  平太阳时S是基本均匀的时间计量系统,与人们的生活息息相关。由于平太阳是假想的,因而无法实际观测它,但它可以间接地从真太阳时S⊙求得,反之,也可以由平太阳时来求真太阳时。为此,需要一个差值来表达二者的关系,这个差值就是时差,以Et表示,即

S⊙=S+Et(6)

  由于真太阳的周年视运动是不均匀的,因此,时差也随时都在变化着,但与地点无关,一年当中有4次为零,并有4次达到极大。时差也可以以式(1)相似的表达式表示:

  Et=0.0028-1.9857sinθ+9.9059sin2θ

-7.0924cosθ-0.6882cos2θ

(7)

  上面,我们给出了3个计算式,从形式上讲,它们与一般书籍中给出的并无不同。我们之所以又重新研究它,是因为以往的公式存在以下的通病:①对平年和闰年不加区分,一方面,这对闰年就不好处理,另一方面,闰年的影响有累计效应,会逐步增长;②即使是从当年天文年历查到的数值,也是格林尼治经度处0点时刻的数值,而我们所需要的数值,会因所在地点的地理经度以及具体时刻与表值有异而不同。具体地讲,一般要进行如下3项订正:
  (1)年度订正:除非我们只用当年的天文年历值,此外均需使用此项订正,引入此项订正的原因就是一回归年的实际长度不是365日,而是365.2422日,但日历上只有整日,不可能有小数日。假定我们选用的是1981年的表值,1982年再用时,就要加上-0.2(-0.2422)日的订正了。这个订正到了1983年为-0.51(-0.4844)日,1984年为-0.7(-0.7266)日,但此年为闰年,多了1日,实际订正应为-0.7+1=0.3(0.2734)日,1985年为0.0(0.0312)日,等等,余类推。
  (2)经度订正:即使我们查阅的是当年的天文年历,也需此项订正。在我国的地理经度范围内,各地的订正值是

≤90°E

-0.2日

>90°E~<128°E

-0.3日

≥128°E

-0.4日

  (3)时刻订正:要求同前一项。即使在格林尼治当地,不同时刻也需加以订正。各时段的订正值是:

时段  336-600  600-824  824-1048  1048-1312

日  +0.2+0.3+0.4+0.5

时段  1312-1536  1536-1800  1800-2024

日  +0.6+0.7+0.8

  由于我国普遍采用的是北京时,它与格林尼治的地方时相差8小时,故具体到我国情况:

时段(北京时)200-424  424-648  648-912  912-1136

订正值(日)-0.2  -0.1  0  0.1

时段  1136-1400  1400-1624  1624-1848  1848-2112

订正值  0.2  0.3  0.4  0.5

  前面3个计算式,项数多计算麻烦,后面多项订正,更显繁琐。为了方便实际应用,特编制如下仅含20句的BASIC语言程序,供使用:

  10  input“经度,经分和年份”,JD,JF,NF

  20  A=NF/4:K=2*3.1415926#/365.2422

  30  N0=79.6764+0.2422*(NF-1985)

      -INT((NF-1985)/4)

  40  input“月,日,时,分(按北京时)”,Y,R,S,F

  50  B=A-INT(A)

  60  C=32.8

  70  ifY≤2thenC=30.6

  80  ifB=0andY>2thenC=31.8

  90  G=INT(30.6*Y-C+0.5)+R

  100  L=(JD+JF/60)/15

  110  H=S-8+F/60

  120  N=G+(H-L)/24

  130=(N-N0)/K

  140  式(1)

  150  式(5)

  160  式(7)

  170  print“Er=”;Er;“Ed=”;Ed,“Et=”;Et

  180  input“是否仍要计算y/n?”,W0

  190  ifW=“Y”orW=“y”then10else200

  200  end

  程序中50-90各句的目的在于计算当天的积日,100句是经度订正,110句是时刻订正,130句包含3年度订正的内容。
  在太阳能利用中,最常见的是要求计算太阳高度和太阳方位。
  太阳高度(h⊙)的计算公式为

sinh⊙=sinδsinφ+cosδcosφcosτ(8)

式中,δ就是太阳赤纬角,即式(5)中的Ed,φ为当地的地理纬度,τ为当时的太阳时角。φ值不难获得,且一旦确定,不会改变。δ值的计算可以从前述程序中得到。唯一需要说明的是太阳时角的计算。其计算式为

(9)

  这里时S和分F的符号均加上了⊙下标,表示是真太阳时,为了从北京时求出真太阳时,需要两个步骤:首先,将北京时换成地方时Sd:

(10)

式中,120°是北京时的标准经度,乘4是将角度转化成时间,即每度相当于4分钟,除60是将分钟化成小时。

  其次,进行时差订正,即

S⊙=Sd+Et/60                     (11)

这里应该指出的是,时角是以太阳正午时刻为0点的,顺时针方向(下午)为正,反之为负。

  太阳方位角的计算式为

  cosA=(sinh⊙sinφ-sinδ)/cosh⊙cosφ

(12)

由此可求出二个A值,第一个A值是午后的太阳方位,

  当cosA≤0时90°≤A≤180°

  当cosA≥0时0≤A≤90°

第2个A值为午前的太阳方位,取360°-A。
  实例:计算东经110°北回归线上1999年6月23日北京时12∶42的太阳高度角及当日的日落时的方位角。

  计算:将JD=110,JF=0,NF=1999,Y=6,R=23,S=12,F=42,各参数输入运行中的程序;屏幕上立即显示:Er=1.0330,Ed=23.438,Et=-1.84

  将北京时12∶42换算成东经110°的地方时,利用式(10),可得Sd=12∶02

  加当日时差Et≈-2,得此时当地的S⊙=12∶00,将其代入式(9)得τ=0°,北回归线处φ=23.442°

  最后根据式(8)求得h⊙=89.966°

  读者可能产生疑问,为何在北回归线上,夏至日的中午时刻的太阳高度不等于90°,大家不妨变换NF的输入值,看一看结果不仅都不等于90°,且各年之间还略有差异。之所以会如此,是因为夏至不仅有日期,还有时刻,很难遇到夏至时刻在正午是12时的。
  在计算日落时的方位角时,由于此时h⊙=0,所以式(12)的形式有所变化:

cosA=-sinδ/cosφ

(13)

将已知参数代入,得cosA=-0.3977
  依照判据90°≤A≤180°,故A=113.44°

地球绕太阳公转,由于地轴的倾斜,地轴与轨道平面始终保持着大概66`34'的夹角,这样,才引起太阳直射点在南北纬23`26’之间往返移动,并决定了太阳可能直射的范围:春,秋分日,太阳直射赤道---即直射点的纬度为0`;冬至日,太阳直射南回归线--即直射点的纬度为23`26’S;夏至日,太阳直射北回归线--即直射点的纬度为23`26’N。如果某地的纬度已经知道,依据下面的公式就可以计算出此地的正午太阳高度角的大小:
H=90`-纬差(纬差是指某地的地理纬度与当日直射点所在纬度之间的差值)
例如:A地的纬度为40`N,求A地夏至日的正午太阳高度。
夏至日太阳直射的纬度为23`26’N,与A地的纬度差=40`-23`26’=16`34’,那么H=90`-16`34’=73`26’。
如果求A地冬至日的正午太阳高度。冬至日太阳直射的纬度为23`26’S,与A地的纬度差=40`+23`26’=63`26’,那么H=90`-63`26’=26`34’。
地球绕太阳公转,由于地轴的倾斜,地轴与轨道平面始终保持着大概66`34'的夹角,这样,才引起太阳直射点在南北纬23`26’之间往返移动,并决定了太阳可能直射的范围:春,秋分日,太阳直射赤道---即直射点的纬度为0`;冬至日,太阳直射南回归线--即直射点的纬度为23`26’S;夏至日,太阳直射北回归线--即直射点的纬度为23`26’N。如果某地的纬度已经知道,依据下面的公式就可以计算出此地的正午太阳高度角的大小:
H=90`-纬差(纬差是指某地的地理纬度与当日直射点所在纬度之间的差值)
例如:A地的纬度为40`N,求A地夏至日的正午太阳高度。
夏至日太阳直射的纬度为23`26’N,与A地的纬度差=40`-23`26’=16`34’,那么H=90`-16`34’=73`26’。
如果求A地冬至日的正午太阳高度。冬至日太阳直射的纬度为23`26’S,与A地的纬度差=40`+23`26’=63`26’,那么H=90`-63`26’=26`34’。
H=90`-|a-b|
http://www.ers.ac.cn/IRSISS/1a5fileformat.htm
文章标题: 高分一号影像的像元的太阳天顶角如何算
文章地址: http://www.xdqxjxc.cn/jingdianwenzhang/115065.html
文章标签:天顶  高分  影像  太阳
Top