# path to osm.pbf file
osm_path = 'C:\planet_osm\tanzania-latest.osm.pbf'

# import through python
driver=ogr.GetDriverByName('OSM')
data = driver.Open(osm_path)
layer = data.GetLayerByName('lines')

features=[x for x in layer]
railway = []
for feature in features:
    if feature.GetField('railway') is not None:
        osm_id = feature.GetField('osm_id')
        shapely_geo = shapely.wkt.loads(feature.geometry().ExportToWkt())
        service = feature.GetField('service')
        maxspeed = feature.GetField('maxspeed')
        highway=feature.GetField('railway')
        railway.append([osm_id,highway,service,maxspeed,shapely_geo])

rail_dir_import = gpd.GeoDataFrame(railway,columns=['osm_id','infra_type','service','maxspeed','geometry'],crs={'init': 'epsg:4326'})

# import through ogr2ogr
shape_in = 'C:\planet_osm\tanzania-railway.shp'

os.system("ogr2ogr -progress -overwrite -f "ESRI Shapefile" -sql
               "SELECT osm_id,railway FROM lines WHERE railway IS NOT NULL"
               -lco ENCODING=UTF-8 -skipfailures -nlt geometry "+shape_in+" "+osm_path)

rail_com_line = gpd.read_file(shape_in)

layer = data.ExecuteSQL("SELECT osm_id,service,maxspeed,railway FROM lines WHERE railway IS NOT NULL")