Convertendo MultiLinestring em Linestring com PostGIS?

16

Eu tenho uma seqüência de linhas múltiplas e quero convertê-lo em uma cadeia de linhas para usar algumas funções que funcionam apenas com cadeias de linhas. Eu usei ST_Dump()e obtive uma coleção das minhas geometrias. Se eu os juntar novamente, ainda recebo a cadeia de linhas múltiplas. Devo pegar os Npoints e combiná-los ou algo assim, para terminar com uma cadeia de linhas simples? ST_LineMerge()também não funciona na minha cadeia de linhas múltiplas, ele retorna o mesmo.

meu exemplo

MULTILINESTRING ((- 3,16420835153456 55,9269166007097, -3,164222 55,926918), (- 3,1642070167833 55,9269296196706, -3,16421351659546 55,9268662214904), (- 3,16421351659546 55,9268662214904, -3,16421636372824 55,9268384509897), (- 3,16421636372824 55,9268384509897, -3,16422182573761 55,9267851753802), (- 3,16422182573761 55,9267851753802, -3,16422870102352 55,926718114886 ), (- 3,16422870102352 55,926718114886, -3,16423309121073 55,926675293667), (- 3,16423309121073 55,926675293667, -3,16423565148822 55,9266503211093), (- 3,16423565148822 55,9266503211093, -3,16424103159897 55,9265978443265), (- 3,16424103159897 55,9265978443265, -3,16424680776317 55,9265415044985), (- 3,16424680776317 55,9265415044985, -3,16425267254583 55,9264843002995 ), (- 3.16425267254583 55.9264843002995, -3.16425541048045 55.9264575949012), ((3.16425541048045 55.9264575949012, -3.16426111146586 55.9264019883556),(-3,16426111146586 55,9264019883556, -3,1642667032531 55,9263474469124), (- 3,1642667032531 55,9263474469124, -3,16426957768543 55,9263194101362), (- 3,16426957768543 55,9263194101362, -3,16427488261739 55,9262676666359), (- 3,16427488261739 55,9262676666359, -3,16428009893088 55,9262167875066), (- 3,16428009893088 55,9262167875066, -3,164282741107 55,9261910161221), (-3,1642875546472 55,9261440655823, -3,164282741107 55,9261910161221), (- 3,1642875546472 55,9261440655823, -3,16429466890915 55,9260746741522), (- 3,16429466890915 55,9260746741522, -3,16430092974527 55,9260136069079), (- 3,16430092974527 55,9260136069079, -3,16430822838418 55,9259424170929), (- 3,16430822838418 55,9259424170929, -3,16431547242401 55,925871759829), (-3.16431547242401 55.925871759829, -3.16431448732505 55.9258328901507), (- 3.16431770120536 55.9257327846001, -3.16431547242401 55.925871759829), (- 3.164339 55.925777, -3.16431770120536 55.9257327846001))

postgis sql Antony
fonte
Dê um exemplo dessas cadeias de linhas múltiplas que você deseja converter para que possamos testá-las.
CaptDragon

Respostas:

6

Você também pode usar o ST_SnapToGrid para corrigir seus dados, mas precisará brincar com os parâmetros para acertar.

SELECT ST_AsText(
    ST_LineMerge(
        ST_SnapToGrid(
            ST_GeomFromText('MULTILINESTRING((-3.16420835153456 55.9269166007097,-3.164222 55.926918),(-3.1642070167833 55.9269296196706,-3.16421351659546 55.9268662214904),(-3.16421351659546 55.9268662214904,-3.16421636372824 55.9268384509897),(-3.16421636372824 55.9268384509897,-3.16422182573761 55.9267851753802),(-3.16422182573761 55.9267851753802,-3.16422870102352 55.926718114886),(-3.16422870102352 55.926718114886,-3.16423309121073 55.926675293667),(-3.16423309121073 55.926675293667,-3.16423565148822 55.9266503211093),(-3.16423565148822 55.9266503211093,-3.16424103159897 55.9265978443265),(-3.16424103159897 55.9265978443265,-3.16424680776317 55.9265415044985),(-3.16424680776317 55.9265415044985,-3.16425267254583 55.9264843002995),(-3.16425267254583 55.9264843002995,-3.16425541048045 55.9264575949012),(-3.16425541048045 55.9264575949012,-3.16426111146586 55.9264019883556),(-3.16426111146586 55.9264019883556,-3.1642667032531 55.9263474469124),(-3.1642667032531 55.9263474469124,-3.16426957768543 55.9263194101362),(-3.16426957768543 55.9263194101362,-3.16427488261739 55.9262676666359),(-3.16427488261739 55.9262676666359,-3.16428009893088 55.9262167875066),(-3.16428009893088 55.9262167875066,-3.164282741107 55.9261910161221),(-3.1642875546472 55.9261440655823,-3.164282741107 55.9261910161221),(-3.1642875546472 55.9261440655823,-3.16429466890915 55.9260746741522),(-3.16429466890915 55.9260746741522,-3.16430092974527 55.9260136069079),(-3.16430092974527 55.9260136069079,-3.16430822838418 55.9259424170929),(-3.16430822838418 55.9259424170929,-3.16431547242401 55.925871759829),(-3.16431547242401 55.925871759829,-3.16431448732505 55.9258328901507),(-3.16431770120536 55.9257327846001,-3.16431547242401 55.925871759829),(-3.164339 55.925777,-3.16431770120536 55.9257327846001))'),
            0.001)
        )
    );
CaptDragon
fonte
25

Você tem certeza de que todos os seus MultiLines que deseja converter podem ser convertidos ?

uma geometria simples é aquela que não possui pontos geométricos anômalos, como auto-interseção ou auto-tangência e refere-se principalmente a geometrias 0 ou 1-dimensionais

Caso contrário, ST_LineMergedeve funcionar:

SELECT ST_AsText(ST_LineMerge(ST_GeomFromText('MULTILINESTRING((-29 -27,-30 -29.7,-36 -31,-45 -33),(-45 -33,-46 -32))')));

ATUALIZAR

Portanto, essa é a cadeia de linhas múltiplas que você forneceu. Parece válido a partir daqui:

insira a descrição da imagem aqui

Porém, ao aproximar o zoom para inspeção, você pode ver claramente que isso não pode ser convertido em uma cadeia de linhas válida.

insira a descrição da imagem aqui insira a descrição da imagem aqui

VÁLIDO:

insira a descrição da imagem aqui insira a descrição da imagem aqui

INVÁLIDO:

insira a descrição da imagem aqui insira a descrição da imagem aqui

Mas pegue um subconjunto VÁLIDO de seus pontos e ele funciona muito bem:

SELECT ST_AsText(ST_LineMerge(ST_GeomFromText('MULTILINESTRING((-3.16422182573761 55.9267851753802,-3.16422870102352 55.926718114886),(-3.16422870102352 55.926718114886,-3.16423309121073 55.926675293667),(-3.16423309121073 55.926675293667,-3.16423565148822 55.9266503211093),(-3.16423565148822 55.9266503211093,-3.16424103159897 55.9265978443265),(-3.16424103159897 55.9265978443265,-3.16424680776317 55.9265415044985),(-3.16424680776317 55.9265415044985,-3.16425267254583 55.9264843002995),(-3.16425267254583 55.9264843002995,-3.16425541048045 55.9264575949012),(-3.16425541048045 55.9264575949012,-3.16426111146586 55.9264019883556),(-3.16426111146586 55.9264019883556,-3.1642667032531 55.9263474469124),(-3.1642667032531 55.9263474469124,-3.16426957768543 55.9263194101362),(-3.16426957768543 55.9263194101362,-3.16427488261739 55.9262676666359),(-3.16427488261739 55.9262676666359,-3.16428009893088 55.9262167875066))')));
CaptDragon
fonte
Isso não está realmente correto. LineStrings com auto-interseções ainda são válidas para OGC. A razão pela qual o MultiLineString original não pode ser convertido em uma única LineString é que ele contém 3 conjuntos separados de vértices conectados.
dr_jts 3/09/19
@dr_jts Certo, não é que o LineString seja inválido, mas a entrada seja inválida. você não pode converter linhas que se cruzam. Portanto, uma entrada válida de cadeias de linhas é uma sem interseções.
CaptDragon 03/09/19
0

Considere usar ST_SubDivide para dividir geometrias grandes em outras menores após a liberação 2.2.0.

Eu validei e testei isso no MultiLineString para quebrá-los no LineString.

https://postgis.net/docs/ST_Subdivide.html

rulhaniam
fonte