ds = Shapefile(file="/../..") for f in ds.features_at_point(Coord(-2,51)): print f TODO: 1.access to Feature properties 2.feature_at_point to accept screen coordinates 3.apply hit_test for geometries.