基于GPS经纬度和当地时间计算日落日出时间实现

最新推荐文章于 2024-04-30 13:32:20 发布
最新推荐文章于 2024-04-30 13:32:20 发布
阅读量1.7k 收藏 9
点赞数 3
分类专栏： NRF52832个人学习笔记 嵌入式 文章标签： c语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mygod2008ok/article/details/120812337
版权
该代码实现了一个基于NRF52MCU平台的时间管理模块，包括UTC时间与本地时间的转换，以及日出、日落、中天、天亮和天黑时间的计算。主要函数有utc_to_local_time()和local_time_to_utc()进行时间转换，calc_sunrise_info()用于计算日出日落等相关时间。此外，还提供了logging功能以方便调试。
摘要由CSDN通过智能技术生成

以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

#define INIT_UTC                     (1609502400) 
#define ONE_DAY_UTC                  86400

typedef struct{
    int tm_sec;   /* seconds after the minute, 0 to 60
                     (0 - 60 allows for the occasional leap second) */
    int tm_min;   /* minutes after the hour, 0 to 59 */
    int tm_hour;  /* hours since midnight, 0 to 23 */
    int tm_mday;  /* day of the month, 1 to 31 */
    int tm_mon;   /* months since January, 0 to 11 */
    int tm_year;  /* years since 1900 */
    int tm_wday;  /* days since Sunday, 0 to 6 */
    int tm_yday;  /* days since January 1, 0 to 365 */
    int tm_isdst; /* Daylight Savings Time flag */
    union {       /* ABI-required extra fields, in a variety of types */
        struct {
            int __extra_1, __extra_2;
        };
        struct {
            long __extra_1_long, __extra_2_long;
        };
        struct {
            char *__extra_1_cptr, *__extra_2_cptr;
        };
        struct {
            void *__extra_1_vptr, *__extra_2_vptr;
        };
    };
}rtc_time_t;



extern int8_t time_zone;
extern uint8_t g_sync_utc_flag;

uint32_t utc_inc_one(void);
uint32_t rtcm_get_sys_utc(void);
void rtcm_set_sys_utc(uint32_t tmp_utc);
void rtcm_read_sys_time(rtc_time_t *tmp_rtc);
void rtcm_secs_to_date_time(uint32_t t, rtc_time_t *p_time);
void utc_to_local_time(uint32_t utc,rtc_time_t *pTime);
void local_time_to_utc(uint32_t *utc,rtc_time_t *pTime);


#ifdef __cplusplus
}
#endif

#endif


#include "rtc_module.h"
#include "main.h"
#include "nrfx.h"

/***********************************************************************
*                               macro define
************************************************************************
*/



#define NRF_LOG_MODULE_NAME rtc_module
#if ANT_HRM_LOG_ENABLED
#define NRF_LOG_LEVEL       3
#define NRF_LOG_INFO_COLOR  3
#else // ANT_HRM_LOG_ENABLED
#define NRF_LOG_LEVEL       0
#endif // ANT_HRM_LOG_ENABLED
#include "nrf_log.h"
NRF_LOG_MODULE_REGISTER();

uint8_t g_sync_utc_flag;
int8_t time_zone = 8;        // 时区

static nrf_atomic_u32_t sys_utc = INIT_UTC;

/*************************************************************************   
**	function name:	rtcm_get_sys_utc 
**	description:	get the system utc                   
**	input para:     
**                  
**	return:         return the system utc        		                                       
**************************************************************************/
uint32_t rtcm_get_sys_utc(void)
{
	 uint32_t utc = nrf_atomic_u32_fetch_or(&sys_utc,0);
	 return utc;
}
/**
* @brieaf utc原子操作加1
*/
uint32_t utc_inc_one(void)
{
	return nrf_atomic_u32_add(&sys_utc,1);
}
/*************************************************************************   
**	function name:	rtcm_set_sys_utc 
**	description:	set the system utc                   
**	input para:     
**                  
**	return:                		                                       
**************************************************************************/
void rtcm_set_sys_utc(uint32_t tmp_utc)
{ 
	 
	  uint32_t utc = nrf_atomic_u32_store(&sys_utc,tmp_utc);
		NRF_LOG_INFO("=========Set_time=%d",sys_utc);
}
/*************************************************************************   
**	function name:	rtcm_read_sys_time 
**	description:	read the system rtc                   
**	input para:     tmp_rtc--the rtc you want to read     
**                  
**	return:                		                                       
**************************************************************************/
void rtcm_read_sys_time(rtc_time_t *tmp_rtc)
{
    //rtc_secs_to_date_time(sys_utc+(u32_t)8*3600, tmp_rtc);
	 
	  uint32_t utc = nrf_atomic_u32_fetch_or(&sys_utc,0);  
    utc_to_local_time(utc, tmp_rtc);
}

