GeographicCellsReference

Reference raster cells to geographic coordinates

Description

For a 2-D raster of cells, a GeographicCellsReference object describes the relationship between the intrinsic coordinate system for the raster and a geographic coordinate system. The object associates each element of the raster with a location in geographic coordinates.

A GeographicCellsReference object assumes that the associated raster is sampled regularly in latitude and longitude, and that the rows and columns of the raster align with meridians and parallels, respectively. The sample spacing from row to row does not need to equal the sample spacing from column to column.

This image illustrates a 8-by-10 raster of cells in geographic coordinates. The raster cells have an extent of 1 degree in latitude and longitude, cover longitude values in the range [40,50] in degrees, and cover latitude values in the range [20,28] in degrees. The boundary of the raster of cells is made up of the outermost boundaries of the outermost cells.

To reference a raster of postings to geographic coordinates, use a GeographicPostingsReference object instead.

Creation

You can use these functions to create GeographicCellsReference objects:

• readgeoraster — Import a raster data file as an array and a reference object.

• wmsread — Import a raster map from a Web Map Service (WMS) server as an array and a reference object.

• georefcells — Create a geographic raster reference object by specifying the latitude and longitude limits of the raster and either the size of the raster or the extent of the cells.

• georasterref — Convert a world file to a geographic raster reference object.

• refmatToGeoRasterReference — Convert a referencing matrix to a geographic raster reference object.

• refvecToGeoRasterReference — Convert a referencing vector to a geographic raster reference object.

Properties

LatitudeLimits — Latitude limits of the geographic quadrangle bounding the georeferenced raster

[0.5 2.5] (default) | two-element vector

Latitude limits of the geographic quadrangle bounding the georeferenced raster, specified as a two-element vector of the form [southern_limit northern_limit].

Example: [-90 90]

Data Types: double

LongitudeLimits — Longitude limits of the geographic quadrangle bounding the georeferenced raster

[0.5 2.5] (default) | two-element vector

Longitude limits of the geographic quadrangle bounding the georeferenced raster, specified as a two-element vector of the form [western_limit eastern_limit].

Example: [-100 180]

Data Types: double

RasterSize — Number of rows and columns of the raster or image associated with the referencing object

[2 2] (default) | two-element vector of positive integers

Number of rows and columns of the raster or image associated with the referencing object, specified as a two-element vector, [m n], where m represents the number of rows and n the number of columns. For convenience, you can assign a size vector having more than two elements. This enables assignments like R.RasterSize = size(RGB), where RGB is m-by-n-by-3. In cases like this, the object stores only the first two elements of the size vector and ignores the higher (nonspatial) dimensions.

Example: [200 300]

Data Types: double

RasterInterpretation — Geometric nature of the raster

Read-only: 'cells' (default)

This property is read-only.

Geometric nature of the raster, specified as 'cells'. The value 'cells' indicates that the raster comprises a grid of quadrangular cells, and is bounded on all sides by cell edges. For an m-by-n raster, points with an intrinsic x-coordinate of 1 or n or an intrinsic y-coordinate of 1 or m fall within the raster, not on its edges.

Data Types: char

AngleUnit — Unit of measurement used for angle-valued properties

'degree' (default)

Unit of measurement used for angle-valued properties, specified as 'degree'.

Cannot be set.

Data Types: char

ColumnsStartFrom — Edge from which column indexing starts

'south' (default) | 'north'

Edge from which column indexing starts, specified as 'south' or 'north'.

Example: 'south'

Data Types: char

RowsStartFrom — Edge from which row indexing starts

'west' (default) | 'east'

Edge from which row indexing starts, specified as 'west' or 'east'.

Example: 'east'

Data Types: char

CellExtentInLatitude — Extent in latitude of individual cells

1 (default) | positive numeric scalar

Extent in latitude of individual cells, specified as a positive numeric scalar. Distance, in units of latitude, between the northern and southern limits of a single raster cell. The value is the same for all cells in the raster.

Example: 2.5

Data Types: double

CellExtentInLongitude — Extent in longitude of individual cells

1 (default) | positive numeric scalar

Extent in longitude of individual cells, specified as a positive numeric scalar. Distance, in units of longitude, between the western and eastern limits of a single raster cell. The value is always positive, and is the same for all cells in the raster.

Example: 2.5

Data Types: double

RasterExtentInLatitude — Latitude extent ("height") of the quadrangle covered by the raster

Read-only: 2 (default) | positive numeric scalar

