#!/usr/bin/env python # coding: utf-8 # In[5]: import geopandas as gpd import pandas as pd import rasterio as rio # In[10]: #set working directory import os os.chdir('C:\\Users\\jtrum\\world_bank\\data') #read in files luanda = gpd.read_file('luanda2clean.geojson') with rio.open('ESA_WorldCover_10m_2021_v200_S12E012_Map.tif') as src: landCover = src.read() profile = src.profile # convert to dataframe landCover = pd.DataFrame(landCover[0]) # In[11]: landCover # In[13]: from sentinelsat import SentinelAPI, read_geojson # In[14]: username = 'jtrummler' password = 'Cf5lvg0503!' # In[15]: api = SentinelAPI(username, password, 'https://scihub.copernicus.eu/dhus') # In[16]: footprint = read_geojson('luanda2clean.geojson') # In[26]: products = api.query(footprint, date=('20/01/31', '20/01/01'), platformname='Sentinel-2', cloudcoverpercentage=(0, 10))