爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2652|回复: 0

python在筛选温度数据时遇到了问题

[复制链接]

新浪微博达人勋

发表于 2022-7-31 12:57:11 | 显示全部楼层 |阅读模式
Python
系统平台:
问题截图:
问题概况: 求问各路大神,导师发给我了一个2014~2019年日最高温度数据,是nc文件,让我先筛选出高温天,这里高温的定义是:区域内的观测站点日最高温度tmax>=35℃的比例占总站点数的30%即为高温,我研究的区域范围是10-35°N,105-125°E,目前我只是将区域定下来了,之后不知道该如何操作,请大家帮我看一下下一步该怎样写代码~

site1是研究区域内的站点
我看过提问的智慧: 看过
自己思考时长(天): 5

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
求问各路大神,导师发给我了一个2014~2019年日最高温度数据,是nc文件,让我先筛选出高温天,这里高温的定义是:区域内的观测站点日最高温度tmax>=35℃的比例占总站点数的30%即为高温,我研究的区域范围是10-35°N,105-125°E,目前我只是将区域定下来了,之后不知道该如何操作,请大家帮我看一下下一步该怎样写代码~
site1是研究区域内的站点
path1 =r'C:\Users\QCN\Desktop\data\T\DailyAirTmax_2014_2019(1).nc'
data1=xr.open_dataset(path1)
lon1=data1.longitude
lat1=data1.latitude
date1=data1.date
tmax1=data1.tmax
site = np.full([len(lon1)],np.nan)
for siteID in range(0,len(lon1)):
           if 10 <= lat1[siteID] <=35 and 105 <= lon1[siteID] <=125:
               site[siteID] = siteID
site1=site[~np.isnan(site)]
site_need=site1[0:1149]
数据的属性是这样的:
<xarray.Dataset>
Dimensions:    (date: 366, site: 2540, year: 6)
Coordinates:
  * year       (year) float64 2.014e+03 2.015e+03 ... 2.018e+03 2.019e+03
  * date       (date) float64 1.0 2.0 3.0 4.0 5.0 ... 363.0 364.0 365.0 366.0
Dimensions without coordinates: site
Data variables:
    tmax       (date, site, year) float64 ...
    longitude  (site) float64 ...
    latitude   (site) float64 ...
Attributes:
    description:  Daily maximum air temperature

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表