geocrop

Обрежьте географический растр

Описание

пример

[B,RB] = geocrop(A,RA,latlim,lonlim) обрезки растр заданы A и растровая ссылка RA и возвращает обрезанный растр B и растровая ссылка RB. Возвращенный растр обрезается к географическим пределам в градусах близко к заданным latlim и lonlim.

Примеры

свернуть все

Обрежьте географический растр и отобразите обрезанный растр на карте.

Во-первых, загрузите растровые данные о вертикальном изменении и географический объект ссылки ячеек для корейского полуострова. Затем обрежьте растр к пределам, заданным latlim и lonlim.

load korea5c
latlim = [34.25 38.72]; 
lonlim = [125.85 129.92];
[B,RB] = geocrop(korea5c,korea5cR,latlim,lonlim);

Отобразите обрезанный растр как поверхность на карте. Примените палитру, подходящую для данных о вертикальном изменении с помощью demcmap функция.

worldmap(latlim,lonlim)
geoshow(B,RB,'DisplayType','surface')
demcmap(B)

Переключите пределы долготы растра с пределами, которые охватывают 360 градусов с помощью geocrop функция. Затем отобразите сдвинутые данные как поверхность на карте.

Во-первых, загрузите топографические растровые данные для мира и создайте связанный географический объект ссылки ячеек.

load('topo','topo','topolatlim','topolonlim')
topoR = georefcells(topolatlim,topolonlim,size(topo));

Затем переключите пределы долготы растра от [0, 360] к [-180, 180].

[B,RB] = geocrop(topo,topoR,topolatlim,[-180 180]);

Сравните растры путем запроса их LongitudeLimits свойства.

topoR.LongitudeLimits
ans = 1×2

     0   360

RB.LongitudeLimits
ans = 1×2

  -180   180

Отобразите сдвинутые данные как поверхность на карте. Переместите метки меридиана в нижнюю часть карты с помощью mlabel функция. Задайте метки меридиана, чтобы отобразить использование MLabelLocation свойство. Затем примените палитру, подходящую для топографических данных с помощью demcmap функция. Обратите внимание на то, что переключенные пределы долготы появляются в ребрах карты.

worldmap(RB.LatitudeLimits,RB.LongitudeLimits)
geoshow(B,RB,'DisplayType','surface')
mlabel('south')
setm(gca,'MLabelLocation',-180:90:180)
demcmap(B)

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

свернуть все

Географический растр в виде M на n или M N P числовым или логическим массивом.

Растровая ссылка для AВ виде GeographicCellsReference объект или GeographicPostingsReference объект.

Широта ограничивает в виде двухэлементного числового вектора формы [nlat slat], где nlat самый северный предел в градусах и slat самый южный предел в градусах.

Долгота ограничивает в виде двухэлементного числового вектора формы [wlon elon], где wlon самый западный предел в градусах и elon самый восточный предел в градусах.

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

свернуть все

Обрезанный географический растр, возвращенный как числовой или логический массив. Тип данных и размер B совпадает с типом данных и размером A.

Если пределы заданы latlim и lonlim не пересекайте растр, заданный A и RA, затем B isempty.

Растровая ссылка для B, возвращенный как GeographicCellsReference объект или GeographicPostingsReference объект. Тип объекта RB совпадает с типом объекта RA.

Точные пределы широты и долготы RB не совпадайте с пределами, заданными latlim и lonlim, если они не совпадают с местоположением регистрации или границей ячейки. В противном случае, пределы RB немного больше, чем latlim и lonlim.

Если пределы заданы latlim и lonlim не пересекайте растр, заданный A и RA, затем RB isempty.

Смотрите также

|

Введенный в R2020a

Для просмотра документации необходимо авторизоваться на сайте