Cari tahu titik yang berada di antara dua garis paralel

Saya menghadapi satu masalah di ArcGIS. Saya bekerja pada basis data navigasi. Dalam basis data kami, jalan jalur tunggal diwakili oleh satu Jalur, sedangkan jalan multi-jalur (jalan dengan pembagi di tengah) diwakili oleh dua garis paralel (garis berwarna merah pada gambar).

Saya memiliki titik shapefile dengan beberapa titik jatuh di dalam jalan multi-Lane dan beberapa di luar.

Saya ingin membuat skrip ArcPy yang akan menemukan titik-titik yang berada di dalam Multi-Lane Streets. yaitu antara garis-garis paralel ini (ditandai dalam gambar).

Saya tidak tahu bagaimana mencapai ini, Adakah yang bisa membantu saya?

Saya melakukan beberapa latihan di atasnya dan saya menemukan bahwa membuat buffer di satu sisi garis dapat membuat di dalam Multi-Lane polygon (Ditampilkan dalam Gambar).

tetapi sekarang masalahnya adalah, poligon sebenarnya melintasi garis (yaitu tumpang tindih batas multi-jalur). sehingga akan menangkap poin yang tidak perlu. apakah ada cara untuk menyelaraskan poligon ini ke garis jalan?

Catatan: integrasi tidak akan berfungsi di sini, karena ia juga memindahkan garis jalan. saya hanya perlu menyelaraskan poligon ke garis jalan.

Saya akan mencoba di bawah ini algoritma arcpy (bahkan manual!)

1. Temukan lebar jalan dua lajur yang tepat - di sini Anda mungkin perlu mengelompokkan jalan dengan lebar yang sama dan ikuti prosedur di bawah ini untuk setiap kluster.
2. Buat penyangga kedua garis ke arah kedua arah (kanan dan kiri) dengan lebar itu (atau sedikit kurang dari itu- untuk memastikan area jalan).
3. Jalankan alat titik-temu untuk mendapatkan wilayah yang tumpang tindih.
4. Jalankan Pilih berdasarkan lokasi untuk memilih titik yang termasuk dalam poligon ini.

Saya akan mengatakan ini adalah latihan geometris.

KODE PSEUDO:

• Untuk setiap titik (titik hitam) temukan jalan terdekat dan temukan proyeksi titik pada jalan ini (titik merah).
• Gambar garis pendek (putus-putus) dalam arah yang berlawanan mulai dari titik hitam
• Cari jika ada persimpangan antara garis pendek dan nama jalan yang sama, bintang biru. Jika ada satu, titik hitam adalah yang kita kejar.

Seperti yang dapat dilihat ada kasus khusus - titik hitam dilingkari:

1. Jalan 1 baris sangat berkelok-kelok. Ini dapat dihilangkan dengan a) bekerja dengan 2 jalur jalan saja atau b) memastikan FID jalan yang memotong titik merah dan bintang berbeda. Namun jika jalan bendy memiliki persimpangan dengan jalan 1 jalur lainnya, ini mungkin tidak berfungsi.
2. Titik hitam terletak pada perpanjangan jalan 1 garis persis tegak lurus. Dalam hal ini ada kemungkinan jalan 1 jalur dapat dipilih sebagai tetangga terdekat.
3. Titik hitam berada di garis.

Semua kasus di atas sangat tidak mungkin, namun tampaknya opsi teraman adalah bekerja dengan 2 jalan saja, yaitu mengekspornya ke kelas fitur yang terpisah. Kasus 3 adalah yang lucu, kita akan membiarkannya kebetulan, karena jarak terpendek ke garis tidak pernah benar nol, dengan demikian arah sinar 'berlawanan' yang menghubungkan 2 titik dapat ditemukan.

Implementasi python:

import arcpy, traceback, os, sys
from arcpy import env
env.overwriteoutput=True

# things to change ---------
maxD=30
mxd = arcpy.mapping.MapDocument("CURRENT")
pointLR = arcpy.mapping.ListLayers(mxd,"NODES")[0]
lineLR = arcpy.mapping.ListLayers(mxd,"LINKS")[0]
sjOneToMany=r'D:\scratch\sj2.shp'
RDNAME='street'
# -------------------------
dDest=arcpy.Describe(lineLR)
SR=dDest.spatialReference

