LineString objects have no attribute Shapes

[Alvin-prc]

I came across this problem when I learn the example in https://flopy.readthedocs.io/en/stable/Notebooks/mfusg_freyberg_example.html.

Code:

from flopy.utils.geometry import LineString # 线状几何形状的类

lines = [
LineString(ls)
for ls in [
[(623000, 3364000), (623000, 3372000)], # 线段的两个端点的坐标
[(623650, 3364000), (623650, 3372000)],
]
]

"""Load the model and retrieve inactive cells."""

gwf = flopy.mfusg.MfUsg.load(
f"{root_name}.nam",
model_ws=str(workspace),
verbose=False,
check=False,
exe_name="mfusg",
)

bas6 = gwf.get_package("bas6")
ibound = bas6.ibound.array

""" Show the map view again with cross-section lines and inactive cells."""

fig = plt.figure(figsize=(8, 8))
ax = fig.add_subplot(1, 1, 1, aspect="equal")
pmv = flopy.plot.PlotMapView(modelgrid=mfgrid, ax=ax)
grid = pmv.plot_grid(alpha=0.2)
shps = pmv.plot_shapes(lines, edgecolor="purple", lw=2, alpha=0.7)
inac = pmv.plot_inactive(ibound=ibound)
plt.show()

Error:

已连接到 pydev 调试器(内部版本号 232.9559.58)Traceback (most recent call last):
File "C:\Program Files\JetBrains\PyCharm 2021.3.1\plugins\python\helpers\pydev\pydevd.py", line 1500, in _exec
pydev_imports.execfile(file, globals, locals) # execute the script
File "C:\Program Files\JetBrains\PyCharm 2021.3.1\plugins\python\helpers\pydev_pydev_imps_pydev_execfile.py", line 18, in execfile
exec(compile(contents+"\n", file, 'exec'), glob, loc)
File "D:\Project_Work\Python_Practice\Flopy-Practice\example_10.py", line 92, in
shps = pmv.plot_shapes(lines, edgecolor="purple", lw=2, alpha=0.7)
File "C:\ProgramData\Anaconda3\lib\site-packages\flopy\plot\map.py", line 597, in plot_shapes
patch_collection = plotutil.plot_shapefile(obj, ax, **kwargs)
File "C:\ProgramData\Anaconda3\lib\site-packages\flopy\plot\plotutil.py", line 2244, in plot_shapefile
pc = shapefile_to_patch_collection(shp, radius=radius, idx=idx)
File "C:\ProgramData\Anaconda3\lib\site-packages\flopy\plot\plotutil.py", line 2109, in shapefile_to_patch_collection
shapes = geofeats.shape
File "C:\ProgramData\Anaconda3\lib\site-packages\flopy\utils\geospatial_utils.py", line 447, in shape
self._shape = self.__shapefile.Shapes()
AttributeError: 'NoneType' object has no attribute 'Shapes'

Thank you for your time, I appreciate it.

[wpbonelli]

@Alvin-prc I suspect you are missing pyshp (available on PyPI, Conda). Many of the notebooks require flopy's optional dependencies. You can install all optional dependencies at once with pip install "flopy[optional]"