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