geoquadpt

Географический четырехугольник, ограничивающий рассеянные точки

Описание

пример

[latlim,lonlim] = geoquadpt(lat,lon) возвращает пределы максимально плотного географического четверика, который ограничивает набор точек географическими координатами lat и lon.

В большинстве случаев tf = ingeoquad(lat, lon, latlim, lonlim) вернёт true, но tf может быть ложным для точек на ребрах четырехугольника, из-за округления tf. будет также ложным для элементов lat которые попадают за пределы интервала [-90 90] и элементов lon которые не являются конечными.

Примеры

свернуть все

В этом случае выход четырёхугольник преодолевает 180-градусный меридиан, отсюда и элементы lonlim находятся в нисходящем численном порядке, хотя и упорядочены с запада на восток.

Считайте набор точек, а затем создайте ограничительный прямоугольник вокруг точек.

cities = shaperead('worldcities.shp','Selector', ...
    {@(name) any(strcmp(name,{'Tokyo','Honolulu'})),'Name'}, ...
    'UseGeoCoords',true);
[latlim,lonlim] = geoquadpt([cities.Lat],[cities.Lon])
latlim =

   21.3178   35.7082


lonlim =

  139.6401 -157.8291

Входные параметры

свернуть все

Широты точек, заданные в виде вектора, матрицы или N-D массива, в единицах степеней.

Типы данных: single | double

Долготы точек, заданные в виде вектора, матрицы или N-D массива, в единицах степеней.

Типы данных: single | double

Выходные аргументы

свернуть все

Пределы широты географического четверика, возвращенные как вектор 1 на 2 вида [southern_limit northern_limit], в единицах степеней. Элементы находятся в порядке возрастания, и оба находятся в закрытом интервале [-90 90].

Пределы долготы географического четверика, возвращенные как вектор 1 на 2 вида [western_limit eastern_limit], в единицах степеней. Пределы обернуты к интервалу [-180 180]. Они не обязательно находятся в численном порядке возрастания.

Введенный в R2012b