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

最新推荐文章于 2020-12-20 07:47:23 发布
最新推荐文章于 2020-12-20 07:47:23 发布
阅读量1.5k 收藏 19
点赞数 2
文章标签： C 日出时间 日落时间
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011983700/article/details/109669320
版权

    代码参考

    利用日期、经纬度求日出日落时间 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
关注
  • 2
    点赞
  • 踩
  • 19
    收藏
    觉得还不错? 一键收藏
  • 5
    评论
日出日落时间计算计算
qq_33169259的博客
04-26 158
【代码】日出日落时间计算计算。
根据经纬度计算日出、日落、中天、天亮、天黑和昼长时间
weixin_30612769的博客
07-04 2434
Jean Meeus的《天文算法》（Astronomical Algorithms，2nd Edition）第二版中第7章第60页内有详细介绍计算儒略日的方法： 设Y为给定年份，M为月份，D为该月日期（带小数，把时：分：秒折算成日的形式）。运算符INT表示为取所给数的整数部分，也即小数点前的部分。 1.若M > 2，Y和M不变。 若 M =1或2，以Y–1代Y，以M+12代M。 换句...
日出日落时间计算程序(C语言).doc
09-30
日出日落时间计算程序(C语言).doc
根据经纬度和日期计算日出日落时间
热门推荐
azoo的专栏
06-28 2万+
计算日出和日落时间会比较复杂,可以通过计算当地的“正午”时间、当地的“白天”时长,用正午时间平分白天时长,得到当地的日出时间和日落时间。 1.计算正午时间 一个地方的正午时间,可以通过当地经度与当地所属时区经度的差计算得到。全世界有24个时区,每个时区相隔1小时,因此每个时区之间相差的经度是 360 / 24 = 15°。已北京为例,北京所属时区为东8区,东8区的中心经度是15 * 8 = 120°(因此东8区的范围是112.5° - 127.5°)。 北京的经度是116.4...
经纬度计算日出日落时间(C语言)
04-12
c语言计算日出日落时间_利用日期、经纬度求日出日落时间 C语言程序代码
weixin_39637256的博客
12-20 1228
展开全部#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语言计算日出日落时间_日出日落时间计算程序（C语言） -
weixin_40008033的博客
12-20 783
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代码
02-25
1，以C语言编写的经纬度日出时刻计算代码 2，已经成功的移植到嵌入式设备中
经纬度计算日出日落时间
04-17
C语言，根据经纬度计算日出、日落、天黑、天亮时间
日出日落时间计算程序(C语言)_C语言_算法_
09-30
通过经纬度以及万年历计算当地的日落日出时间
日出日落时间计算(sunset)
12-06
C#项目，根据经纬度计算当地日出时间和日落时间。
C# 计算日出时间日落时间太阳倾角仰角
11-16
一份计算日出时间日落时间太阳倾角仰角的代码工具类。
C语言获取系统时间的几种方式
01-31
核心提示：C语言中如何获取时间？精度如何？1使用time_ttime(time_t*timer)精确到秒2使用clock_tclock()得到的是CPU时间精确到1/CLOCKS_PER_SEC秒3计算时间差使用doubledifftime(...1使用time_ttime(time_t*timer)...
如何根据纬度，经度和日期判断太阳的日出和日落的时间
smilestone322的专栏
01-30 6355
如题： 日出时间计算公式： 24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360 我国时区为东8区，时区=8 经度、纬度采用角度制，东经、北纬为正，西经、南纬为负 日期序列数为当天在这一年中的序列，如2月11日就是42 因为计算机一般采用弧度制，上面公式可变化为： 24*(180+时区*15
计算机语言发展史
weixin_45989987的博客
09-11 532
计算机语言发展史 第一代语言 机器语言 计算机的基本计算方式都是基于二进制的方式。 二进制：010111001010110010110100 这种代码是直接输入给计算机使用的，不经过任何的转换！ 1 1 2 10 3 11 4 100 5 101 6 110 7 111 8 1000 16 10000 32 100000 第二代语言 汇编语言 解决人类无法读懂机器语言的问题 指令代替二进制 codesg segment start: mov ax, 0123h
C语言时间计算
qq_25677349的博客
12-01 4385
时间计算 1. clock函数 精确度为秒 clock函数返回从程序开始执行起处理器始终滴答的次数，为了把这个值转化为秒，应该除以CLOCKS_PER_SEC #include<time.h> #include<stdio.h> #include<unistd.h> int main() { clock_t t_begin,t_end; t_...
c语言计算日出日落时间的代码
最新发布
06-11
以下是使用NOAA的Solar Calculator库计算日出日落时间的完整代码示例： ``` #include <stdio.h> #include "solar_calculator.h" int main() { double lat = 37.7749; // 纬度 double lon = -122.4194; // 经度 int year = 2021; // 年份 int month = 8; // 月份 int day = 25; // 日期 int tz = -7; // 时区，如美国西部时区为-7 double sunrise = calculate_sunrise(lat, lon, year, month, day, tz); double sunset = calculate_sunset(lat, lon, year, month, day, tz); printf("日出时间：%.2f\n", sunrise); printf("日落时间：%.2f\n", sunset); return 0; } ``` 需要注意的是，计算日出日落时间的精度与数据来源和计算方法有关，可能存在一定误差。如果需要更高精度的计算，可以使用更复杂的算法和数据源。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值