# vec2mtx

Convert latitude-longitude vectors to regular data grid

## Syntax

`[Z,R] = vec2mtx(lat,lon,density)`

[Z,R] = vec2mtx(lat,lon,density,latlim,lonlim)

[Z,R] = vec2mtx(lat,lon,Z1,R1)

[Z,R] = vec2mtx(...,'filled')

## Description

`[Z,R] = vec2mtx(lat,lon,density)`

creates a regular data grid
`Z`

from vector data, placing ones in grid cells intersected by a
vector and zeroes elsewhere. `R`

is the raster reference object for the
computed grid. `lat`

and `lon`

are vectors of equal
length containing geographic locations in units of degrees. `density`

indicates the number of grid cells per unit of latitude and longitude (a value of 10
indicates 10 cells per degree, for example), and must be scalar-valued. Whenever there
is space, a buffer of two grid cells is included on each of the four sides of the grid.
The buffer is reduced as needed to keep the latitudinal limits within [-90 90] and to
keep the difference in longitude limits from exceeding 360 degrees.

`[Z,R] = vec2mtx(lat,lon,density,latlim,lonlim)`

uses
the two-element vectors `latlim`

and `lonlim`

to
define the latitude and longitude limits of the grid.

`[Z,R] = vec2mtx(lat,lon,Z1,R1)`

uses a pre-existing
data grid `Z1`

, georeferenced by `R1`

, to define the
limits and density of the output grid. Specify `R1`

as a `GeographicCellsReference`

object. The `RasterSize`

property of `R1`

must be consistent with `size(Z)`

.
With this syntax, the output `R`

is equal to
`R1`

.

`[Z,R] = vec2mtx(...,'filled')`

, where
`lat`

and `lon`

form one or more closed polygons
(with `NaN`

-separators), fills the area outside the polygons with the
value two instead of the value zero.

## Notes

Empty `lat,lon`

vertex arrays will result in
an error unless the grid limits are explicitly provided (via `latlim,lonlim`

or `Z1,R1`

).
In the case of explicit limits, `Z`

will be filled
entirely with 0s if the `'filled'`

parameter is omitted,
and 2s if it is included.

It's possible to apply `vec2mtx`

to sets of
polygons that tile without overlap to cover an area, as in Example
1 below, but using `'filled'`

with polygons that
actually overlap may lead to confusion as to which areas are inside
and which are outside.

## Examples

## Version History

**Introduced before R2006a**