Python 底图"Assertion failed"

Question

我正在使用 Anaconda 2.1.0 并且底图安装成功

conda install basemap
conda install gdal

我有 geos 版本 3.3.3 和 proj4 4.8.0-0。但是，当我尝试将（纬度，纬度）度投影到平面时，我不断收到以下错误：

 Assertion failed: (0), function query, file AbstractSTRtree.cpp, line 286.

我的代码如下所示：

import shapely.geometry as sg
from pymongo import MongoClient
from mpl_toolkits.basemap import Basemap
p = Basemap(projection='sinu',lon_0=0,resolution='c')
projected_xy = p(24.4,45.1)

Answer 1

我在与地理熊猫一起工作时遇到了同样的问题。同样，进口顺序是相关的。在这种情况下，我必须在底图之后导入地理熊猫。我的导入列表如下所示：

from mpl_toolkits.basemap import Basemap
import geopandas as gp
from shapely.geometry import Point, Polygon, MultiPoint, MultiPolygon
from shapely.prepared import prep

请注意，如果您在 Jupyter 笔记本中工作，您的内核连接将崩溃，并且不会提供错误信息。使用 python shell 后，给出了以下错误：

Assertion failed: (0), function query, file AbstractSTRtree.cpp, line 285.
Abort trap: 6

Answer 2

出于某种原因，底图必须在形状之前导入。

from mpl_toolkits.basemap import Basemap
import shapely.geometry as sg

相反，执行导入将导致错误。