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