| import numpy as np | |
| from matplotlib import colors | |
| import matplotlib.pyplot as plt | |
| def lch_to_lab(LCh): | |
| a = np.cos(LCh[2] * np.pi / 180) * LCh[1] | |
| b = np.sin(LCh[2] * np.pi / 180) * LCh[1] | |
| return np.array((LCh[0], a, b)) |
Steps
gdalbuildvrt g.vrt greenland-modis-median-julyaugust-100m-*.tif
| import numpy as np | |
| import pandas as pd | |
| import geopandas as gpd | |
| from shapely.geometry import Point, Polygon | |
| zwally = np.loadtxt('grndrainagesystems_ekholm.txt',skiprows=7) | |
| regions = np.unique(zwally[:,0]) | |
| for region in regions: | |
| z = zwally[np.abs(zwally[:,0]-region)<0.001,1:] |
| import rioxarray as rio | |
| from rasterio.enums import Resampling | |
| import xarray as xr | |
| import numpy as np | |
| ds = rio.open_rasterio('BedMachineGreenland-v5.nc') | |
| props_to_save = ["thickness", "surface", "bed", "errbed"] | |
| ds.thickness.rio.write_nodata(0, inplace=True) |
| figure.facecolor: white | |
| figure.dpi: 200 | |
| savefig.dpi: 600 | |
| figure.figsize: 3.3,2.2 | |
| axes.facecolor: white | |
| figure.subplot.bottom: 0.01 | |
| figure.subplot.left: 0.01 | |
| figure.subplot.right: 0.99 | |
| figure.subplot.wspace: 0.2 |
| # | |
| # Snippet to empirically determine credible intervals for return period plots | |
| # | |
| # aslak grinsted 2021 | |
| # | |
| import numpy as np | |
| from scipy.stats import beta |
| # -*- coding: utf-8 -*- | |
| """ | |
| Created on Thu Feb 6 15:01:40 2020 | |
| @author: aslak grinsted | |
| """ | |
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| from numpy import nan |
| import numpy as np | |
| import heapq | |
| import scipy.ndimage | |
| # lebrocq et al. | |
| # https://www.sciencedirect.com/science/article/pii/S0098300406000781?via%3Dihub | |
| def balanceflux(Z, P, n=1): | |
| """ |