Геоплот не будет наносить на карту объект геопанды - PullRequest
Новая
       31

Геоплот не будет наносить на карту объект геопанды

0 голосов
/ 10 декабря 2018

Я импортировал шейп-файл с округами штата Висконсин, используя geopandas и избирательный запас каждого округа для предыдущих выборов.Сейчас я пытаюсь разделить штат Висконсин по районам, используя тесселяцию Вороного доли голосов каждого округа.Я хочу использовать geoplot.voronoi, чтобы сделать это, но geoplot не может построить мою карту и вылетает, когда я вызываю geoplot.voronoi для данных моего шейп-файла.

Следующий код успешно генерирует карту доли голосования

import pandas as pd
import matplotlib.pyplot as plt
import geopandas as gp

wi = gp.GeoDataFrame.from_file('./data/shape/Wards_111312.shp')
elec_dat = pd.read_csv('./data/wi_2014_heda.csv')

wi['vrat'] = elec_dat['g2014_SOS_rv']
wi['vdat'] = elec_dat['g2014_SOS_dv']
wi['mr'] = (wi['vdat'] - wi['vrat']) / ((wi['vdat'] + wi['vrat']))
wi['mr'] = wi['mr'].fillna(0)

wi.plot(column = 'mr', cmap='OrRd')

imagegeopandas plotted margin map">

Но, когда я пытаюсь построить то же самое, используя geoplot, я получаю пустое изображение (видимый пробел ниже - это изображение). 

import geoplot
geoplot.polyplot(wi)

imagegeoplot plotted map (blank)">

Более того, если я попытаюсь построить тесселяцию Вороного: 

ax = geoplot.voronoi(
    wi.sample(1000),
    hue='mr', cmap='Reds', scheme='fisher_jenks',
    clip=wi.geometry,
    linewidth=0)
geoplot.polyplot(wi, ax=ax)

Я получаю следующую ошибку: 

/home/rtse/anaconda3/envs/gerry/lib/python3.6/site-packages/pysal/__init__.py:65: VisibleDeprecationWarning: PySAL's API will be changed on 2018-12-31. The last release made with this API is version 1.14.4. A preview of the next API version is provided in the `pysal` 2.0 prelease candidate. The API changes and a guide on how to change imports is provided at https://pysal.org/about
  ), VisibleDeprecationWarning)

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-2-5ed5e333efcd> in <module>
      7     hue='mr', cmap='Reds', scheme='fisher_jenks',
      8     clip=wi.geometry,
----> 9     linewidth=0)
     10 geoplot.polyplot(wi, ax=ax)
     11 

~/anaconda3/envs/gerry/lib/python3.6/site-packages/geoplot/geoplot.py in voronoi(df, projection, edgecolor, clip, hue, scheme, k, cmap, categorical, vmin, vmax, legend, legend_kwargs, legend_labels, extent, figsize, ax, **kwargs)
   2126 
   2127     # Finally we draw the features.
-> 2128     geoms = _build_voronoi_polygons(df)
   2129     if projection:
   2130         for color, geom in zip(colors, geoms):

~/anaconda3/envs/gerry/lib/python3.6/site-packages/geoplot/geoplot.py in _build_voronoi_polygons(df)
   2628     """
   2629     from scipy.spatial import Voronoi
-> 2630     geom = np.array(df.geometry.map(lambda p: [p.x, p.y]).tolist())
   2631     vor = Voronoi(geom)
   2632 

~/anaconda3/envs/gerry/lib/python3.6/site-packages/pandas/core/series.py in map(self, arg, na_action)
   2996         """
   2997         new_values = super(Series, self)._map_values(
-> 2998             arg, na_action=na_action)
   2999         return self._constructor(new_values,
   3000                                  index=self.index).__finalize__(self)

~/anaconda3/envs/gerry/lib/python3.6/site-packages/pandas/core/base.py in _map_values(self, mapper, na_action)
   1002 
   1003         # mapper is a function
-> 1004         new_values = map_f(values, mapper)
   1005 
   1006         return new_values

pandas/_libs/src/inference.pyx in pandas._libs.lib.map_infer()

~/anaconda3/envs/gerry/lib/python3.6/site-packages/geoplot/geoplot.py in <lambda>(p)
   2628     """
   2629     from scipy.spatial import Voronoi
-> 2630     geom = np.array(df.geometry.map(lambda p: [p.x, p.y]).tolist())
   2631     vor = Voronoi(geom)
   2632 

AttributeError: 'MultiPolygon' object has no attribute 'x'

Я пропустил шаг, необходимый для построения геометрии в geoplot?Я следовал этому примеру кода , но они загружают свои данные из geopandas.datasets, а не предоставляют свои собственные шейп-файлы.Что я делаю не так?

Заранее спасибо за любую помощь!

1 Ответ

0 голосов
/ 26 июля 2019

У меня просто была похожая проблема, и я обнаружил, что моя проблема заключалась в том, что мои значения LAT и LON были задом наперед.shapely.Point tuples должно быть (LON, LAT).

Другая возможная причина для пустого графика состоит в том, что график вообще не является пустым, это просто карта очень большого региона.Я думаю, что это маловероятно для вашего графика, но вы можете проверить это, передав кортеж со значениями min и max LON и LAT в переменную extent в gplt.polyplot ().

...