Use df.overlay instead of homegrown method to compute interstections?

[fmaussion]

See http://geopandas.org/set_operations.html#set-operations-with-overlay