C语言 根据经纬度如何获取日落日出时间

最新推荐文章于 2024-08-30 08:23:32 发布
最新推荐文章于 2024-08-30 08:23:32 发布
阅读量1.7k 收藏 22
点赞数 3
文章标签： C 日出时间 日落时间
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011983700/article/details/109669320
版权
这是一个使用C语言编写的程序，能够根据输入的日期、经纬度计算日出和日落时间。程序中包含了闰年的判断、日期处理、地理坐标转换等关键算法，并通过不断迭代调整结果以减小误差。最终，程序会输出格林威治标准时间的日出和日落时间，并转换为当地时区的时间。
摘要由CSDN通过智能技术生成

    代码参考

    利用日期、经纬度求日出日落时间 C语言程序代码_百度知道

    源代码有几次错误，已经改成可以使用的C版本，计算结果和查询网站的结果会有1分钟左右的误差。

    测试网站地址日出日落时间查询 - 国旗标准升降时间

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define PI 3.1415926

//定义全局变量
int days_of_month_1[] ={ 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 };
int days_of_month_2[] ={ 31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31 };
long double h = -0.833;

//输入日期
void input_date(int *c)
{
	printf("Enter the date (form: 2009 03 10):\n");
	scanf("%d %d %d",&c[0],&c[1],&c[2]);
}
//输入纬度
void input_glat(int *c)
{
	printf("Enter the degree of latitude(range: 0°- 60°,form: 40 40 40 (means 40°40′40″)):\n");
	scanf("%d %d %d",&c[0],&c[1],&c[2]);
}
//输入经度
void input_glong(int *c)
{
	printf("Enter the degree of longitude(west is negativ,form: 40 40 40 (means 40°40′40″)):\n");
	scanf("%d %d %d",&c[0],&c[1],&c[2]);
}
//判断是否为闰年:若为闰年，返回1；若非闰年，返回0
int leap_year(int year)
{
	if ((year % 400 == 0) || ((year % 100 != 0) && (year % 4 == 0)))
		return 1;
	else
		return 0;
}
//求从格林威治时间公元2000年1月1日到计算日天数days
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
 double t_century(int days,  double UTo)
{
	return (( double) days + UTo / 360) / 36525;
}
//求太阳的平黄径
 double L_sun( double t_century)
{
	return (280.460 + 36000.770 * t_century);
}
//求太阳的平近点角
 double G_sun( double t_century)
{
	return (357.528 + 35999.050 * t_century);
}
//求黄道经度
 double ecliptic_longitude( double L_sun,  double G_sun)
{
	return (L_sun + 1.915 * sin(G_sun * PI / 180)
			+ 0.02 * sin(2 * G_sun * PI / 180));
}
//求地球倾角
 double earth_tilt( double t_century)
{
	return (23.4393 - 0.0130 * t_century);
}
//求太阳偏差
 double sun_deviation( double earth_tilt,
		 double ecliptic_longitude)
{
	return (180 / PI* asin(sin(PI / 180 * earth_tilt)* sin(PI / 180 * ecliptic_longitude)));
}
//求格林威治时间的太阳时间角GHA
 double GHA( double UTo,  double G_sun,
		 double ecliptic_longitude)
{
	return (UTo - 180 - 1.915 * sin(G_sun * PI / 180)
			- 0.02 * sin(2 * G_sun * PI / 180)
			+ 2.466 * sin(2 * ecliptic_longitude * PI / 180)
			- 0.053 * sin(4 * ecliptic_longitude * PI / 180));
}
//求修正值e
 double e( double h,  double glat,  double sun_deviation)
{
	return 180 / PI* acos((sin(h * PI / 180)- sin(glat * PI / 180)* sin(sun_deviation * PI / 180))/ (cos(glat * PI / 180)* cos(sun_deviation * PI / 180)));
}
//求日出时间
 double UT_rise( double UTo,  double GHA,  double glong,
		 double e)
{
	return (UTo - (GHA + glong + e));
}
//求日落时间
 double UT_set( double UTo,  double GHA,  double glong,
		 double e)
{
	return (UTo - (GHA + glong - e));
}
//判断并返回结果（日出）
 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;
}
//判断并返回结果（日落）
 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;
}
//求时区
int Zone( double glong)
{
    int timeZone ;
    int shangValue = (int)(glong / 15);
    double yushuValue =  glong - shangValue*15;
//    printf("%lf\n",yushuValue);
    if (yushuValue <= 7.5) {
        timeZone = shangValue;
    } else {
        timeZone = shangValue +(glong > 0 ?  1 :-1);
    }
    return timeZone;
}

void output( double rise,  double set,  double glong)
{
	if ((int) (60 * (rise / 15 + Zone(glong) - (int) (rise / 15 + Zone(glong))))< 10)
	{
		printf("The time at which the sun rises is %d:0%d \n",(int) (rise / 15 + Zone(glong)),
				(int) (60* (rise / 15 + Zone(glong)- (int) (rise / 15 + Zone(glong)))));
	}
	else
	{
		printf("The time at which the sun rises is %d:%d \n",(int) (rise / 15 + Zone(glong)),
				(int) (60* (rise / 15 + Zone(glong)- (int) (rise / 15 + Zone(glong)))));
	}

	if ((int) (60 * (set / 15 + Zone(glong) - (int) (set / 15 + Zone(glong))))< 10)
	{
		printf("The time at which the sun sets is %d:%d \n",(int) (set / 15 + Zone(glong)),
				(int) (60* (set / 15 + Zone(glong)- (int) (set / 15 + Zone(glong)))));
	}
	else
	{
		printf("The time at which the sun sets is %d:%d \n",(int) (set / 15 + Zone(glong)),
				(int) (60* (set / 15 + Zone(glong)- (int) (set / 15 + Zone(glong)))));
	}
}
//打印结果
int main()
{
	double UTo = 180.0;
	int year, month, date;
	double glat, glong;
	int c[3];
//	input_date(c);
	year = 2020;
	month =11;
	date = 13;
	printf("year=%d month=%d date=%d\n",year,month,date);
//	input_glat(c);
//	glat = c[0] + c[1] / 60 + c[2] / 3600;
	glat=29.026046;
	printf("glat=%lf\n",glat);
//	input_glong(c);
//	glong = c[0] + c[1] / 60 + c[2] / 3600;
	glong=121.126078;
	printf("glong=%lf\n",glong);

	double rise, set;
	rise = 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);
	set = 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);
	output(rise, set, glong);
