Java 根据经纬度和日期计算出日出和日落时间

最新推荐文章于 2024-06-11 13:52:03 发布
最新推荐文章于 2024-06-11 13:52:03 发布
阅读量5.1k 收藏 12
点赞数 4
分类专栏： 随心笔记 文章标签： java html js jquery 
package com.rtmc.cloud.commons.utils.calculationSunriseAndSunsetTime;

import java.math.BigDecimal;
import java.text.SimpleDateFormat;
import java.util.Date;

public class SunRiseSet {

    private static int[] days_of_month_1 = {31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};

    private static int[] days_of_month_2 = {31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};

    private final static double h = -0.833;//日出日落时太阳的位置

    private final static double UTo = 180.0;//上次计算的日落日出时间，初始迭代值180.0

    //输入日期

//输入经纬度

//判断是否为闰年：若为闰年，返回1；若不是闰年,返回0

    public static boolean leap_year(int year) {

        if (((year % 400 == 0) || (year % 100 != 0) && (year % 4 == 0))) return true;

        else return false;

    }

//求从格林威治时间公元2000年1月1日到计算日天数days

    public static int days(int year, int month, int date) {

        int i, a = 0;

        for (i = 2000; i < year; i++) {

            if (leap_year(i)) a = a + 366;

            else a = a + 365;

        }

        if (leap_year(year)) {

            for (i = 0; i < month - 1; i++) {

                a = a + days_of_month_2[i];

            }

        } else {

            for (i = 0; i < month - 1; i++) {

                a = a + days_of_month_1[i];

            }

        }

        a = a + date;

        return a;

    }

//求格林威治时间公元2000年1月1日到计算日的世纪数t

    public static double t_century(int days, double UTo) {

        return ((double) days + UTo / 360) / 36525;

    }

//求太阳的平黄径

    public static double L_sun(double t_century) {

        return (280.460 + 36000.770 * t_century);

    }

//求太阳的平近点角

    public static double G_sun(double t_century) {

        return (357.528 + 35999.050 * t_century);

    }

//求黄道经度

    public static double ecliptic_longitude(double L_sun, double G_sun) {

        return (L_sun + 1.915 * Math.sin(G_sun * Math.PI / 180) + 0.02 * Math.sin(2 * G_sun * Math.PI / 180));

    }

//求地球倾角

    public static double earth_tilt(double t_century) {

        return (23.4393 - 0.0130 * t_century);

    }

//求太阳偏差

    public static double sun_deviation(double earth_tilt, double ecliptic_longitude) {

        return (180 / Math.PI * Math.asin(Math.sin(Math.PI / 180 * earth_tilt) * Math.sin(Math.PI / 180 * ecliptic_longitude)));

    }

//求格林威治时间的太阳时间角GHA

    public static double GHA(double UTo, double G_sun, double ecliptic_longitude) {

        return (UTo - 180 - 1.915 * Math.sin(G_sun * Math.PI / 180) - 0.02 * Math.sin(2 * G_sun * Math.PI / 180) + 2.466 * Math.sin(2 * ecliptic_longitude * Math.PI / 180) - 0.053 * Math.sin(4 * ecliptic_longitude * Math.PI / 180));

    }

//求修正值e

    public static double e(double h, double glat, double sun_deviation) {

        return 180 / Math.PI * Math.acos((Math.sin(h * Math.PI / 180) - Math.sin(glat * Math.PI / 180) * Math.sin(sun_deviation * Math.PI / 180)) / (Math.cos(glat * Math.PI / 180) * Math.cos(sun_deviation * Math.PI / 180)));

    }

//求日出时间

    public static double UT_rise(double UTo, double GHA, double glong, double e) {

        return (UTo - (GHA + glong + e));

    }

//求日落时间

    public static double UT_set(double UTo, double GHA, double glong, double e) {

        return (UTo - (GHA + glong - e));

    }

//判断并返回结果（日出）

    public static double result_rise(double UT, double UTo, double glong, double glat, int year, int month, int date) {

        double d;

        if (UT >= UTo) d = UT - UTo;

        else d = UTo - UT;

        if (d >= 0.1) {

            UTo = UT;

            UT = UT_rise(UTo,

                    GHA(UTo, G_sun(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))),

                    glong,

                    e(h, glat, sun_deviation(earth_tilt(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo))))));

            result_rise(UT, UTo, glong, glat, year, month, date);


        }

        return UT;

    }

//判断并返回结果（日落）

    public static double result_set(double UT, double UTo, double glong, double glat, int year, int month, int date) {

        double d;

        if (UT >= UTo) d = UT - UTo;

        else d = UTo - UT;

        if (d >= 0.1) {

            UTo = UT;

            UT = UT_set(UTo,

                    GHA(UTo, G_sun(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))),

                    glong,

                    e(h, glat, sun_deviation(earth_tilt(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo))))));

            result_set(UT, UTo, glong, glat, year, month, date);

        }

        return UT;

    }

//求时区

    public static int Zone(double glong) {

        if (glong >= 0) return (int) ((int) (glong / 15.0) + 1);

        else return (int) ((int) (glong / 15.0) - 1);

    }

