Aslak Grinsted grinsted

Steps

use earth engine to export a geotiff of the median modis in higher resolution than you need. https://code.earthengine.google.com/5ce94569601c58a14bcaf28c1852e006
download them
build a vrt file

gdalbuildvrt g.vrt greenland-modis-median-julyaugust-100m-*.tif


	import requests
	_host= 'https://eu5.fusionsolar.huawei.com'
	_username = 'test'
	_password = 'test'

	url = _host + '/thirdData/login'
	headers = {
	'accept': 'application/json',
	}

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

	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)

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