<div>Hello list!</div>
<div> </div>
<div>I&#39;m trying to get the projection from a shapefile using a python script and gdal, but it didn&#39;t work. I want to show the EPSG and UTM projection, my code is:</div>
<div> </div>
<div>from osgeo import ogr<br>from osgeo import osr<br>shapeData = ogr.Open(filepath)</div>
<div>layer = shapeData.GetLayer()<br>feature = layer.GetNextFeature()</div>
<div>spatial = layer.GetSpatialRef()</div>
<div>print spatial</div>
<div> </div>
<div>The printing that shows to me, doesn&#39;t please me because it prints many times &quot;WGS_1984_UTM_Zone_13N&quot;, some of mt shapes don&#39;t work with this projection, any idea what I&#39;m doing wrong?<br></div>