I have two questions :
I have some lines (eg: road network). I want to find coordinates of each junction.
Should i use crosses or intersects function from python's Shapely module? They are giving different counts.
Moreover these functions return boolean value. How to get the junction coordinates ?
with open('./lines.geojson') as f: lines = json.load(f) line =  for feature in lines['features']: line.append(feature['geometry']['coordinates']) count = 0 for i in range(len(line)-1): j =i+1 while (j < len(line)): l1 = LineString(line[i]) l2 = LineString(line[j]) if l1.crosses(l2)== True : count =count+1 j = j+1