1. 程式人生 > >NDVI 資料處理,及估算植被覆蓋度(二)arcgis 柵格轉ascii

NDVI 資料處理,及估算植被覆蓋度(二)arcgis 柵格轉ascii

經過上一步的處理,得到了11年*12月 共132各tif影象。

接下來把tif影象轉成我們習慣使用的文字格式,即經度,緯度,ndvi。共三列。當然也可以根據個人偏好製作為nc檔案等(製作nc檔案可參見我的部落格)。

一般在 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檔案。後續處理請見下篇。