fix to pyart (fields) bugs

kaomiao · Sep 23, 2019 · c497888 · c497888
1 parent 4e06129
commit c497888
Show file tree

Hide file tree

Showing 4 changed files with 8 additions and 4 deletions.
diff --git a/libs/io/CCFile.py b/libs/io/CCFile.py
@@ -367,7 +367,8 @@ def ToPyartRadar(self):
  if field_name is None:
  continue
  field_dic = get_metadata(field_name)
- field_dic['data'] = self.fields[field_name_abbr]
+ field_dic['data'] = np.ma.masked_array(self.fields[field_name_abbr],\
+ mask=np.isnan(self.fields[field_name_abbr]), fill_value=get_fillvalue())
  field_dic['_FillValue'] = get_fillvalue()
  fields[field_name] = field_dic
  return Radar(time, _range, fields, metadata, scan_type,

diff --git a/libs/io/SABFile.py b/libs/io/SABFile.py
@@ -433,7 +433,8 @@ def ToPyartRadar(self):
  if field_name is None:
  continue
  field_dic = get_metadata(field_name)
- field_dic['data'] = self.fields[field_name_abbr]
+ field_dic['data'] = np.ma.masked_array(self.fields[field_name_abbr],\
+ mask=np.isnan(self.fields[field_name_abbr]), fill_value=get_fillvalue())
  field_dic['_FillValue'] = get_fillvalue()
  fields[field_name] = field_dic
 

diff --git a/libs/io/SCFile.py b/libs/io/SCFile.py
@@ -374,7 +374,8 @@ def ToPyartRadar(self):
  if field_name is None:
  continue
  field_dic = get_metadata(field_name)
- field_dic['data'] = self.fields[field_name_abbr]
+ field_dic['data'] = np.ma.masked_array(self.fields[field_name_abbr],\
+ mask=np.isnan(self.fields[field_name_abbr]), fill_value=get_fillvalue())
  field_dic['_FillValue'] = get_fillvalue()
  fields[field_name] = field_dic
  return Radar(time, _range, fields, metadata, scan_type,

diff --git a/libs/io/WSR98DFile.py b/libs/io/WSR98DFile.py
@@ -441,7 +441,8 @@ def ToPyartRadar(self):
  if field_name is None:
  continue
  field_dic = get_metadata(field_name)
- field_dic['data'] = self.fields[field_name_abbr]
+ field_dic['data'] = np.ma.masked_array(self.fields[field_name_abbr],\
+ mask=np.isnan(self.fields[field_name_abbr]), fill_value=get_fillvalue())
  field_dic['_FillValue'] = get_fillvalue()
  fields[field_name] = field_dic
  return Radar(time, _range, fields, metadata, scan_type,