ValueError on setting a geometry to a multipolygon ("equal len keys and value")

I have 2 geodataframes with a geometry column and I copy some geometries from 1 to the other. This works well with Polygons but returns a ValueError for any (valid) MultiPolygon.

geodata01.loc[index, 'geometry'] = geodata02.loc[index, 'geometry']

This works well with Polygons. Only with MultiPolygons I get:

ValueError 
if len(labels) != len(value):
--> 611      raise ValueError('Must have equal len keys and value '
    612      'when setting with an iterable')

I also cannot do assignments of a buffered or simplified MultiPolygon either (the MultiPolygon is valid and I can plot, buffer, simplify but I cannot assign it):

geodata01.loc[index, 'geometry'] = geodata01.loc[index, 'geometry'].buffer(0)
#or
geodata01.loc[index, 'geometry'] = geodata01.loc[index, 'geometry'].simplify(tolerance=0)

This returns the same ValueError. I have tried many geometries and (only) encountered this with (all) multipolygons that I tried. Attached an example of a multipolygon that yields the error:

MULTIPOLYGON.txt

[elizabethswkim]

Thank you so much for raising this, @wbijster! And for your responses, @jorisvandenbossche! I couldn’t make this work by assigning to a list of the scalar element, nor could I by indexing on a mask matching the scalar element. I was only able to get this assignment to work using @koshy1123’s solution at @jorisvandenbossche’s follow-up pandas issue, specifically: gdf.loc[scalar_index_loc, 'geometry'] = geopandas.GeoDataFrame(geometry=[shapely_geo]).geometry.values

[martinfleis]

Just posting another workaround here. Go through single element GeoSeries. Might be helpful for some cases.

geom = df.loc[1, 'geometry']
df.loc[[0], 'geometry'] = geopandas.GeoSeries([geom])

~If I am correct, this bug could be in future resolved by Shapely 2.0, which makes geometries non-iterable by default - https://github.com/shapely/shapely-rfc/pull/1.~