Hi. Newbie here. I’m trying to create a workflow for processing ODIM OPERA structured radar data, and export rain depth maps. I was thinking this, and would appreciate any comment/correction/mockery (won’t take personally, i’m completely green…). couldn’t find any notebooks that do exactly this. most notebooks examples use the german input, which is quite different from OPERA.
A general workflow for processing radar and exporting rain maps to GIS:
- Date selector using range.
the simple and efficient method is to use file names, which contain the dates. it might not be suitable for selecting files from different sources.
optionally: download these files directly from API using the date selector - Load the database in the correct database format
wradlib uses a plenty of data structures: multidimencional numpy arrays, dictionaries, xarrays, datatrees and datasets. each algorithm recieves a different in put. too confusing… - Georef the dbz
I think I read somewhere it is essential to georef the data prior to cleaning and aggregating - Correct for attenuation
- aggregate the volume do a single 2d image (cappi)
use wradlib.vpr.CAPPI. Not sure about the product this method, but I’d like to have a 2d image of aggregated sweeps. - Convert to rain rate
using palmer or different - convert to rain depth
multiply by time difference (300 seconds) - Correct for rain gauges
use Adjusting radar-base rainfall estimates by rain gauge observations — wradlib - Export to gis
export to asc files. naming using the date and product.