用arcgis和maxent适应性区划中的一些点

2023年8月16日#python #arcgis #maxent #适应性区划 #课题17

AI 生成的摘要

最近的课题使用了ArcGIS和MaxEnt软件。在使用MaxEnt时，需要气象因子，但从WorldClim下载的未来生物气象因子只有一个包含19个波段的tif文件，需要手动拆分。由于找不到用ArcGIS拆分的方法，转而使用Python进行处理。以下是用于拆分tif文件的Python代码： ```python import gdal import os input_file = "G:/4.其他/wc2.1_30s_bioc_ACCESS-CM2_ssp126_2041-2060.tif" output_folder = "G:/4.其他/feature_bio/" if not os.path.exists(output_folder): os.makedirs(output_folder) dataset = gdal.Open(input_file) num_bands = dataset.RasterCount for i in range(num_bands): band = dataset.GetRasterBand(i + 1) output_file = os.path.join(output_folder, f"bio_{i + 1}.tif") driver = gdal.GetDriverByName("GTiff") new_dataset = driver.Create(output_file, dataset.RasterXSize, dataset.RasterYSize, 1, gdal.GDT_Float32) new_dataset.SetProjection(dataset.GetProjection()) new_dataset.SetGeoTransform(dataset.GetGeoTransform()) new_dataset.GetRasterBand(1).WriteArray(band.ReadAsArray()) print(f"Creating {output_file}...") new_dataset.FlushCache() new_dataset = None dataset = None ``` 在安装GDAL时可能会遇到错误，建议使用conda安装。如果安装失败，可以从GitHub下载对应的.whl文件并手动安装。如果无法调用GDAL，可以在Python的site-packages文件夹下新建gdal.py文件，并复制相关代码以解决问题。

最近的课题有用到这两个软件，其中做 maxent 的时候需要气象因子，而在 worldclim 下载的未来生物气象因子只有一张包含 19 个波段的 tif 文件，需要手动拆分成 19 个文件。实在是没有找到用 arcgis 拆分的方法，所以转向了 python

代码如下：

import gdal
import os

input_file = "G:/4.其他/wc2.1_30s_bioc_ACCESS-CM2_ssp126_2041-2060.tif"#下载的tif文件
output_folder = "G:/4.其他/feature_bio/"#输出文件夹

if not os.path.exists(output_folder):
   os.makedirs(output_folder)

dataset = gdal.Open(input_file)
num_bands = dataset.RasterCount

for i in range(num_bands):
   band = dataset.GetRasterBand(i + 1)
   output_file = os.path.join(output_folder, f"bio_{i + 1}.tif")
   driver = gdal.GetDriverByName("GTiff")
   new_dataset = driver.Create(output_file, dataset.RasterXSize, dataset.RasterYSize, 1, gdal.GDT_Float32)
   new_dataset.SetProjection(dataset.GetProjection())
   new_dataset.SetGeoTransform(dataset.GetGeoTransform())
   new_dataset.GetRasterBand(1).WriteArray(band.ReadAsArray())
   print(f"Creating {output_file}...")
   new_dataset.FlushCache()
   new_dataset = None

dataset = None

但是有个问题，安装 gdal 的时候大概率会报错，具体原因我也没查到。网上的推荐办法是使用 conda：

conda install gdal

实际情况可能会出现无法安装或者安装了没法调用，解决办法：

如果安装失败：
到 Github 仓库下载对应的.whl 文件，比如 python 3.11 就下载 GDAL-3.7.1-cp311-cp311-win_amd64.whl
，然后手动安装pip install GDAL-3.7.1-cp311-cp311-win_amd64.whl
如果无法调用：
在 python 安装各种环境包的文件夹下，如：D:Path\To\python3.9.12\Lib\site-packages\文件夹下，新建 gdal.py 文件，将以下代码复制进去：

# import osgeo.gdal as a convenience
from osgeo.gdal import deprecation_warn
deprecation_warn('gdal')
from osgeo.gdal import *