globalrelief

[x, y, z, xyz, xgrid, ygrid, zgrid] = globalrelief(filepath, xmin, xmax, ymin, ymax, dispout)

Description

Return x (longitude), y (latitude) and z (elevation) data from ETOPO1 Global Relief Model (Amante & Eakins, 2009) interpolated on 0.08 degree grid
ETOPO1 is 1 arc-minute global relief, however, this data are interpolated on 4.8 arc-minute (0.08 degree)
This data are obtained from ETOPO1 Global Relief Bedrock (grid-registered)
ETOPO1 horizontal datum: WGS 84 geographic
ETOPO1 vertical datum: sea level
https://www.ngdc.noaa.gov/mgg/global/global.html

Inputs

filepath
Path of the folder that contains ‘ETOPO1_XYZ_0_08_Deg_Matlab.mat’ file
Download ETOPO1_XYZ_0_08_Deg_Matlab.mat file from
https://github.com/akarimp/ScientiMate/releases/download/1.1/ETOPO1_XYZ_0_08_Deg_Matlab.mat
Example: filepath = ‘C:’
xmin=-180
Minimum x (longitude) of domain to be returned in degree
It should be between -180 and 180 degree
xmax=180
Maximum x (longitude) of domain to be returned in degree
It should be between -180 and 180 degree
ymin=-90
Minimum y (latitude) of domain to be returned in degree
It should be between -90 and 90 degree
ymax=90
Maximum y (latitude) of domain to be returned in degree
It should be between -90 and 90 degree
dispout=’no’
Define to display outputs or not (‘yes’: display, ‘no’: not display)

Outputs

x
Interpolated x (longitude) data in degree
y
Interpolated y (latitude) data in degree
z
Interpolated z (elevation) data in degree
xyz
Interpolated xyz data
xyz is a 3-column array
1st column contains longitude (x) data in (Degree)
2nd column contains latitude (y) data in (Degree)
3rd column contains elevation (z) data in (m)
xgrid
Interpolated x (longitude) data on 2d mesh in degree
ygrid
Interpolated y (latitude) data on 2d mesh in degree
zgrid
Interpolated z (elevation) data on 2d mesh in degree

Examples

%Path of ETOPO1_XYZ_0_08_Deg_Matlab.mat
filepath = 'C:'

%Globe
[x,y,z,xyz,xgrid,ygrid,zgrid]=globalrelief(filepath,-180,180,-90,90,'yes');

%Middle East
[x,y,z,xyz,xgrid,ygrid,zgrid]=globalrelief(filepath,24,64,9,43,'yes');

%North America
[x,y,z,xyz,xgrid,ygrid,zgrid]=globalrelief(filepath,-169,-8,5,90,'yes');