關(guān)于Python 的簡單柵格圖像邊界提取方法
在GIS中,柵格屬性里有關(guān)于柵格自身的信息,背景(nodata value)對于識別一張圖像的邊界像元尤為重要,我們目的只要把每行每列中的第一次出現(xiàn)不是nodata的像元和最后一次出現(xiàn)nodata的前一個(gè)像元就可以了。
對于柵格,可以用ArcPy中的RasterToNumpyArray函數(shù)將將柵格轉(zhuǎn)成numpy數(shù)組,然后就可以按照所想讀取出每行列中首尾像元。
以下是部分代碼提取邊界像元的核心算法,其實(shí)是很簡單的一個(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 的簡單柵格圖像邊界提取方法就是小編分享給大家的全部內(nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。
相關(guān)文章
Pycharm學(xué)習(xí)教程(7)虛擬機(jī)VM的配置教程
這篇文章主要為大家詳細(xì)介紹了最全的Pycharm學(xué)習(xí)教程第七篇,Python快捷鍵相關(guān)設(shè)置,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2017-05-05使用Python爬取Json數(shù)據(jù)的示例代碼
這篇文章主要介紹了使用Python爬取Json數(shù)據(jù)的示例代碼,本文給大家介紹的非常詳細(xì),對大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2020-12-12pycharm遠(yuǎn)程開發(fā)項(xiàng)目的實(shí)現(xiàn)步驟
這篇文章主要介紹了pycharm遠(yuǎn)程開發(fā)項(xiàng)目的實(shí)現(xiàn)步驟,小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,也給大家做個(gè)參考。一起跟隨小編過來看看吧2019-01-01Python selenium實(shí)現(xiàn)斷言3種方法解析
這篇文章主要介紹了Python selenium實(shí)現(xiàn)斷言3種方法解析,文中通過示例代碼介紹的非常詳細(xì),對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-09-09python解析html提取數(shù)據(jù),并生成word文檔實(shí)例解析
這篇文章主要介紹了python解析html提取數(shù)據(jù),并生成word文檔實(shí)例解析,小編覺得還是挺不錯(cuò)的,具有一定借鑒價(jià)值,需要的朋友可以參考下2018-01-01python發(fā)送郵件示例(支持中文郵件標(biāo)題)
python發(fā)送中文郵件示例,支持中文郵件標(biāo)題和中文郵件內(nèi)容。支持多附件。根據(jù)用戶名推測郵件服務(wù)器提供商2014-02-02Python合并2個(gè)字典成1個(gè)新字典的方法(9種)
這篇文章主要介紹了Python合并2個(gè)字典成1個(gè)新字典的方法,本文通過實(shí)例代碼給大家分享9中方法,需要的朋友可以參考下2019-12-12