Верно для точек внутри или на лат-лонном четырёхугольнике
tf = ingeoquad(lat, lon, latlim, lonlim)
tf = ingeoquad(lat, lon, latlim, lonlim) возвращает массив tf который имеет тот же размер, что и lat и lon. tf(k) true тогда и только тогда, когда точка lat(k), lon(k) попадает в пределы или на край географического четырехугольника, определяемого latlim и lonlim. latlim - вектор формы, [southern-limit northern-limit], и lonlim - вектор формы, [western-limit eastern-limit]. Все углы в единицах градусов.
Загрузить данные о высоте и географические ячейки опорного объекта для Корейского полуострова. Отображение данных на карте мира. Применение карты цветов, соответствующей данным отметки, с помощью demcmap.
load korea5c
figure('Color','white')
worldmap([20 50],[90 150])
geoshow(korea5c,korea5cR,'DisplayType','texturemap');
demcmap(korea5c)
Очертить четырехугольник, содержащий данные отметки:
[outlineLat, outlineLon] = outlinegeoquad(korea5cR.LatitudeLimits, ... korea5cR.LongitudeLimits,90,5); geoshow(outlineLat,outlineLon,'DisplayType','line', ... 'Color','k')
Создайте дорожку, пересекающую данные отметки:
[lat,lon] = track2(23,110,48,149,[1 0],'degrees',20); geoshow(lat,lon,'DisplayType','line')
Определите и отметьте точки на дорожке, которые попадают в четырехугольник с указанием данных отметки:
tf = ingeoquad(lat,lon,korea5cR.LatitudeLimits, ...
korea5cR.LongitudeLimits);
geoshow(lat(tf),lon(tf),'DisplayType','point')