maprefcells

R2026b

Reference raster cells to map coordinates

Syntax

R = maprefcells(xlimits,ylimits,rasterSize)

R = maprefcells(xlimits,ylimits,xcellextent,ycellextent)

R = maprefcells(___,Name=Value)

Description

R = maprefcells(xlimits,ylimits,rasterSize) creates a reference object for a regular raster of cells in planar coordinates. Specify the xy-limits of the raster using xlimits and ylimits. Specify the number of rows and columns in the raster using rasterSize.

R = maprefcells(xlimits,ylimits,xcellextent,ycellextent) specifies the width and height of each cell using xcellextent and ycellextent. If necessary, the maprefcells function adjusts the xy-limits to ensure that each dimension of the raster has an integer number of cells.

R = maprefcells(___,Name=Value) specifies properties of the reference object using one or more name-value arguments, in addition to the input arguments in the previous syntaxes.

Examples

Create Reference Object for Raster of Cells

Create a reference object for a raster of cells, where the first cell is in the northwest corner. Specify the xy-limits using world coordinates. Specify the projected CRS as NAD83 / Massachusetts Mainland, which has the EPSG authority code 26986.

xlimits = [207000 208000];
ylimits = [912500 913000];
rasterSize = [1000 2000];
pcrs = projcrs(26986);

R1 = maprefcells(xlimits,ylimits,rasterSize, ...
    ColumnsStartFrom="north",ProjectedCRS=pcrs)

R1 = 
  MapCellsReference with properties:

            XWorldLimits: [207000 208000]
            YWorldLimits: [912500 913000]
              RasterSize: [1000 2000]
    RasterInterpretation: 'cells'
        ColumnsStartFrom: 'north'
           RowsStartFrom: 'west'
      CellExtentInWorldX: 1/2
      CellExtentInWorldY: 1/2
    RasterExtentInWorldX: 1000
    RasterExtentInWorldY: 500
        XIntrinsicLimits: [0.5 2000.5]
        YIntrinsicLimits: [0.5 1000.5]
      TransformationType: 'rectilinear'
    CoordinateSystemType: 'planar'
            ProjectedCRS: [1×1 projcrs]

Alternatively, specify the cell extents instead of the raster size.

extent = 1/2;
R2 = maprefcells(xlimits,ylimits,extent,extent, ...
    ColumnsStartFrom="north",ProjectedCRS=pcrs)

R2 = 
  MapCellsReference with properties:

            XWorldLimits: [207000 208000]
            YWorldLimits: [912500 913000]
              RasterSize: [1000 2000]
    RasterInterpretation: 'cells'
        ColumnsStartFrom: 'north'
           RowsStartFrom: 'west'
      CellExtentInWorldX: 1/2
      CellExtentInWorldY: 1/2
    RasterExtentInWorldX: 1000
    RasterExtentInWorldY: 500
        XIntrinsicLimits: [0.5 2000.5]
        YIntrinsicLimits: [0.5 1000.5]
      TransformationType: 'rectilinear'
    CoordinateSystemType: 'planar'
            ProjectedCRS: [1×1 projcrs]

Input Arguments

xlimits — x-limits, in world coordinates

two-element vector

x-limits of the rectangle that bounds the raster, in world coordinates, specified as a two-element vector of the form [xmin xmax], where xmax is greater than xmin.

ylimits — y-limits, in world coordinates

two-element vector

y-limits of the rectangle that bounds the raster, in world coordinates, specified as a two-element vector of the form [ymin ymax], where ymax is greater than ymin.

rasterSize — Number of rows and columns

two-element vector

Number of rows and columns in the raster, specified as a two-element vector of the form [m n], where m is the number of rows and n is the number of columns. This argument sets the RasterSize property of R.

xcellextent — Width of each cell, in world coordinates

positive scalar

Width of each cell, in world coordinates, specified as a positive scalar. This argument sets the CellExtentInWorldX property of R.

ycellextent — Height of each cell, in world coordinates

positive scalar

Height of each cell, in world coordinates, specified as a positive scalar. This argument sets the CellExtentInWorldY property of R.

Name-Value Arguments

Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

Example: R = maprefcells(xlimits,ylimits,rasterSize,ColumnsStartFrom="north") creates a map cells reference object with columns that start from north.

ColumnsStartFrom — Edge from which column indexing starts

"south" (default) | "north"

Edge from which column indexing starts, specified as one of these options:

• "south" — Column indexing starts from south.

• "north" — Column indexing starts from north.

RowsStartFrom — Edge from which row indexing starts

"west" (default) | "east"

Edge from which row indexing starts, specified as one of these options:

• "west" — Row indexing starts from west.

• "east" — Row indexing starts from east.

ProjectedCRS — Projected coordinate reference system

[] (default) | projcrs object

Projected coordinate reference system (CRS), specified as a projcrs object.

Output Arguments

R — Map cells reference object

MapCellsReference object

Map cells reference object, returned as a MapCellsReference object.

Tips

To create a map raster reference object from a world file matrix, use the maprasterref function.

Version History

Introduced in R2015b

R2026b: Specify projected CRS using name-value argument

Specify the projected CRS of the raster reference object during creation by using the ProjectedCRS name-value argument.