python 计算日出、中午、日落

最新推荐文章于 2024-04-26 11:58:59 发布
最新推荐文章于 2024-04-26 11:58:59 发布
阅读量1.3k 收藏 5
点赞数
分类专栏： Python基础 文章标签： python 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42102546/article/details/125802868
版权
这篇博客介绍了一段JavaScript代码，用于根据给定的经纬度、时区和日期计算日出和日落时间。它包含了多个函数，如mjd、hrsmin等，用于进行日期转换和坐标计算。通过调用calculate_rise_set_center函数，可以获取到日出、日中和日落的具体时间。
摘要由CSDN通过智能技术生成

原版出处：https://www.jianshu.com/p/05d9b2dbb3aa

import execjs

# pip install PyExecJS

# 经度、维度、时区
# longitude = 115.81567872975157
# dimensionality = 28.900022800650184
# timezone = 8

# 百度查询：拾取坐标系统  复制经纬度坐标替换
longitude = 119.711112
dimensionality = 29.788806
timezone = 8

# 以下是js用于计算日出日落时间
mjd_def = """
function mjd(day, month, year, hour) {
    var a, b;
    if (month <= 2) {
        month = month + 12;
        year = year - 1;
        }
    a = 10000.0 * year + 100.0 * month + day;
    if (a <= 15821004.1) {
        b = -2 * Math.floor((year + 4716)/4) - 1179;
        }
    else {
        b = Math.floor(year/400) - Math.floor(year/100) + Math.floor(year/4);
        }
    a = 365.0 * year - 679004.0;
    return (a + b + Math.floor(30.6001 * (month + 1)) + day + hour/24.0);
}
"""
cal_def = """
function hrsmin(hours) {
    var hrs, h, m, dum;
    hrs = Math.floor(hours * 60 + 0.5)/ 60.0;
    h = Math.floor(hrs);
    m = Math.floor(60 * (hrs - h) + 0.5);
    if(h<10)h="0"+h;
    if(m<10)m="0"+m;
    dum = h +":"+ m;

    if (dum < 1000) dum = "0" + dum;
    if (dum <100) dum = "0" + dum;
    if (dum < 10) dum = "0" + dum;
    return dum;
    }


function ipart(x) {
    var a;
    if (x> 0) {
        a = Math.floor(x);
        }
    else {
        a = Math.ceil(x);
        }
    return a;
    }


function frac(x) {
    var a;
    a = x - Math.floor(x);
    if (a < 0) a += 1;
    return a;
    }


function round(num, dp) {
   return Math.round (num * Math.pow(10, dp)) / Math.pow(10, dp);
 }


function range(x) {
    var a, b;
    b = x / 360;
    a = 360 * (b - ipart(b));
    if (a  < 0 ) {
        a = a + 360
    }
    return a
}


function mjd(day, month, year, hour) {
    var a, b;
    if (month <= 2) {
        month = month + 12;
        year = year - 1;
        }
    a = 10000.0 * year + 100.0 * month + day;
    if (a <= 15821004.1) {
        b = -2 * Math.floor((year + 4716)/4) - 1179;
        }
    else {
        b = Math.floor(year/400) - Math.floor(year/100) + Math.floor(year/4);
        }
    a = 365.0 * year - 679004.0;
    return (a + b + Math.floor(30.6001 * (month + 1)) + day + hour/24.0);
}   



function quad(ym, yz, yp) {
    var nz, a, b, c, dis, dx, xe, ye, z1, z2, nz;
    var quadout = new Array;

    nz = 0;
    a = 0.5 * (ym + yp) - yz;
    b = 0.5 * (yp - ym);
    c = yz;
    xe = -b / (2 * a);
    ye = (a * xe + b) * xe + c;
    dis = b * b - 4.0 * a * c;
    if (dis > 0)    {
        dx = 0.5 * Math.sqrt(dis) / Math.abs(a);
        z1 = xe - dx;
        z2 = xe + dx;
        if (Math.abs(z1) <= 1.0) nz += 1;
        if (Math.abs(z2) <= 1.0) nz += 1;
        if (z1 < -1.0) z1 = z2;
        }
    quadout[0] = nz;
    quadout[1] = z1;
    quadout[2] = z2;
    quadout[3] = xe;
    quadout[4] = ye;
    return quadout;
    }


function lmst(mjday, glong) {
    var lst, t, d;
    d = mjday - 51544.5
    t = d / 36525.0;
    lst = range(280.46061837 + 360.98564736629 * d + 0.000387933 *t*t - t*t*t / 38710000);
    return (lst/15.0 + glong/15);
    }


function minisun(t) {
    var p2 = 6.283185307, coseps = 0.91748, sineps = 0.39778;
    var L, M, DL, SL, X, Y, Z, RHO, ra, dec;
    var suneq = new Array;

    M = p2 * frac(0.993133 + 99.997361 * t);
    DL = 6893.0 * Math.sin(M) + 72.0 * Math.sin(2 * M);
    L = p2 * frac(0.7859453 + M / p2 + (6191.2 * t + DL)/1296000);
    SL = Math.sin(L);
    X = Math.cos(L);
    Y = coseps * SL;
    Z = sineps * SL;
    RHO = Math.sqrt(1 - Z * Z);
    dec = (360.0 / p2) * Math.atan(Z / RHO);
    ra = (48.0 / p2) * Math.atan(Y / (X + RHO));
    if (ra <0 ) ra += 24;
    suneq[1] = dec;
    suneq[2] = ra;
    return suneq;
    }

function sin_alt(iobj, mjd0, hour, glong, cglat, sglat) {
    var mjday, t, ra, dec, tau, salt, rads = 0.0174532925;
    var objpos = new Array;
    mjday = mjd0 + hour/24.0;
    t = (mjday - 51544.5) / 36525.0;
    if (iobj == 1) {
        objpos = minimoon(t);
                }
    else {
        objpos = minisun(t);
        }
    ra = objpos[2];
    dec = objpos[1];
    tau = 15.0 * (lmst(mjday, glong) - ra);
    salt = sglat * Math.sin(rads*dec) + cglat * Math.cos(rads*dec) * Math.cos(rads*tau);
    return salt;
    }

function getzttime(mjday, tz, glong) {
    var sglong, sglat, date, ym, yz, utrise, utset, j;
    var yp, nz, hour, z1, z2, iobj, rads = 0.0174532925;
    var quadout = new Array;

    sinho = Math.sin(rads * -0.833);        
    date = mjday - tz/24;
    hour = 1.0;
    ym = sin_alt(2, date, hour - 1.0, glong, 1, 0) - sinho;

    while(hour < 25) {
        yz = sin_alt(2, date, hour, glong, 1, 0) - sinho;
        yp = sin_alt(2, date, hour + 1.0, glong, 1, 0) - sinho;
        quadout = quad(ym, yz, yp);
        nz = quadout[0];
        z1 = quadout[1];
        z2 = quadout[2];
        xe = quadout[3];
        ye = quadout[4];

        if (nz == 1) {
            if (ym < 0.0) 
                utrise = hour + z1;         
            else 
                utset = hour + z1;

        } 
        if (nz == 2) {
            if (ye < 0.0) {
                utrise = hour + z2;
                utset = hour + z1;
            }
            else {
                utrise = hour + z1;
                utset = hour + z2;
            }
        } 
        ym = yp;
        hour += 2.0;
    } 
    var zt=(utrise+utset)/2;
    if(zt<utrise)
        zt=(zt+12)%24;
    return zt;
}

function Cal(mjday, tz, glong, glat) {

    var sglong, sglat, date, ym, yz, above, utrise, utset, j;
    var yp, nz, rise, sett, hour, z1, z2, iobj, rads = 0.0174532925;
    var quadout = new Array;
    var   always_up = "日不落";
    var always_down = "日不出";
    var outstring = "";
    var resobj = {};

    sinho = Math.sin(rads * -0.833);        
    sglat = Math.sin(rads * glat);
    cglat = Math.cos(rads * glat);
    date = mjday - tz/24;

    rise = false;
    sett = false;
    above = false;
    hour = 1.0;
    ym = sin_alt(2, date, hour - 1.0, glong, cglat, sglat) - sinho;
    if (ym > 0.0) above = true;
    while(hour < 25 && (sett == false || rise == false)) {
        yz = sin_alt(2, date, hour, glong, cglat, sglat) - sinho;
        yp = sin_alt(2, date, hour + 1.0, glong, cglat, sglat) - sinho;
        quadout = quad(ym, yz, yp);
        nz = quadout[0];
        z1 = quadout[1];
        z2 = quadout[2];
        xe = quadout[3];
        ye = quadout[4];

        if (nz == 1) {
            if (ym < 0.0) {
                utrise = hour + z1;
                rise = true;
            }
            else {
                utset = hour + z1;
                sett = true;
            }
        } 

        if (nz == 2) {
            if (ye < 0.0) {
                utrise = hour + z2;
                utset = hour + z1;
            }
            else {
                utrise = hour + z1;
                utset = hour + z2;
            }
        } 

        ym = yp;
        hour += 2.0;

    } 

    if (rise == true || sett == true ) {
        if (rise == true) {
            resobj["rise"] = hrsmin(utrise);
            outstring += "&nbsp;&nbsp;日出时刻:" + hrsmin(utrise)+"<BR>";
        }
        else{
            outstring += "&nbsp;&nbsp;日不出或日不落"+"<BR>";  
            resobj["pole"] = "日不出或日不落";
        }
        var zt=getzttime(mjday, tz, glong);
        resobj["center"] = hrsmin(zt);
        outstring+= "&nbsp;&nbsp;日中时刻:" + hrsmin(zt)+"<BR>";
        if (sett == true){ 
            outstring += "&nbsp;&nbsp;日没时刻:" + hrsmin(utset);
            resobj["set"] = hrsmin(utset);
        }else {
            outstring += "&nbsp;&nbsp;日不出或日不落";
            resobj["pole"] = "日不出或日不落";
        }
    }
    else {
        if (above == true){
            outstring += always_up;
            var zt=getzttime(mjday, tz, glong);
            resobj["center"] = hrsmin(zt);
            resobj["pole"] = "极昼";
            outstring+="<BR>"+"&nbsp;&nbsp;日中时刻:"+hrsmin(zt);
        } else {
            outstring += always_down;
            resobj["pole"] = "极夜";
        }
    }       
    return resobj;
}

"""


