Извлечение и соединение полигонов или отрезков

В этом примере показано, как идентифицировать сегменты линий или сегментов сегментов после их объединения в большие сегменты. NaN- отсеченные векторы. Эти векторы многоугольников или линий можно разделить на сегменты компонентов с помощью команды polysplit функция, которая принимает векторы столбцов в качестве входных данных. Чтобы объединить отдельные векторы многоугольников или линий, используйте polyjoin.

Создать два NaN-различенные массивы в виде векторов столбцов.

lat = [45.6 -23.47 78 NaN 43.9 -67.14 90 -89]';
lon = [13 -97.45 165 NaN 0 -114.2 -18 0]';

Разбейте векторы столбцов на отдельные массивы ячеек линейных сегментов в NaN сепараторы с использованием polysplit.

[latc,lonc] = polysplit(lat,lon)

latc=2×1 cell array
    {3x1 double}
    {4x1 double}

lonc=2×1 cell array
    {3x1 double}
    {4x1 double}

Проверьте содержимое массивов ячеек. Обратите внимание, что каждый элемент массива ячеек содержит сегмент исходной строки.

[latc{1} lonc{1}]

ans = 3×2

   45.6000   13.0000
  -23.4700  -97.4500
   78.0000  165.0000

[latc{2} lonc{2}]

ans = 4×2

   43.9000         0
  -67.1400 -114.2000
   90.0000  -18.0000
  -89.0000         0

Для сторнирования процесса используйте polyjoin.

[lat2,lon2] = polyjoin(latc,lonc);

Выполните логическое сравнение соединенных сегментов. Обратите внимание, что они идентичны начальным lat и lon массивы. Логическое сравнение имеет значение false для NaN разделители, по определению.

[lat lon] == [lat2 lon2]

ans = 8x2 logical array

   1   1
   1   1
   1   1
   0   0
   1   1
   1   1
   1   1
   1   1

Проверка глобального равенства, в том числе NaN значения.

isequaln(lat,lat2) & isequaln(lon,lon2)

ans = logical
   1

См. также

polyjoin | polysplit

Связанные темы

• Создание и отображение полигонов