Преобразуйте определенный геолокацию массив данных в обычную сетку данных
[Z, refvec] = geoloc2grid(lat, lon,
A, cellsize)
[Z, refvec] = geoloc2grid(lat, lon,
A, cellsize) преобразовывает определенный геолокацию массив данных A, учитывая точки геолокации в lat и lon, чтобы произвести обычную сетку данных, Z и соответствующий трехэлементный ссылающийся векторный refvec. cellsize является скаляром, который задает ширину и высоту ячеек данных в обычной сетке данных, с помощью тех же угловых модулей в качестве lat и lon. Ячейки данных в Z, выходящем за пределы области, покрытой A, установлены в NaN.
% Load the geolocated data array 'map1'
% and grid it to 1/2-degree cells.
load mapmtx
cellsize = 0.5;
[Z, refvec] = geoloc2grid(lt1, lg1, map1, cellsize);
% Create a figure
f = figure;
[cmap,clim] = demcmap(map1);
set(f,'Colormap',cmap,'Color','w')
% Define map limits
latlim = [-35 70];
lonlim = [0 100];
% Display 'map1' as a geolocated data array in subplot 1
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');
% Display 'Z' as a regular data grid in subplot 2
subplot(1,2,2)
ax = axesm('mercator','MapLatLimit',latlim,...
'MapLonLimit',lonlim,'Grid','on',...
'MeridianLabel','on','ParallelLabel','on');
set(ax,'Visible','off')
geoshow(Z, refvec, 'DisplayType', 'texturemap');
geoloc2grid обеспечивает, простая в использовании альтернатива gridding определила геолокацию массивов данных с imbedm. Нет никакой потребности предварительно выделить выходную карту; нет никаких разрывов данных в выводе (даже если cellsize выбран, чтобы быть очень маленьким), и выходная карта более сглаженна.