gpx_project/test15.py at master · swjmy/gpx_project · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
import gpxpy
import gpxpy.gpx
from geopy.distance import vincenty


gpx_file=open('data\CAStd_116.gpx','r+')
gpx=gpxpy.parse(gpx_file)
gpx_file.close()

print('time_data/n',gpx.get_moving_data())

gpx_file1=open('data\CAStd_116_out.gpx', "w") #打开文件，进行写操作
out=gpxpy.gpx.GPX()


#分割
def splite(start,avg,seg,new_trk):
    i_p=0
    pre_point=None
    new_seg=gpxpy.gpx.GPXTrackSegment()

    for point in seg.points[start:]:
        i_p+=1
        #访问到segment.points中第几个point
        index = i_p+start
        #判断distance
        if point.time:
            if pre_point:
                newport_ri = (point.latitude, point.longitude)
                distance = vincenty(pre_point, newport_ri).meters
                if distance > avg:
                    # 出口1：查出不合格的distance，将当前new_seg加入
                    new_trk.segments.append(new_seg)
                    # 继续后面的检查
                    splite(index - 1, avg=avg, seg=seg, new_trk=new_trk)
                    break
                else:
                    new_seg.points.append(seg.points[index - 1])
            else:
                new_seg.points.append(seg.points[index - 1])
            if index == len(seg.points):
                # 出口2：已检查完全部point，将最后一个new_seg加入
                new_trk.segments.append(new_seg)
            pre_point = (point.latitude, point.longitude)


#计算平均值的三倍
i_trk=0
for track in gpx.tracks:
    i_trk+=1
    i_seg=0
    new_trk=gpxpy.gpx.GPXTrack()

    for segment in track.segments:
        print(type(segment.points))
        i_seg+=1
        pre_point=None
        i_p=0
        # 记录量：总距离
        total_dis = 0
        #判断是否有时间的flag
        #第一次循环计算出avg*3
        for point in segment.points:
            if point.time:
                i_p += 1
                if pre_point:
                    newport_ri = (point.latitude, point.longitude)
                    distance = vincenty(pre_point, newport_ri).meters
                    total_dis += distance
                pre_point = (point.latitude, point.longitude)

        if i_p != 0:
            avg = (total_dis / (i_p - 1)) * 3
            #splite(start=0, avg=avg, seg=segment, new_trk=new_trk)
        else:
            print('Empty Segment')
        #分割
        splite(start=0, avg=avg, seg=segment, new_trk=new_trk)
    if new_trk.segments:
        out.tracks.append(new_trk)


gpx_file1.write(out.to_xml())
#print(avg_list)