gpxread
读取 GPX 文件
说明
___ = gpxread(___, 从 GPX 文件中读取数据，并通过一个或多个 Name,Value)Name,Value 参量对来指定附加选项，这些选项用于控制导入过程的各项特性。Name 是参量名称，Value 是对应的值。Name 必须位于单引号 ('') 之内，且不区分大小写。您可以按任意顺序指定多个名称-值参量。
示例
从 boston_placenames.gpx 文件中读取并显示航点，并将这些点叠加到 boston.tif 图像上。
首先，导入航点和 GeoTIFF 文件。
p = gpxread("boston_placenames.gpx"); [A,R] = readgeoraster("boston.tif");
获取 GeoTIFF 图像的投影结构。将 X 和 Y 限制的长度单位转换为米，以便与投影结构配合使用。
proj = geotiffinfo("boston.tif");
mstruct = geotiff2mstruct(proj);
R.XWorldLimits = R.XWorldLimits * proj.UOMLengthInMeters;
R.YWorldLimits = R.YWorldLimits * proj.UOMLengthInMeters;显示地图。
axesm(mstruct) mapshow(A,R)
显示每个点的名称和位置。
for k=1:length(p) textm(p(k).Latitude, p(k).Longitude, p(k).Name, ... Color="k",BackgroundColor=[0.9 0.9 0],... Interpreter="none"); end geoshow(p) xlim(R.XWorldLimits) ylim(R.YWorldLimits)
读取并显示从波士顿洛根国际机场到马萨诸塞州纳蒂克市 MathWorks 的路线。
从 GPX 文件中读取路线信息。
route = gpxread('sample_route.gpx');使用 0.05 的缓冲量计算 latlim 和 lonlim。
[latlim, lonlim] = geoquadline(route.Latitude, route.Longitude); [latlim, lonlim] = bufgeoquad(latlim, lonlim, .05, .05);
显示路线。
fig = figure;
pos = fig.Position;
fig.Position = [300 300 1.25*pos(3) 1.25*pos(4)];
ax = usamap(latlim, lonlim);
setm(ax, 'MLabelParallel', 43.5)
geoshow(route.Latitude, route.Longitude)
提取路线中包含转弯描述的元素，在地图上为每个转弯标注标记和颜色代码，并制作一个显示这些描述的图例。将顺序颠倒过来，使图例中第一圈显示在顶部，最后一圈显示在底部。
turns = route(~cellfun(@isempty, route.Description)); turns = turns(end:-1:1); n = length(turns); colors = cool(n); for k=1:n geoshow(turns(k).Latitude, turns(k).Longitude, ... 'DisplayType','point','MarkerEdgeColor',colors(k,:),... 'Tag','turn','DisplayName',turns(k).Description) end legend(findobj(ax,'Tag','turn'),'Location','SouthOutside')
读取 GPX 文件中的轨迹日志。gpxread 函数返回 geoshape 对象中的数据。
idx = 1:2;
tracks = gpxread("sample_tracks.gpx",Index=idx);在地图上显示轨迹。为每条轨迹日志使用不同的颜色。
figure geobasemap streets-light hold on for k = idx track = tracks(k); lat = track.Latitude; lon = track.Longitude; geoplot(lat,lon) end
放大到围绕第一条轨迹的区域。
[latlim,lonlim] = geoquadline(tracks(1).Latitude,tracks(1).Longitude); [latlim,lonlim] = bufgeoquad(latlim,lonlim,0.0005,0.0005); geolimits(latlim,lonlim)
从 GPX 文件中读取轨迹日志和航点。
trk = gpxread("sample_mixed.gpx",FeatureType="track"); wpt = gpxread("sample_mixed.gpx");
从 geopoint 向量中提取坐标。然后，将坐标显示在地图上。使用一条线表示轨迹，用图标表示航点。
figure geobasemap topographic hold on trkLat = trk.Latitude; trkLon = trk.Longitude; geoplot(trkLat,trkLon,DisplayName="Track Log") wptLat = wpt.Latitude; wptLon = wpt.Longitude; geoiconchart(wptLat,wptLon,DisplayName="Waypoints") legend
本示例演示了如何显示海拔和时间区域图，以及如何利用轨迹日志计算距离。
读取 sample_mixed.gpx 文件中的轨迹日志。
trk = gpxread('sample_mixed.gpx','FeatureType','track');
通过去除秒数的小数部分和时区偏移量，简化时间值字符串的格式。将字符串转换为 datetime 数组。
timeStr = strrep(trk.Time,'.000Z',''); t = datetime(timeStr);
绘制时间和海拔值的值图。
figure area(t,trk.Elevation) xtickformat('HH:mm:ss') xlabel('time (hours:minutes:seconds)') ylabel('elevation (meters)') title('Elevation Area Plot');
计算并显示地面轨迹距离。将以米为单位的距离转换为以美国测绘英里为单位的距离。
e = wgs84Ellipsoid; lat = trk.Latitude; lon = trk.Longitude; d = distance(lat(1:end-1), lon(1:end-1), lat(2:end), lon(2:end), e); d = d * unitsratio('sm', 'meter');
显示累计地面轨迹距离和经过时间。
elapsedTime = t - t(1); figure line(elapsedTime(2:end),cumsum(d)) ylabel('cumulative ground track distance (statute mile)') xlabel('elapsed time (hours:minutes:seconds)') title({'Cumulative Ground Track Distance in Miles', ... ['Total Distance in Miles: ' num2str(sum(d))]});
输入参数
名称-值参数
输出参量
详细信息
提示
不包括扩展名，GPX 1.1 版本得到全面支持。如果检测到其他版本，系统将发出警告。不过，在大多数情况下，1.0 版本的 GPX 文件都能被成功读取，除非其中包含某些元数据标签。有关详细信息，请参阅 GPX 1.1 架构文档。
版本历史记录
在 R2012a 中推出
