I ran into an error when trying to read in shapes from SpatialData in the _wk_polygon call in _read_shapes_parquet. I figured out that, at lines 554 and 555, the geometry and index columns (that are read into geoms_wkb and raw_ids) were swapped around in my data from how they are expected to be. The function expects geometry to be first and index to be second, but when I generated mine, the index column was first and geometry was second. I am not sure if that is always the case or if this is fixable when constructing the SpatialData object but if this is variable across spatial platforms, would it be possible to check whether each column is an index and assume the other is the geometry (or vice versa)?
I ran into an error when trying to read in shapes from SpatialData in the _wk_polygon call in _read_shapes_parquet. I figured out that, at lines 554 and 555, the geometry and index columns (that are read into geoms_wkb and raw_ids) were swapped around in my data from how they are expected to be. The function expects geometry to be first and index to be second, but when I generated mine, the index column was first and geometry was second. I am not sure if that is always the case or if this is fixable when constructing the SpatialData object but if this is variable across spatial platforms, would it be possible to check whether each column is an index and assume the other is the geometry (or vice versa)?