/**
* @brieaf 本地时间转UTC
*/
void local_time_to_utc(uint32_t *utc,rtc_time_t *pTime)
{
  struct tm p;    
	memcpy(&p,pTime,sizeof(struct tm));
	p.tm_year -= 1900;
	p.tm_mon -= 1;
	*utc =mktime(&p);
}

/**
* @brieaf utc转本地时间
*/
void utc_to_local_time(uint32_t utc,rtc_time_t *pTime)
{
		struct tm *p;  
		p = localtime(&utc); 
		p->tm_year += 1900;
		p->tm_mon++;
    memcpy(pTime,p,sizeof(struct tm));
}

计算日出日落相关函数

#ifndef __BSP_SUNRISE_CALC_H
#define __BSP_SUNRISE_CALC_H
#include <math.h>
#include <stdio.h>
#include "rtc_module.h"

#define M_PI 3.1415926

#ifdef __cplusplus
extern "C" {
#endif

typedef struct TAG_SUNRISE
{
	double sunrise_time;   // 日出
	double sunset_time;    // 日落
	double midDayTime;  //中天时间
	double dawnTime;   //天亮时间
	double darkness_time;   //天黑时间
}SUNRISE_TIME_T;


void calc_sunrise_info(uint32_t utc,double lo,double la,int32_t zone);

 /**
 * @brieaf 获取日出时间
 */
 void get_sunrise_string(char *timestr);

	 
 /**
 * @brieaf 获取日落时间
 */
 void get_sunset_string(char *timestr);

 /**
 * @brieaf 获取中天时间
 */
 void get_midDayTime_string(char *timestr);

 
 /**
 * @brieaf 获取天亮时间
 */
 void get_dawnTime_string(char *timestr);

 /**
 * @brieaf 获取天黑时间
 */
 void get_darkness_time_string(char *timestr);

 

#ifdef __cplusplus
}
#endif

#endif

#include "BSP_sunrise_calc.h"
#include "nrf_log_ctrl.h"
#include "nrf_log_default_backends.h"

#define NRF_LOG_MODULE_NAME sunrise_calc
#define NRF_LOG_LEVEL       3
#define NRF_LOG_INFO_COLOR  4
#define NRF_LOG_DEBUG_COLOR 3
#include "nrf_log.h"
NRF_LOG_MODULE_REGISTER();


static double RAD = 180.0 * 3600 /M_PI;



static SUNRISE_TIME_T s_sun_time;






/*************************
     * 儒略日的计算
     *
     * @param y 年
     *
     * @param M 月
     *
     * @param d 日
     *
     * @param h 小时
     *
     * @param m 分
     *
     * @param s秒
     *
     * @return int
***************************/
static double timeToDouble(int y, int M, double d)
{
//        double A=0;
        double B=0;
        double jd=0;

        //设Y为给定年份，M为月份，D为该月日期（可以带小数）。
        //若M > 2，Y和M不变，若 M =1或2，以Y–1代Y，以M+12代M，换句话说，如果日期在1月或2月，则被看作是在前一年的13月或14月。
        //对格里高利历有 ：A = INT（Y/100）   B = 2 - A + INT(A/4)
        //对儒略历，取 B = 0
        //JD = INT(365.25(Y+4716))+INT(30.6001(M+1))+D+B-1524.5 （7.1）
        B=-13;
        jd=floor(365.25 * (y + 4716))+ floor(30.60001 * (M + 1)) +B+ d- 1524.5;
        return jd;
 }

 


