Преобразовать геоориентированный массив данных в обычную сетку данных
[Z,R] = geoloc2grid(lat,lon,A,cellsize)
[Z,R] = geoloc2grid(lat,lon,A,cellsize) преобразует геоориентированный массив данных A, заданные точки геолокации в lat и lon, для создания регулярной сетки данных, Zи соответствующий объект растровой ссылки R. cellsize скаляр, задающий ширину и высоту ячеек данных в регулярной сетке данных, используя те же угловые единицы, что и lat и lon. Ячейки данных в Z выпадение за пределы зоны, охватываемой A имеют значение NaN.
Загрузить геоориентированный массив данных map1 и распределить его по ½-градусным ячейкам.
load mapmtx
cellsize = 0.5;
[Z,R] = geoloc2grid(lt1,lg1,map1,cellsize);Создайте фигуру.
f = figure; [cmap,clim] = demcmap(map1); set(f,'Colormap',cmap,'Color','w')
Определите границы карты.
latlim = [-35 70]; lonlim = [0 100];
Показ map1 как геоориентированный массив данных и Z в качестве обычной сетки данных с использованием вложенных графиков.
subplot(1,2,1) ax = axesm('mercator','MapLatLimit',latlim,... 'MapLonLimit',lonlim,'Grid','on',... 'MeridianLabel','on','ParallelLabel','on'); set(ax,'Visible','off') geoshow(lt1,lg1,map1,'DisplayType','texturemap'); subplot(1,2,2) ax = axesm('mercator','MapLatLimit',latlim,... 'MapLonLimit',lonlim,'Grid','on',... 'MeridianLabel','on','ParallelLabel','on'); set(ax,'Visible','off') geoshow(Z,R,'DisplayType','texturemap');
geoloc2grid обеспечивает простую в использовании альтернативу сетке геоориентированных массивов данных с imbedm. Нет необходимости предварительно распределять карту вывода; в выходных данных отсутствуют пробелы в данных (даже если cellsize выбрано очень маленьким), а выходная карта более плавная.