Tag: 几何学

太阳的位置，一天的时间，经度和纬度

这个问题在三年前已经被问过了。 有一个答案，但我发现在解决scheme中的一个小故障。 下面的代码在R.我已经将它移植到另一种语言，但是已经在R中直接testing了原始代码，以确保问题不在我的移植中。 sunPosition <- function(year, month, day, hour=12, min=0, sec=0, lat=46.5, long=6.5) { twopi <- 2 * pi deg2rad <- pi / 180 # Get day of the year, eg Feb 1 = 32, Mar 1 = 61 on leap years month.days <- c(0,31,28,31,30,31,30,31,31,30,31,30) day <- day + cumsum(month.days)[month] leapdays <- year %% 4 […]