Обрезка многоугольников до четырехугольника широты-долготы
[latTrimmed,lonTrimmed] = maptrimp(lat,lon,latlim,lonlim)
[latTrimmed,lonTrimmed] = maptrimp(lat,lon,latlim,lonlim) обрезка полигонов в lat и lon в четырехугольник, указанный latlim и lonlim. latlim и lonlim - двухэлементные векторы, определяющие пределы широты и долготы соответственно. lat и lon должны быть векторами, представляющими допустимые многоугольники.
Считывайте данные береговой линии и отображайте их на карте.
load coastlines figure mapshow(coastlon,coastlat,'DisplayType','polygon');
Обрезать набор данных.
latlim = [-50 50];
lonlim = [-100 50];
[latTrimmed,lonTrimmed] = maptrimp(coastlat,coastlon, ...
latlim, lonlim);Отображение обрезанного набора данных.
figure mapshow(lonTrimmed,latTrimmed, 'DisplayType', 'polygon');
maptrimp определяет пределы долготы таким образом, что:
lonlim(2) всегда превышает lonlim(1)
lonlim(2) никогда не превышает lonlim(1) более чем на 360
lonlim(1) < 180 или lonlim(2) > -180
Убедитесь, что если четырехугольник охватывает Гринвичский меридиан, то этот меридиан появляется на долготе 0.