当前位置: 首页 > 编程笔记 >

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

易京
2023-03-14
本文向大家介绍python获取栅格点和面值的实现,包括了python获取栅格点和面值的实现的使用技巧和注意事项,需要的朋友参考一下

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获取栅格点和面值的实现就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持小牛知识库。

 类似资料:
  • Framework7 有灵活的布局网格,允许你按需求放置内容: <!-- Each "cells" row should be wrapped with div class="row" --> <div class="row"> <!-- Each "cell" has col-[widht in percents] class --> <div class="col-50">50

  • 一个节点就象征着一个数据,现在的或者未来的,它是一种预期,所以我们也对应的有获取即时值和获取未来值两种方式。 获取即时值 对于一个节点来说,访问T value属性是最简单有效的形式,但是由于空值,我们可能需要特殊注意一下: EZRMutableNode<NSNumber *> *node = [EZRMutableNode new]; NSNumber *value = node.value;

  • 栅格系统简介: MUI 提供了非常简单实用的12列响应式栅格系统。使用时只需在外围容器上添加.mui-row,在列上添加 .mui-col-[sm|xs]-[1-12],即可 栅格参数: 尺寸 超小屏幕(<400px)(Extrasmall) 小屏幕(≥400px) Small 类前缀 .mui-col-xs-[1-12] .mui-col-sm-[1-12] 列(column)数 12 可嵌套

  • 定义 栅格布局组件。 图片展示 代码演示 import Grid from 'pile/dist/components/grid' <Grid flexCells> <div></div> <div></div> </Grid> 属性 参数 描述 数据类型 默认值 flexCells 是否flex布局 true

  • 本文向大家介绍Python 获取numpy.array索引值的实例,包括了Python 获取numpy.array索引值的实例的使用技巧和注意事项,需要的朋友参考一下 举个例子: 我想获取其中值等于7的那个值的下标,以便于用于其他计算。 如果使用np.where,如: 运行结果是: 显然(array([7]),)中的数字7我是没法提取出来做运算的,这是一个tuple。 处理方法是: 运行结果为:

  • 在我的PostgreSQL数据库中,有一个栅格表和一个具有相同投影的点表。点表具有诸如 id、地址和几何等列。栅格表具有 id、r_proj4 和 rast 等列。栅格表的每一行描述一个栅格切片。如何检索每个点的栅格值? 我希望粗略的指导方针如何解决这个问题和PostGIS代码示例。