GIS遥感方面的数据处理、arcpy及批处理脚本的咨询、交流

应部分博友的需求，欢迎arcpy及批处理脚本的咨询、交流。博客中文章也有相关的应用案例，如果有特殊定制化需要的可以私信联系。
（1）arcpy功能开发及批处理；
（2）arctool box小工具制作；
（3）栅格、矢量数据批处理；
（4）影像数据下载；
（5）python批处理实现

例如：GDAL读取、保存栅格

 imgarr = read_as_bsq(fint, bands, rows, col)

    #将提取的数组存储为tif格式图像.
    #注意这里未设置地理坐标和仿射变换信息，所以不能用ENVI等软件打开。
    #tif_savedir = r"C:\Users\HP\Desktop\tifFolder\NEP_2019_364.tif"
    datatype = gdal.GDT_Int16   #gdal.GDT_UInt16
    bands, high, width = imgarr.shape
    driver = gdal.GetDriverByName("GTiff")
    datas = driver.Create(tif_savedir, col, rows, bands, datatype)

    #设置地理坐标和仿射变换信息
    # datas.SetGeoTransform(image_geotrans)
    # datas.SetProjection(image_projetion)
    for i in range(bands):
        datas.GetRasterBand(i + 1).WriteArray(imgarr[i])
    del datas