In Python, calls on a Point geometry cause crash

[lakshmanok]

In Python, calls on a Point geometry cause the interpreter to crash. The Windows error reads "calls on a pure virtual function".

I have attached a bare-bones Python script that will cause the crash. The shapefile that I was using is attached, but the bug is not limited to this shapefile. I believe that the crash happens on any shapefile that contains Points. I was able to reproduce the bug on the point shapefile available at: ​http://invisibleroads.com/tutorials/_downloads/gdal-shapefile-points.zip

thanks Lak

import sys import os import shutil from osgeo import ogr, osr

def get_location(feature):

geometry = feature.GetGeometryRef() x = geometry.GetX() #### Causes Crash y = geometry.GetY() return x,y

### my workaround ... to get it to work ##def get_location(feature): ## s = feature.ExportToJson() ## c = s[s.find(')+1 : s.find(')].split(',') ## x,y = float(c[0]),float(c[1]) ## return x,y

# read input inputshp = 'IOW_Fire_Hydrants.shp' datasource = ogr.GetDriverByName('ESRI Shapefile').Open(inputshp) if datasource is None:

print 'Could not open ' + inputshp sys.exit(-1)

inlayer = datasource.GetLayer() get_location(inlayer.GetFeature(1))

[lakshmanok]

Example shapefile with Points that illustrates the bug