Считайте файл 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');Вычислите 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 и отображения, наложенного на webmap.
Считайте журналы дорожки из файла GPX. gpxread возвращает данные в геообъекте фигуры.
tracks = gpxread('sample_tracks', 'Index', 1:2);
Отобразите вход в систему дорожки webmap с различным цветом для каждого журнала дорожки.
webmap('openstreetmap') colors = {'cyan', 'red'}; wmline(tracks, 'Color', colors)
Масштабируйте webmap, чтобы просмотреть первую дорожку около кампуса MathWorks в Натике.
[latlim, lonlim] = geoquadline(tracks(1).Latitude, tracks(1).Longitude); wmlimits(latlim, lonlim)
Считайте waypoints и отследите журнал от sample_mixed.gpx файл.
wpt = gpxread('sample_mixed'); trk = gpxread('sample_mixed', 'FeatureType', 'track');
Отобразите waypoints и вход в систему дорожки webmap.
webmap('worldimagery') wmline(trk, 'OverlayName', 'Track Logs');
Добавьте веб-маркеры, чтобы отметить положения каждого пути точка.
wmmarker(wpt, 'FeatureName', wpt.Name, 'OverlayName', 'Waypoints')
В этом примере показано, как отобразить вертикальное изменение и карты области времени и вычислить расстояние с помощью журналов дорожки.
Считайте журнал дорожки из sample_mixed.gpx файл.
trk = gpxread('sample_mixed', 'FeatureType', 'track');
Преобразуйте строки временной стоимости в последовательные числа даты с помощью datenum, и затем вычислите время суток в мелких секундах часов.
timeStr = strrep(trk.Time, 'T', ' '); timeStr = strrep(timeStr, '.000Z', ''); trk.DateNumber = datenum(timeStr, 31); day = fix(trk.DateNumber(1)); trk.TimeOfDay = trk.DateNumber - day;
Отобразите график области вертикального изменения и временных стоимостей.
figure area(trk.TimeOfDay, trk.Elevation) datetick('x', 13, 'keepticks', 'keeplimits') ylabel('elevation (meters)') xlabel('time(Z) hours:minutes:seconds') title({'Elevation Area Plot', datestr(day)});
Вычислите и отобразите наземное расстояние дорожки. Преобразуйте расстояние в метрах к расстоянию в американских милях обзора.
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');
Отобразите совокупное наземное расстояние дорожки и прошедшее время.
trk.ElapsedTime = trk.TimeOfDay - trk.TimeOfDay(1); figure line(trk.ElapsedTime(2:end), cumsum(d)) datetick('x', 13) ylabel('cumulative ground track distance (statute mile)') xlabel('elapsed time (hours:minutes:seconds)') title({'Cumulative Ground Track Distance in Miles', datestr(day), ... ['Total Distance in Miles: ' num2str(sum(d))]});
Исключая расширения полностью поддержана версия 1.1 GPX. Если какая-либо другая версия обнаруживается, предупреждение выдано. Однако в большинстве случаев, версия 1.0, файлы GPX могут быть считаны успешно, если они не содержат определенные теги метаданных. Для получения дополнительной информации см. документацию GPX 1.1 Схемы.