exponenta event banner

геоурожай

Географический растр подрезки

Описание

пример

[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 функция. Затем отобразите сдвинутые данные как поверхность на карте.

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

load topo60c
latlim = topo60cR.LatitudeLimits;
[B,RB] = geocrop(topo60c,topo60cR,latlim,[-180 180]);

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

topo60cR.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-by-N или M-by-N-by-P числовой или логический массив.

Ссылка на растр для A, указано как GeographicCellsReference объект или GeographicPostingsReference объект.

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

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

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

свернуть все

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

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

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

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

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

См. также

|

Представлен в R2020a