-
Notifications
You must be signed in to change notification settings - Fork 2
/
ConvertPAR.py
executable file
·37 lines (24 loc) · 1.12 KB
/
ConvertPAR.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
#!/usr/bin/env python
import os
import glob
from PARDataUtils import *
import numpy # for .meshgrid()
import scipy.signal # for .medfilt2d()
from RadarRastify import *
dataDirname = "/tint/PAR_wdssii/"
rastDirname = "/tint/PAR_wdssii/"
filelist = glob.glob(os.path.join(dataDirname, '*.netcdf'))
filelist.sort()
for filename in filelist :
print filename
# Loads the data from a file.
parData = LoadPAR_wdssii(filename)
parData['vals'] = scipy.signal.medfilt2d(parData['vals'], [1, 25])
# Rastify the data onto a Lat-Lon grid
[rastData, latAxis, lonAxis] = Rastify(parData['stat_lat'], parData['stat_lon'],
parData['vals'], parData['azimuth'], parData['range_gate'],
parData['elev_angle'], parData['beam_width'], parData['gate_length'],
0.005)
outfile = os.sep.join([rastDirname, datetime.datetime.utcfromtimestamp(parData['scan_time']).strftime('PAR_%Y%m%d_%H%M%S_rast.nc')])
SavePAR_NetCDF(outfile, rastData, latAxis, lonAxis, parData['scan_time'], parData['var_name'])
print outfile