Using wradlib to obtain PBB and CBB, testing the terrain display of GeoTIFF file.

Code has been adapted by code from a notebook by Scott Collis: https://github.com/EVS-ATMOS/ARM-Azores/blob/master/notebooks/2.%20Terrain%20of%20ENA.ipynb

Code for projection conversion by Kai Muehlbauer and Wradlib notebook by Kai: http://wradlib.org/wradlib-docs/latest/notebooks/beamblockage/wradlib_beamblock.html

And Code by Nick Guy from PyRadarMet: https://github.com/nguy/PyRadarMet


In [1]:
from matplotlib import pyplot as plt
import gdal
import osr
import numpy as np
from mpl_toolkits.basemap import Basemap
import pickle
import netCDF4
import re
import pyart
import wradlib as wrl
import plotly.plotly as py
import plotly.graph_objs as go
import matplotlib
from matplotlib import cm
from mpl_toolkits.mplot3d import Axes3D
from plotly.graph_objs import Surface, Layout, Scene
from plotly.offline import(
    download_plotlyjs, init_notebook_mode, plot, iplot)

#warnings.filterwarnings("ignore")
init_notebook_mode()
%matplotlib inline


## You are using the Python ARM Radar Toolkit (Py-ART), an open source
## library for working with weather radar data.
##
## If you use this software to prepare a publication, please cite:
##
##     JJ Helmus and SM Collis, JORS 2016, doi: 10.5334/jors.119