(最近在做数学建模的大报告,自己记录一下)
下面是使用python将地球表面一点的经纬度转换为位置矢量:
import math
def degrees_to_radians(degrees):
return degrees * math.pi / 180.0
# 经纬度 -> 位置矢量
def convert_to_cartesian(longitude, latitude):
# 将经纬度转换为弧度
lat_radians = degrees_to_radians(latitude)
lon_radians = degrees_to_radians(longitude)
earth_radius = 6371.0 # 地球半径(单位:km)
# 计算直角坐标系中的坐标
x = earth_radius * math.cos(lat_radians) * math.cos(lon_radians)
y = earth_radius * math.cos(lat_radians) * math.sin(lon_radians)
z = earth_radius * math.sin(lat_radians)
return x, y, z