Documentation Center

  • 評価版
  • 製品アップデート

satbath

Read 2-minute terrain/bathymetry from Smith and Sandwell

Syntax

[latgrat,longrat,z] = satbath
[latgrat,longrat,z] = satbath(scalefactor)
[latgrat,longrat,z] = satbath(scalefactor,latlim,lonlim)
[latgrat,longrat,z] = satbath(scalefactor,latlim,lonlim,gsize)

Description

[latgrat,longrat,z] = satbath reads the global topography file for the entire world (topo_8.2.img), returning every 50th point. The result is returned as a geolocated data grid. If you use a different version of the global topography file, you need to rename it to "topo_8.2.img". If the file is not found on the MATLAB® path, a dialog opens to request the file.

[latgrat,longrat,z] = satbath(scalefactor) returns the data for the entire world, subsampled by the integer scalefactor. A scalefactor of 10 returns every 10th point. The matrix at full resolution has 6336 by 10800 points.

[latgrat,longrat,z] = satbath(scalefactor,latlim,lonlim) returns data for the specified region. The returned data extends slightly beyond the requested area. If omitted, the entire area covered by the data file is returned. The limits are two-element vectors in units of degrees, with latlim in the range [-90 90] and lonlim in the range [-180 180].

[latgrat,longrat,z] = satbath(scalefactor,latlim,lonlim,gsize) controls the size of the graticule matrices. gsize is a two-element vector containing the number of rows and columns desired. If omitted, a graticule the size of the data grid is returned.

Background

This is a global bathymetric model derived from ship soundings and satellite altimetry by W.H.F. Smith and D.T. Sandwell. The model was developed by iteratively adjusting gravity anomaly data from Geosat and ERS-1 against historical track line soundings. This technique takes advantage of the fact that gravity mirrors the large variations in the ocean floor as small variations in the height of the ocean's surface. The computational procedure uses the ship track line data to calibrate the scaling between the observed surface undulations and the inferred bathymetry. Land elevations are reduced-resolution versions of GTOPO30 data.

Examples

Read the data for the Falklands Islands (Islas Malvinas) at full resolution.

[latgrat,longrat,mat] = satbath(1,[-55 -50],[-65 -55]);
whos

  Name          Size         Bytes  Class

  latgrat     247x301       594776  double array
  longrat     247x301       594776  double array
  mat         247x301       594776  double array

More About

expand all

Tips

Land elevations are given in meters above mean sea level. The data is stored in a Mercator projection grid. As a result, spatial resolution varies with latitude. The grid spacing is 2 minutes (about 4 kilometers) at the equator.

This data is available over the Internet, but subject to copyright. The data file is binary, and should be transferred with no line-ending conversion or byte swapping. This function carries out any byte swapping that might be required. The data requires about 133 MB uncompressed.

The data and documentation are available over the Internet via http and anonymous ftp. Download the latest version of file topo_x.2.img, where x is the version number, and rename it topo_8.w.img for compatibility with the satbath function.

satbath returns a geolocated data grid rather than a regular data grid and a referencing vector or matrix. This is because the data is in a Mercator projection, with columns evenly spaced in longitude, but with decreasing spacing for rows at higher latitudes. Referencing vectors and matrices assume that the number of cells per degrees of latitude and longitude are both constant across a data grid.

    Note   For details on locating map data for download over the Internet, see the following documentation at the MathWorks Web site: Finding Geospatial Data.

See Also

| |

Was this topic helpful?