時刻、緯度、経度を指定した太陽の位置
この質問は、3年ちょっと前に尋ねられました。答えがありましたが、解決策に不具合が見つかりました。 以下のコードは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 …