/****************************
     * @param t 儒略世纪数
     *
     * @return 太阳黄经
*****************************/
static double sunHJ(double t)
{
        double j;
        t = t + (32.0 * (t + 1.8) * (t + 1.8) - 20) / 86400.0 / 36525.0;
        // 儒略世纪年数,力学时
        j = 48950621.66 + 6283319653.318 * t + 53 * t * t - 994 + 334166 *cos(4.669257 + 628.307585 * t) + 3489 * cos(4.6261 + 1256.61517 * t) + 2060.6 * cos(2.67823 + 628.307585 * t) * t;
        return (j / 10000000);
}

static double mod(double num1, double num2)
{
        num2 = fabs(num2);
        // 只是取决于Num1的符号
        return num1 >= 0 ?(num1 - (floor(num1 / num2)) * num2 ): ((floor(fabs(num1) / num2)) * num2 - fabs(num1));
}
/********************************
     * 保证角度∈(-π,π)
     *
     * @param ag
     * @return ag
***********************************/
static double degree(double ag)
{
        ag = mod(ag, 2 * M_PI);
        if(ag<=-M_PI){
            ag=ag+2*M_PI;
        }
        else if(ag>M_PI){
            ag=ag-2*M_PI;
        }
        return ag;
}

/***********************************
     *
     * @param date  儒略日平午
     *
     * @param lo    地理经度
     *
     * @param la    地理纬度
     *
     * @param tz    时区
     *
     * @return 太阳升起时间
*************************************/
 static double sunRiseTime(double date, double lo, double la, double tz)
 {
        double t,j,sinJ,cosJ,gst,E,a,D,cosH0,cosH1,H0,H1,H;
        date = date - tz;
        // 太阳黄经以及它的正余弦值
        t = date / 36525;
        j = sunHJ(t);
        // 太阳黄经以及它的正余弦值
        sinJ = sin(j);
        cosJ = cos(j);
        // 其中2*M_PI*(0.7790572732640 + 1.00273781191135448*jd)恒星时（子午圈位置）
        gst = 2 * M_PI * (0.779057273264 + 1.00273781191135 * date) + (0.014506 + 4612.15739966 * t + 1.39667721 * t * t) / RAD;
        E = (84381.406 - 46.836769 * t) / RAD; // 黄赤交角
        a = atan2(sinJ * cos(E), cosJ);// '太阳赤经
        D = asin(sin(E) * sinJ); // 太阳赤纬
        cosH0 = (sin(-50 * 60 / RAD) - sin(la) * sin(D)) / (cos(la) * cos(D)); // 日出的时角计算，地平线下50分
        cosH1 = (sin(-6 * 3600 / RAD) - sin(la) * sin(D)) / (cos(la) * cos(D)); // 天亮的时角计算，地平线下6度，若为航海请改为地平线下12度
        // 严格应当区分极昼极夜，本程序不算
        if (cosH0 >= 1 || cosH0 <= -1){
            return -0.5;// 极昼
        }
        H0 = -acos(cosH0); // 升点时角（日出）若去掉负号 就是降点时角，也可以利用中天和升点计算
        H1 = -acos(cosH1);
        H = gst - lo - a; // 太阳时角
        s_sun_time.midDayTime = date - degree(H) / M_PI / 2 + tz; // 中天时间
        s_sun_time.dawnTime = date - degree(H - H1) / M_PI / 2 + tz;// 天亮时间
        return date - degree(H - H0) / M_PI / 2 + tz; // 日出时间，函数返回值
}

 

static void doubleToStr(double time,char *str)
{
       double t;
       int h,m,s;

        t = time + 0.5;
        t = (t - (int) t) * 24;
        h = (int) t;
        t = (t - h) * 60;
        m = (int) t;
        t = (t - m) * 60;
        s = (int) t;
        sprintf(str,"%02d:%02d:%02d",h,m,s);
}
 

//########################################################################################################################


