python - 如何使用 python 查找 ESRI shapefile 中特定纬度/经度的信息？

Question

我有一个 ESRI shapefile(来自这里:http://pubs.usgs.gov/ds/425/)。我希望使用 python 在给定的纬度/经度下从形状文件(在本例中为表面 Material )中查找信息。

解决这个问题的最佳方法是什么？

谢谢。

最终解决方案:

#!/usr/bin/python

from osgeo import ogr, osr

dataset = ogr.Open('./USGS_DS_425_SHAPES/Surficial_materials.shp')
layer = dataset.GetLayerByIndex(0)
layer.ResetReading()

# Location for New Orleans: 29.98 N, -90.25 E
point = ogr.CreateGeometryFromWkt("POINT(-90.25 29.98)")

# Transform the point into the specified coordinate system from WGS84
spatialRef = osr.SpatialReference()
spatialRef.ImportFromEPSG(4326)
coordTransform = osr.CoordinateTransformation(
        spatialRef, layer.GetSpatialRef())

point.Transform(coordTransform)

for feature in layer:
    if feature.GetGeometryRef().Contains(point):
        break

for i in range(feature.GetFieldCount()):
    print feature.GetField(i)

Answer 1

检查 Python Shapefile Library

这应该给你几何和不同的信息。