I'm trying to retrieve a database with geometry information from a SQL server in Python using `sqlalchemy` and `geopandas`. There are 2 columns in geometry format and I want to retain both. The code I'm using is

```
import geopandas as gpd
from sqlalchemy import create_engine

db_connection_url = "postgresql://username:password@host:5432/database"
con = create_engine(db_connection_url)  
sql = 'SELECT osm_id, way, tags, way_centroid FROM osm.bldg WHERE height IS NOT NULL;'
df = gpd.read_postgis(sql, con, geom_col='way')

```

I want to achieve something similar in R which is `geom_col = c('way', 'way_centroid')` for the function`gpd.real_postgis`, but I'm aware that it's not going to work this way in python. How can I achieve this in python?