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 相关文章推荐
Saltstack快速入门简单汇总
Mar 01 Python
python3 pillow生成简单验证码图片的示例
Sep 19 Python
python 实现语音聊天机器人的示例代码
Dec 02 Python
10 分钟快速入门 Python3的教程
Jan 29 Python
解决Pycharm调用Turtle时 窗口一闪而过的问题
Feb 16 Python
使用python爬取抖音视频列表信息
Jul 15 Python
结合OpenCV与TensorFlow进行人脸识别的实现
Oct 10 Python
python超时重新请求解决方案
Oct 21 Python
关于tf.matmul() 和tf.multiply() 的区别说明
Jun 18 Python
python开发一个解析protobuf文件的简单编译器
Nov 17 Python
利用Python将多张图片合成视频的实现
Nov 23 Python
Python爬虫新手入门之初学lxml库
Dec 20 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
很温暖很温暖的Lester Young
2021/03/03 冲泡冲煮
测试您的 PHP 水平的题目
2007/05/30 PHP
帝国cms目录结构分享
2015/07/06 PHP
使用PHP+AJAX让WordPress动态加载文章的教程
2015/12/11 PHP
php实现贪吃蛇小游戏
2016/07/26 PHP
javascript开发中因空格引发的错误
2010/11/08 Javascript
基于Jquery的开发个代阴影的对话框效果代码
2011/07/28 Javascript
JS去除右边逗号的简单方法
2013/07/03 Javascript
JS的参数传递示例介绍
2014/02/08 Javascript
js实现数组转换成json
2015/06/26 Javascript
深入探讨Vue.js组件和组件通信
2016/09/12 Javascript
Bootstrap中点击按钮后变灰并显示加载中实例代码
2016/09/23 Javascript
KnockoutJS 3.X API 第四章之数据控制流foreach绑定
2016/10/10 Javascript
微信小程序 开发之顶部导航栏实例代码
2017/02/23 Javascript
AngularJS实现表格的增删改查(仅限前端)
2017/07/04 Javascript
AngularJS2 与 D3.js集成实现自定义可视化的方法
2017/12/01 Javascript
Vue父子组建的简单通信之控制开关Switch的实现
2018/06/04 Javascript
jQuery 函数实例分析【函数声明、函数表达式、匿名函数等】
2020/05/19 jQuery
JS+canvas五子棋人机对战实现步骤详解
2020/06/04 Javascript
使用Vue实现一个树组件的示例
2020/11/06 Javascript
js+html+css实现手动轮播和自动轮播
2020/12/30 Javascript
python修改FTP服务器上的文件名
2019/09/11 Python
Selenium元素定位的30种方式(史上最全)
2020/05/11 Python
使用keras和tensorflow保存为可部署的pb格式
2020/05/25 Python
python日志通过不同的等级打印不同的颜色(示例代码)
2021/01/13 Python
解决pip安装tensorflow中出现的no module named tensorflow.python 问题方法
2021/02/20 Python
TripAdvisor越南:全球领先的旅游网站
2017/09/21 全球购物
澳大利亚领先的内衣店:Bendon Lingerie澳大利亚
2020/05/15 全球购物
简述数据库的设计过程
2015/06/22 面试题
什么是TCP/IP
2014/07/27 面试题
初中生个人学习的自我评价
2013/12/04 职场文书
行政人员工作职责
2013/12/05 职场文书
办公室秘书岗位职责范本
2014/02/11 职场文书
咖啡店创业计划书
2014/08/15 职场文书
离婚协议书格式
2015/01/26 职场文书
Nginx报错104:Connection reset by peer问题的解决及分析
2022/07/23 Servers