//打印结果

// public static void output(double rise, double set, double glong){

//     if((int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))<10)

//         System.out.println("The time at which the sunrise is: "+(int)(rise/15+Zone(glong))+":"+(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))+" .\n");

//     else System.out.println("The time at which the sunrise is: "+(int)(rise/15+Zone(glong))+":"+(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))+" .\n");

//

//     if((int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))<10)

//         System.out.println("The time at which the sunset is: "+(int)(set/15+Zone(glong))+": "+(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))+" .\n");

//     else System.out.println("The time at which the sunset is: "+(int)(set/15+Zone(glong))+":"+(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))+" .\n");

// }

    public static String getSunrise(BigDecimal longitude, BigDecimal latitude, Date sunTime) {
        if (sunTime != null && longitude != null && latitude != null) {
            double sunrise, glong, glat;
            int year, month, date;
            SimpleDateFormat sdf = new SimpleDateFormat("yyyy-MM-dd");
            String dateTime = sdf.format(sunTime);
            String[] rq = dateTime.split("-");
            String y = rq[0];
            String m = rq[1];
            String d = rq[2];
            year = Integer.parseInt(y);
            if (m != null && m != "" && m.indexOf("0") == -1) {
                m = m.replaceAll("0", "");
            }
            month = Integer.parseInt(m);

            date = Integer.parseInt(d);

            glong = longitude.doubleValue();

            glat = latitude.doubleValue();

            sunrise = result_rise(UT_rise(UTo,

                    GHA(UTo, G_sun(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))),

                    glong,

                    e(h, glat, sun_deviation(earth_tilt(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))))), UTo, glong, glat, year, month, date);

//System.out.println("Sunrise is: "+(int)(sunrise/15+Zone(glong))+":"+(int)(60*(sunrise/15+Zone(glong)-(int)(sunrise/15+Zone(glong))))+" .\n");

//        Log.d("Sunrise", "Sunrise is: "+(int)(sunrise/15+8)+":"+(int)(60*(sunrise/15+8-(int)(sunrise/15+8)))+" .\n");

            //return "Sunrise is: "+(int)(sunrise/15+Zone(glong))+":"+(int)(60*(sunrise/15+Zone(glong)-(int)(sunrise/15+Zone(glong))))+" .\n";

            return (int) (sunrise / 15 + 8) + ":" + (int) (60 * (sunrise / 15 + 8 - (int) (sunrise / 15 + 8)));
        }
        return null;
    }


    public static String getSunset(BigDecimal longitude, BigDecimal latitude, Date sunTime) {
        if (sunTime != null && latitude != null && longitude != null) {
            double sunset, glong, glat;
            int year, month, date;
            SimpleDateFormat sdf = new SimpleDateFormat("yyyy-MM-dd");
            String dateTime = sdf.format(sunTime);
            String[] rq = dateTime.split("-");
            String y = rq[0];
            String m = rq[1];
            String d = rq[2];
            year = Integer.parseInt(y);
            if (m != null && m != ""  && m.indexOf("0") == -1) {
                m = m.replaceAll("0", "");
            }
            month = Integer.parseInt(m);

            date = Integer.parseInt(d);

            glong = longitude.doubleValue();

            glat = latitude.doubleValue();

            sunset = result_set(UT_set(UTo,

                    GHA(UTo, G_sun(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))),

                    glong,

                    e(h, glat, sun_deviation(earth_tilt(t_century(days(year, month, date), UTo)),

                            ecliptic_longitude(L_sun(t_century(days(year, month, date), UTo)),

                                    G_sun(t_century(days(year, month, date), UTo)))))), UTo, glong, glat, year, month, date);

//System.out.println("The time at which the sunset is: "+(int)(sunset/15+Zone(glong))+":"+(int)(60*(sunset/15+Zone(glong)-(int)(sunset/15+Zone(glong))))+" .\n");

//        Log.d("Sunset", "Sunset is: "+(int)(sunset/15+8)+":"+(int)(60*(sunset/15+8-(int)(sunset/15+8)))+" .\n");

            //return "Sunset is: "+(int)(sunset/15+Zone(glong))+":"+(int)(60*(sunset/15+Zone(glong)-(int)(sunset/15+Zone(glong))))+" .\n";

            return (int) (sunset / 15 + 8) + ":" + (int) (60 * (sunset / 15 + 8 - (int) (sunset / 15 + 8)));
        }
        return null;
    }

    public static void main(String[] args) {
        String str1 = SunRiseSet.getSunrise(new BigDecimal(104.066507),new BigDecimal(30.625124),new Date());
        String str2 = SunRiseSet.getSunset(new BigDecimal(104.066507),new BigDecimal(30.625124),new Date());
        System.out.println("日出时间：" + str1);
        System.out.println("日落时间：" + str2);
    }
}

 

