JimShady opened a new issue, #1523:
URL: https://github.com/apache/sedona/issues/1523

   Two (almost identical) queries below. One in PostGIS returns 78 features, 
the same one in apache-sedona returns 1 feature. I don't understand why - can 
anyone help/advise please?
   
   ### PostGIS query - returns 78 features
   
   ```
   WITH a_table AS (
   SELECT st_boundary(geom) geom
   FROM    emea_eu_seu.albania_prefectures_2021
   UNION ALL
   SELECT st_boundary(geom) geom
   FROM    emea_eu_seu.albania_municipalities_2021
   UNION ALL
   SELECT st_boundary(geom) geom
   FROM    emea_eu_seu.albania_2_digit_postcode_areas_2021
   ),
   b_table AS (
   SELECT   geom
   FROM     a_table
   GROUP BY geom
   ),
   c_table AS (
   SELECT   ST_Union(geom) geom
   FROM     b_table
   ),
   d_table AS (
   SELECT ST_Polygonize(geom) AS geom 
   FROM c_table
   ),
   e_table AS (
   SELECT st_dump(geom) geom
   FROM   d_table
   )
   SELECT  COUNT(*) FROM e_table
   ```
   
   ### Sedona query - returns 1 feature
   
   %sql
   
   ```
   WITH a_table AS (
   SELECT st_boundary(st_geomfromwkb(geom_wkb)) geom
   FROM   hive_metastore.gfk_bronze.albania_prefectures_2021
   UNION ALL
   SELECT st_boundary(st_geomfromwkb(geom_wkb)) geom
   FROM   hive_metastore.gfk_bronze.albania_municipalities_2021
   UNION ALL
   SELECT st_boundary(st_geomfromwkb(geom_wkb)) geom
   FROM   hive_metastore.gfk_bronze.albania_2_digit_postcode_areas_2021
   ),
   b_table AS (
   SELECT   geom
   FROM     a_table
   GROUP BY geom
   ),
   c_table AS (
   SELECT   ST_Union_Aggr(geom) geom
   FROM     b_table
   ),
   d_table AS (
   SELECT ST_Polygonize(geom) AS geom 
   FROM c_table
   ),
   e_table AS (
   SELECT st_dump(geom) geom
   FROM   d_table
   )
   SELECT  COUNT(*) FROM e_table
   ```


-- 
This is an automated message from the Apache Git Service.
To respond to the message, please log on to GitHub and use the
URL above to go to the specific comment.

To unsubscribe, e-mail: dev-unsubscr...@sedona.apache.org.apache.org

For queries about this service, please contact Infrastructure at:
us...@infra.apache.org

Reply via email to