I am using geopandas to import a shapefile and then check if a point is within that shape. I am using the shapefile found in http://opendata.atlantaregional.com/datasets/cities-georgia?geometry=-87.766%2C33.036%2C-81.448%2C34.633 to get the shape of Atlanta. However, whenever I "ask" if a point is within that shape (the point is surely within the shape), the function returns False. I have tried creating a Polygon and then asking if a point is within that polygon and it works as expected.
I don't know what else to do to troubleshot or debug this issue.
I am doing all of this using a jupyter notebook. Here is the MWE: https://github.com/vcaballero-salle/mwe-shapefile
The code looks like (full code with imports can be found at the github repo, in the README.md):
plt.clf() shape = gp.read_file(r'./Cities_Georgia/Cities_Georgia.shp') atlanta_shape = shape[shape["Name"] == "Atlanta"] atlanta_shape.plot() plt.show() plt.ion() # In: def within_shape(df, shapes): in_shape =  for sh in shapes.geometry: within = df.within(sh) in_shape.append(within) return in_shape print(within_shape(Point(33.751079, -84.443274), atlanta_shape)) # This returns [False]