GIS: How to use OGR Simplify in Python?

GIS: How to use OGR Simplify in Python?

[manuelpopp]

I wondered how to overwrite the entire shapefile with simplified features. I think I've figured it out, maybe this helps someone else:
def simplify_poly(path, tolerance = 0.05):
shp = ogr.Open(path, 1)
layer = shp.GetLayer()
feature = layer.GetNextFeature()
while feature:
geom = feature.geometry()
simple = geom.Simplify(tolerance)
feature.SetGeometry(simple)
layer.SetFeature(feature)
feature = layer.GetNextFeature()
layer = None
shp = None