/**
* @brieaf 计算日落,日出,中天,天亮,天黑时间
* @param[in] utc - 当前时间utc
* @param[in] lo -  经度
* @param[in] la - 纬度
* @param[in] zone - 时区
*/
 void calc_sunrise_info(uint32_t utc,double lo,double la,int32_t zone)
 {
	  rtc_time_t time;
	  
	  utc_to_local_time(utc,&time);
	  double jd = -lo/180*M_PI;
    double wd = la/180*M_PI;  
	  
	  double richu = timeToDouble(time.tm_year,time.tm_mon,time.tm_mday) - 2451544.5;
	 
	  for (uint32_t i = 0; i < 10; i++){
        richu = sunRiseTime(richu, jd, wd, zone/24.0);// 逐步逼近法算10次
    }
	 
		s_sun_time.sunrise_time = richu;
		s_sun_time.sunset_time = s_sun_time.midDayTime + s_sun_time.midDayTime - richu;
		s_sun_time.darkness_time = s_sun_time.midDayTime + s_sun_time.midDayTime-s_sun_time.dawnTime;
		
 }
 
 
 /**
 * @brieaf 获取日出时间
 */
 void get_sunrise_string(char *timestr)
 {
	 doubleToStr(s_sun_time.sunrise_time,timestr);
   NRF_LOG_INFO("------> sunrise %s",timestr);
 }
	 
 /**
 * @brieaf 获取日落时间
 */
 void get_sunset_string(char *timestr)
 {
	 doubleToStr(s_sun_time.sunset_time,timestr);
   NRF_LOG_INFO("------> sunset_time %s",timestr);
 }
 
 /**
 * @brieaf 获取中天时间
 */
 void get_midDayTime_string(char *timestr)
 {
	 doubleToStr(s_sun_time.midDayTime,timestr);
   NRF_LOG_INFO("------> midDayTime %s",timestr);
 }
 
 
 /**
 * @brieaf 获取天亮时间
 */
 void get_dawnTime_string(char *timestr)
 {
	 doubleToStr(s_sun_time.dawnTime,timestr);
   NRF_LOG_INFO("------> dawnTime %s",timestr);
 }
 
 /**
 * @brieaf 获取天黑时间
 */
 void get_darkness_time_string(char *timestr)
 {
	 doubleToStr(s_sun_time.darkness_time,timestr);
   NRF_LOG_INFO("------> darkness_time %s",timestr);
 }

调用calc_sunrise_info函数计算

calc_sunrise_info(1633750052,113.9167f,22.6446,8);
  • 调用获取结果函数
  char buf[40];
  get_sunrise_string(buf);
  get_sunset_string(buf);
  get_midDayTime_string(buf);
  get_dawnTime_string(buf);
  get_darkness_time_string(buf);
  • 运行结果
    <00> info> sunrise_calc: ------> sunrise 06:18:30
    00>
    <00> info> sunrise_calc: ------> sunset_time 18:04:52
    00>
    <00> info> sunrise_calc: ------> midDayTime 12:11:41
    00>
    <00> info> sunrise_calc: ------> dawnTime 05:56:01
    00>
    <00> info> sunrise_calc: ------> darkness_time 18:27:22
