Extract data grid values for specified locations

`val = ltln2val(Z, R, lat, lon)`

val = ltln2val(Z, R, lat, lon, method)

`val = ltln2val(Z, R, lat, lon)`

interpolates
a regular data grid `Z`

with referencing vector `R`

at
the points specified by vectors of latitude and longitude, `lat`

and `lon`

. `R`

can
be a geographic raster reference object, a referencing vector, or
a referencing matrix.

If `R`

is a geographic raster reference object,
its `RasterSize`

property must be consistent with `size(Z)`

.

If `R`

is a referencing vector, it must be
a 1-by-3 with elements:

[cells/degree northern_latitude_limit western_longitude_limit]

If `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

If `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 `lat`

or `lon`

contain
NaN. All angles are in units of degrees.

`val = ltln2val(Z, R, lat, lon, method)`

accepts
a `method`

string to specify the type of interpolation: `'bilinear'`

for
linear interpolation, `'bicubic'`

for cubic interpolation,
or `'nearest'`

for nearest neighbor interpolation.

Find the elevations in `topo`

associated with
three European cities—Milan, Bern, and Prague (`topo`

elevations
are in meters):

load topo % The city locations, [Milan Bern Prague] lats = [45.45; 46.95; 50.1]; longs = [9.2; 7.4; 14.45]; elevations = ltln2val(topo,topolegend,lats,longs) elevations = 313 1660 297

Was this topic helpful?