I have an array of points (in latitude, longitude coordinates) which I would like to buffer and dissolve (into a polygon) and then intersect the Ordnance Survey Open Roads (projected in BNG) with that polygon and export the results.
Note: I am not interested in preserving any of the shape-file information -> I just want to effectively clip the road-network to show a line only within my polygon as efficiently as possible.
I have had some luck producing with the first step (producing my polygon)
path_to_data = 'Raw_Points_Iso_Mined.csv' # Read in pts = [] with open(path_to_data) as f: for x in csv.reader(f): # My data is: latitude, longitude, minutes # I only want those points which are within 10 minutes if float(x[2])
Note: I am not interested in preserving any of the shape-file information -> I just want to effectively clip the road-network to show a line only within my polygon as efficiently as possible.
I have had some luck producing with the first step (producing my polygon)
path_to_data = 'Raw_Points_Iso_Mined.csv' # Read in pts = [] with open(path_to_data) as f: for x in csv.reader(f): # My data is: latitude, longitude, minutes # I only want those points which are within 10 minutes if float(x[2])