第1个回答 2012-09-27
有公式，根据纬度可以计算出来
计算原理
一、已知天体赤道坐标，求黄道坐标。
赤经：α；
赤纬：δ；
黄赤交角：ε；
黄经：λ；
黄纬：β；
计算公式如下：
sinβ=cosε*sinδ-sinε*cosδ*sinα
cosβ*cosλ=cosδ*cosα
cosβ*sinλ=sinε*sinδ+cosε*cosδ*sinα
二、根据黄道坐标求空间直角坐标。
空间直角坐标的定义：以天球球心为原点，以春分点的方向为X轴的正向（黄经0度，黄纬90度），以夏至点的方向为Y轴的正向（黄经0度，黄纬90度），以北黄极的方向为Z轴的正向（黄纬90度）。
假设天球半径为R，计算公式如下：
X=R*cosλ*cosβ
Y=R*sinλ*cosβ
Z=R*sinβ
三、计算天体在椭圆轨道上运行时的任意时刻的坐标。
在下图中，红色椭圆为天体的运行轨迹（位于XY坐标平面内），箭头所指为天体的运行方向（即作逆时针方向的运行）。蓝色点P为任意时刻的天体，它的坐标为（x，y）。C点为椭圆的一个焦点，中心天体位于椭圆的这个焦点上，则A点为近点，B点为远点。O点为椭圆的中心。将坐标的原点定在中心天体的位置上，天体P从近点A开始，沿椭圆轨道围绕中心天体做逆时针方向的公转，在任意时刻t走到了P位置。
已知：
椭圆半长轴：a（图中OA）；
椭圆偏心率：e（e=c/a，c为焦点到中心距离，即OC）；
公转周期：T；
任意时刻：t（应小于T，>=0，从A点开始计时）；
在时刻t处的天体坐标：x，y；
算法如下：
设n=2*π/T；（π为圆周率）
设M=n*t；
设E（0）=M；
然后进行递归迭代运算
E（1）=M+e*sinE（0）；
E（2）=M+e*sinE（1）；
E（3）=M+e*sinE（2）；
……
E（N）=M+e*sinE（N-1）；
N越大，表明迭代的次数越多，精度越高，越接近理论真实值，但永远也达不到理论真实值，取适当的次数即可，然后计算x，y的值。
四、计算太阳的高度与方位角。
太阳高度：h；
太阳方位：A；
地理纬度：φ；
太阳赤纬：δ；
真太阳时：t；
计算公式如下：（单位为度）
h=arcsin（sinφ*sinδ-cosφ*cosδ*cos（15*t））；
A=arccos（tgh*tgφ-sinδ/（cosh*cosφ））；（午前取负，午后取正）
当在两极时：φ=±90度，h=±δ，A不存在；
当在赤道上，且太阳赤纬为0：φ=0度，午前h=15*t-90，A=-90；午后h=90-（t-12）*15，A=90；
当在正午时：t=12，h=90-│φ-δ│，A=0或180；
当在子夜时：t=0，h=│φ+δ│-90，A=0或180；
五、计算真太阳时与地方方标准时。
真太阳时：T⊙；
平太阳时：TM；
地方标准时：T；
地理经度：λ；
时区中心经度：λ。；
时差：η
计算公式如下：
T⊙=TM+η；
TM=T+（λ-λ。）/15；
T⊙=T+（λ-λ。）/15+η；
T=T⊙-η-（λ-λ。）/15；
六、计算晨光始、日出、日落、昏影终、昼长、晨昏蒙影长、黑夜长、日出日落太阳方位角。
地理纬度：φ；
太阳赤纬：δ；
日出日落太阳方位角：A；
晨光始时刻：TM；
日出时刻：TR；
日落时刻：TD；
昏影终时刻：TE；
昼长：LD；
晨昏蒙影长：LS；
黑夜长：LN；
太阳上边沿与地平线相切时的太阳高度：hr；
晨昏蒙影开始或结束时的太阳高度：hb；
计算公式如下：（所有时间均为真太阳进，单位为小时）
TM=arccos（tgφ*tgδ-sinhb/（cosφ*cosδ））/15；
TR=arccos（tgφ*tgδ）；（不考虑空气折射和太阳视半径）
TR=arccos（tgφ*tgδ-sinhr/（cosφ*cosδ））/15；（考虑空气折射和太阳视半径）
以上三式中，括号中的值大于1时，视为等于1；小于-1时，视为等于-1。
TD=24-TR；
TE=24-TM；
LD=24-2*TR；
LS=(TR-TM)*2；
LN=TM*2；
A=arccos（-sinδ/cosφ）；（不考虑空气折射和太阳视半径）
A=arccos（tghr*tgφ-sinδ/（coshr*cosφ））；（考虑空气折射和太阳视半径）
日出取负，日落出正。