python获取栅格点和面值的实现


Posted in Python onMarch 10, 2020

1、获取指定位置的点值:

OutputFile = open(statisticResultTXT, 'w') 

cellvalue=arcpy.GetCellValue_management(inputfilepath+filenname+".tif",staionXY, "1")

OutputFile.write(stationID+"_"+filenname+""+str(cellvalue)+'\n')

OutputFile.close()

2、获取区域的均值极值

utputFile = open(statisticResult, 'w') 
env.workspace = ('E:\loessplateau\ET_Spline\extract_year')
rlist = arcpy.ListRasters('*','tif')
for file in rlist:         
minValueInfo= arcpy.GetRasterProperties_management(file,'MINIMUM')
minValue = minValueInfo.getOutput(0)
maxValueInfo = arcpy.GetRasterProperties_management(file,'MAXIMUM')
maxValue = maxValueInfo.getOutput(0)
meanValueInfo = arcpy.GetRasterProperties_management(file,'MEAN')
meanValue = meanValueInfo.getOutput(0)
print file + ',' + str(minValue) + ',' + str(maxValue) + ',' + str(meanValue) +'\n'
OutputFile.write(file + ',' + str(minValue) + ',' + str(maxValue) + ',' + str(meanValue) +'\n')
OutputFile.close()

补充知识:python_计算格点数据降水变化斜率

我就废话不多说了,还是直接看代码吧!

import numpy as np
import pandas as pd
import codecs
 
path=r'C:\try.csv' #读取格点数据CSV
data=pd.read_csv(path,encoding='gbk')
outpath=r'C:\rate.csv'#新建存放斜率的CSV
file_csv = codecs.open(outpath,'w+','utf-8') #打开方式
data1=data.T #转置
n=80 #格点数
m=17 #年数
sl=[] #空列表
for i in list(range(0,80)):#循环计算80个格点的变化斜率
  print (i)
  sp=[]
  y=data1[i]
  y=y.values
 
  for k in range(0,m-1):
    for j in range(k+1,m):
      print(k,j)
     
      s=(y[j]-y[k])/(j-k)  #计算斜率
      sp.append(s)
  sl.append(np.median(sp))
  
slop= pd.DataFrame(sl)
slop.to_csv(outpath,header=['slop']) #将SLOP写入CSV

以上这篇python获取栅格点和面值的实现就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python并发2之使用asyncio处理并发
Dec 21 Python
python字符串string的内置方法实例详解
May 14 Python
python查看列的唯一值方法
Jul 17 Python
Python通过调用有道翻译api实现翻译功能示例
Jul 19 Python
python自带tkinter库实现棋盘覆盖图形界面
Jul 17 Python
django数据关系一对多、多对多模型、自关联的建立
Jul 24 Python
python文档字符串(函数使用说明)使用详解
Jul 30 Python
Python自动生成代码 使用tkinter图形化操作并生成代码框架
Sep 18 Python
Python接口测试get请求过程详解
Feb 28 Python
Python实现从N个数中找到最大的K个数
Apr 02 Python
Django 解决阿里云部署同步数据库报错的问题
May 14 Python
TensorFlow保存TensorBoard图像操作
Jun 23 Python
Python列表切片常用操作实例解析
Mar 10 #Python
Python numpy多维数组实现原理详解
Mar 10 #Python
python中使用you-get库批量在线下载bilibili视频的教程
Mar 10 #Python
Python字符串hashlib加密模块使用案例
Mar 10 #Python
Python中求对数方法总结
Mar 10 #Python
Python标准库shutil模块使用方法解析
Mar 10 #Python
python实现提取COCO,VOC数据集中特定的类
Mar 10 #Python
You might like
PHP PDO fetch 模式各种参数的输出结果一览
2015/01/07 PHP
php源码分析之DZX1.5字符串截断函数cutstr用法
2015/06/17 PHP
php生成静态html页面的方法(2种方法)
2015/09/14 PHP
PHP生成可点击刷新的验证码简单示例
2016/05/13 PHP
修改yii2.0用户登录使用的user表为其它的表实现方法(推荐)
2017/08/01 PHP
PHP基于DateTime类解决Unix时间戳与日期互转问题【针对1970年前及2038年后时间戳】
2018/06/13 PHP
php输出反斜杠的实例方法
2019/09/19 PHP
Yii2.0框架behaviors方法使用实例分析
2019/09/30 PHP
Prototype Selector对象学习
2009/07/23 Javascript
用js解决数字不能换行问题
2010/08/10 Javascript
Extjs TimeField 显示正常时间格式的代码
2011/06/28 Javascript
JQuery中的事件及动画用法实例
2015/01/26 Javascript
jquery实现上下左右滑动的方法
2015/02/09 Javascript
jQuery实现按键盘方向键翻页特效
2015/03/18 Javascript
JavaScript随机生成信用卡卡号的方法
2015/04/07 Javascript
jQuery实现控制文字内容溢出用省略号(…)表示的方法
2016/02/26 Javascript
JavaScript体验异步更好的解决办法
2018/01/08 Javascript
微信小程序支付前端源码
2018/08/29 Javascript
vue实现pdf导出解决生成canvas模糊等问题(推荐)
2018/10/18 Javascript
three.js利用gpu选取物体并计算交点位置的方法示例
2019/11/25 Javascript
typescript编写微信小程序创建项目的方法
2021/01/29 Javascript
Python设计模式编程中解释器模式的简单程序示例分享
2016/03/02 Python
Python系统监控模块psutil功能与经典用法分析
2018/05/24 Python
Python turtle绘画象棋棋盘
2019/08/21 Python
Python使用psutil获取进程信息的例子
2019/12/17 Python
Django Haystack 全文检索与关键词高亮的实现
2020/02/17 Python
什么是Python中的顺序表
2020/06/02 Python
非常漂亮的CSS3百叶窗焦点图动画
2016/02/24 HTML / CSS
HTML5 层的叠加的实现
2020/07/07 HTML / CSS
美国乡村商店:Plow & Hearth
2016/09/12 全球购物
精彩的大学生自我评价
2013/11/17 职场文书
研究生毕业鉴定
2014/01/29 职场文书
网上快餐厅创业计划书
2014/02/01 职场文书
聚美优品陈欧广告词
2014/03/14 职场文书
《青海高原一株柳》教学反思
2014/04/25 职场文书
解决pytorch-gpu 安装失败的记录
2021/05/24 Python