快捷導(dǎo)航

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

更新時(shí)間：2020年03月10日 14:24:42 作者：忠言睿長(zhǎng)

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

1、獲取指定位置的點(diǎn)值：

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()

補(bǔ)充知識(shí)：python_計(jì)算格點(diǎn)數(shù)據(jù)降水變化斜率

我就廢話(huà)不多說(shuō)了，還是直接看代碼吧！

import numpy as np
import pandas as pd
import codecs
 
path=r'C:\try.csv' #讀取格點(diǎn)數(shù)據(jù)CSV
data=pd.read_csv(path,encoding='gbk')
outpath=r'C:\rate.csv'#新建存放斜率的CSV
file_csv = codecs.open(outpath,'w+','utf-8') #打開(kāi)方式
data1=data.T #轉(zhuǎn)置
n=80 #格點(diǎn)數(shù)
m=17 #年數(shù)
sl=[] #空列表
for i in list(range(0,80)):#循環(huán)計(jì)算80個(gè)格點(diǎn)的變化斜率
  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)  #計(jì)算斜率
      sp.append(s)
  sl.append(np.median(sp))
  
slop= pd.DataFrame(sl)
slop.to_csv(outpath,header=['slop']) #將SLOP寫(xiě)入CSV

以上這篇python獲取柵格點(diǎn)和面值的實(shí)現(xiàn)就是小編分享給大家的全部?jī)?nèi)容了，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。

您可能感興趣的文章:

相關(guān)文章

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