电脑知识|欧美黑人一区二区三区|软件|欧美黑人一级爽快片淫片高清|系统|欧美黑人狂野猛交老妇|数据库|服务器|编程开发|网络运营|知识问答|技术教程文章 - 好吧啦网

您的位置:首頁技術文章
文章詳情頁

用Python進行柵格數據的分區統計和批量提取

瀏覽:125日期:2022-06-18 10:10:35

有時候我們會有這樣的想法,就是針對某個區域的柵格數據,要提取它的平均值或者其他統計指標,比如在一個省內提取多年的降雨數據,最后分區域地計算一些統計值,或者從多個柵格數據中提取某個區域的數值形成一個序列。為了方便,畫一個示意圖看看,比如就像提取這個區域中的某一個市的區域,然后形成一個序列數據,這就可以使用rasterstats庫了,此外的分區統計也可以用這個庫

這個實驗使用的數據格式分別是柵格(*.tif)和矢量(.shp),之后的分區統計操作和柵格數據的提取都是源于這兩類數據。為了能使用上這個rasterstats庫,選擇了在google colab平臺運行腳本,因為安裝庫實在是太方便了,在win上老是安裝不上的,在google notebook立馬就搞定了,而且可以把數據存儲到谷歌云盤,直接在notebook中就是可以鏈接使用的

用Python進行柵格數據的分區統計和批量提取

那么現在就開始做測試,使用的數據就是左側的柵格和矢量數據集導入相關的模塊

import geopandas as gpdimport pandas as pdimport numpy as npimport matplotlib.pyplot as pltimport rasterioimport rasterstatsfrom rasterio.plot import show# show()方法用來展示柵格圖形from rasterio.plot import show_hist# 用來展示直方圖import cartopy.crs as ccrsimport cartopy.feature as cfeaturefrom cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

使用geopandas和rasterio分別讀取矢量和柵格數據

# 使用geopandas讀取矢量數據districts = gpd.read_file(’/content/drive/MyDrive/Datashpraster/Data/Districts/districts.shp’)# 使用rasterio讀取柵格數據,柵格數據和矢量數據的坐標投影需要一致raster = rasterio.open(’/content/drive/MyDrive/Datashpraster/Data/Rainfall Data Rasters/2020-4-1.tif’)

# 把矢量數據和柵格數據繪制到一個axis上,這個axis不是坐標軸,而是圖形plt.rcParams[’font.family’] = ’Times New Roman’plt.rcParams[’font.size’] = 20fig, (ax1,ax2) = plt.subplots(1,2,figsize=(15,6))show(raster, ax=ax1,title=’Rainfall’)# 讀取進來的矢量數據可以直接調用gpd的plot()方法繪制districts.plot(ax=ax1, facecolor=’None’, edgecolor=’red’)show_hist(raster,ax=ax2,title=’hist’)plt.show()

先繪制一下結果看看

用Python進行柵格數據的分區統計和批量提取

讀取柵格數據:

# 提取雨量柵格值到numpy數組# 遵循GDAL規則從第一波段讀取rainfall_data = raster.read(1)rainfall_data

用Python進行柵格數據的分區統計和批量提取

開始分區統計:

# 設置坐標變換信息affine = raster.transform# 準備開始進行空間分區計算# 第一個參數是矢量分區,第二個是柵格,第三個是坐標變換信息,第四個是統計均值avg_rallrain = rasterstats.zonal_stats(districts,rainfall_data,affine=affine,stats=[’mean’],geojson_out=True)# avg_rallrain# 除了統計平均值之外,還有最大最小值那些

繪制一下,只是一個簡單的圖形而已

用Python進行柵格數據的分區統計和批量提取

當然第二部分更有意思,就是從多個分散的柵格數據中提取數據形成一個序列

,就是這些tif數據

用Python進行柵格數據的分區統計和批量提取

loop這些柵格數據集:

用Python進行柵格數據的分區統計和批量提取

獲得提取到的結果,沒錯,就是這么一個序列數據,然后就是繪圖了

用Python進行柵格數據的分區統計和批量提取

轉換數據格式

# 將Date列轉為時間型data[’Date’] = pd.to_datetime(data[’Date’], infer_datetime_format=True)# print(data)data[’Date’] = data[’Date’].dt.dateprint(data)

用Python進行柵格數據的分區統計和批量提取

繪圖結果就是簡單的圖形而已

# 準備繪制圖形fig,(ax1,ax2)= plt.subplots(2,1,figsize=(18,6))plt.rcParams[’font.size’] = 15data.plot(x=’Date’, y=’Average_RF_Porto’, ax=ax1, kind=’bar’, title=’Avg_Rail_Porto’)data.plot(x=’Date’, y=’Average_RF_Faro’, ax=ax2, kind=’bar’, title=’Avg_Rail_Faro’,color=’red’)#自動調整圖形的分布plt.tight_layout()plt.show()

用Python進行柵格數據的分區統計和批量提取

結果就這樣一個序列圖,目的就是從柵格提取指定的研究區,然后提取柵格的值,再來繪圖

