Обрезка многоугольников по четырехугольнику широта-долгота
[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.