ポリゴン フィーチャを含むレイヤーがあります。各機能には属性と値があります。座標のリストもあり、座標がどのフィーチャ (またはポリゴン) にあるかを知りたいです。
誰かがこれを行う方法について私を案内してもらえますか? 目標を達成するのに役立つ API の機能はありますか、それとも計算幾何学アルゴリズムを使用して自分で行う必要がありますか? 後者の方法は知っていますが、組み込み関数が既にある場合は時間を節約できます。
ありがとう。
ポリゴン フィーチャを含むレイヤーがあります。各機能には属性と値があります。座標のリストもあり、座標がどのフィーチャ (またはポリゴン) にあるかを知りたいです。
誰かがこれを行う方法について私を案内してもらえますか? 目標を達成するのに役立つ API の機能はありますか、それとも計算幾何学アルゴリズムを使用して自分で行う必要がありますか? 後者の方法は知っていますが、組み込み関数が既にある場合は時間を節約できます。
ありがとう。
while provider.nextFeature(feature):
if (feature.geometry().contains(QgsGeometry.fromPoint(QgsPoint(lon, lat)))):
print 'Contained in feature %d' % feature.id()
私は最終的にそれを自分で行うことができました。
import sys
import os
from qgis.core import *
import matplotlib.pyplot as plt
from matplotlib.path import Path
import matplotlib.patches as patches
LATITUDE = 1.29306
LONGITUDE = 103.856
QgsApplication.setPrefixPath("/usr", True)
QgsApplication.initQgis()
layer=QgsVectorLayer("/home/shubham/SMART/zones/mtz1092p.shp", "mtz1092p", "ogr")
if not layer.isValid():
print "Layer failed to load!"
provider = layer.dataProvider()
def findFeatureId(point):
feat = QgsFeature()
allAttrs = provider.attributeIndexes()
provider.select(allAttrs)
while provider.nextFeature(feat):
geom = feat.geometry()
x = geom.asPolygon()
if len(x) == 0:
print "Feature ID %d has no ring" % feat.id()
else:
codes = []
codes.append(Path.MOVETO)
for i in range (0, len(x[0]) - 2):
codes.append(Path.LINETO)
codes.append(Path.CLOSEPOLY)
path = Path(x[0], codes)
if (path.contains_point(point, None, 0.0)):
print "Point contained in feature ID %d" %feat.id()
if __name__ == "__main__":
crsSrc = QgsCoordinateReferenceSystem(4326) # WGS84
crsDest = QgsCoordinateReferenceSystem(3414)# SVY21
xform = QgsCoordinateTransform(crsSrc, crsDest)
pt = xform.transform(QgsPoint(LONGITUDE, LATITUDE))
findFeatureId(pt)