"Given a latitude and a distance west, return the change in longitude."

# Find the radius of a circle around the earth at given latitude.

r = earth_radius*math.cos(latitude*degrees_to_radians)

slat, nlat = latitude+change_in_latitude(-3.75), latitude+change_in_latitude(3.75)

wlon = longitude+change_in_longitude(latitude,-3.75)

elon = longitude+change_in_longitude(latitude, 3.75)

'''First argument degrees longitude (E is positive, W negative)

of the landslide location,

second argument latitude (N positive, S negative),

in decimal format(not minutes etc.)'''

nlat, wlon, slat, elon = ten_km_square(lat,lon)

#print("(NLat:{:.4f},WLon:{:.4f}),(SLat:{:.4f},ELon:{:.4f});".format(nlat, wlon, slat, elon))

region = ('[{:.4f},{:.4f}], [{:.4f},{:.4f}], [{:.4f},{:.4f}], [{:.4f},{:.4f}]'.format(wlon,nlat,elon,nlat,wlon,slat,elon,slat))

rectangle = ('[{:.4f},{:.4f},{:.4f},{:.4f}]'.format(wlon,slat,elon,nlat))

lister=region.split('/')

finallist=[ast.literal_eval(i) for i in lister] # to remove apostrophe

final = [list(value) for value in finallist] # final list

print(final)

import ee

import ee.mapclient

import ee

ee.Initialize()

import ee.mapclient

start = ee.Date('2017-01-01')

finish = ee.Date('2017-03-20')

rectangle = ee.Geometry.Polygon(final) #passing the final value as the polygon boundaries defined 10km square area

collection = ee.ImageCollection('COPERNICUS/S2')

filteredCollection = collection.filterBounds(rectangle).filterDate(start,finish)

first = filteredCollection.first()

path = first.getDownloadURL({

'region':final,

'scale': 10

})

print(path)

urls = [path]