def calculate_rise_set_center(datatime, timezone, longitude_latitude):
    """
    parameter:
    ----------
    datatime: tuple(day, month, year, hour)
    timezone: int 8
    longitude_latitude: tuple(longitude, latitude)
    ====================
    return:
    dict(rise, center, set)

    example:
    calculate_rise_set_center((12, 10, 2018, 0.0),8,(116.3833, 39.9))
    return: {'rise': '06:21', 'center': '12:01', 'set': '17:41'}
    """
    # 计算mjd
    ctx = execjs.compile(mjd_def)
    mjd = ctx.call("mjd", datatime[0], datatime[1], datatime[2], datatime[3])  # 传入：天、月、年、小时

    # 计算日出日落
    ctx = execjs.compile(cal_def)
    result = ctx.call("Cal", mjd, timezone, longitude_latitude[0], longitude_latitude[1])  # 传入：时区、经度、维度
    return result


import datetime

day = datetime.datetime.now().day
month = datetime.datetime.now().month
year = datetime.datetime.now().year
# 日 月 年 0.0 时区：8 经纬度
print(day, month, year)
# a = calculate_rise_set_center((day, month, year, 0.0), timezone, (longitude, dimensionality))
a = calculate_rise_set_center((31, 8, 2022, 0.0), timezone, (longitude, dimensionality))
print(a["set"])
默执_
关注
专栏目录
根据经纬度计算日出日落时间
10-10
根据经纬度计算日出日落时间。只需要经度和纬度。
Python 超级简单精准计算地点日出日落时间
12-20
原文来自Python实用宝典：Python 超级简单精准计算地点日出日落时间 本来想浪漫一把去看个日出，过去后看见天上挂着个大太阳的尴尬情景你经历过吗？ 本来想去海边看个日落，结果车到了，太阳也落了的悲伤情景你经历过吗？ 我们今天要用Python解决的，就是这种尴尬情况：想到某个地点看日出日落，却不知道那个地方准确的日出日落时间。 需要使用的包是 Astral ，这个包能通过你提供的经纬度并利用相应的数学知识计算日出日落时间。下面就来以西涌三号沙滩为例，告诉大家怎么样超级简单地计算某个地点日出日落的时间。 一、准备 进行这个实验，你需要安装 python3. 如果你还没有安装，可以看这个教程
python日出日落时间实现和详解
Boy_Martin的博客
04-19 1115
python实现日出日落得.exe程序教程
利用Python包Astral计算日出日落时间
最新发布
南山老沙
04-26 368
使用astral计算日落日出时间
sunrisesunset:Python脚本确定日出和日落
05-10
日出日落 确定日出和日落的Python脚本。
Python计算日出日落时间/Astral工具包
weixin_46250368的博客
08-19 1867
Python 日出日落时间计算 astral 2.2
根据经纬度和日期计算日出日落时间
06-28
在IT领域，根据经纬度和日期计算日出日落时间是一项常见的地理信息处理任务，它涉及到天文学、地球物理学和编程技术。这个程序的核心是利用天文学公式来预测太阳在地平线上的位置，从而确定一天中的日出和日落时刻。...
suncycle.rar_日出 日落_日出、日落时间 matlab_日出日落_日出时间_日出时间计算
07-13
在IT领域，特别是数据分析和科学计算中，经常需要处理与天文事件相关的问题，例如计算日出和日落的时间。这个名为“suncycle.rar”的压缩包文件提供了相关的资源，包括一个MATLAB脚本（suncycle.m）和一个文档...
怎么使用Python计算日出日落
05-27
Python中有一些库可以帮助我们计算日出日落时间，比如`astral`和`ephem`。我这里给你一个使用`astral`库计算日出日落时间的示例代码： ```python from astral.sun import sun from astral import LocationInfo ...
批量抓取日落时刻数据python
11-14
该段代码可以实现批量抓取日落时刻数据。方便提取网站的数据。
astral:Python计算太阳和月亮的位置
05-03
星体 这是“ astral” Python模块，该模块计算 太阳各个位置的时间：黎明，日出，太阳正午，日落，黄昏，太阳高度，太阳方位角和rahukaalam。 月相。 有关文档，请参见
pysolar:Pysolar是Python库的集合，用于模拟太阳对地球上任何点的照射。 它包括用于非常精确的星历表计算的代码
05-14
烈火 Pysolar是Python库的集合，用于模拟太阳对地球上任何点的照射。 它包含用于极其精确的星历表计算的代码，等等。 注意：目前，Pysolar的最新提交不适用于Python 2.x 版本0.6适用于2.x： ： 但0.7和更高版本有很多更改。 它们已经针对Python 3.4进行了验证，但是3.2版或更早的版本缺少更改所需的功能。 * 此外，API略有变化： Pysolar现在希望您提供时区感知的datetime ，而不是UTC中的原始日期时间。 如果您的结果看起来很疯狂，这可能就是原因。 现在，函数名称为lowercase_separated_by_underscores ，符合 。 安装 假设您已安装Python 3.4或更高版本，则可以使用pip安装Pysolar： sudo pip install pysolar 文档现在出现在。 会费 所有贡献都通过Gi
日出日落时间计算.html （js）
05-20
日出日落时间计算.html （js） js脚本
日出日落时间计算
11-09
计算某地某天的日出日落时间、太阳正午时间，以及正午太阳高度角
sun.js:Calculate sunrise and sunset times in Javascript.(JS计算日出日落时间)
05-19
sun.js 新增内容：由Hollow Man修改，因此指定的时间始终在日出和日落时间之间。 兰州大学（2020）版权所有:copyright:2020hollowman6。 用Javascript计算日出和日落时间。 该库扩展了Javascript Date对象，添加了计算日出和日落时间的方法。 它还扩展了Math对象，以添加一些有用的支持功能。 如果您不喜欢扩展核心JS类型，则可能应该重写它！ 松散且间接地基于Kevin Boone对美国海军天文台算法的SunTimes Java实现。 用法很简单： //Sunset tonight at the Triggertrap office var sunset = new Date ( ) . sunset ( 51.4541 , - 2.5920 ) ; //Sunrise at Stonehenge on midsummer's day 2
python通过astral,matpotlin库实现计算日出日落并绘制折线图
sumido的博客
06-16 1664
python通过astral,matpotlin库实现计算日出日落并绘制折线图 from astral import LocationInfo #引用库 import datetime from astral.sun import sun import numpy as np import matplotlib.pyplot as plt #import r
Python实现天文计算
09-13 701
迷途小书童读完需要2分钟速读仅需 1 分钟1 简介ephem 模块为 Python 提供了精确的天文计算能力，可以预测星球、卫星的轨道信息，计算日出日落、经星时间等数据，它的算法准确可靠。最初由 Brandon Craig Rhodes 在 20 世纪 90 年代开发，经过多年扩展补充，应用非常广泛。2 安装首先，确保您的系统上安装了 Python。您可以从 Python 的官方网站 ( ...
计算日落日出时间算法与代码
热门推荐
AndyLizh的专栏
10-20 2万+
算法： 怎样用经纬度计算日出日落的时间 下面是一种随经纬度变化的日出日落时间计算方法,我成功运用在一智能路灯控制器中,希望对需要的朋友有帮助。    已知：日出日落时太阳的位置h＝-0.833°,要计算地的地理位置,经度Long,纬度G1at,时区zone,UTo为上次计算的日出日落时间,第一次计算时UTo＝180°。  (1)先计算出从格林威治时间公元2000年1
JS查询日出日落时间和太阳月亮位置
weixin_46147602的博客
03-29 1642
基于suncalc实现计算太阳日出日落时间
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

默执_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值