给定时间,纬度和经度的太阳位置
三年多以前才问过这个问题。给出了一个答案,但是我在解决方案中发现了一个小故障。 下面的代码在R中。我已将其移植到另一种语言,但是已经在R中直接测试了原始代码,以确保问题不在于我的移植。 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, e.g. 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 …