This property is read-only.

Latitude extent ("height") of the quadrangle covered by the raster, specified as a positive numeric scalar.

Example: 2

Data Types: double

RasterExtentInLongitude — Longitude extent ("width") of the quadrangle covered by the raster

Read-only: 2 (default) | positive numeric scalar

This property is read-only.

Longitude extent ("width") of the quadrangle covered by the raster, specified as a positive numeric scalar.

Data Types: double

XIntrinsicLimits — Raster limits in intrinsic x coordinates

Read-only: [0.5 2.5] (default) | two-element row vector of positive integers

This property is read-only.

Raster limits in intrinsic x coordinates, specified as a two-element row vector of positive integers, [xMin xMax]. For an m-by-n raster, XIntrinsicLimits equals [0.5, m+0.5], because the RasterInterpretation is 'cells'.

Example: [0.5 2.5]

Data Types: double

YIntrinsicLimits — Raster limits in intrinsic y coordinates

Read-only: [0.5 2.5] (default) | two-element row vector of positive integers

This property is read-only.

Raster limits in intrinsic y coordinates, specified as a two-element row vector of positive integers, [yMin yMax]. For an m-by-n raster, YIntrinsicLimits equals [0.5, n+0.5], because the RasterInterpretation is 'cells'.

Data Types: double

CoordinateSystemType — Type of coordinate system to which the image or raster is referenced

Read-only: 'geographic' (default)

This property is read-only.

Type of coordinate system to which the image or raster is referenced, specified as 'geographic'.

Data Types: char

GeographicCRS — Geographic coordinate reference system

[] (default) | geocrs object

Geographic coordinate reference system (CRS), specified as a geocrs object. A geographic CRS consists of a datum (including its ellipsoid), prime meridian, and angular unit of measurement.

Object Functions

contains	Determine if geographic or map raster contains points
cropToBlock	Crop raster reference object to row and column limits of raster
geographicGrid	Geographic coordinates of raster elements
geographicToDiscrete	Convert geographic coordinates to discrete subscripts
geographicToIntrinsic	Transform geographic to intrinsic coordinates
intrinsicToGeographic	Transform intrinsic to geographic coordinates
intrinsicXToLongitude	Convert from intrinsic x to longitude coordinates
intrinsicYToLatitude	Convert from intrinsic y to latitude coordinates
latitudeToIntrinsicY	Convert from latitude to intrinsic y coordinates
longitudeToIntrinsicX	Convert from longitude to intrinsic x coordinates
sizesMatch	Determine if geographic or map raster object and image or raster are size-compatible
worldFileMatrix	Return world file parameters for transformation

Examples

Display Geographic Raster with Regularly Spaced Cell Centers

Import a GeoTIFF image as an array and a GeographicCellsReference object by using the readgeoraster function.

[Z,R] = readgeoraster('katrina.tif');

View the properties of the GeographicCellsReference object.

R

R = 
  GeographicCellsReference with properties:

             LatitudeLimits: [10 40]
            LongitudeLimits: [-100 -70]
                 RasterSize: [1024 1024]
       RasterInterpretation: 'cells'
           ColumnsStartFrom: 'north'
              RowsStartFrom: 'west'
       CellExtentInLatitude: 0.029296875
      CellExtentInLongitude: 0.029296875
     RasterExtentInLatitude: 30
    RasterExtentInLongitude: 30
           XIntrinsicLimits: [0.5 1024.5]
           YIntrinsicLimits: [0.5 1024.5]
       CoordinateSystemType: 'geographic'
              GeographicCRS: [1×1 geocrs]
                  AngleUnit: 'degree'

Create a world map with limits that match the limits of the reference object. Then, display the image.

worldmap(Z,R)
geoshow(Z,R)

More About

Intrinsic Coordinate System

An intrinsic coordinate system is a 2-D Cartesian coordinate system where:

• The x-axis is parallel to the rows of the raster. From column to column, the value of x increases by 1. At the center of a cell, the value of x has an integer value that is equal to the column index.

• The y-axis is parallel to the columns of the raster. From row to row, the value of y increases by 1. At the center of a cell, the value of y has an integer value equal to the row index.

Coordinates in an intrinsic coordinate system are related to the indices of the raster. For example, when column indexing starts from the north edge of the raster and row indexing starts from the west edge of the raster, the intrinsic coordinates x = 3.0, y = 5.0 correspond to the center point of the cell in row 5 and column 3.

Version History

Introduced in R2013b