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

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

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

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

根据基础的地理知识，

纬度：太阳高度角为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