java获取日出日落时间
import java.math.BigDecimal;
import java.text.ParseException;
import java.text.SimpleDateFormat;
import java.util.Date; /**
* (1)先计算出从格林威治时间公元2000年1月1日到计算日天数days;
* (2)计算从格林威治时间公元2000年1月1日到计算日的世纪数t, 则t＝(days+UTo／360)／36525;
* (3)计算太阳的平黄径 ： L=280.460+36000.770×t;
* (4)计算太阳的平近点角 ：G＝357.528+35999.050×t
* (5)计算太阳的黄道经度 ：λ＝L+1.915×sinG+0.020xsin(2G);
* (6)计算地球的倾角 ε＝23.4393-0.0130×t;
* (7)计算太阳的偏差 δ＝arcsin(sinε×sinλ);
* (8)计算格林威治时间的太阳时间角GHA： GHA=UTo-180-1.915×sinG-0.020×sin(2G) +2.466×sin(2λ)-0.053×sin(4λ)
* (9)计算修正值e： e=arcos{[ sinh-sin(Glat)sin(δ)]/cos(Glat)cos(δ)}
* (10)计算新的日出日落时间 ：UT＝UTo-(GHA+Long±e); 其中“+”表示计算日出时间,“-”表示计算日落时间;
* (11)比较UTo和UT之差的绝对值,如果大于0.1°即0.007小时,把UT作为新的日出日落时间值,重新从第(2)步开始进行迭代计算,如果UTo和UT之差的绝对值小于0.007小时,则UT即为所求的格林威治日出日落时间;
* (12)上面的计算以度为单位,即180°=12小时,因此需要转化为以小时表示的时间,再加上所在的时区数Zone,即要计算地的日出日落时间为 ：T=UT/15+Zone
* 上面的计算日出日落时间方法适用于小于北纬60°和南纬60°之间的区域,如果计算位置为西半球时,经度Long为负数。
*/
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、判断是否为闰年：若为闰年，返回1；若不是闰年,返回0
public static boolean leap_year(int year) {
if (((year % 400 == 0) || (year % 100 != 0) && (year % 4 == 0))) return true;
else return false;
} // 1、求从格林威治时间公元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 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); 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); return (int) (sunset / 15 + 8) + ":" + (int) (60 * (sunset / 15 + 8 - (int) (sunset / 15 + 8)));
}
return null;
} /**
* 将当前时间转换为16进制
*
* @return
*/
public static String getTimeTo16(String time) {
Date date = null;
SimpleDateFormat formatter = new SimpleDateFormat("yyyy-MM-dd HH:mm:ss");
try {
date = formatter.parse(time);
} catch (ParseException e) {
e.printStackTrace();
}
//Date格式
Long t = date.getTime() / 1000;
String hexString = Long.toHexString(t);
System.out.println("十六进制：" + hexString); return hexString;
} //将指定时间转换成 date 格式
public static Date getTime(String time) {
Date date = null;
SimpleDateFormat formatter = new SimpleDateFormat("yyyy-MM-dd"); //日期
try {
date = formatter.parse(time);
} catch (ParseException e) {
e.printStackTrace();
}
return date;
} public static void main(String[] args) {
System.out.println(SunRiseSet.getSunrise(new BigDecimal(108.94359),new BigDecimal(34.352276),new Date()));
System.out.println(SunRiseSet.getSunset(new BigDecimal(108.94359),new BigDecimal(34.352276),new Date()));
}
}
java获取日出日落时间的更多相关文章
- java学习第13天( java获取当前时间,有关大数据的运算及精确数字运算,Date类)
一 java获取当前时间 学习一个函数,得到当前时间的准确值 System.currectTimeMillis(). 可以得到以毫秒为单位的当前时间.它主要用于计算程序运行时间,long start= ...
- java获取前一天时间SimpleDateFormat，java判断某个时间段
java获取前一天时间SimpleDateFormat SimpleDateFormat predf = new SimpleDateFormat("yyyy-MM-dd"); D ...
- Java获取系统时间少了八个小时
Java获取系统时间少了八个小时 今天忽然遇到需要获取当前时间的问题,我向来谨慎,先测试获取到的系统时间是否正确,结果竟然发现少了八个小时,晕死了,记得之前在页面用javascript获取过当前时间, ...
- java 获取当前时间，前一天时间
java获取当前时间,并按一定格式输出 1.用Calendar获取Date Calendar calendar=Calendar.getInstance(); SimpleDateFormat for ...
- android java获取当前时间的总结
import java.text.SimpleDateFormat; SimpleDateFormat formatter = new SimpleDateFormat (&q ...
- java获取当前日期时间代码总结
1.获取当前时间,和某个时间进行比较.此时主要拿long型的时间值. 方法如下: 要使用 java.util.Date .获取当前时间的代码如下 代码如下 复制代码 Date date = new ...
- Java获取指定时间的毫秒值的方法
有以下两种方法获取指定时间的毫秒值: 1.Calendar类 先由getInstance获取Calendar对象,然后用clear方法将时间重置为(1970.1.1 00:00:00),接下来用set ...
- Java 获取当前时间及实现时间倒计时功能
引言 在一些项目中或是一些特殊的业务场景中,需要用到显示系统的当前时间,以及一些固定的时间倒计时,时间到后做一些什么事情的业务 .接下来咱们就具体看看代码是怎么实现的: <%@ page lan ...
- JAVA获取当前时间加一天
01.获取当前时间 SimpleDateFormat df = new SimpleDateFormat("yyyy-MM-dd HH:mm:ss"); return df.for ...
随机推荐
- 一次性搞清Java中的类加载问题
摘要:很多时候提到类加载,大家总是没法马上回忆起顺序,这篇文章会用一个例子为你把类加载的诸多问题一次性澄清. 本文分享自华为云社区<用1个例子加5个问题,一次性搞清java中的类加载问题[奔跑吧 ...
- Vue指令实现原理
前言 自定义指令是vue中使用频率仅次于组件,其包含bind.inserted.update.componentUpdated.unbind五个生命周期钩子.本文将对vue指令的工作原理进行相应介绍, ...
- 二、JavaSE语言基础之常量与变量
1.常量 所谓常量值的是数据处理过程中值不能更改的数据. 2.变量 所谓变量值的是运算过程中值可以改变的数据,类似于代数中的未知数. 在Java语言中,使用变量时必须遵循先定义,而后赋值, ...
- 17、lnmp_php编译安装
17.1.FastCGI介绍: 1.什么是CGI: CGI的全称为"通用网关接口",为http服务器与其他机器上的程序服务通信交流的一种工具,CGI程序 必须运行在网络服务器上:传 ...
- LCA总结
作为一名合格的 OIer ,一定要有自我总结的意识,一定要通过写博客的方式来验证自己的掌握程度 ----沃.茨基硕德 目录 作为一名合格的 OIer ,一定要有自我总结的意识,一定要通过写博客的方式来 ...
- 互联网时代CRM软件帮助企业销售升级
随着信息技术的发展,互联网+的浪潮逐渐改变了我们的生活.对于企业来说,他们的管理模式和服务模式也需要作出改变,企业不再满足只进行内部业务的优化和管理,传统CRM开始不再符合企业的需求.由此可见,在网络 ...
- CentOS-自定义SFTP用户及目录
ftp功能说明:通过SSH启动CentOS的sftp功能 创建用户组及用户(sftp可变) $ groupadd sftp $ useradd -g sftp -s /sbin/nologin -d ...
- php 经典的算法题-偷苹果
有5个人偷了一堆苹果,准备在第二天分赃.晚上,有一人遛出来,把所有菜果分成5份,但是多了一个,顺手把这个扔给树上的猴了,自己先拿1/5藏了.没想到其他四人也都是这么想的,都如第一个人一样分成5份把多的 ...
- ESP32S2获取sht30温湿度
static const char *TAG = "i2c-temp"; static unsigned char sht30_buf[6]={0}; static float g ...
- POJ 尺取法
poj3061 Subsequence 题目链接: http://poj.org/problem?id=3061 挑战P146.题意:给定长度为n的数列整数a0,a1,...,a(n-1)以及整数S, ...