В этом примере показано, как упростить многоугольник и данные о линии с помощью reducem функция. Упрощение многоугольника и данных о линии может ускорить определенные вычисления, не оказывая значимого влияния на данные. Один способ приблизиться к упрощению состоит в том, чтобы использовать reducem со значением допуска по умолчанию сначала и представлением выход. Если результаты не удовлетворяют ваши требования, повторите операцию, увеличившись или уменьшив значение допуска, чтобы достигнуть результата, которого вы желаете. reducem возвращает используемое значение допуска. Для получения дополнительной информации о векторном упрощении данных, смотрите, Упрощают Векторные Координатные Данные.
Извлеките береговые линии Массачусетса и государственные границы от usastatehi.shp файл форм.
ma = shaperead('usastatehi.shp','UseGeoCoords',true, ... 'Selector', {@(name)strcmpi(name,'Massachusetts'),'Name'}); masslat = ma.Lat; masslon = ma.Lon;
Результатом является пара векторов, которые обрисовывают в общих чертах состояние с 957 вершинами.
numel(masslat)
ans = 957
Упростите схему с помощью reducem функция со значением допуска по умолчанию.
[masslat1, masslon1, cerr, tol] = reducem(masslat', masslon');
Смотрите результаты. Количество вершин было сокращено к 252. Векторы уменьшались приблизительно до четверти своих исходных длин.
numel(masslat1)
ans = 252
numel(masslat1)/numel(masslat)
ans = 0.2633
Исследуйте ошибку и значения допуска, возвращенные reducem. cerr значение указывает на тот reducem уменьшал общую длину схемы приблизительно на 3,3 процента. Допуск, что reducem используемый, чтобы достигнуть этого сокращения были 0,006 градуса, или приблизительно 660 метров.
[cerr tol]
ans = 1×2
0.0331 0.0060
Постройте уменьшаемую схему красного цвета по исходной схеме в синем, использующем geoshow. В этом разрешении трудно видеть любое различие между исходной схемой и уменьшаемой схемой.
figure axesm('MapProjection', 'eqdcyl', 'MapLatLim', [41.1 43.0],... 'MapLonLim', [-73.6, -69.8], 'Frame', 'off', 'Grid', 'off'); geoshow(masslat, masslon, 'DisplayType', 'line', 'color', 'blue') geoshow(masslat1, masslon1, 'DisplayType', 'line', 'color', 'red')
Чтобы получить лучший взгляд на два основ, используйте xlim и ylim увеличить масштаб фрагмента карты. Заметьте, как уменьшаемая схема соответствует общим контурам исходной карты, но теряет много детали.
axesm('MapProjection', 'eqdcyl', 'MapLatLim', [41.1 43.0],... 'MapLonLim', [-73.6, -69.8], 'Frame', 'off', 'Grid', 'off'); xlim([0.0104 0.0198]) ylim([0.7202 0.7264]) geoshow(masslat, masslon, 'DisplayType', 'line', 'color', 'blue') geoshow(masslat1, masslon1, 'DisplayType', 'line', 'color', 'red')
Выполните операцию снова, на этот раз удвоив значение допуска.
[masslat2,masslon2,cerr2,tol2] = reducem(masslat', masslon', 0.012); numel(masslat2)
ans = 157
numel(masslat2)/numel(masslat)
ans = 0.1641
Исследуйте ошибку и значения допуска, возвращенные reducem. На этот раз, cerr значение указывает на тот reducem уменьшал общую длину схемы приблизительно на 5,2 процентов. Допуск, что reducem используемый, чтобы достигнуть этого сокращения были 0,012 градуса.
[cerr2 tol2]
ans = 1×2
0.0517 0.0120
Постройте эту уменьшаемую схему темно-зеленого цвета по исходной схеме синего цвета. Отметьте, как эта уменьшаемая схема обеспечивает общую форму исходной карты, но теряет большую часть мелких деталей.
geoshow(masslat2, masslon2, 'DisplayType', 'line', 'color', [0 .6 0])
Увеличьте допуск до 0,1 градусов.
[masslat3, masslon3, cerr3, tol3] = reducem(masslat', masslon', 0.1);
Постройте эту уменьшаемую схему черного цвета цвета. Отметьте, как эта уменьшаемая карта только сохраняет самые широкие элементы исходной формы и теряет большую часть детали.
geoshow(masslat3, masslon3, 'DisplayType', 'line', 'color', 'black')
