Error code, SF-files

[Johannes Paul]

Hello,

I just started using the wradlib package for the reading of SF-files from the DWD. I use the following code which is pretty much the same as presented in the tutorials. I am not really familiar with coding. So I cannot find the reason for some errors that accur. Could you help me?

Thank you in advance.

Johannes

import wradlib as wrl
import numpy as np
import matplotlib.pyplot as pl
pl.interactive(True)
import os

sf_filename = 'C:\Users\jmoellm\Desktop\raa01-sf_10000-0610300750-dwd---bin'
print(sf_filename)
sfdata, sfattrs = wrl.io.read_RADOLAN_composite(sf_filename)
   
print("SF Attributes:", sfattrs)

sec = sfattrs['secondary']
sfdata.flat[sec] = -9999
sfdata = np.ma.masked_equal(sfdata, -9999)

radolan_grid_xy = wrl.georef.get_radolan_grid(900,900)
x = radolan_grid_xy[:,:,0]
y = radolan_grid_xy[:,:,1]

pl.pcolormesh(x, y, sfdata, cmap="spectral")
cb = pl.colorbar(shrink=0.75)
cb.set_label("mm\h")
pl.title('RADOLAN sf Product Polar Stereo \n' + sfattrs['datetime'].isoformat())
pl.grid(color='r')

pl.show()


>>>
>>> import wradl>>> import wradlib as wrl
>>> import numpy as np
>>> import matplotlib.pyplot as pl
>>> pl.interactive(True)
>>> import os
>>>
>>> sf_filename = 'C:\Users\jmoellm\Desktop\raa01-sf_10000-0610300750-dwd---bin'

>>> print(sf_filename)
aa01-sf_10000-0610300750-dwd---bin
>>> sfdata, sfattrs = wrl.io.read_RADOLAN_composite(sf_filename)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\wradlib\io.py", line 818, in read_RADOLAN_composite
    f = get_radolan_filehandle(fname)
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\wradlib\io.py", line 746, in get_radolan_filehandle
    f = open(fname, 'rb')
IOError: [Errno 22] invalid mode ('rb') or filename: 'C:\\Users\\jmoellm\\Deskto
p\raa01-sf_10000-0610300750-dwd---bin'
>>>
... print("SF Attributes:", sfattrs)
('SF Attributes:', {'maxrange': '128 km', 'radarlocations': ['bln', 'drs', 'eis'
, 'emd', 'ess', 'fbg', 'fld', 'fra', 'ham', 'han', 'muc', 'neu', 'nhb', 'ros', '
tur', 'umd'], 'nrow': 900, 'intervalseconds': 86400, 'cluttermask': array([], dt
ype=int64), 'precision': 0.1, 'datetime': datetime.datetime(2006, 10, 30, 18, 50
), 'ncol': 900, 'radolanversion': '1.5.0.1', 'producttype': 'SF', 'nodataflag':
-9999, 'datasize': 1620000, 'radarid': '10000', 'secondary': array([], dtype=int
64)})
>>>
>>> sec = sfattrs['secondary']
>>> sfdata.flat[sec] = -9999
>>> sfdata = np.ma.masked_equal(sfdata, -9999)
>>>
>>> radolan_grid_xy = wrl.georef.get_radolan_grid(900,900)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\wradlib\georef.py", line 1442, in get_radolan_grid
    x_0, y_0 = get_radolan_coords(9.0, 51.0, trig=trig)
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\wradlib\georef.py", line 1331, in get_radolan_coords
    x, y = reproject(lon, lat, projection_source=proj_wgs, projection_target=pro
j_stereo)
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\wradlib\georef.py", line 1188, in reproject
    trans = np.array(ct.TransformPoints(C))
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\osgeo\osr.py", line 800, in TransformPoints
    return _osr.CoordinateTransformation_TransformPoints(self, *args)
TypeError: in method 'CoordinateTransformation_TransformPoints', argument 1 of t
ype 'OSRCoordinateTransformationShadow *'
>>> x = radolan_grid_xy[:,:,0]
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
NameError: name 'radolan_grid_xy' is not defined
>>> y = radolan_grid_xy[:,:,1]
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
NameError: name 'radolan_grid_xy' is not defined
>>>
>>> pl.pcolormesh(x, y, sfdata, cmap="spectral")
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
NameError: name 'x' is not defined
>>> cb = pl.colorbar(shrink=0.75)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "C:\Users\jmoellm\AppData\Local\Continuum\Anaconda2\envs\wradlib\lib\site
-packages\matplotlib\pyplot.py", line 2230, in colorbar
    raise RuntimeError('No mappable was found to use for colorbar '
RuntimeError: No mappable was found to use for colorbar creation. First define a
 mappable such as an image (with imshow) or a contour set (with contourf).
>>> cb.set_label("mm\h")
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
NameError: name 'cb' is not defined
>>> pl.title('RADOLAN sf Product Polar Stereo \n' + sfattrs['datetime'].isoforma
t())
<matplotlib.text.Text object at 0x00000000127A6BE0>
>>> pl.grid(color='r')
>>>
>>> pl.show()ib as wrl

[Johannes Paul]

Esspecally this error is a problem:

TypeError: in method 'CoordinateTransformation_TransformPoints', argument 1 of t
ype 'OSRCoordinateTransformationShadow *'

The error occurs, when this line is executed:

>>> radolan_grid_xy = wrl.georef.get_radolan_grid(900,900)

Could it be possible that the GDAL_DATA environment variable isn't set?