根據你的資料型別選用不同的資料生成方法
若你的資料是
raster
格式:則有如下方法
1 erdas
imagine----interpreter---gis analysis---matrix
,輸入兩個時相的
raster
資料即可
做這一步之前記得先對兩時相的資料進行重編碼(
nterpreter---gis analysis---recode
)一般執行如果出現錯誤
肯定是重編碼沒做好,請繼續查證。2先在
erdas
中利用modeler
計算如下公式
nc(i,j)=nc(i)*10+nc(j),(j>i)
其中:nc(i,j)表示i
,j兩年份的土地利用變化圖;nc(
i)表示
i年份遙感分類影像;nc(
j)表示
j年份的遙感分類影像。
在此計算的基礎上,將以上變化影像圖轉化為
bil格式,再利用
arc/info grid
模組將影像轉為
grid
格式,然後利用
grid
模組中的屬性表(
vat)檢視命令對影像灰度值進行統計,最後得出土地利用轉化舉證。(注
:此方法本人尚未實現過,不知可行否)。
若資料是
vector
格式1 erdas
imagine----interpreter---gis analysis---matrix
,輸入兩個時相的
vector
資料即可
此時注意
輸出柵格大小不應設的太小
要不一執行就會提示你的空間不足
做這一步之前,請做好前期的地理編碼。
2 arcview3.3
載入spatial analysis
模組把兩時相的
vector
圖轉成grid
格式(當然中間有一些單位的設定根據你做的圖的解析度來設定即可)
analysis---mapcaculate
直接計算即可。
3把兩期解譯完的
vector
檔案在arctoolbox
——overlay
——union
中疊加,注意:兩個檔案不能用同乙個欄位名,比如乙個用
93type
,另乙個時相則用
00type
疊加後的檔案在
arcmap
中開啟,選中檔案,然後點右鍵——
property
——空間查詢,輸入條件語句,比如:
93type=『1
』and 00type=『2
』;查詢結果即為第一種型別轉化為第二種型別的圖形,可以另建一圖層比如:
12,把查詢結果複製到
12圖層上。統計出面積,依進行,就可以得到土地利用型別轉移矩陣
。
基於ENVI下的土地利用資訊提取(一)
基於envi下的土地利用資訊提取 1.前言 對於土地利用資訊的提取的研究,國內外各類文獻資料中介紹的方法十分繁多,可大致分為 基於監督分類法 基於非監督分類法 基於決策樹分類法 基於物件導向分類法 以及基於人工神經網路分類法5種。對於不同區域以及不同時相的資料,選取的方法不盡相同。提取土地利用資訊的...
土地利用覆被變化的概念 欠壓實概念極其荒謬
李傳亮 地表疏鬆的沉積物,孔隙度極高,有時可高達80 以上,但隨著埋藏深度的加大,在上覆壓力的作用下,其排列方式不斷由鬆散排列趨於緻密排列,同時伴隨有孔隙度的大幅度減小,此即岩石的壓實階段 i 當埋深到達成岩深度之後,岩石完成了排列方式的轉變,即完全達到了緻密排列,再加上膠結作用,岩石的排列方式基本...
定義矩陣的幾種方法
本文主要討論的是在矩陣大小未知的情況。在陣列大小未知時,一般常用的定義方法是動態記憶體分配,對於乙個高維陣列,以2d的int 型陣列為例 1.將陣列定義為一維的,如 int array new int row colum 2.另一種方法是仍然定義為高維的形式,如 int array new int ...