ingeoquad
Верный для точек внутри или на lat-lon четырехугольнике
Синтаксис
tf = ingeoquad(lat, lon, latlim, lonlim)
Описание
tf = ingeoquad(lat, lon, latlim, lonlim) возвращает массив tf это имеет тот же размер как lat и lon. tf(k) верно если и только если точка lat(k), lon(k) находится в пределах или на ребре географического четырехугольника, заданного latlim и lonlimlatlim вектор из формы [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')