Question

我正在尝试实现两个功能，1）将日期转换为朱利安日数和2）将朱利安日数转换回日期（年，月和日）。我的代码返回的结果是一天。我不熟悉算法是如何工作的但是我可以看到，如果我使用Julian数的Math.Ceiling它可以工作，但我不确定它是修复代码的最佳方法。任何帮助将不胜感激。谢谢。

public class Program
{
    /// <summary>
    /// Converts year, month and day to a Julian number
    /// </summary>
    /// <param name="year"></param>
    /// <param name="month"></param>
    /// <param name="day"></param>
    /// <returns></returns>
    public static decimal JulianNumber(int year, int month, int day)
    {
        decimal a, b, c, e, f;

        if (month == 1 || month == 2)
        {
            year -= 1;
            month += 12;
        }
        a = Math.Truncate((decimal) year / 100);
        b = Math.Truncate(a / 4);
        c = 2 - a + b;
        e = Math.Truncate((365.25m * (year + 4716)));
        f = Math.Truncate((30.6001m * (month + 1)));

        return (c + day + e + f - 1524.5m);
        // return Math.Ceiling(c + day + e + f - 1524.5m);
    }

    /// <summary>
    /// Converts Julian number to year, month and day
    /// </summary>
    /// <param name="julianNumber"></param>
    /// <returns></returns>
    public static (int year, int month, int day) GregorianDate(decimal julianNumber)
    {
        int l, n, i, j, k;

        l = (int)julianNumber + 68569;
        n = 4 * l / 146097;
        l = l - (146097 * n + 3) / 4;
        i = 4000 * (l + 1) / 1461001;
        l = l - 1461 * i / 4 + 31;
        j = 80 * l / 2447;
        k = l - 2447 * j / 80;
        l = j / 11;
        j = j + 2 - 12 * l;
        i = 100 * (n - 49) + i + l;

        return (i, j, k);
    }

    public static void Main(string[] args)
    {
        var (year1, month1, day1) = (2010, 1, 2);

        var (year2, month2, day2) = GregorianDate(JulianNumber(year1, month1, day1));

        Console.WriteLine(year1 == year2);     // True
        Console.WriteLine(month1 == month2);   // True
        Console.WriteLine(day1 == day2);       // False!
    }
}

Answer 1

问题在于朱利安日的定义。取自您链接的页面：

继Herschel的首席天文学家采用该系统后，中午GMT -4712-01-01 JC（公元前4713年1月1日）为零点

所以2018-05-28 00.00是2458266.5,201-05-28,12.00是2458267.如果你看JulianNumber，Math.Ceiling实际上回报2458266.5。现在，the page you linked（您从中使用第二个方法GregorianDate仅使用整数，因此它适用于中午（12.00）的日期。所以通过四舍五入（天花板）JulianNumber的结果，您将日期移至12.00小时，并使其与GregorianDate兼容，并且＃34;

可能的解决方案：用于JulianNumber same page中存在的算法，并且在任何地方仅使用int（以显示您忽略小时，分钟，秒），或搜索GregorianDate的其他算法。

将日期（年，月，日）转换为Julian Day Number并重新更新

1 个答案: