PyCINRAD模块处理雷达基数据指南（发布1.7.0版本）(修复经纬度标注问题)

Eegle

eeeee 发表于 2020-1-10 11:14
这个问题在发布1.5版本之后才发现，你可以去GitHub更新一下源码

好的 感谢老师

李山山0103

老师你好，我用的是win10的子系统ubuntu，在PPI画图的过程中有如下提示：老师可曾遇到过？
脚本：import cinrad
f = cinrad.io.CinradReader('/mnt/e/ubuntu/PyCINRAD/shuju/Z_RADR_I_Z9551_20190620234900_O_DOR_SA_CAP.bin.BZ2')
r = f.get_data(0, 230, 'REF') #获取数据
print(r)
fig = cinrad.visualize.PPI(r)
fig('/mnt/e/ubuntu/PyCINRAD/') #Pass the path to save the fig
报错：Geometry must be a Point or LineString
python3: geos_ts_c.cpp:3991: int GEOSCoordSeq_getSize_r(GEOSContextHandle_t, const geos::geom::CoordinateSequence*, unsigned int*): Assertion `0 != cs' failed.
Aborted (core dumped)

eeeee

李山山0103 发表于 2020-2-2 22:04
老师你好，我用的是win10的子系统ubuntu，在PPI画图的过程中有如下提示：老师可曾遇到过？
脚本：import c ...

有没有更详细一点的错误报告

李山山0103

eeeee 发表于 2020-2-3 05:01
有没有更详细一点的错误报告

<Datatype: REF Station name: 合肥 Scan time: 2019-06-20 23:49:21.375000 Elevation angle: 0.61 Range: 230>
Geometry must be a Point or LineString
python3: geos_ts_c.cpp:3991: int GEOSCoordSeq_getSize_r(GEOSContextHandle_t, const geos::geom::CoordinateSequence*, unsigned int*): Assertion `0 != cs' failed.
Aborted (core dumped)
以上是全部报错信息，数据没问题，画图有问题

eeeee

李山山0103 发表于 2020-2-3 09:46
Geometry must be a Point or LineString
python3: geos_ts_c.cpp:3991: int GEOSCoordSeq_getSize_r( ...

我感觉是cartopy或者shapely的问题，你画其他数据有没有问题

李山山0103

eeeee 发表于 2020-2-5 03:24
我感觉是cartopy或者shapely的问题，你画其他数据有没有问题

谢谢老师，再重装试一下

李山山0103

eeeee 发表于 2020-2-5 03:24
我感觉是cartopy或者shapely的问题，你画其他数据有没有问题

Traceback (most recent call last):
  File "test.py", line 16, in <module>
    fig = PPI(cr) #Plot PPI
  File "/mnt/d/PyCINRAD/cinrad/visualize/ppi.py", line 102, in __init__
    self._plot(**kwargs)
  File "/mnt/d/PyCINRAD/cinrad/visualize/ppi.py", line 190, in _plot
    extent=self.geoax.get_extent(self.data_crs),
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 209, in add_shp
    shps = get_shp()
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 197, in get_shp
    shps = [list(ShpReader(i).geometries()) for i in flist]
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 197, in <listcomp>
    shps = [list(ShpReader(i).geometries()) for i in flist]
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 135, in __init__
    self._geometry_factory = shapereader.GEOMETRY_FACTORIES.get(shapeType)
AttributeError: module 'cartopy.io.shapereader' has no attribute 'GEOMETRY_FACTORIES'
你好，我重新安装了cartopy又报上面错误

eeeee

李山山0103 发表于 2020-2-5 11:51
Traceback (most recent call last):
  File "test.py", line 16, in
    fig = PPI(cr) #Plot PPI

你试一下from cartopy.io.shapereader import GEOMETRY_FACTORIES，看看会不会报错

李山山0103

eeeee 发表于 2020-2-6 11:14
你试一下from cartopy.io.shapereader import GEOMETRY_FACTORIES，看看会不会报错

按照您的要求输入仍然会报错
ImportError: cannot import name 'GEOMETRY_FACTORIES' from 'cartopy.io.shapereader' (/home/lss/anaconda3/lib/python3.7/site-packages/Cartopy-0.17.1.dev237_-py3.7-linux-x86_64.egg/cartopy/io/shapereader.py)
我进入utils.py后注释掉下面语句就可以画图了
##        self._geometry_factory = shapereader.GEOMETRY_FACTORIES.get(shapeType)
##        if self._geometry_factory is None:
##            raise ValueError("Unsupported shape type: %s" % shapeType)

李山山0103

eeeee 发表于 2020-2-6 11:14
你试一下from cartopy.io.shapereader import GEOMETRY_FACTORIES，看看会不会报错

Traceback (most recent call last):
  File "test.py", line 18, in <module>
    fig = PPI(cr,style='white', highlight='合肥市') #Plot PPI
  File "/mnt/d/PyCINRAD/cinrad/visualize/ppi.py", line 102, in __init__
    self._plot(**kwargs)
  File "/mnt/d/PyCINRAD/cinrad/visualize/ppi.py", line 193, in _plot
    draw_highlight_area(self.settings["highlight"])
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 273, in draw_highlight_area
    lines = highlight_area(area)
  File "/mnt/d/PyCINRAD/cinrad/visualize/utils.py", line 263, in highlight_area
    target = np.array(rec)[(name == i).nonzero()[0]]
AttributeError: 'bool' object has no attribute 'nonzero'
我画区域高亮时出现以上问题，不知道您遇到过没？