在GIS中，柵格屬性里有關于柵格自身的信息，背景（nodata value）對于識別一張圖像的邊界像元尤為重要，我們目的只要把每行每列中的第一次出現不是nodata的像元和最后一次出現nodata的前一個像元就可以了。
對于柵格，可以用ArcPy中的RasterToNumpyArray函數將將柵格轉成numpy數組，然后就可以按照所想讀取出每行列中首尾像元。
以下是部分代碼提取邊界像元的核心算法，其實是很簡單的一個思路(假設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 的簡單柵格圖像邊界提取方法就是小編分享給大家的全部內容了，希望能給大家一個參考，也希望大家多多支持腳本之家。
更多文章、技術交流、商務合作、聯系博主
微信掃碼或搜索：z360901061
微信掃一掃加我為好友
QQ號聯系： 360901061
您的支持是博主寫作最大的動力，如果您喜歡我的文章，感覺我的文章對您有幫助，請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧，狠狠點擊下面給點支持吧，站長非常感激您！手機微信長按不能支付解決辦法：請將微信支付二維碼保存到相冊，切換到微信，然后點擊微信右上角掃一掃功能，選擇支付二維碼完成支付。
【本文對您有幫助就好】元
