점 (x, y)가 ArcMap으로 작성된 다중 다각형의 경계 (및 shapefile 형식) 내에 있는지 확인하기 위해 shapefile, fiona 및 ogr과 같은 라이브러리를 사용하여 몇 가지 코드 예제를 시도했습니다. 그러나 멀티 폴리곤에서는 잘 작동하지 않지만 일반적인 단일 폴리곤 쉐이프 파일에서는 잘 작동합니다. 내가 시도한 일부 스 니펫은 다음과 같습니다.
# First example using shapefile and shapely:
from shapely.geometry import Polygon, Point, MultiPolygon
import shapefile
polygon = shapefile.Reader('shapefile.shp')
polygon = polygon.shapes()
shpfilePoints = []
for shape in polygon:
shpfilePoints = shape.points
polygon = shpfilePoints
poly = Polygon(poly)
point = Point(x, y)
# point in polygon test
if polygon.contains(point):
print 'inside'
else:
print 'OUT'
# Second example using ogr and shapely:
from shapely.geometry import Polygon, Point, MultiPolygon
from osgeo import ogr, gdal
driver = ogr.GetDriverByName('ESRI Shapefile')
dataset = driver.Open("shapefile.shp", 0)
layer = dataset.GetLayer()
for index in xrange(layer.GetFeatureCount()):
feature = layer.GetFeature(index)
geometry = feature.GetGeometryRef()
polygon = Polygon(geometry)
print 'polygon points =', polygon # this prints 'multipoint' + all the points fine
point = Point(x, y)
# point in polygon test
if polygon.contains(point):
print 'inside'
else:
print 'OUT'
첫 번째 예제는 한 번에 하나의 다각형으로 잘 작동하지만 다중 다각형 모양 파일의 모양 중 하나에 점을 입력하면 여러 부분 중 하나에 속하더라도 "out"을 반환합니다.
두 번째 예제에서는 지오메트리 필드를 색인화 된 일반 목록 / 배열로 읽을 수 없기 때문에 " 'Geometry'유형의 객체에 len ()이 없습니다"라는 오류가 발생합니다.
추가로 작동하지 않는 코드의 일부가 아닌지 확인하기 위해 여기 에 제안 된대로 다각형 코드의 실제 지점을 바꾸려고 했습니다. 링크의 예제는 간단한 다각형 모양 파일로 잘 작동하지만 복잡한 다중 다각형을 올바르게 테스트 할 수는 없습니다.
그래서 파이썬을 통해 점이 다중 다각형 모양 파일 내에 있는지 여부를 테스트하는 다른 방법을 생각할 수 없습니다 ... 아마도 내가 누락 된 다른 라이브러리가 있습니까?
두 번째 예는 다중 다각형을 다각형으로 강제 변환하는 것 같습니다. 다중 다각형의 첫 번째 부분에 대해서만 점을 확인하는 것일 수 있습니다. 포인트를 다른 부분으로 옮기고 검사가 성공했는지 확인하십시오.
—
obrl_soil
@obrl_soil 제안 해 주셔서 감사합니다. 그러나 두 번째 예는 위에서 설명한 오류 메시지 ( 'Geometry'유형의 객체에는 len ()이 없습니다) "때문에 작동하지 않습니다." . 첫 번째 예와는 메인 다각형 일 이내에 그 희망이 설명이 도움이됩니다.
—
spartmar
예, 오류가 발생하면
—
obrl_soil
polygon = Polygon(geometry)스위치가있는 try 루프 로 교체해야한다고 생각 polygon = MultiPolygon(geometry)합니다.
첫 번째 예제의 문제는 첫 번째 루프에 있습니다.
—
xunilk
