关于Python 的简单栅格图像边界提取方法


Posted in Python onJuly 05, 2019

在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 的简单栅格图像边界提取方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
Pyramid Mako模板引入helper对象的步骤方法
Nov 27 Python
列举Python中吸引人的一些特性
Apr 09 Python
举例讲解Python的Tornado框架实现数据可视化的教程
May 02 Python
Python使用turtule画五角星的方法
Jul 09 Python
Windows平台Python连接sqlite3数据库的方法分析
Jul 12 Python
windows下添加Python环境变量的方法汇总
May 14 Python
jupyter notebook 中输出pyecharts图实例
Apr 23 Python
python join方法使用详解
Jul 30 Python
Django rstful登陆认证并检查session是否过期代码实例
Aug 13 Python
python 操作hive pyhs2方式
Dec 21 Python
基于Django OneToOneField和ForeignKey的区别详解
Mar 30 Python
Python应用实现双指数函数及拟合代码实例
Jun 19 Python
Python3+Appium安装使用教程
Jul 05 #Python
Python叠加两幅栅格图像的实现方法
Jul 05 #Python
用vue.js组件模拟v-model指令实例方法
Jul 05 #Python
python买卖股票的最佳时机(基于贪心/蛮力算法)
Jul 05 #Python
Python 计算任意两向量之间的夹角方法
Jul 05 #Python
python实现两个经纬度点之间的距离和方位角的方法
Jul 05 #Python
Python3+Appium实现多台移动设备操作的方法
Jul 05 #Python
You might like
php下通过POST还是GET来传值
2008/06/05 PHP
神盾加密解密教程(一)PHP变量可用字符
2014/05/28 PHP
js/jQuery简单实现选项卡功能
2014/01/02 Javascript
浅谈Javascript数组索引
2015/07/29 Javascript
javascript跨域的方法汇总
2015/10/23 Javascript
详解javascript中的事件处理
2015/11/06 Javascript
通过jsonp获取json数据实现AJAX跨域请求
2017/01/22 Javascript
详解微信开发中snsapi_base和snsapi_userinfo及静默授权的实现
2017/03/11 Javascript
form表单序列化详解(推荐)
2017/08/15 Javascript
JS实现前端页面的搜索功能
2018/06/12 Javascript
vue读取本地的excel文件并显示在网页上方法示例
2019/05/29 Javascript
layui实现数据分页功能
2019/07/27 Javascript
[04:54]DOTA2 2017国际邀请赛:上届冠军WINGS采访短片
2017/08/09 DOTA
基于python3 类的属性、方法、封装、继承实例讲解
2017/09/19 Python
利用Pandas读取文件路径或文件名称包含中文的csv文件方法
2018/07/04 Python
python实现批量注册网站用户的示例
2019/02/22 Python
tensorflow从ckpt和从.pb文件读取变量的值方式
2020/05/26 Python
python 实现超级玛丽游戏
2020/11/25 Python
html5小技巧之通过document.head获取head元素
2014/06/04 HTML / CSS
详解html2canvas截图不能截取圆角图片的解决方案
2018/01/30 HTML / CSS
英国网上香水店:Fragrance Direct
2016/07/20 全球购物
极简鞋类,赤脚的感觉:Lems Shoes
2019/08/06 全球购物
Kickers鞋英国官网:男士、女士和儿童鞋
2021/03/08 全球购物
外贸公司实习自我鉴定
2013/09/24 职场文书
机电专业大学生求职信
2013/10/04 职场文书
投资合作协议书范本
2014/04/17 职场文书
环保建议书500字
2014/05/14 职场文书
学习之星事迹材料
2014/05/17 职场文书
会计专业毕业生自荐书
2014/06/25 职场文书
关于读书的活动方案
2014/08/14 职场文书
捐款感谢信
2015/01/20 职场文书
2016年秋季运动会加油稿
2015/12/21 职场文书
《田忌赛马》教学反思
2016/02/19 职场文书
社交电商模式的兴起:这些新的商机千万别错过
2019/07/26 职场文书
基于HTML十秒做出淘宝页面
2021/10/24 HTML / CSS
十大最强电系宝可梦,阿尔宙斯电系之一,第七被称为雷神
2022/03/18 日漫