scientimate.globalrelief#

x, y, z, xyz, xgrid, ygrid, zgrid = scientimate.globalrelief(filepath, xmin=-180, xmax=180, ymin=-90, ymax=90, dispout='no')

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_Python.npz’ file
Download ETOPO1_XYZ_0_08_Deg_Python.npz file from
akarimp/ScientiMate
Example: filepath = r’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#

import scientimate as sm

#Path of ETOPO1_XYZ_0_08_Deg_Python.npz
filepath = r'C:'

#Globe
x,y,z,xyz,xgrid,ygrid,zgrid=sm.globalrelief(filepath,-180,180,-90,90,'yes')

#Middle East
x,y,z,xyz,xgrid,ygrid,zgrid=sm.globalrelief(filepath,24,64,9,43,'yes')

#North America
x,y,z,xyz,xgrid,ygrid,zgrid=sm.globalrelief(filepath,-169,-8,5,90,'yes')

References#

ETOPO1 Global Relief Model

Amante, C. and B.W. Eakins, 2009. ETOPO1 1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis. NOAA Technical Memorandum NESDIS NGDC-24. National Geophysical Data Center, NOAA. doi:10.7289/V5C8276M

GEBCO Global ocean & land terrain models

Natural Earth 1:10m Raster Data

Geospatial data