Skip to content

TypeError When Plotting Mesh with Cell Boundary on Antimeridian #61

@gecoombs

Description

@gecoombs

If I run plot_mesh on a mesh with just two cells, each one exactly touching the anti-meridian I get the below error:

(venv) gecoomb@BSL-GECOOMB-L1:~/PycharmProjects/PolarRoute$ plot_mesh test_meshes/horizontal_0lat_antimeridian_mesh.json
[15-05-24 12:00:10 :INFO    ] - plot_mesh 0.1.0
Traceback (most recent call last):
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/bin/plot_mesh", line 8, in <module>
    sys.exit(plot_mesh_cli())
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/lib/python3.8/site-packages/bas_geoplot/utils.py", line 31, in wrapper
    res = func(*args, **kwargs)
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/lib/python3.8/site-packages/bas_geoplot/cli.py", line 93, in plot_mesh_cli
    mp.Maps(mesh, 'MeshGrid', predefined='cx')
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/lib/python3.8/site-packages/bas_geoplot/interactive.py", line 538, in Maps
    dataframe_pandas = split_at_antimeridian(dataframe_pandas)
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/lib/python3.8/site-packages/bas_geoplot/utils.py", line 187, in split_at_antimeridian
    new_polygon = MultiPolygon([poly_a, poly_b])
  File "/home/gecoomb/PycharmProjects/PolarRoute/venv/lib/python3.8/site-packages/shapely-2.0.1-py3.8-linux-x86_64.egg/shapely/geometry/multipolygon.py", line 79, in __new__
    shell = ob[0]
TypeError: 'MultiLineString' object is not subscriptable

Mesh file attached:
horizontal_0lat_antimeridian_mesh.json

Metadata

Metadata

Assignees

Labels

bugSomething isn't working

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions