Read Shapefile with GDAL

python

from osgeo import ogr
shapefile = ogr.Open("shapefile.shp")
layer = shapefile.GetLayer(0)
crs = layer.GetSpatialRef()
print(crs)
print(crs.GetAttrValue('AUTHORITY',1))
for i in range(layer.GetFeatureCount()):
    feature = layer.GetFeature(i)
    print(feature.GetGeometryRef())
    geometry = feature.GetGeometryRef()