¿Dibujando el límite del conjunto de puntos usando QGIS?

Tengo un conjunto de puntos que me gustaría cambiar a un polígono.

Intenté el casco cóncavo pero no obtengo el polígono previsto como puedes ver a continuación,

Y cuando uso la triangulación Denaulay, necesito eliminar los triángulos exteriores.

¿Hay una manera eficiente de hacer lo mismo?

También probé el casco convexo.

Supongo que sabe qué puntos pertenecen a la "izquierda" o "derecha", porque de lo contrario hay muchas soluciones. Si es así, puede utilizar la triangulación delaunay seguida de una selección de los triángulos "centrales". Con este método, no es necesario tener los puntos en un orden específico. Todo lo que necesita considerar es el hecho de que los triángulos que necesita deben tocar al menos un punto de cada lado.

abrir tabla de atributos> seleccionar por expresión> "código" = 1

vector> herramienta de investigación> seleccionar por ubicación (nueva selección, intersección)

abrir tabla de atributos> invertir la selección de los puntos (ctrl + R)

vector> herramienta de investigación> seleccionar por ubicación (eliminar de la selección, disjunto)

vector> herramienta de geoprocesamiento> disolver

Tengo el mismo resultado que @radouxju.

Operaciones:

1. Seleccione 1 polilínea [seleccionar al hacer clic]
2. Actualice su tabla de polilínea y afecte un número a cada entidad [calculadora de campo]
3. Seleccione los puntos que intersecan una polilínea [herramienta de consulta espacial]
4. Actualice los puntos seleccionados por el valor de polilínea afectado previamente [calculadora de campo]

5. Mantenga el punto de selección y seleccione los triángulos delaunay que intersecan estos puntos seleccionados [herramienta de consulta espacial]
   
   

6. Invierta la selección de puntos o seleccione los puntos que se cruzan con la otra polilínea [selección por expresión] o [selección invertida]

7. Elimine de la selección actual los triángulos delaunay que son disjuntos con los otros puntos [herramienta de consulta espacial]
   
   

8. Combinar las características por selección [herramienta de geoprocesamiento> disolver]
   
   

IMPORTANTE : ¡DEBE tener las polilíneas para obtener este resultado!

Si tiene experiencia con Python, puede usar la biblioteca Shapely y crear un Polígono a partir de los puntos en las dos líneas. Deberá decirle a Python cuáles son los puntos de inicio y finalización de ambas líneas.

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)

Aún mejor: también puedes usar geopandas para hacer esto. geopandas le permite guardar fácilmente en múltiples formatos, incluidos los archivos de forma

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')

Puede leer las geometrías de línea usando la función geopandas gpd.read_file ().

Orden de vértices poligonales en general SIG: en sentido horario o antihorario

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