Description:
bomber 0.0.5
# Read data from the Bureau of Meterology (BoM)author: Jess Robertson (@jesserobertson)The BoM, in it's wisdom, have made getting their geospatial data fairly difficult. This package aims to fix thatTo see this readme as an ipython notebook, check out https://github.com/jesserobertson/bomber/blob/master/examples/request_bom_data.ipynb## Example usageCheck out the measurement/observation datasets that we have available:```python>>> import bomber>>> print(bomber.measurements.DATASETS)``````{'ndvi': ['ndviave'], 'rainfall': ['totals'], 'solar': ['solarave'], 'temperature': ['maxave', 'minave'], 'vprp': ['vprp09', 'vprp15']}```Then we can get the bit that we want as a geotiff:```python>>> geotiff = bomber.get_measurements(dataset='rainfall', year=2015, month=1)``````Downloaded data to rainfall_totals_month_2015010120150131.geotiff```and then plot it using rasterio```python>>> import rasterio, numpy>>> import matplotlib.pyplot as plt>>> with rasterio.drivers():... with rasterio.open(geotiff) as src:... fig = plt.figure(figsize=(11, 11))... data = numpy.ma.MaskedArray(... data=src.read(1), ... mask=src.read_masks(1))... ax = fig.gca()... ax.imshow(data, cmap=plt.get_cmap('coolwarm'))... ax.set_aspect('equal')... ax.set_axis_off()... ax.set_title('Rainfall, January 2015')```You can also get the climatic average datasets as well:```python>>> print(bomber.climate.DATASETS)``````{'decadal-rainfall': ['r'], 'decadal-temperature': ['mx', 'mn']}```There's also a borked version for the climate regions but I'm too lazy to finish that one. Pull requests welcome!
License
For personal and professional use. You cannot resell or redistribute these repositories in their original state.
Customer Reviews
There are no reviews.
