def cut(lat,lon):
	import numpy as np
	root, is_new = nc.open(localpath+"/dem.nc")
	dem = nc.getvar(root, 'dem', 'i2', lat.dimensions)
	x, y = p.pixels_from_coordinates(lat[:], lon[:], dem.shape[0], dem.shape[1])
	# In the transformation the 'y' dimension is twisted because the map is inverted.
	result = p.transform_data(dem,x,dem.shape[0] - y)
	nc.close(root)
	return result, x, y
def project_coordinates(lat, lon):
	shape = get_month(1).shape
	return p.pixels_from_coordinates(lat, lon, shape[0], shape[1])
def project_coordinates(lat, lon):
    ds = gdal.Open(path + "/tifs/01_longlat_wgs84.tif")
    return p.pixels_from_coordinates(lat, lon, ds.RasterYSize, ds.RasterXSize)