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

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
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')