风雨依依
关注
专栏目录
【雕爷学编程】Arduino月球基地之使用经纬度和时间库计算日出和日落时间
雕爷学编程
12-30 984
因此，在设计和实施Arduino月球基地场景时，需要根据实际需求和目标权衡这些局限性，并结合其他适当的教学工具和资源，以提供更全面和深入的学习体验。灵活应用：经纬度和时间库的使用不仅限于计算日出和日落时间，还可以根据需要进行其他时间相关计算，如黄昏时间、黎明时间、日中时间等。通过使用经纬度和时间库，可以计算出每天的日出和日落时间，并将其用于太阳能电池板的定时控制，实现最佳的能源捕捉。这些库通常包含了地球的几何模型和天文计算方法，可以根据给定的经纬度和日期，计算出对应地点的日出和日落时间。
c语言计算日出日落时间_利用日期、经纬度求日出日落时间 C语言程序代码
weixin_39637256的博客
12-20 1296
展开全部#definePI3.1415926#include#includeusingnamespacestd;intdays_of_month_1[]={31,28,31,30,31,30,31,31,30,31,30,31};intdays_of_month_2[]={31,29,31,30,31,30,31,31,30,31,30,31};longdoubleh=-0.833...
根据经纬度和日期计算日出日落时间
06-28
可以自行输入经纬度和时间，计算出日出日落时间，亲测误差很小，也就正负2到3分钟
根据经纬度计算日出日落时间
10-10
根据经纬度计算日出日落时间。只需要经度和纬度。
经纬度计算日出日落时间
04-17
C语言，根据经纬度计算日出、日落、天黑、天亮时间
如何根据纬度，经度和日期判断太阳的日出和日落的时间
smilestone322的专栏
01-30 6599
如题： 日出时间计算公式： 24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360 我国时区为东8区，时区=8 经度、纬度采用角度制，东经、北纬为正，西经、南纬为负 日期序列数为当天在这一年中的序列，如2月11日就是42 因为计算机一般采用弧度制，上面公式可变化为： 24*(180+时区*15
经纬度计算日出日落C代码
02-25
【标题】"经纬度计算日出日落C代码"涉及的核心知识点主要集中在地理坐标系统、时间计算以及嵌入式系统编程。以下是对这些知识点的详细解释： 1. **地理坐标系统**：地球上的每一个位置都可以用经度和纬度来表示。...
输入经纬度和当前时间自动求太阳高度角、方位角、天顶角
04-20
在IT行业中，尤其是在数据分析、地理信息系统（GIS）和天文计算领域，了解如何根据地理位置和时间计算太阳的位置是非常重要的。这个程序使用Python编程语言实现，它能够自动计算出太阳的高度角、方位角和天顶角，...
经纬度计算日出日落时间(C语言)
04-12
iOS根据日期和经纬度计算日出日落时间
08-31
根据当前的日期和经纬度，计算出日出和日落时间，OC写的，需要的可以下载看一下，已经验证，误差两分钟之内
android 根据坐标日期，计算日出日落时间表
03-31
根据坐标，时间段算出，日出日落作息时间表，并转换为16进制，存入本地。经纬度，时间段可自行配置，实测时间相差一分钟内
根据当地的经纬度与时间，计算当前太阳高度、角度
09-11
给定当地的经纬度，可以自动计算当前时间下太阳的高度，角度。该算法经过实际验证，与网上提供的软件计算结果相同。
android获得日升日落时间
最新发布
2401_84132544的博客
04-30 223
【代码】android获得日升日落时间。
C语言 根据经纬度如何获取日落日出时间
u011983700的专栏
11-13 1639
代码参考 https://zhidao.baidu.com/question/574747751.html 源代码有几次错误，已经改成可以使用的C版本，计算结果和查询网站的结果会有1分钟左右的误差。 测试网站地址https://richurimo.51240.com/121.126078__jw__29.026046__richurimo/ #include <stdio.h> #include <stdlib.h> #include ...
已知经纬度，基于matlab计算太阳和月亮的升起和落下时间以及黄昏时间
m0_60703264的博客
10-08 596
随着科技的不断发展，我们现在可以通过多种方式获取到地球上任意一个地点的经纬度。这为我们了解当地的天文现象提供了便利。其中，计算太阳和月亮的升起和落下时间以及黄昏时间是非常有用的。在本文中，我们将探讨如何根据已知的经纬度来计算这些时间。首先，让我们来了解一下天文学中的一些基本概念。在地球上观测天文现象时，我们通常使用地平坐标系来描述天体的位置。在这个坐标系中，天体的位置由两个参数确定：方位角和高度角。方位角是从正北方向开始逆时针测量的角度，范围是0到360度。高度角是天体与地平面的夹角，范围是0到90度。
Java 根据经纬度和日期计算出日出和日落时间
qq_40999938的博客
09-01 5190
package com.rtmc.cloud.commons.utils.calculationSunriseAndSunsetTime; import java.math.BigDecimal; import java.text.SimpleDateFormat; import java.util.Date; public class SunRiseSet { private s...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

风雨依依

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

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

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

打赏作者

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

抵扣说明：

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

余额充值