GeoSeries.touches(other) and GeoSeries.intersects(other) unexpected behavior

So using these two datasets, we ran into some unexpected behavior of the intersect and touches functions.

assigned_bins = cc.read('assigned_hexbins')
unassigned_bins = cc.read('unassigned_hexbins')

assigned_hexbins.csv.zip unassigned_hexbins.csv.zip

We converted the dataframes into geodataframes

unassigned_bins =cc.query('SELECT *, ST_AsText(the_geom) as geom FROM mspichigercarto.unassigned_hexbins').reset_index()
assigned_bins =cc.query('SELECT *, ST_AsText(the_geom) as geom FROM mspichigercarto.assigned_hexbins').reset_index()

import shapely.wkt

geometry = unassigned_bins['geom'].map(shapely.wkt.loads)
unassigned_bins = unassigned_bins.drop('geom', axis=1)
crs = {'init': 'epsg:4326'}
unassigned_bins_geo = gpd.GeoDataFrame(unassigned_bins, crs=crs, geometry=geometry)`
import shapely.wkt

geometry = assigned_bins['geom'].map(shapely.wkt.loads)
assigned_bins = assigned_bins.drop('geom', axis=1)
crs = {'init': 'epsg:4326'}
assigned_bins_geo = gpd.GeoDataFrame(assigned_bins, crs=crs, geometry=geometry)

Then we tried running touches and intersects on them, but instead of returning expected values it would return false for everything and none of the geometries intersected.

assigned_bins_geo.touches(unassigned_bins_geo)

We then ran this query in cartoframes

st_touch_query = '''
SELECT assigned_hexbins.*,
      unassigned_hexbins.cartodb_id as index_unassigned
FROM assigned_hexbins, unassigned_hexbins
WHERE ST_Touches(assigned_hexbins.the_geom, unassigned_hexbins.the_geom)
'''

according to this query, the hexbins are touching so touches and intersects are not working as expected.

[michellemho]

Yes, it doesn’t seem to be related to cartoframes. It has to do with Geopandas intersection methods.

On Tue, Jun 19, 2018 at 9:58 AM Andy Eschbacher [e-mail] wrote:

Hey all, how’s this going? Overall it looks like it doesn’t have to do with cartoframes

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/CartoDB/cartoframes/issues/449#issuecomment-398408594, or mute the thread https://github.com/notifications/unsubscribe-auth/AL-ECYlWt9REZNTieNA-Du3YqaFA4stkks5t-QOTgaJpZM4UjUwr .

– Michelle Ho Data Scientist [e-mail]

[andy-esch]

Closing this in favor of the geopandas one