Верный для точек внутри или на 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]. Все углы находятся в модулях степеней.
Загрузите и отобразите цифровую модель вертикального изменения (DEM) включая корейский Полуостров:
korea = load('korea');
R = refvecToGeoRasterReference(korea.refvec, size(korea.map));
figure('Color','white')
worldmap([20 50],[90 150])
geoshow(korea.map, R, 'DisplayType', 'texturemap');
demcmap(korea.map)
Обрисуйте в общих чертах четырехугольник, содержащий DEM:
[outlineLat, outlineLon] = outlinegeoquad(R.Latlim, ... R.Lonlim, 90, 5); geoshow(outlineLat,outlineLon,'DisplayType','line', ... 'Color','black')
Сгенерируйте дорожку, которая пересекает DEM:
[lat, lon] = track2(23, 110, 48, 149, [1 0], 'degrees', 20); geoshow(lat, lon, 'DisplayType', 'line')
Идентифицируйте и отметьте точки на дорожке, которые находятся в пределах четырехугольника, обрисовывающего в общих чертах DEM:
tf = ingeoquad(lat, lon, R.Latlim, R.Lonlim); geoshow(lat(tf), lon(tf), 'DisplayType', 'point')
