Shapely: 기능 요청: STRtree.nearest 기능 추가

@FuriousRococo님 , 저는 가장 가까운 이웃 검색 기능을 사용하는 데 익숙하지 않지만 대략 올바른 방향으로 가고 있는 것 같습니다. 그러나 https://github.com/libgeos/geos/blob/master/capi/geos_c.h.in#L1824 에서 트리의 항목이 Python 개체이기 때문에 일반 메서드를 사용해야 한다는 표시를 봅니다. . https://github.com/libgeos/geos/blob/master/capi/geos_c.h.in#L1850 을 시도해야 한다고 생각합니다.

조언 @sgillies 에 감사드립니다. 지금 작동하는 것 같습니다.
strtree.py

class STRtree:
    def nearest(self, geom):
        if self._n_geoms == 0:
            return None

        envelope = geom.envelope

        def callback(item1, item2, distance, userdata):
            try:
                geom1 = ctypes.cast(item1, ctypes.py_object).value
                geom2 = ctypes.cast(item2, ctypes.py_object).value
                dist = ctypes.cast(distance, ctypes.POINTER(ctypes.c_double))
                lgeos.GEOSDistance(geom1._geom, geom2._geom, dist)
                return 1
            except:
                return 0 

        item = lgeos.GEOSSTRtree_nearest_generic(self._tree_handle, ctypes.py_object(geom), envelope._geom, \
            lgeos.GEOSDistanceCallback(callback), None)
        geom = ctypes.cast(item, ctypes.py_object).value

        return geom

ctypes_declarations.py

def prototype(lgeos, geos_version):
    if geos_version >= (3, 6, 0):
        lgeos.GEOSDistanceCallback = CFUNCTYPE(c_int, c_void_p, c_void_p, c_void_p, c_void_p)

        lgeos.GEOSSTRtree_nearest_generic.argtypes = [
            c_void_p, py_object, c_void_p, lgeos.GEOSDistanceCallback, py_object]
        lgeos.GEOSSTRtree_nearest_generic.restype = c_void_p

test.py

from shapely.geometry import Point, Polygon
from shapely.strtree import STRtree

a = []
a.append(Polygon([(1,0),(2,0),(2,1),(1,1)]))
a.append(Polygon([(0,2),(1,2),(1,3),(0,3)]))
a.append(Polygon([(-1.5,0),(-2.5,0),(-2.5,-1),(-1.5,-1)]))
a.append(Point(0,0.5))
tree = STRtree(a)
pt = tree.nearest(Point(0,0))
print(pt.wkt)

산출
POINT (0 0.5)

@FuriousRococo 이 새로운 기능으로 풀 리퀘스트를 제출하시겠습니까?

@FuriousRococo 기능을 추가해주셔서 감사합니다. 이 가장 가까운 기능을 이해하려고 시도하는 동안 질문입니다. 한 지점에 대해 가장 가까운 도로를 찾으려고 한다고 가정해 보겠습니다. 따라서 STR 트리는 MultiLineString 도로의 트리입니다. 따라서 이 경우 가장 가까운 도로는 점에서 MultiLineStrings까지의 거리 투영 거리를 사용하여 지정됩니까? 아니면 도로를 나타내는 MultiLineStrings의 중심까지의 가장 작은 거리를 기반으로 합니까? 감사합니다!

테스트를 작성했고 가장 가까운 함수가 MultiLineString의 모든 줄을 확인하여 가장 가까운 줄이 MultiLineString 안에 있는지 확인한다고 생각합니다. 여기 내 스크립트가 있습니다 @asif-rehan
test.py

from shapely.strtree import STRtree
from shapely.geometry import Point, LineString, MultiLineString

tree_list = [LineString([[0,1],[1,1]])]
tree_list.append(MultiLineString([LineString([[0,0],[1,0]]),LineString([[0,3],[1,3]])]))
tree = STRtree(tree_list)

tree.nearest(Point(0.5,0.25)).wkt

산출
MULTILINESTRING ((0 0, 1 0), (0 3, 1 3))