gpxファイルのグラフ化プログラム~ver.0~

<trk>
<name>冬山大山　その2</name>
<trkseg>
<trkpt lat="35.394167" lon="133.529682">
<ele>746.000000</ele>
<time>2019-03-01T22:52:09Z</time>
<speed>0.000000</speed>
<trueHeading>0.000000</trueHeading>
<magneticHeading>0.000000</magneticHeading>
<gpsHeading>0.000000</gpsHeading>
<haccuracy>17.312000</haccuracy>
<vaccuracy>17.312000</vaccuracy>
<internetconnect>2</internetconnect>
</trkpt>
．．．

# simple.py

import gpxpy
import matplotlib.pyplot as plt

gpx_file = open('20190302_daisen_winter_2.gpx', 'r')
gpx = gpxpy.parse(gpx_file)

ele = []

for track in gpx.tracks:
    for segment in track.segments:
        for point in segment.points:
            ele.append(point.elevation)

plt.figure()
plt.plot(ele, color = 'blue', lw = 1, alpha = 1)
plt.show()

ヒュベニの公式

D=sqrt((M*dP)*(M*dP)+(N*cos(P)*dR)*(N*cos(P)*dR))

       D:  ２点間の距離(m)
       P:  ２点の平均緯度
       dP: ２点の緯度差
       dR: ２点の経度差
       M:  子午線曲率半径
       N:  卯酉線曲率半径
       M=6334834/sqrt((1-0.006674*sin(P)*sin(P))^3)
       N=6377397/sqrt(1-0.006674*sin(P)*sin(P))

# simple_2.py

import gpxpy
import matplotlib.pyplot as plt
import numpy as np

# gpxファイルの指定
gpx_file = open('20190302_daisen_winter_2.gpx', 'r')
gpx = gpxpy.parse(gpx_file)

# ヒュベニの公式
def hyubeni(lat_old,lat_new,lon_old,lon_new):
	lat_old = lat_old * np.pi / 180
	lat_new = lat_new * np.pi / 180
	lon_old = lon_old * np.pi / 180
	lon_new = lon_new * np.pi / 180
	elat = lat_new - lat_old
	elon = lon_new - lon_old
	mlat = elat/2
	
	M = 6334834/ np.sqrt((1- 0.006674 * np.sin(mlat) * np.sin(mlat)**3))
	N = 6377397/ np.sqrt(1-0.006674 * np.sin(mlat) * np.sin(mlat))
	
	dis = np.sqrt((M * elon) * (M * elon) + (N * np.cos(mlat) * elon) * (N * np.cos(mlat) * elon))
	return(dis)

lat = []
lon = []
ele = []
dis = []
dis_all = []

# 緯度，経度，距離データの格納
for track in gpx.tracks:
    for segment in track.segments:
        for point in segment.points:
            lat.append(point.latitude)
            lon.append(point.longitude)
            ele.append(point.elevation)

# 二点間の距離を計算
for count in range(len(lon)):
	if count == 0:
		dis.append(0)
	else:
		lon_old = lon[count-1]
		lon_new = lon[count]
		lat_old = lat[count-1]
		lat_new = lat[count]
		dis.append(hyubeni(lat_old,lat_new,lon_old,lon_new))
# 距離を積算し，走行距離の計算と格納
for count in range(len(lon)):
	if count == 0:
		dis_all.append(dis[count])
	else:
		dis_all.append(dis_all[count-1]+dis[count])
# 距離の単位をkmに
for count in range(len(lon)):
		dis_all[count] = dis_all[count]/1000

# 数値を確認
print(dis_all)

# グラフの出力  
plt.figure()
#標高，データ数
plt.subplot(211)
plt.plot(ele, color = 'blue', lw = 1, alpha = 1)
#標高，距離
plt.subplot(212)
plt.plot(dis_all, ele, color = 'red', lw = 1, alpha = 1)
plt.show()