解决GDAL写投影问题——proj.db报错

处理高分影像时，因为原数据没有存储投影信息，所以需要自己定义。使用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'