GeoPandas - данные с разбросом по сетке и перепроектирование

Мне нужно сопоставить разбросанные данные в кадре данных GeoPandas с регулярной сеткой (например, 1 градус) и получить средние значения отдельных ячеек сетки, а затем построить эти данные с различными проекциями.

Первое, чего мне удалось достичь, используя gpd_lite_toolbox.

Этот результат я могу изобразить на простой карте широты и долготы, однако попытка преобразовать его в любую другую проекцию не удалась.

Вот небольшой пример с искусственными данными, показывающими мою проблему:

import gpd_lite_toolbox as glt
import geopandas as gpd
import matplotlib.pyplot as plt
import pandas as pd
from shapely import wkt

# creating the artificial df
df = pd.DataFrame(
{'data': [20, 15, 17.5, 11.25, 16],
 'Coordinates': ['POINT(-58.66 -34.58)', 'POINT(-47.91 -15.78)',
                 'POINT(-70.66 -33.45)', 'POINT(-74.08 4.60)',
                 'POINT(-66.86 10.48)']})

# converting the df to a gdf with projection
df['Coordinates'] = df['Coordinates'].apply(wkt.loads)
crs = {'init': 'epsg:4326'}
gdf = gpd.GeoDataFrame(df, crs=crs, geometry='Coordinates')

# gridding the data using the gridify_data function from the toolbox and setting grids without data to nan
g1 = glt.gridify_data(gdf, 1, 'data', cut=False)
g1 = g1.where(g1['data'] > 1)

# simple plot of the gridded data
fig, ax = plt.subplots(ncols=1, figsize=(20, 10))
g1.plot(ax=ax, column='data', cmap='jet')

# trying to convert to (any) other projection
g2 = g1.to_crs({'init': 'epsg:3395'})

# I get the following error
---------------------------------------------------------------------------

AttributeError: 'float' object has no attribute 'is_empty'

Я также был бы рад использовать другую функцию сетки, если это решит проблему.