go-sunrise 库使用指南：计算指定位置日出日落时间

go-sunrise 是一个用于计算地球上任意位置日出日落时间的 Go 语言库。它基于精确的天文学算法，可以计算出特定日期和地理坐标下的日出、日落、正午等时间。

安装

go get github.com/nathan-osman/go-sunrise

基本用法

1. 计算日出日落时间

package main

import (
	"fmt"
	"time"
	
	"github.com/nathan-osman/go-sunrise/sunrise"
)

func main() {
	// 设置地理位置（纬度，经度）
	lat := 39.9042  // 北京纬度
	long := 116.4074 // 北京经度

	// 设置日期（年，月，日）
	year := 2023
	month := time.November
	day := 15

	// 计算日出日落时间
	sunriseTime, sunsetTime := sunrise.SunriseSunset(
		lat, long,
		year, month, day,
	)

	fmt.Printf("北京 %d年%d月%d日:\n", year, month, day)
	fmt.Printf("日出: %02d:%02d\n", sunriseTime.Hour(), sunriseTime.Minute())
	fmt.Printf("日落: %02d:%02d\n", sunsetTime.Hour(), sunsetTime.Minute())
}

2. 计算更详细的天文时间

func main() {
	lat := 31.2304  // 上海纬度
	long := 121.4737 // 上海经度
	year := 2023
	month := time.December
	day := 21 // 冬至日

	// 计算更详细的天文时间
	sunTimes := sunrise.Calculate(
		lat, long,
		year, month, day,
	)

	fmt.Printf("上海 %d年%d月%d日:\n", year, month, day)
	fmt.Printf("日出: %02d:%02d\n", sunTimes.Sunrise.Hour(), sunTimes.Sunrise.Minute())
	fmt.Printf("日落: %02d:%02d\n", sunTimes.Sunset.Hour(), sunTimes.Sunset.Minute())
	fmt.Printf("正午: %02d:%02d\n", sunTimes.SolarNoon.Hour(), sunTimes.SolarNoon.Minute())
	fmt.Printf("日照时长: %.1f 小时\n", sunTimes.DayLength.Hours())
}

3. 计算特定时区的时间

func main() {
	lat := 35.6762  // 东京纬度
	long := 139.6503 // 东京经度
	year := 2023
	month := time.March
	day := 20 // 春分日

	// 计算时间
	sunriseTime, sunsetTime := sunrise.SunriseSunset(
		lat, long,
		year, month, day,
	)

	// 东京时区 (UTC+9)
	loc, _ := time.LoadLocation("Asia/Tokyo")
	sunriseInTokyo := sunriseTime.In(loc)
	sunsetInTokyo := sunsetTime.In(loc)

	fmt.Printf("东京 %d年%d月%d日:\n", year, month, day)
	fmt.Printf("日出: %02d:%02d (JST)\n", sunriseInTokyo.Hour(), sunriseInTokyo.Minute())
	fmt.Printf("日落: %02d:%02d (JST)\n", sunsetInTokyo.Hour(), sunsetInTokyo.Minute())
}

4. 计算多日的日出日落时间

func main() {
	lat := 40.7128  // 纽约纬度
	long := -74.0060 // 纽约经度
	year := 2023
	month := time.July

	// 计算7月每天的日出日落时间
	for day := 1; day <= 31; day++ {
		sunriseTime, sunsetTime := sunrise.SunriseSunset(
			lat, long,
			year, month, day,
		)
		
		// 纽约时区 (UTC-4 夏令时)
		loc, _ := time.LoadLocation("America/New_York")
		sunriseNY := sunriseTime.In(loc)
		sunsetNY := sunsetTime.In(loc)
		
		fmt.Printf("%d-%02d-%02d: 日出 %02d:%02d, 日落 %02d:%02d\n",
			year, month, day,
			sunriseNY.Hour(), sunriseNY.Minute(),
			sunsetNY.Hour(), sunsetNY.Minute(),
		)
	}
}

高级功能

1. 计算民用晨昏时间

func main() {
	lat := 48.8566  // 巴黎纬度
	long := 2.3522   // 巴黎经度
	year := 2023
	month := time.April
	day := 1

	// 计算民用晨昏时间（太阳在地平线下6度时）
	sunTimes := sunrise.CalculateWithElevation(
		lat, long,
		year, month, day,
		-6, // 太阳在地平线下6度（民用晨昏）
	)

	fmt.Printf("巴黎 %d年%d月%d日:\n", year, month, day)
	fmt.Printf("民用晨光始: %02d:%02d\n", sunTimes.Sunrise.Hour(), sunTimes.Sunrise.Minute())
	fmt.Printf("民用昏影终: %02d:%02d\n", sunTimes.Sunset.Hour(), sunTimes.Sunset.Minute())
}

2. 计算日出日落方位角

func main() {
	lat := 34.0522  // 洛杉矶纬度
	long := -118.2437 // 洛杉矶经度
	year := 2023
	month := time.June
	day := 21 // 夏至日

	// 计算日出日落方位角
	_, _, sunriseAzimuth, sunsetAzimuth := sunrise.SunriseSunsetAzimuth(
		lat, long,
		year, month, day,
	)

	fmt.Printf("洛杉矶 %d年%d月%d日:\n", year, month, day)
	fmt.Printf("日出方位角: %.1f° (正北为0°，顺时针增加)\n", sunriseAzimuth)
	fmt.Printf("日落方位角: %.1f°\n", sunsetAzimuth)
}

注意事项

1. 所有时间默认返回UTC时间，使用时需要根据实际时区进行转换
2. 高纬度地区在夏季可能出现极昼（无日落），冬季可能出现极夜（无日出）
3. 计算结果与实际观测可能有几分钟差异，受大气折射、地形等因素影响
4. 对于精确应用场景，建议使用专业天文计算库或API

go-sunrise 是一个轻量级的库，适合大多数基本日出日落计算需求。对于更复杂的天文计算，可以考虑使用更专业的库如 github.com/soniakeys/meeus。