verbose = True print(sys.argv, year, country_code) ### get land cover data fname = f'data/landcover_{year}_{country_code}.npz' ofile = Path('data') / fname done = False if ofile.exists(): done = True # try to download it from server done = procure_dataset(fname, verbose=True) from geog0111.get_modis_files import get_modis_files from geog0111.process_timeseries import mosaic_and_clip ''' Get the MODIS LC files from the server to store in data ''' if not done: try: url = 'https://e4ftl01.cr.usgs.gov//MODV6_Cmp_C/MOTA/' filename = get_modis_files(1,year,[tiles],base_url=url,\ version=6,verbose=True,\ destination_folder='data',\ product='MCD12Q1')[0] print(filename)
x = np.arange(-3*sigma,3*sigma+1) gaussian = np.exp((-(x/sigma)**2)/2.0) FIPS = country_code dates, lai_array, weights_array = lai['dates'],lai['lai'],lai['weights'] print(lai_array.shape, weights_array.shape) #Check the output array shapes from geog0111.geog_data import procure_dataset import numpy as np from pathlib import Path ''' LAI data ''' # read in the LAI data for given country code lai_filename = f'data/lai_data_{year}_{country_code}.npz' # get the dataset in case its not here procure_dataset(Path(lai_filename).name,verbose=False) lai = np.load(lai_filename) print(lai_filename,list(lai.keys())) ''' T 2m data ''' t2_filename = f'data/europe_data_{year}_{country_code}.npz' # get the dataset in case its not here procure_dataset(Path(t2_filename).name,verbose=False) t2data = np.load(t2_filename) print(t2_filename,list(t2data.keys()))
import sys year = 2017 if len(sys.argv) == 3: country_code = sys.argv[1] year = int(sys.argv[2]) elif len(sys.argv) == 2: country_code = sys.argv[1] else: country_code = 'UK' verbose = True print(sys.argv,year,country_code) files = list(Path('data').glob('MCD15A3H.A201[6-7]*h1[7-8]v0[3-4].006*hdf')) if len(files) < 732: _ = procure_dataset("lai_files",verbose=False) # Get the shapefile for country borders import requests import shutil from pathlib import Path force = False # zip file zipfile = 'TM_WORLD_BORDERS-0.3.zip' # URL tm_borders_url = f"http://thematicmapping.org/downloads/{zipfile}" # destibnation folder destination_folder = Path('data') # set up some filenames