try:
    def showPyMessage():
        arcpy.AddMessage(str(time.ctime()) + " - " + message)
    g = arcpy.Geometry()
    geometryList=arcpy.CopyFeatures_management(pointLR,g)
    n=len(geometryList)
    endPoint=arcpy.Point()

    arcpy.SpatialJoin_analysis(pointLR, lineLR,sjOneToMany,"JOIN_ONE_TO_MANY","KEEP_COMMON","","WITHIN_A_DISTANCE",maxD)
    initFidList=(-1,)
    for fid in range(n):
        query='"TARGET_FID" = %s' %str(fid)
        nearTable=arcpy.da.TableToNumPyArray(sjOneToMany,("TARGET_FID","JOIN_FID"),query)
        if len(nearTable)<2:continue
        fidLines=[int(row[1]) for row in nearTable]
        query='"FID" in %s' %str(tuple(fidLines))
        listOfLines={}
        blackPoint=geometryList[fid]
        with arcpy.da.SearchCursor(lineLR,("FID", "Shape@","STREET"),query) as rows:
            dMin=100000
            for row in rows:
                shp=row[1];dCur=blackPoint.distanceTo(shp)
                listOfLines[row[0]]=row[-2:]
                if dCur<dMin:
                    fidNear,lineNear, roadNear=row
                    dMin=dCur
            chainage=lineNear.measureOnLine(blackPoint)
            redPoint=lineNear.positionAlongLine (chainage).firstPoint
            smallD=blackPoint.distanceTo(redPoint)
            fp=blackPoint.firstPoint
            dX=(redPoint.X-fp.X)*(maxD-smallD)/smallD
            dY=(redPoint.Y-fp.Y)*(maxD-smallD)/smallD
            endPoint.X=fp.X-dX;endPoint.Y=fp.Y-dY
            dashLine=arcpy.Polyline(arcpy.Array([fp,endPoint]),SR)

            for n in listOfLines:
                if n==fidNear:continue
                line, road=listOfLines[n]
                if road!=roadNear:continue
                blueStars=dashLine.intersect(line,1)
                if blueStars.partCount==0:continue
                initFidList+=(fid,); break
    query='"FID" in %s' %str(initFidList)
    arcpy.SelectLayerByAttribute_management(pointLR, "NEW_SELECTION", query)
    arcpy.AddMessage ('\n %i point(s) found' %(len(initFidList)-1))
except:
    message = "\n*** PYTHON ERRORS *** "; showPyMessage()
    message = "Python Traceback Info: " + traceback.format_tb(sys.exc_info()[2])[0]; showPyMessage()
    message = "Python Error Info: " +  str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"; showPyMessage()            

Ada solusi lain yang mungkin lebih elegan. Ini melibatkan triangulasi. Beri tahu saya jika ini menarik dan saya akan memperbarui jawaban saya

Karena jalanannya paralel, saya berasumsi mereka dibuat dengan Copy Parallelalat di bilah alat Edit sehingga membuat sepasang garis memiliki arah yang sama. Kita kemudian dapat beralih di atas koordinat baris pertama dan menambahkannya ke poligon dan kemudian beralih ke kebalikan dari baris kedua. Pasti ada cara yang lebih baik untuk mendekati meraih pasangan garis; pendekatan OID berfungsi, tetapi tidak terlalu cantik.

import collections
import arcpy

FC = "fc"
points = "points"
pgons = "pgons"
arcpy.env.overwriteOutput = True

def buildpoly(oid_coords):
    #create ddict of the form OID:<x1y1, x2y2, ..., xn-1yn-1, xnyn>
    ddict = collections.defaultdict(list)    
    for k,v in oid_coords:
        ddict[k].append(v)

    line1,line2 = ddict.keys()    

    #Assume that the parallel lines have same direction, so reverse the second
    arr = arcpy.Array()
    arr.extend(arcpy.Point(*pt) for pt in ddict[line1])    
    arr.extend(arcpy.Point(*pt) for pt in ddict[line2][::-1])

    return arcpy.Polygon(arr)

#id is an integer field that pairs parallel lines together
unique = list(set(t[0] for t in arcpy.da.SearchCursor(FC, "id")))
polygons = []
for uni in unique:
    polygons.append(buildpoly([r for r in row] for row in arcpy.da.SearchCursor(FC,
                                                                                ["OID@", "SHAPE@XY"],
                                                                                "id={}".format(uni),
                                                                                explode_to_points=True)))


arcpy.CopyFeatures_management(polygons, pgons)

Dari sana, itu panggilan ke Intersect / Select Layer dengan lokasi / apa pun. Perhatikan bahwa Spoligon yang berbentuk tidak sempurna karena saya menggambar secara bebas dan ada beberapa busur yang explode_to_pointstidak menanganinya dengan baik. Jalankan saja Densifyatau sederajat.