快捷導(dǎo)航

python獲取柵格點和面值的實現(xiàn)

更新時間：2020年03月10日 14:24:42 作者：忠言睿長

這篇文章主要介紹了python獲取柵格點和面值的實現(xiàn)，具有很好的參考價值，希望對大家有多幫助。一起跟隨小編過來看看吧

1、獲取指定位置的點值：

OutputFile = open(statisticResultTXT, 'w') 

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

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

OutputFile.close()

2、獲取區(qū)域的均值極值

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_計算格點數(shù)據(jù)降水變化斜率

我就廢話不多說了，還是直接看代碼吧！

import numpy as np
import pandas as pd
import codecs
 
path=r'C:\try.csv' #讀取格點數(shù)據(jù)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 #轉(zhuǎn)置
n=80 #格點數(shù)
m=17 #年數(shù)
sl=[] #空列表
for i in list(range(0,80)):#循環(huán)計算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獲取柵格點和面值的實現(xiàn)就是小編分享給大家的全部內(nèi)容了，希望能給大家一個參考，也希望大家多多支持腳本之家。

您可能感興趣的文章:

相關(guān)文章

YOLOv5目標(biāo)檢測之a(chǎn)nchor設(shè)定
在訓(xùn)練yolo網(wǎng)絡(luò)檢測目標(biāo)時,需要根據(jù)待檢測目標(biāo)的位置大小分布情況對anchor進行調(diào)整,使其檢測效果盡可能提高,下面這篇文章主要給大家介紹了關(guān)于YOLOv5目標(biāo)檢測之a(chǎn)nchor設(shè)定的相關(guān)資料,需要的朋友可以參考下
2022-05-05
python+Word2Vec實現(xiàn)中文聊天機器人的示例代碼
本文主要介紹了python+Word2Vec實現(xiàn)中文聊天機器人，文中通過示例代碼介紹的非常詳細，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
2023-03-03
Python matplotlib畫圖與中文設(shè)置操作實例分析
這篇文章主要介紹了Python matplotlib畫圖與中文設(shè)置操作,結(jié)合實例形式分析了Python使用matplotlib進行圖形繪制及中文設(shè)置相關(guān)操作技巧,需要的朋友可以參考下
2019-04-04
python如何查看微信消息撤回
這篇文章主要為大家詳細介紹了python實現(xiàn)查看微信消息撤回的方法，具有一定的參考價值，感興趣的小伙伴們可以參考一下
2018-11-11
淺談Python從全局與局部變量到裝飾器的相關(guān)知識
今天給大家?guī)淼氖顷P(guān)于Python的相關(guān)知識,文章圍繞著Python從全局與局部變量到裝飾器的相關(guān)知識展開,文中有非常詳細的介紹及代碼示例,需要的朋友可以參考下
2021-06-06
Python+Opencv實現(xiàn)把圖片、視頻互轉(zhuǎn)的示例
這篇文章主要介紹了Python+Opencv實現(xiàn)把圖片、視頻互轉(zhuǎn)的示例，幫助大家更好的理解和實用python，感興趣的朋友可以了解下
2020-12-12
Python使用pip安裝pySerial串口通訊模塊
這篇文章主要為大家詳細介紹了Python使用pip安裝pySerial串口通訊模塊，具有一定的參考價值，感興趣的小伙伴們可以參考一下
2018-04-04
詳解Python核心編程中的淺拷貝與深拷貝
本片文章給大家講述了Python核心編程中的淺拷貝與深拷貝的相關(guān)知識點，有需要的朋友跟著學(xué)習(xí)下吧。
2018-01-01
簡單了解Python中的幾種函數(shù)
這篇文章主要介紹了簡單了解Python中的幾種函數(shù)，具有一定參考價值。需要的朋友可以了解下。
2017-11-11
基于Python實現(xiàn)Excel轉(zhuǎn)Markdown表格
Markdown（也簡稱md）作為一種輕量級標(biāo)記語言，因其易寫易讀，效果美觀大方，不僅被眾多網(wǎng)站使用，也是程序員們做筆記、寫文檔的首選。本文將利用Python實現(xiàn)Excel轉(zhuǎn)Markdown表格，感興趣的可以了解一下
2022-04-04