function [data_GPS,str]=read_GPS(gpxname);
% by Longqi Yang
% 读取GPX文件时,先将头三行给删掉
% gpxname为文件名称
fid=fopen(gpxname,'r');
str='<name>';
data_GPS=[];
while (~strncmp(str,'</gpx>',6))
str=fgetl(fid);
if (strncmp(str,' <trkpt',9))
lat=sscanf(str(16:24),'%f');
lon=sscanf(str(32:41),'%f');
elseif (strncmp(str,' <time>',10))
month=sscanf(str(16:17),'%d');
day=sscanf(str(19:20),'%d');
hour=sscanf(str(22:23),'%d');
minute=sscanf(str(25:26),'%d');
second=sscanf(str(28:29),'%d');
time=datenum(0000,month,day,hour,minute,second)+8/24;
data_1=[time lon lat];
data_GPS=[data_GPS;data_1];
end
end
return