True для точек внутри или на четверике лат-лона
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')