gpxread
Считайте файл GPX
Синтаксис
Описание
P = gpxread(filename)filename, и возвращает n-by-1 вектор геоточки, P, где n является количеством waypoints или точками, которые задают маршрут или дорожку.
gpxread ищет файл сначала waypoints, затем направляет, и затем отслеживает, и это возвращает первый тип данных, которые это находит. Поле метаданных P идентифицирует тип функции ('waypoint', 'track', или 'route') и любые дополнительные метаданные, сопоставленные с waypoint, маршрутом или дорожкой. Если файл содержит несколько дорожек или маршрутов, P содержит точки, которые задают первую дорожку или маршрут в файле. Если gpxread не может найти функции в файле, он возвращает пустой вектор геоточки.
S = gpxread(___,'Index',V)'Index' как вектор, V. Используйте этот синтаксис, когда это необходимо, чтобы работать с данными как линия, а не как набор точек.
___ = gpxread(___, считывает данные из файла GPX с дополнительными опциями, заданными одним или несколькими Name,Value)Name,Value парные аргументы, то управление различные характеристики импорта. Name имя аргумента и Value соответствующее значение. Name должен появиться в одинарных кавычках ('') и является нечувствительным к регистру. Можно задать несколько аргументов пары "имя-значение" в любом порядке.
Примеры
Рид Вейпойнтс и отображает их по изображению
Считайте и отобразите waypoints от boston_placenames.gpx файл и наложение точки на boston.tif изображение.
Во-первых, импортируйте файл GeoTIFF и waypoints.
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 в Натике, MA.
Считайте информации маршрута из файла GPX.
route = gpxread('sample_route.gpx');Вычислите latlim и lonlim с 0,05 буферами.
[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 и отображения, наложенного на веб-карте.
Считайте журналы дорожки из файла GPX. gpxread возвращает данные в геообъекте фигуры.
tracks = gpxread('sample_tracks', 'Index', 1:2);
Отобразите вход в систему дорожки веб-карта с различным цветом для каждого журнала дорожки.
webmap('openstreetmap') colors = {'cyan', 'red'}; wmline(tracks, 'Color', colors)
Масштабируйте веб-карту, чтобы просмотреть первую дорожку около кампуса MathWorks в Натике.
[latlim, lonlim] = geoquadline(tracks(1).Latitude, tracks(1).Longitude); wmlimits(latlim, lonlim)
Считайте и отобразите Waypoints и отследите вход в систему веб-карта
Считайте waypoints и отследите журнал от sample_mixed.gpx файл.
wpt = gpxread('sample_mixed'); trk = gpxread('sample_mixed', 'FeatureType', 'track');
Отобразите waypoints и вход в систему дорожки веб-карта.
webmap('worldimagery') wmline(trk, 'OverlayName', 'Track Logs');
Добавьте веб-маркеры, чтобы отметить положения каждого пути точка.
wmmarker(wpt, 'FeatureName', wpt.Name, 'OverlayName', 'Waypoints')
Отобразите карты области времени и вертикальное изменение
В этом примере показано, как отобразить вертикальное изменение и карты области времени и вычислить расстояние с помощью журналов дорожки.
Считайте журнал дорожки из 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))]});
Входные параметры
Выходные аргументы
Больше о
Советы
Исключая расширения полностью поддерживается версия 1.1 GPX. Если какая-либо другая версия обнаруживается, предупреждение выдано. Однако в большинстве случаев, версия 1.0, файлы GPX могут быть считаны успешно, если они не содержат определенные теги метаданных. Для получения дополнительной информации см. документацию GPX 1.1 Схемы.