############ 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,