NDVI 数据处理，及估算植被覆盖度（二）arcgis 栅格转ascii

公告：最近忙其他事情，暂时停止更新。

经过上一步的处理，得到了11年*12月 共132各tif图像。

接下来把tif图像转成我们习惯使用的文本格式，即经度，纬度，ndvi。共三列。当然也可以根据个人偏好制作为nc文件等（制作nc文件可参见我的博客）。

第一步：把tif文件转为ascii格式的文本文件。所谓ascii格式的文本文件可参见博客：https://blog.csdn.net/kikitamoon/article/details/45174145 。如下：

一般在 ASCII 栅格的文件头中我们会看到这几行：

ncols 和 nrows 表示这份数据的行列数，这份数据是10×10。

xllcorner 和 yllcorner 表示的就是栅格的左下角（low left）角点的坐标。

cellsize 也就是像元大小，表示每个正方形单元的尺寸。

通过这几个值不就确定了栅格数据的范围了吗？嗯，确定了。

 

看图说话：

 

文本文件下面的值就对应到每个格子里的像元值：

 

那么如何把这些tif批量转为这种格式？若数据不多则建议手动即可。

同上篇博客的使用方法一样，打开arcgis，粘贴上下面的代码，把目录改成自己的，敲回车即可。代码的意思就是把这些tif全部转换成这样的ascii格式。

import arcpy
import os
raster_path = 'F:/ZPATR-TIME-JOB/ndvi-china-/xinjiang/'   # tif 文件的目录
out_path = 'F:/ZPATR-TIME-JOB/ndvi-china-/xinjiang-ascii/'  # 输出目录

files = [f for f in os.listdir(raster_path) if f.endswith('TIF')]
for file in files:
    arcpy.RasterToASCII_conversion(raster_path+file, out_path+file+'.txt')

运行完成后，则得到一堆这样的ascii文件。后续处理请见下篇。