Чтение файла GPX
P = gpxread(filename)filenameи возвращает вектор географической точки n-by-1, P, где n - количество ППМ или точек, определяющих маршрут или трассу.
gpxread ищет в файле сначала ППМ, затем маршруты, а затем дорожки и возвращает первый тип найденных данных. Поле метаданных P определяет тип элемента ('waypoint', 'track', или 'route') и любые дополнительные метаданные, связанные с ППМ, маршрутом или дорожкой. Если файл содержит несколько путей или маршрутов, P содержит точки, определяющие первую дорожку или маршрут в файле. Если gpxread не удается найти элементы в файле, он возвращает пустой вектор географической точки.
___ = 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');Вычислить 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 возвращает данные в объекте geoshape.
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)
Прочитайте ППМ и журнал отслеживания из sample_mixed.gpx файл.
wpt = gpxread('sample_mixed'); trk = gpxread('sample_mixed', 'FeatureType', 'track');
Отображение ППМ и журнала дорожек на веб-карте.
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))]});
За исключением расширений, GPX версии 1.1 полностью поддерживается. При обнаружении другой версии выдается предупреждение. Однако в большинстве случаев файлы GPX версии 1.0 можно считать успешно, если они не содержат определенных тегов метаданных. Дополнительные сведения см. в документации по схеме GPX 1.1.