############
x, y, zb, dists = glaclib.load_flowline(glacier,
shapefilename='flowline_flightline',
filt_len=2.0e3,
verticaldatum='geoid',
bedsource='cresis')
############
# Get DEMs #
############
xdem, ydem, zdem, timedem, errordem = zslib.dem_grid(glacier,
xmin,
xmax,
ymin,
ymax,
years='all',
verticaldatum='geoid',
return_error=True)
####################################
# Get bed elevations near flowline #
####################################
# Get radar thicknesses close to flightline
cresis = bedlib.cresis('all', glacier)
if glacier == 'Helheim':
cresis2001 = bedlib.cresis('2001', glacier)
cresis = np.row_stack([cresis, cresis2001])
cutoff = 200.
ymax = -2545300.
ximage, yimage, image = geotifflib.readrgb(
os.path.join(
os.getenv("DATA_HOME"),
"Imagery/Landsat/Helheim/TIF/20140704140535_LC82330132014185LGN00.tif"
))
dx = 300.
x = np.arange(xmin - 1e3, xmax + 1e3, dx)
y = np.arange(ymin - 1e3, ymax + 1e3, dx)
xdem, ydem, zdem, timedem, errordem = zslib.dem_grid(glacier,
xmin - 2e3,
xmax + 2e3,
ymin - 2e3,
ymax + 2e3,
years='all',
verticaldatum='ellipsoid',
return_error=True,
data='TDM')
xbed, ybed, zbed = bedlib.morlighem_grid(xmin - 2e3,
xmax + 2e3,
ymin - 2e3,
ymax + 2e3,
verticaldatum='ellipsoid')
# Smooth surface DEMs
#zdem_blur = np.zeros_like(zdem)
#for k in range(0,len(timedem)):
# zdem_blur[:,:,k] = scipy.ndimage.filters.gaussian_filter(zdem[:,:,k],sigma=2,truncate=4)
# # "Blurring DEM over 17 pixels (roughly 500 m in each direction)..."
'Helheim',
shapefilename='flowline_flightline',
filt_len=2.0e3,
bedsource='cresis')
x_K, y_K, zb_K, dists_K = glaclib.load_flowline(
'Kanger',
shapefilename='flowline_flightline',
filt_len=2.0e3,
bedsource='cresis')
# Load DEMs
xdem_H, ydem_H, zdem_H, timedem_H, errordem_H = zslib.dem_grid(
'Helheim',
285000.0,
320000.0,
-2588000.0,
-2566000.0,
years='all',
verticaldatum='ellipsoid',
method='nearest',
return_error=True)
xdem_K, ydem_K, zdem_K, timedem_K, errordem_K = zslib.dem_grid(
'Kanger',
449800.0,
503000.0,
-2302000.0,
-2266000.0,
years='all',
verticaldatum='ellipsoid',
method='nearest',
return_error=True)
glacier = args[1][:] # Options: Kanger, Helheim
# Plot extent
if glacier == 'Helheim':
xmin = 304000.0
xmax = 314000.0
ymin = -2582500.0
ymax = -2572500.0
elif glacier == 'Kanger':
xmin = 485000.0
xmax = 498000.0
ymin = -2298000.0
ymax = -2285000.0
# Load worldview DEMs
xwv,ywv,zwv,timewv = zslib.dem_grid(glacier,xmin,xmax,ymin,ymax,years='all',resolution=32,verticaldatum='geoid')
# Calculate height above flotation at CreSIS radar picks
xf,yf,zabovefloat = floatlib.extent(xwv,ywv,zwv,timewv,glacier,rho_i=917.0,rho_sw=1020.0,bedsource='cresis',verticaldatum='geoid')
# Length of plot
N = len(timewv)
years = np.arange(np.floor(np.min(timewv)),np.floor(np.max(timewv))+1)
# Find dimensions for plot
ncol = 0
for year in years:
ndem = len(np.where(np.floor(timewv)==year)[0])
if ndem > ncol:
ncol = ndem
nrow = len(years)
# Flowlines
x_H, y_H, zb_H, dists_H = glaclib.load_flowline(
'Helheim',
shapefilename='flowline_flightline',
filt_len=2.0e3,
bedsource='cresis')
# Flux gates
xgate_H, ygate_H = fluxlib.fluxbox_geometry('Helheim', "fluxgate1")
# Load DEMs
xdem_H, ydem_H, zdem_H, timedem_H, errordem_H = zslib.dem_grid(
'Helheim',
285000.0,
320000.0,
-2588000.0,
-2566000.0,
years='all',
verticaldatum='ellipsoid',
method='nearest',
return_error=True)
# And flotation condition for Helheim
xmin = 283000.
xmax = 313000.
ymin = -2587000.
ymax = -2552000.
xwv_H, ywv_H, zwv_H, timewv_H = zslib.dem_grid('Helheim',
xmin,
xmax,
ymin,
ymax,