Преобразуйте определенный геолокацию массив данных в обычную сетку данных
[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
выбран, чтобы очень быть малыми), и выходная карта более сглаженна.