王根基
关注
  • 4
    点赞
  • 踩
  • 12
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
计算日出日落算法实现
从一点一滴做起
03-01 7417
首先感谢帮我算题的小玉玉~ 算了一周的日出日落，倍觉痛苦，写下来防止下次再算。 根据经纬度和UTC计算日出日落的思路、算法及scala代码如下~ 日出和日落的计算思路： 一天的日出和日落可以根据正午时刻平分昼长得到。 正午时刻可以根据所在时区和经度得到。 根据经度可以知道所在时区。 昼长可以根据太阳赤角和所在纬度计算得到。 符号说明： rc：日出 rl：日落 zw：正午...
根据经纬度和日期计算日出日落时间
热门推荐
azoo的专栏
06-28 2万+
计算日出和日落时间会比较复杂,可以通过计算当地的“正午”时间、当地的“白天”时长,用正午时间平分白天时长,得到当地的日出时间和日落时间。 1.计算正午时间 一个地方的正午时间,可以通过当地经度与当地所属时区经度的差计算得到。全世界有24个时区,每个时区相隔1小时,因此每个时区之间相差的经度是 360 / 24 = 15°。已北京为例,北京所属时区为东8区,东8区的中心经度是15 * 8 = 120°(因此东8区的范围是112.5° - 127.5°)。 北京的经度是116.4...
热门api接口：开发者的高效指南
最新发布
海碗吃饭的博客
06-11 382
热门api接口：开发者的高效指南
java版根据日期和经纬度计算日出日落
07-10
main函数需要三个参数，经纬度和距离2000-01-01的天数
Java根据地址计算日出日落时间（百度地图API）
qq_31350875的博客
05-22 2883
获取百度地图API返回值（基于IP大致定位），并转化为JSON Object import cn.hutool.core.lang.Console; import cn.hutool.http.HttpRequest; import cn.hutool.json.JSONObject; import cn.hutool.json.JSONUtil; /* * 基于IP大致定位 * ...
java 日出日落时间计算
04-27
通过坐标及时区计算日出日落工具类，内含 java及C++两个版本
根据经纬度计算日出日落时间
10-10
根据经纬度计算日出日落时间。只需要经度和纬度。
获取日出日落时间的Java工具类
04-27
获取日出日落时间的Java工具类，如果要获取国外的日出日落时间，需要减去对应时差哦
根据本地时间和经纬度来计算真太阳时
01-15
3. **纬度**：虽然纬度对计算真太阳时的影响较小，但在考虑日出日落和太阳高度角时，它会变得重要。 计算步骤如下： 1. **确定格林尼治标准时间（GMT）**：将本地时间转换为UTC，考虑到时区和夏令时的影响。 2. **...
太阳高度角与太阳方位角计算器
06-24
根据经纬度计算日出时间、日落时间 根据经纬度计算日出方位、日落方位 根据经纬度、北京时间计算地方时间 根据地方时间，计算时角 根据经纬度、北京时间，计算太阳高度角、太阳方位角 根据隧道高度、宽度、深度，...
iOS根据日期和经纬度计算日出日落时间
08-31
根据当前的日期和经纬度，计算出日出和日落时间，OC写的，需要的可以下载看一下，已经验证，误差两分钟之内
经纬度计算日出日落时间
04-17
C语言，根据经纬度计算日出、日落、天黑、天亮时间
android 根据坐标日期，计算日出日落时间表
03-31
根据坐标，时间段算出，日出日落作息时间表，并转换为16进制，存入本地。经纬度，时间段可自行配置，实测时间相差一分钟内
Java获取当地的日出日落时间代码分享
09-04
主要介绍了Java获取当地的日出日落时间代码分享,国外猿友写的一个类,需要的朋友可以参考下
如何根据纬度，经度和日期判断太阳的日出和日落的时间
smilestone322的专栏
01-30 6506
如题： 日出时间计算公式： 24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360 我国时区为东8区，时区=8 经度、纬度采用角度制，东经、北纬为正，西经、南纬为负 日期序列数为当天在这一年中的序列，如2月11日就是42 因为计算机一般采用弧度制，上面公式可变化为： 24*(180+时区*15
日出日落时间和年均光照时长计算 java
Empty_mountain的博客
12-13 3785
Java计算日出日落时间，求算年均光照时长。
日出日落知多少？Java一招调用API轻松查询
java专栏
05-10 326
查询日出日落时间涉及到复杂的天文计算，一般不直接通过Java标准库实现，而是采用外部API或者天文计算库。一个常用的API是Sunrise-Sunset API，它允许你通过经纬度和日期查询日出日落时间。然而，如果你需要在Java中自行实现，可以参考一些天文算法，如 NOAA's Solar Calculator Algorithm 或者使用第三方库如 JPARSEC（Java Package for Astronomical Research and Space Calculations）。
基于GPS经纬度和当地时间计算日落日出时间实现
mygod2008ok的专栏
10-17 1720
以NRF52 MCU平台实现，先实现UTC及要本时间等转换函数 #ifndef __RTC_MODULE_H__ #define __RTC_MODULE_H__ #include "time.h" #include "stdbool.h" #include "stdint.h" #include "sdk_config.h" #include "nrf_atomic.h" #pragma anon_unions #ifdef __cplusplus extern "C" { #endif #defi
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值