太阳直射点经纬度计算方法

在一个三维系统中，要根据时间确定太阳位置，并将此位置设置为具体的光源

以产生真实光照及阴影效果

太阳的直射点经纬度计算方法如下：

根据基础的地理知识，

纬度：太阳高度角为90度的地方的纬度

经度：当地时间为正午12点的地方的经度

经度比较好算，一天24小时，地球360度，

经度每差一度，时间差4分钟

经度几乎是每分钟都会变化的，纬度则是每天变化

先来看看经度的计算方法

EXPORT double CALLBACK CalcLongitude(double currLong,double currTime) { double res = 0.0; if(currLong>=0) { res = (12.0 - currTime)*15.0+currLong; if(res>180.0) { res -= 360.0; } } else { res = currLong - (currTime - 12.0)*15.0; if(res<-180.0) res += 360.0; } return res; } 

纬度的计算方法，参考网上的python代码改编过来的

具体的参数也是不懂，不过实际测试了一下，与实际情况一致

int CalcDay(int year,int month,int day) { int leap=0; if((year%4==0 && year%100 !=0) || year%400 ==0) { leap = 1; } else { leap = 0; } int mm[12]= {0,31,leap+59,leap+90,leap+120,leap+151,leap+181,leap+212,leap+243,leap+273,leap+304,leap+334}; return mm[month-1]+day; } EXPORT double CALLBACK CalcLatitude(int year,int month,int day) { double tick = CalcDay(year,month,day) - 79.6764-0.2422*(year-1985)+int((year-1985)/4); double sita = 2*3.14159265*tick/365.2422; double delta = 0.3723+23.2567*sin(sita) +0.1149*sin(2*sita)-0.1712*sin(3*sita) -0.758*cos(sita)+0.3656*cos(2*sita) +0.0201*cos(3*sita); return delta; } 

参考：http://blog.chinaunix.net/u/21370/showart_390554.html

最后

以上就是精明导师为你收集整理的太阳直射点经纬度计算方法的全部内容，希望文章能够帮你解决太阳直射点经纬度计算方法所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。