//	system("pause");
	return 0;
}

juelianhuayao
关注
c语言计算日出日落时间_计算日出日落算法实现
weixin_39765280的博客
12-20 1132
首先感谢帮我算题的小玉玉~算了一周的日出日落，倍觉痛苦，写下来防止下次再算。根据经纬度和UTC计算日出日落的思路、算法及scala代码如下~日出和日落的计算思路：一天的日出和日落可以根据正午时刻平分昼长得到。正午时刻可以根据所在时区和经度得到。根据经度可以知道所在时区。昼长可以根据太阳赤角和所在纬度计算得到。符号说明：rc：日出rl：日落zw：正午day：昼长zone：时区δ：太阳赤纬lon：经度...
c语言计算日出日落时间_利用日期、经纬度求日出日落时间 C语言程序代码
weixin_39637256的博客
12-20 1349
展开全部#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...
经纬度计算日出日落时间(C语言)
04-12
计算日出日落时间：基于经纬度的精准算法库
最新发布
gitblog_00914的博客
08-30 721
计算日出日落时间：基于经纬度的精准算法库 项目地址:https://gitcode.com/gh_mirrors/su/sunset 在寻找如何精确判断天色是否转暗的过程中，我们发现了这样一项杰出的开源项目——一个基于 Mike Chirico 在2004年发布的 sunrise.c 的进化版本。这个库自那以后经历了多次迭代，现今已成为能够准确计算多种日出日落时刻的强大工具，为开发者们提供了前所未...
经纬度计算日出日落C代码
02-25
1，以C语言编写的经纬度日出时刻计算代码 2，已经成功的移植到嵌入式设备中
根据经纬度计算日出日落时间
10-10
根据经纬度计算日出日落时间。只需要经度和纬度。
c语言计算日出日落时间_日出日落时间计算程序（C语言） -
weixin_40008033的博客
12-20 922
long 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),UT...
日出日落时间计算程序(C语言)_C语言_算法_
09-30
在这个特定的案例中，"日出日落时间计算程序(C语言)"是一个利用C语言编写的程序，其目标是根据地理位置的经纬度以及万年历数据来计算特定地点的日出和日落时间。 日出日落时间的计算涉及到天文学知识，主要是地球...
根据经纬度和日期计算日出日落时间
热门推荐
liu877260630的博客
05-28 3万+
接计算日出和日落时间会比较复杂,可以通过计算当地的“正午”时间、当地的“白天”时长,用正午时间平分白天时长,得到当地的日出时间和日落时间。计算正午时间一个地方的正午时间,可以通过当地经度与当地所属时区经度的差计算得到。全世界有24个时区,每个时区相隔1小时,因此每个时区之间相差的经度是 360 / 24 = 15°。已北京为例,北京所属时区为东8区,东8区的中心经度是15 * 8 = 120°(因...
根据经纬度计算日出日落时间算法实现(C#)
05-09
本代码已经在项目中应用。 压缩包里附送 日升日落时间计算(2.0版).xls
经纬度计算日出日落时间
04-17
C语言，根据经纬度计算日出、日落、天黑、天亮时间
日出日落时间计算程序(C语言).doc
09-30
日出日落时间计算程序(C语言).doc
日出日落时间算法实现(C#)
06-12
C#完整实现日出日落时间计算，代码可直接使用。欢迎交流指正。
c语言计算日出日落时间_在C日出日落时间
weixin_39678089的博客
12-20 252
In my C application, I want to calculate sunrise/sunset times for a given date, latitude and longitude. i have been searching on the net but i can not find a working sample.But this sample didnt work ...
天黑时间跟经度还是纬度有关_根据经纬度计算日出、日落、中天、天亮、天黑和昼长时间...
weixin_39870413的博客
12-21 501
importjava.io.IOException;importjava.util.Date;importjava.util.HashMap;importjava.util.Map;importjava.util.Properties;public classDayTime {private staticProperties propt;private static double RAD = 18...
C#可用的日出日落时间类
weixin_34059951的博客
09-05 419
  一个现成代码的公共类库，复制下来作为一个类文件就可以调用了。一般不需要了解实现过程，各种数学公式太麻烦。   调用方法： SunTimeResult result = SunTimes.GetSunTime(DateTime.Now, double.Parse(Longitude), double.Parse(Latitude)); MessageBox.Show(string....
基于GPS经纬度和当地时间计算日落日出时间实现
mygod2008ok的专栏
10-17 1855
以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
C语言实现：根据日期和经纬度计算日落时间程序
本资源是一份名为"利用日期,经纬度求日落时间的C语言程序文件.pdf"的文档，它详细介绍了如何使用C语言编写一个程序来计算给定日期和地理位置的日出和日落时间。程序的关键知识点包括： 1. **C语言编程基础**：文件...
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值