您好,登錄后才能下訂單哦!
這篇文章將為大家詳細講解有關Python中如何實現(xiàn)柵格圖像邊界提取方法,小編覺得挺實用的,因此分享給大家做個參考,希望大家閱讀完這篇文章后可以有所收獲。
在GIS中,柵格屬性里有關于柵格自身的信息,背景(nodata value)對于識別一張圖像的邊界像元尤為重要,我們目的只要把每行每列中的第一次出現(xiàn)不是nodata的像元和最后一次出現(xiàn)nodata的前一個像元就可以了。
對于柵格,可以用ArcPy中的RasterToNumpyArray函數(shù)將將柵格轉成numpy數(shù)組,然后就可以按照所想讀取出每行列中首尾像元。
以下是部分代碼提取邊界像元的核心算法,其實是很簡單的一個思路(假設0是nodata value)。
a=[[0 for col in range(Raster.width)]for row in range(Raster.height)] for i in range(0,Raster.width): ... for j in range(0,Raster.height): ... if(myRaster[j][i]!=0 and myRaster[j-1][i]==0): ... a[j][i]=myRaster[j][i] ... if(myRaster[j][i]==0 and myRaster[j-1][i]!=0): ... a[j-1][i]=myRaster[j-1][i] ... for i in range(0,myRaster.height): ... for j in range(0,myRaster.width): ... if(arr[i][j]!=0 and arr[i][j-1]==0): ... a[i][j]=arr[i][j] ... if(arr[i][j]==0 and arr[i][j-1]!=0): ... a[i][j-1]=arr[i][j-1]
關于“Python中如何實現(xiàn)柵格圖像邊界提取方法”這篇文章就分享到這里了,希望以上內(nèi)容可以對大家有一定的幫助,使各位可以學到更多知識,如果覺得文章不錯,請把它分享出去讓更多的人看到。
免責聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉載和分享為主,文章觀點不代表本網(wǎng)站立場,如果涉及侵權請聯(lián)系站長郵箱:is@yisu.com進行舉報,并提供相關證據(jù),一經(jīng)查實,將立刻刪除涉嫌侵權內(nèi)容。