Encode data points into regular data grid
Z = imbedm(lat, lon, value, Z, R)
Z = imbedm(lat, lon, value, Z, R, units)
[Z, indxPointOutsideGrid] = imbedm(...)
Z = imbedm(lat, lon, value, Z, R) resets
certain entries of a regular data grid,
be a geographic raster reference object, a referencing vector, or
a referencing matrix.
R is a geographic raster reference object,
RasterSize property must be consistent with
R is a referencing vector, it must be
a 1-by-3 with elements:
[cells/degree northern_latitude_limit western_longitude_limit]
R is a referencing matrix, it must be
3-by-2 and transform raster row and column indices to or from geographic
coordinates according to:
[lon lat] = [row col 1] * R
R is a referencing matrix, it must define
a (non-rotational, non-skewed) relationship in which each column of
the data grid falls along a meridian and each row falls along a parallel.
Nearest-neighbor interpolation is used by default. NaN is returned
for points outside the grid limits or for which
NaN. All angles are in units of degrees.
Z = imbedm(lat, lon, value, Z, R, units) specifies the units of the
units is any valid angle units character vector
'degrees' by default).
[Z, indxPointOutsideGrid] = imbedm(...) returns
the indices of
to points outside the grid in the variable
Create a simple grid map and embed new values in it:
Z = ones(3,6) Z = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 refvec = [1/60 90 -180] refvec = 0.0167 90.0000 -180.0000 newgrid = imbedm([23 -23], [45 -45],[5 5],Z,refvec) newgrid = 1 1 1 1 1 1 1 1 5 5 1 1 1 1 1 1 1 1