[GIS] TypeError when trying to plot a combined geopandas geodataframe

geopandaspythonshapefile

I have GeoDataFrame that I am trying to plot in relation to 4 reference points. I'm creating the reference points as part of a separate GeoDataFrame, combining them and then trying to plot the result. I am using the following code:

import pandas
import geopandas as gpd
from geopandas import GeoDataFrame
from shapely.geometry import Point

def dfTogpdf(df):
    geometry = [Point(xy) for xy in zip(df.Lon, df.Lat)]
    df = df.drop(['Lon', 'Lat'], axis=1)
    crs = {'init': 'epsg:4326'}

    return GeoDataFrame(df, crs=crs, geometry=geometry)


habitat = gpd.read_file('desktop/species_19377/species_19377.shp')

latitude = [90, -90, 0, 0]
longitude = [0, 0, -179, 180]
frame = dfTogpdf(pandas.DataFrame({'Lat': latitude, 'Lon': longitude}))

framed = gpd.sjoin(frame, habitat,  how="inner", op='intersects')
framed.geometry.plot()

However this gives me the following error: TypeError: invalid type comparison

I am able to plot habitat and frame without any issue, but after they have been combined, I am left with an empty GeoDataFrame and not able to plot the result. Is there an extra step I need to take to plot my combined GeoDataFrames?

species_19377.shp

Traceback:

C:\Users\Alec\Anaconda3\lib\site-packages\pandas\core\ops.py:798: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  result = getattr(x, name)(y)
Traceback (most recent call last):

  File "<ipython-input-255-31bfbe140734>", line 1, in <module>
    runfile('C:/Users/Alec/test.py', wdir='C:/Users/Alec')

  File "C:\Users\Alec\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 705, in runfile
    execfile(filename, namespace)

  File "C:\Users\Alec\Anaconda3\lib\site-packages\spyder\utils\site\sitecustomize.py", line 102, in execfile
    exec(compile(f.read(), filename, 'exec'), namespace)

  File "C:/Users/Alec/test.py", line 29, in <module>
    framed.geometry.plot()

  File "C:\Users\Alec\Anaconda3\lib\site-packages\geopandas\geoseries.py", line 253, in plot
    return plot_series(self, *args, **kwargs)

  File "C:\Users\Alec\Anaconda3\lib\site-packages\geopandas\plotting.py", line 273, in plot_series
    poly_idx = np.asarray((geom_types == 'Polygon')

  File "C:\Users\Alec\Anaconda3\lib\site-packages\pandas\core\ops.py", line 861, in wrapper
    res = na_op(values, other)

  File "C:\Users\Alec\Anaconda3\lib\site-packages\pandas\core\ops.py", line 800, in na_op
    raise TypeError("invalid type comparison")

TypeError: invalid type comparison

habitat.head():

             ID_NO        BINOMIAL  PRESENCE  ORIGIN COMPILER  YEAR  \
index_right                                                           
0            19377  Ratufa bicolor         1       1     IUCN  2008   
1            19377  Ratufa bicolor         1       1     IUCN  2008   
2            19377  Ratufa bicolor         1       1     IUCN  2008   
3            19377  Ratufa bicolor         1       1     IUCN  2008   
4            19377  Ratufa bicolor         1       1     IUCN  2008   

                                                      CITATION SOURCE  \
index_right                                                             
0            IUCN (International Union for Conservation of ...   None   
1            IUCN (International Union for Conservation of ...   None   
2            IUCN (International Union for Conservation of ...   None   
3            IUCN (International Union for Conservation of ...   None   
4            IUCN (International Union for Conservation of ...   None   

            DIST_COMM    ISLAND SUBSPECIES SUBPOP TAX_COMMEN  \
index_right                                                    
0                None   Con Son       None   None       None   
1                None      None       None   None       None   
2                None      Bali       None   None       None   
3                None      None       None   None       None   
4                None  Billiton       None   None       None   

                        LEGEND  SEASONAL  \
index_right                                
0            Extant (resident)         1   
1            Extant (resident)         1   
2            Extant (resident)         1   
3            Extant (resident)         1   
4            Extant (resident)         1   

                                                      geometry  
index_right                                                     
0            POLYGON ((106.6621000000001 8.752335593000055,...  
1            POLYGON ((113.1460000000001 -8.175164407999944...  
2            POLYGON ((115.2932000000001 -8.102764407999928...  
3            POLYGON ((97.17740000000003 2.226135593000038,...  
4            POLYGON ((107.8088000000001 -2.535364406999975...

framed.head():

Empty GeoDataFrame
Columns: [geometry, index_right, ID_NO, BINOMIAL, PRESENCE, ORIGIN, COMPILER, YEAR, CITATION, SOURCE, DIST_COMM, ISLAND, SUBSPECIES, SUBPOP, TAX_COMMEN, LEGEND, SEASONAL]
Index: []

Best Answer

You are trying to plot an empty GeoSeries, as evidenced by your post:

>>> framed.head():
Empty GeoDataFrame
Columns: [geometry, index_right, ID_NO, BINOMIAL, PRESENCE, ORIGIN, COMPILER, YEAR, CITATION, SOURCE, DIST_COMM, ISLAND, SUBSPECIES, SUBPOP, TAX_COMMEN, LEGEND, SEASONAL]
Index: []

The TypeError is being triggered by the fact that this is an empty GeoDataFrame; the plotting code did not expect this edge case. This was a known bug and has been fixed, but the geopandas project hasn't had a new release since then. If you really need it fixed, you will have to install bleeding-edge geopandas from the source on GitHub.

But on the other hand, maybe you're looking in the wrong place. Perhaps you didn't expect an empty result from the spatial join?

Related Question