Search code examples
pythoncartopy

Geometry must be a Point or LineString error using Cartopy


I'm trying to run a simple Cartopy example:

   import cartopy.crs as ccrs
   import matplotlib.pyplot as plt
   ax = plt.axes(projection=ccrs.PlateCarree())
   ax.coastlines()
   plt.show()

But I'm getting this error:

Geometry must be a Point or LineString
python: geos_ts_c.cpp:4179: int GEOSCoordSeq_getSize_r(GEOSContextHandle_t, const geos::geom::CoordinateSequence*, unsigned int*): Assertion0 != cs' failed`.

I installed Cartopy using miniconda3: conda install -c conda-forge cartopy I also tried to install Cartopy with pip (in a virtual environment), but I get the same error. My OS is Debian Buster.

Any idea?


Solution

  • The problem is a wrong version of shapely, with Cartopy the binary package shoudn't be used, it should be built from source instead. This is explained here and here.

    So I did:

    pip uninstall shapely
    pip install shapely --no-binary shapely