Skip to content

Commit

Permalink
fix pyart start dtype problems
Browse files Browse the repository at this point in the history
  • Loading branch information
zy committed Oct 10, 2019
1 parent be3fe30 commit d8032c7
Show file tree
Hide file tree
Showing 3 changed files with 8 additions and 9 deletions.
4 changes: 2 additions & 2 deletions NuistRadar/__init__.py
Original file line number Diff line number Diff line change
@@ -1,2 +1,2 @@
from . import core, draw, io, interp, GraphicalInterface
__all__ = ["core", "draw", "io", "interp", "GraphicalInterface"]
from . import core, draw, io, interp
__all__ = ["core", "draw", "io", "interp"]
7 changes: 4 additions & 3 deletions NuistRadar/io/SCFile.py
Original file line number Diff line number Diff line change
Expand Up @@ -56,9 +56,10 @@ def _parse_BaseDataHeader(self, buf_header):
##解码其余一些观测参数
BaseDataHeader_dict['RadarObserationParam_2'], _ = _unpack_from_buf(buf_header,\
dtype_sc.RadarObserationParamPos_2, dtype_sc.BaseDataHeader['RadarObserationParam_2'])
self.nrays = np.sum(BaseDataHeader_dict['LayerParam']['recordnumber'])
self.sweep_end_ray_index_add1 = np.cumsum(BaseDataHeader_dict['LayerParam']['recordnumber']) ##python格式的结束
self.sweep_start_ray_index = self.sweep_end_ray_index_add1 - BaseDataHeader_dict['LayerParam']['recordnumber']
self.nrays = np.sum((BaseDataHeader_dict['LayerParam']['recordnumber']).astype(np.int64))
self.sweep_end_ray_index_add1 = np.cumsum((BaseDataHeader_dict['LayerParam']['recordnumber']).astype(np.int64)) ##python格式的结束
self.sweep_start_ray_index = self.sweep_end_ray_index_add1 - \
(BaseDataHeader_dict['LayerParam']['recordnumber']).astype(np.int64)
return BaseDataHeader_dict

def _parse_radial(self):
Expand Down
6 changes: 2 additions & 4 deletions examples/test.py
Original file line number Diff line number Diff line change
Expand Up @@ -5,9 +5,7 @@
import matplotlib.pyplot as plt

file = r"E:\RadarBaseData\CINRAD-SA\温州\2015080816.59A"
file = r"C:\Users\zy\Desktop\Z_RADR_I_Z9831_20190427200200_O_DOR_SC_CAP.bin"
NRadar = radar_io(file).ToNuistRadar()

sec = VerticalSection(NRadar)
dat = sec.section((0, 0), (-80000, 80000), "dBZ")
plt.show()

Radar = radar_io(file).ToPyartRadar()

0 comments on commit d8032c7

Please sign in to comment.