1

Bom dia a Todos. Estou fazendo um código e me deparei com um problema.

Estou usando o GeoPandas e preciso encontrar a intersecção entre duas Linhas(Linestring Z) e os respectivos valores de elevação (Z).

o que tem acontecido é que ao executar a função intersection esta me retorna corretamento o X e o Y, mas não o Z para cada linha.

Exemplo: Intersecção na linha_Laranja = (2,1,3) Intersecção na linha_Verde = (2,1,4)

Sendo assim, X e Y seriam iguais, mas as elevações não (Z1=3 e Z2=4)

Código usado e resultado

linha_Laranja.intersection.(linha_Verde).z

Resultado = 3

linha_Verde.intersection.(linha_Laranja).z

Resultado = 3 (deveria ser 4)

Abaixo ilustrei o problema.

inserir a descrição da imagem aqui

inserir a descrição da imagem aqui

Como os pontos das minhas linhas não são iguais, este XYZ da intersecção deve ser provenientes de uma interpolação.

Abaixo segue o código que estou utilizando:

import numpy as np
import pandas as pd
import geopandas as gpd
from shapely.geometry import LineString, Point


def load_xyz(file_path, file_name, spatial_reference, sep=';', header=0):
    df = pd.read_csv(file_path+'\\'+file_name, header=header, sep=sep,
                     names=['Line','X','Y','Z'])

    gdf = gpd.GeoDataFrame(df.values,
                           crs=spatial_reference,
                           geometry=[Point(xyz) for xyz in zip(df['X'], df['Y'], df['Z'])],
                           columns=['Line','X','Y','Z'])

    gdf = gdf.groupby(['Line'])['geometry'].apply(lambda x: LineString(x.tolist()))

    return gdf, df



spatial_reference = {'init': 'epsg:32724'}

file_path = r'C:\Users\VN\Desktop'
file_name = r'xyz.csv'
df_lines, df_lines_raw = load_xyz(file_path, file_name, spatial_reference)

file_path = r'C:\Users\VN\Desktop'
file_name = r'xyz_cross.csv'
df_cross, df_cross_raw = load_xyz(file_path, file_name, spatial_reference)


cross_validation = gpd.GeoDataFrame([],
                                    columns=['Line','Cross_Line','X','Y','Z1','Z2'],
                                    crs=spatial_reference)


lines_names = df_lines.index.tolist()
cross_names  = df_cross.index.tolist()
idx = 0


for i in range(len(df_lines)):
    for j in range(len(df_cross)):
        if(df_lines[i].crosses(df_cross[j])):
            cross_validation.loc[idx, 'Line'] = lines_names[i]
            cross_validation.loc[idx, 'Cross_Line'] = cross_names[j]
            cross_validation.loc[idx, 'X'] = df_lines[i].intersection(df_cross[j]).x
            cross_validation.loc[idx, 'Y'] = df_lines[i].intersection(df_cross[j]).y
            cross_validation.loc[idx, 'Z1'] = df_lines[i].intersection(df_cross[j]).z
            cross_validation.loc[idx, 'Z2'] = df_cross[j].intersection(df_lines[i]).z
            idx = idx + 1


print(cross_validation)

Obrigado!

1
  • Você colocou todo o código, mas ele depende dos dados que estão no arquivoCSV. sem saber como estão esses dados é muito mais difícil pra alguém reproduzir o seu problema e conseguir te ajudar. Se tiver todos os elementos para as pessoas poderem reproduzir exatamente o que está acontecendo é a única forma de responder a pergutna, com exceção do reduzidíssimo número de pessoas que já estão familiarizados com GeoPandas - esses são os únicos que poderiam bater o olho e dizer o que tem de errado sem poder reproduzir o problema localmente.
    – jsbueno
    2/06/2019 às 15:29

1 Resposta 1

1

Tente esta função:

   def exists_intersection(line1, line2):

    def ccw(A,B,C):
        return (C[1]-A[1])*(B[0]-A[0]) > (B[1]-A[1])*(C[0]-A[0])

    def intersect(line1, line2):
        A,B = line1
        C,D = line2
        bol1 = ccw(A,C,D) != ccw(B,C,D) and ccw(A,B,C) != ccw(A,B,D)
        bol2 = (A == C or A == D) or (B == C or B == D)
        return bol1 or bol2

    resultado = []

    x1, y1 = line1[0]
    x2, y2 = line1[1]

    x3, y3 = line2[0]
    x4, y4 = line2[1]

    px_a, px_b = ((x1 * y2 - y1 * x2) * (x3 - x4) - (x1 - x2) * (x3 * y4 - y3 * x4)), ((x1 - x2) * (y3 - y4) - (y1 - y2) * (x3 - x4))
    py_a, py_b = ((x1 * y2 - y1 * x2) * (y3 - y4) - (y1 - y2) * (x3 * y4 - y3 * x4)), ((x1 - x2) * (y3 - y4) - (y1 - y2) * (x3 - x4))

    if ([px_a, px_b] != [0, 0] and [px_b, py_b] != [0, 0]):
        if intersect(line1, line2):
            px, py = px_a/px_b , py_a/py_b
            resultado = [px, py]

    return resultado

Eu usei para linhas retas, não tenho certeza do resultado em curvas, mas talvez te ajude.

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .