处理高分影像时,因为原数据没有存储投影信息,所以需要自己定义。使用GDAL(osr)对遥感影像定义投影。使用以下代码

from osgeo import gdal, osr
’‘’
其他读取影像数据的程序
‘’‘
# 获取地理坐标系统信息,用于选取需要的地理坐标系统
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)  # 定义输出的坐标系为"WGS 84"
out_IRS.SetProjection(srs.ExportToWkt())  # 给新建图层赋予投影信息

出现问题:

ERROR 1: PROJ: proj_create_from_database: D:\Anaconda\pkgs\proj-6.2.1-h9f7ef89_0\Library\share\proj\proj.db lacks DATABASE.LAYOUT.VERSION.MAJOR / DATABASE.LAYOUT.VERSION.MINOR metadata. It comes from another PROJ installation.

这个报错是说,使用的proj.db是由另一个程序创建的。
原因分析
使用osgeo.osr能够直接从安装的GDAL包读取指定投影类型的信息,在proj.db中存储了常用的投影系统的参数,因此指定EPSG号就能从proj.db中取到投影信息。
经过排查,我的情况是在anaconda的虚拟环境中同时安装了GDAL和其他地理空间库,此时同一环境路径下可能会产生多个proj.db文件,访问它的时候程序会无法判断访问哪个。因此只需要去查找gdal安装的那个proj.db在哪里,在程序中指定环境变量即可解决
资源管理器中虚拟环境的路径下查找proj.db,找到两个结果:
在这里插入图片描述
显然,应该使用osgeo目录下的那个。

解决方法
在程序开头指定环境变量,指定上一步查找的路径即可。

os.environ['PROJ_LIB'] = r'D:\Anaconda\envs\geoenv\Lib\site-packages\osgeo\data\proj'
Logo

为开发者提供学习成长、分享交流、生态实践、资源工具等服务,帮助开发者快速成长。

更多推荐