¿Cómo dividir las características de ST_Split en una tabla por las características en otra?

Question

Necesito dividir los polígonos (capa 'pol') por líneas cerradas y no cerradas (capa 'lin').

Desafortunadamente no obtengo resultados adecuados al realizar la siguiente consulta.

CREATE VIEW splitted_pol AS
SELECT 
    g.path[1] as gid, 
    g.geom::geometry(polygon, SRID) as geom 
FROM
    (SELECT 
    (ST_Dump(ST_Split(pol.geom, lin.geom))).* 
    FROM pol, lin
) as g;

En mi ejemplo ST_Split debería crear seis polígonos (capa 'splitted_pol').

¿Alguien sabe cómo usar ST_Split con QGIS/PostGIS?

Answer 1

Podrías crear una función como ésta:

create or replace function ST_MultiSplit(geom Geometry, blades Geometry) RETURNS Geometry AS $$
BEGIN
  FOR i IN 1..ST_NumGeometries(blades)
     LOOP
        geom = ST_Split(geom, ST_GeometryN(blades, i));
     END LOOP;
  RETURN geom;
END;
$$ LANGUAGE plpgsql;

Entonces úsalo como:

Select ST_AsEWKT(a.geom) from (
   select (ST_Dump(ST_MultiSplit(pol.geom, (
      select ST_MemUnion(lin.geom) from lin where ST_Intersects(pol.geom, lin.geom) = 't')
 ))).geom geom from pol) a;

Así se obtienen los seis registros que se esperan. Es posible que quieras añadir alguna comprobación/manejo de errores y no estoy seguro de la escalabilidad.

Answer 2

Estoy usando postGIS sql para dividir característica por línea en JAVA, y mi código funcionó. mi código es:

public List splitGeometry(String geom1, String geom2) {

    List<String> result=new ArrayList<String>() ;
    try {

        Statement s = connection.createStatement();
        String sql_stat = null;
        sql_stat = "select st_astext (a.geom)from (select ( st_dump(p_geom)).geom as geom FROM (SELECT  st_split(ST_GeomFromEWKT('"+geom1+"') ,ST_GeomFromEWKT('"+geom2+"'))AS p_geom) AS b) AS a;";
        System.out.println(sql_stat);
        ResultSet rs=  s.executeQuery(sql_stat);
        while (rs.next()){
            result.add(rs.getString(1)) ;
        }
    } catch (Exception e) {
        e.printStackTrace();
    }
    return result ;
}

Espero que te ayude.