这个问题在三年前已经被问过了。 有一个答案,但我发现在解决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 […]
我想知道什么时候太阳在地平线以下X度以上。 例如,我想find太阳在地平线以下19.75度的时间。 我认为这与函数date_sunrise / date_sunset但我不确定。 提前致谢!
我如何画这样的星图? 我有: 具有坐标(赤经和赤纬) 观察者坐标(纬度,经度) 观察的时间 我需要一个考虑这些参数的公式。