雖然感覺不是那么花里胡哨的圖,但這個應該還是比較實用的,特別是大批量提取柵格值的時候。由于在google colab里面操作的步驟比較多,中間可能有省略的地方,但重要的應該都在文中了,當然也可以遷移運用到其他地方,也可以查看一下這個第三方庫的教程,比如read(1)是什么意思,官網的docs就寫得有,實在是很方便的

以上就是用Python進行柵格數據的分區統計和批量提取的詳細內容,更多關于Python 柵格數據的分區統計和批量提取 的資料請關注好吧啦網其它相關文章!

標簽: Python 編程
相關文章:
主站蜘蛛池模板: 飞象网 - 通信人每天必上的网站| 广东护栏厂家-广州护栏网厂家-广东省安麦斯交通设施有限公司 | 聚丙烯酰胺_阴离子_阳离子「用量少」巩义亿腾厂家直销,售后无忧 聚合甘油__盐城市飞龙油脂有限公司 | 沈阳网站建设_沈阳网站制作_沈阳网页设计-做网站就找示剑新零售 沈阳缠绕膜价格_沈阳拉伸膜厂家_沈阳缠绕膜厂家直销 | 运动木地板_体育木地板_篮球馆木地板_舞台木地板-实木运动地板厂家 | 江苏皓越真空设备有限公司| 天津货架厂_穿梭车货架_重型仓储货架_阁楼货架定制-天津钢力仓储货架生产厂家_天津钢力智能仓储装备 | 伸缩器_伸缩接头_传力接头-巩义市润达管道设备制造有限公司 | 液压压力机,液压折弯机,液压剪板机,模锻液压机-鲁南新力机床有限公司 | 双杰天平-国产双杰电子天平-美国双杰-常熟双杰仪器 | 天一线缆邯郸有限公司_煤矿用电缆厂家_矿用光缆厂家_矿用控制电缆_矿用通信电缆-天一线缆邯郸有限公司 | 卓能JOINTLEAN端子连接器厂家-专业提供PCB接线端子|轨道式端子|重载连接器|欧式连接器等电气连接产品和服务 | 汽车润滑油厂家-机油/润滑油代理-高性能机油-领驰慧润滑科技(河北)有限公司 | 伺服电机维修、驱动器维修「安川|三菱|松下」伺服维修公司-深圳华创益 | 家德利门业,家居安全门,别墅大门 - 安徽家德利门业有限公司 | 在线钠离子分析仪-硅酸根离子浓度测定仪-油液水分测定仪价格-北京时代新维测控设备有限公司 | 番茄畅听邀请码怎么输入 - Dianw8.com | 东风体检车厂家_公共卫生体检车_医院体检车_移动体检车-锦沅科贸 | 直线模组_滚珠丝杆滑台_模组滑台厂家_万里疆科技 | 金属切削液-脱水防锈油-电火花机油-抗磨液压油-深圳市雨辰宏业科技发展有限公司 | 韦伯电梯有限公司| 硫化罐-胶管硫化罐-山东鑫泰鑫智能装备有限公司 | 驾驶人在线_专业学车门户网站 | 东莞韩创-专业绝缘骨架|马达塑胶零件|塑胶电机配件|塑封电机骨架厂家 | 深圳快餐店设计-餐饮设计公司-餐饮空间品牌全案设计-深圳市勤蜂装饰工程 | 深圳激光打标机_激光打标机_激光焊接机_激光切割机_同体激光打标机-深圳市创想激光科技有限公司 深圳快餐店设计-餐饮设计公司-餐饮空间品牌全案设计-深圳市勤蜂装饰工程 | 家用净水器代理批发加盟_净水机招商代理_全屋净水器定制品牌_【劳伦斯官网】 | 医用酒精_84消毒液_碘伏消毒液等医用消毒液-漓峰消毒官网 | 南京泽朗生物科技有限公司-液体饮料代加工_果汁饮料代加工_固体饮料代加工 | 冷藏车厂家|冷藏车价格|小型冷藏车|散装饲料车厂家|程力专用汽车股份有限公司销售十二分公司 | 中高频感应加热设备|高频淬火设备|超音频感应加热电源|不锈钢管光亮退火机|真空管烤消设备 - 郑州蓝硕工业炉设备有限公司 | 必胜高考网_全国高考备考和志愿填报信息平台| NM-02立式吸污机_ZHCS-02软轴刷_二合一吸刷软轴刷-厦门地坤科技有限公司 | 双能x射线骨密度检测仪_dxa骨密度仪_双能x线骨密度仪_品牌厂家【品源医疗】 | 警用|治安|保安|不锈钢岗亭-售货亭价格-垃圾分类亭-移动厕所厂家-苏州灿宇建材 | 消防设施操作员考试报名时间,报名入口,报考条件 | 济南品牌设计-济南品牌策划-即合品牌策划设计-山东即合官网 | 河南卓美创业科技有限公司-河南卓美防雷公司-防雷接地-防雷工程-重庆避雷针-避雷器-防雷检测-避雷带-避雷针-避雷塔、机房防雷、古建筑防雷等-山西防雷公司 | 3dmax渲染-效果图渲染-影视动画渲染-北京快渲科技有限公司 | 渗透仪-直剪仪-三轴仪|苏州昱创百科 | 合肥注册公司|合肥代办营业执照、2024注册公司流程 |