Desenhando o limite do conjunto de pontos usando o QGIS?

Eu tenho um conjunto de pontos que gostaria de transformar em um polígono.

Tentei o casco côncavo, mas não obtive o polígono pretendido, como você pode ver abaixo,

E quando uso a triangulação de Denaulay, preciso excluir os triângulos externos.

Existe uma maneira eficiente de fazer o mesmo?

Eu também tentei o casco convexo.

Presumo que você saiba quais pontos pertencem à "esquerda" ou "direita", porque, caso contrário, existem muitas soluções. Nesse caso, você pode usar a triangulação delaunay seguida de uma seleção dos triângulos "centrais". Com esse método, não é necessário ter os pontos em uma ordem específica. Tudo o que você precisa considerar é o fato de que os triângulos necessários devem tocar pelo menos um ponto de cada lado.

abrir tabela de atributos> selecionar por expressão> "código" = 1

vetor> ferramenta de pesquisa> selecione por local (nova seleção, interseção)

abrir tabela de atributos> inverter a seleção dos pontos (ctrl + R)

vetor> ferramenta de pesquisa> selecionar por local (remover da seleção, separar)

vetor> ferramenta de geoprocessamento> dissolver

Eu tenho o mesmo resultado que @radouxju.

Operações :

1. Selecione 1 polilinha [selecione ao clicar]
2. Atualize sua tabela de polilinha e afeta um número para cada recurso [calculadora de campo]
3. Selecione os pontos que cruzam uma polilinha [ferramenta de consulta espacial]
4. Atualize os pontos selecionados pelo valor da polilinha afetado anteriormente [calculadora de campo]

5. Mantenha o ponto de seleção e selecione os triângulos delaunay que cruzam esses pontos selecionados [ferramenta de consulta espacial]
   
   

6. Inverta a seleção de pontos ou selecione os pontos que cruzam a outra polilinha [seleção por expressão] ou [inverter seleção]

7. Remova da seleção atual os triângulos delaunay que são disjuntos com os outros pontos [ferramenta de consulta espacial]
   
   

8. Mesclar os recursos pela seleção [ferramenta de geoprocessamento> dissolver]
   
   

IMPORTANTE : DEVE ter as polilinhas para obter este resultado!

Se você tem experiência com python, pode usar a biblioteca Shapely e criar um polígono a partir dos pontos nas duas linhas. Você precisará dizer ao python quais são os pontos inicial e final de ambas as linhas.

from shapely.geometry import Point, Polygon, LineString
import geopandas as gpd
import pandas as pd
line1 = [(1,1),(2,1.2),(3,1)]
line2 = [(1,2),(2,2.2),(3,2)]
# you need to reverse the order of one line to make it a polygon
line2reverse = list(reversed(line2))
polgonList2 = line1 + line2reverse
Polygon(polgonList2)

Ainda melhor: você também pode usar geopandas para fazer isso. geopandas permite salvar facilmente em vários formatos, incluindo shapefiles

d = {'identifier' : [1, 2],
 'name' : ["Netherlands", "Germany"],
 "line1": [[(1,1),(2,1.2),(3,1)], [(1,1),(2,1.2),(3,1)]],
 "line2": [[(1.1,2.1),(2.1,2.3),(3.1,2.2)],[(1,2),(2,2.2),(3,2)]]
}

df = pd.DataFrame(d)
def makePolygon(row):
    line2reverse = list(reversed(row["line2"]))
    return Polygon(line1+line2reverse)    
geometries = []

for index, row in df.iterrows():
    geometries.append(makePolygon(row))
crs = {'init': 'epsg:4326'} 
gdf = gpd.GeoDataFrame(df, crs=crs, geometry=geometries)
gdf.to_file('MyGeometries.shp', driver='ESRI Shapefile')

você pode ler as geometrias da linha usando a função geopandas gpd.read_file ().

Ordem dos vértices poligonais em geral GIS: no sentido horário ou anti-horário

https://nbviewer.jupyter.org/gist/rutgerhofste/b01c17aa6851ea577f10c21a4c3717bc