id summary reporter owner description type status priority milestone component version severity resolution keywords cc 4487 Using SetGeom() for a PostGIS feature does not set geometry spatial reference system (Python bindings) hairymapper warmerdam "When doing: {{{ out_srs = ogr.osr.SpatialReference() out_srs.ImportFromEPSG(27700) geom.AssignSpatialReference(out_srs) feature.SetGeom(geom) layer.CreateFeature(feature) # where layer is a PostGIS table. }}} The resulting EKWT string passed to PostGIS contains an SRID value of -1 and not the SRS applied from out_srs. You can see my related question on GIS StackExchange at: http://gis.stackexchange.com/questions/19614/how-to-assign-ogr-geometry-a-spatial-reference-when-writing-to-a-new-layer-using" defect closed normal closed_because_of_github_migration default 1.8.1 normal wontfix Python, Spatial Reference, PostGIS