Is "polygon_query" function working properly?

[yt723]

Hi,

I'm using napari_spatialdata to segment regions of interest.
It seems the filtered table doesn't include the correct cell IDs.
I noticed all polygon shapes changed into rectangles.
I attached an example below. The triangle is the ROI I placed and the blue rectangle area shows cells obtained by polygon_query.
I'm using Seurat for clustering etc. So, I export the filtered table and add annotations as metadata.
This happens no matter what shapes of polygon I create.
However, "sdata.pl.render_images().pl.render_shapes(color='annotation').pl.show(coordinate_systems=["x"])" displays the actual polygon.

Why does this discrepancy happen (is this happening only to me?) and is there any way to solve this problem?
Thank you so much for your help!

[LucaMarconato]

Hi @yt723 thanks for reporting and happy to look into this but I will need some help in reproducing the bug please.

1. Does the bug appear also with the latest main of spatialdata and napari-spatialdata?

2. Could you please try reproducing the bug with the following dataset:

from spatialdata.datasets import blobs_annotating_element

sdata = blobs_annotating_element('blobs_polygons')

3. If you can't reproduce with the dataset above, is there any change you could share the dataset you are using in your screenshot?

Thanks!