Question

我跟着一本书然后输入示例代码，但是当我运行它时，它给出了这些错误。我使用Enthought Canopy以及所有必需的包。怎么能解决这个问题呢？我不想使用另一个包，因为我需要使用其他一些步骤ogr。在Enthought Canopy上，我更新了ogr，但没有帮助。

ERROR 6: Unable to load PROJ.4 library (proj.dll), creation of
OGRCoordinateTransformation failed.

这是示例代码：

from    __future__  import  print_function
import ogr
import osr

def open_shape_file(file_path):
    #Open   the shapefile,  get the first   layer   and returns
    #the    ogr datasource.

    datasource=ogr.Open(file_path)
    layer=datasource.GetLayerByIndex(0)
    print ("opening {}".format(file_path))
    print ("Number of feature:{}".format(layer.GetFeatureCount()))
    return datasource
def transform_geometries(datasource,    src_epsg,   dst_epsg):
#Transform  the coordinates of  all geometries  in  the
#first  layer.

#   Part    1
    src_srs =   osr.SpatialReference()
    src_srs.ImportFromEPSG(src_epsg)
    dst_srs =   osr.SpatialReference()
    dst_srs.ImportFromEPSG(dst_epsg)
    transformation  =   osr.CoordinateTransformation(src_srs,   dst_srs)
    layer   =   datasource.GetLayerByIndex(0)
#   Part    2
    geoms   =   []
    layer.ResetReading()
    for feature in  layer:
        geom    =   feature.GetGeometryRef().Clone()
        geom.Transform(transformation)
        geoms.append(geom)
    return  geoms

datasource=open_shape_file("D:/python/python_geospe/exampledata/TM_WORLD_BORDERS/TM_WORLD_BORDERS-0.3.shp")
layer   =   datasource.GetLayerByIndex(0)
feature =   layer.GetFeature(0)
print("Before   transformation:")
print(feature.GetGeometryRef())
transformed_geoms   =   transform_geometries(datasource,    4326,   3395)
print("After    transformation:")
print(transformed_geoms[0])

open_shape_file("D:/python/python_geospe/exampledata/TM_WORLD_BORDERS/TM_WORLD_BORDERS-0.3.shp")

Answer 1

您是否正确设置了环境变量？ proj.dll通常位于C:\Program Files (x86)\GDAL。您需要使用此路径设置环境变量。

我建议你按照this installation guide解释在Windows操作系统中正确安装GDAL / OGR的过程。

另一个指南：here。

无法加载PROJ.4库

1 个答案: