prl900

prl900

prl900

from pylab import *
from mpl_toolkits.mplot3d import Axes3D
from matplotlib.cbook import get_sample_data
from matplotlib._png import read_png
import imageio
from PIL import Image

img1 = imageio.imread("image1.jpg")
img1 = img1/255

prl900

from osgeo import gdal
from osgeo import ogr
from osgeo import gdalconst

#This raster is the model for our output (CRS, extent)
ndsm = 'geotif_output_extent_proj.tif'
#This shapefile contains the features we want to burn
shp = 'shapefile_features_fields.shp'

data = gdal.Open(ndsm, gdalconst.GA_ReadOnly)

prl900

prl900

import os
import requests # get the requsts library from https://github.com/requests/requests
import datetime
import urllib.request
from datetime import timedelta, datetime

def daterange(start_date, end_date):
    n = 0

    while True:

prl900

prl900

prl900