near/earth-observation-for-journalism
1
0
Fork 0
mirror of https://github.com/heyarne/earth-observation-for-journalism.git synced 2026-05-06 19:13:40 +02:00
Code Issues Projects Releases Wiki Activity

Add notebook containing researching the tiles needed to cover Brandenburg within the month of June

Browse source
This commit is contained in:
heyarne 2020-08-14 12:01:30 +00:00
commit a48cf96b14
15 changed files with 6830 additions and 1797 deletions
Download patch file Download diff file Expand all files Collapse all files

24
True-Color Image/sentinel_helpers.py
View file

@ -1,9 +1,12 @@
import zipfile
import glob
from Pathlib import Path
import urllib.parse
from pathlib import Path
import zipfile
import geopandas as gp
def band_paths(p, bands, resolution=None):
'''
Given a zip file or folder at `p`, returns the paths inside p to the raster files containing
@ -18,8 +21,7 @@ def band_paths(p, bands, resolution=None):
files = f.namelist()
rasters = [f for f in files if f.endswith('.jp2')]
else:
rasters = glob.glob(Path(p) / '**/*.jp2'
rasters = glob.glob(Path(p) / '**/*.jp2')
rasters = [r for r in rasters for b in bands if b in r]
if resolution:
rasters = [r for r in rasters if resolution in r]
@ -29,4 +31,14 @@ def band_paths(p, bands, resolution=None):
def rgb_paths(zip_file, resolution='10m'):
return band_paths(zip_file, ['B02', 'B03', 'B04'], resolution)
return band_paths(zip_file, ['B02', 'B03', 'B04'], resolution)
def search_osm(place):
'''
Returns a GeoDataFrame with results from OpenStreetMap Nominatim for the given search string.
This allows us to fetch detailed geometries for virtually any place on earth.
'''
urlescaped_place = urllib.parse.quote(place)
search_url = 'https://nominatim.openstreetmap.org/search/?q={}&format=geojson&polygon_geojson=1'.format(urlescaped_place)
return gp.read_file(search_url)