首页
学习
活动
专区
圈层
工具
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

利用 pandas 和 xarray 整理气象站点数据

作者:石异 (南京大学大气科学学院,硕士生) 利用 pandas 和 xarray 整理气象站点数据 平时用 xarray 库在处理 nc 格式的数据非常方便,但偶尔还是要用到一些站点数据来辅助分析,而站点数据一般都是用文本文件存储的...,此外,其中有不少特征值比如30XXX代表缺测/微量的情况,用Fortran处理也有不小的麻烦。...一、 目标和步骤 将上图示例的文件处理为(站点,时间)坐标的 nc 格式数据,方便以后直接读取,主要有以下几个步骤: 将文本文件读取为 DataFrame 并将无效值替换为 Nan 将时间信息处理为...\s 代表空白字符,+ 表示前面的字符至少重复一次(具体查看正则表达式的用法) na_values 选项将把指定的值替换为 Nan parse_dates=False 防止将某些字符解析为日期 StaDir...lat = stainfo.loc[ind, '纬度'].apply(LatLng_Rad2Dec) # 转换为十进制小数 lon = stainfo.loc[ind, '经度'].apply(LatLng_Rad2Dec

5.4K13

利用 pandas 和 xarray 整理气象站点数据

利用 pandas 和 xarray 整理气象站点数据 平时用 xarray 库在处理 nc 格式的数据非常方便,但偶尔还是要用到一些站点数据来辅助分析,而站点数据一般都是用文本文件存储的,比如下图这种格式...一、 目标和步骤 将上图示例的文件处理为(站点,时间)坐标的 nc 格式数据,方便以后直接读取,主要有以下几个步骤: 将文本文件读取为 DataFrame 并将无效值替换为 Nan 将时间信息处理为...df_t.loc[df_t['20-20时降水量'] >= 29999, '20-20时降水量'] = np.nan # 替换掉所有特征值 df_t.insert( # 插入日期列,此时并不以此为索引...\s 代表空白字符,+ 表示前面的字符至少重复一次(具体查看正则表达式的用法) na_values 选项将把指定的值替换为 Nan parse_dates=False 防止将某些字符解析为日期 StaDir...lat = stainfo.loc[ind, '纬度'].apply(LatLng_Rad2Dec) # 转换为十进制小数 lon = stainfo.loc[ind, '经度'].apply(LatLng_Rad2Dec

10.2K41
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    带你走近AngularJS - 体验指令实例

    一旦元素拥有了ID值,方法将通过jQuery来选择具有"accordion-toggle"类的子元素并且设置它的 "data-parent" 和 "href" 属性。...transclude 属性为true表明选项卡包含HTML标签。scope 下的 "title" 属性将会被实例所替代。 这个例子中的模板比较复杂。...模板中"{{title}}" 属性将会显示标签名称。目前我们仅仅实现了纯文本显示,没有定义其样式。我们使用link 方法可以替换标题为HTML源码从而得到更丰富的样式。...就这样,我们完成了第一个具有实用价值的指令。...Google Maps 指令 下一个例子是创建Google地图的指令: ? Google Maps 指令 在我们创建指令之前,我们需要添加Google APIs 引用到页面中: <!

    2.4K50

    基于NB-IoT的智慧路灯监控系统(NB-IoT专栏—实战篇5:手机应用开发)

    画像分析也可对某地方或某用户进行大数据AI分析得到监测数据,并且能够实时向用户推送用电情况,并为其用户进行合理的用电安排及方案。 ?...此时需要对路灯上的二维码进行扫码,通过扫码得到该路灯的设备信息,从而在主界面中可查看该路灯的其他信息。 效果如下所示: ?...Y值 // post是因为在前面得到当前view是再post中完成,如果不这样直接获取位置信息,会出现老的数据,不能获取正确的数据...childLayout.isWeekState()) { //不是周状态也不是月状态时,自动滑动 autoScroll(); }...周状态下 月日历的getY 是个负值 * 用于在 周状态下日期改变设置正确的y值 * * @return */ protected abstract float

    80430

    如何解决 Render Props 的回调地狱

    ---- 术语 “render prop” 是指一种在 React 组件之间使用一个值为函数的 prop 共享代码的简单技术 简而言之,只要一个组件中某个属性的值是函数,那么就可以说该组件使用了 Render...那有什么方法可以避免这种方式的修改呢? 在原生 JS 中,如果咱们调用函数后,还要做些骚操作,咱们一般使用回调函数来处理这种情况。...,并不是只有在 render 属性中传入函数才能叫 Render Props,实际上任何属性只要它的值是函数,都可称之为 Render Props,比如上面这个例子把 render 属性名改成 children...这样的组件设计更容易理解,因为渲染逻辑封装在一个单独的方法中。 如果需要更多嵌套,类的方式是垂直增加(通过添加新方法),而不是水平(通过相互嵌套函数),回调地狱问题消失。...实用的方法 如果想要在如何处理render props回调方面具有更大的灵活性,那么使用React-adopt是一个不错的选择。

    92220

    气象绘图——复杂的三维图

    修改普通的三维图固定设置 在普通的matplotlib的三维投影中,我们似乎并不能获得我们想要的结果,尤其是视觉上的,虽然倾斜了图形,但是文字等标注仍然是二维的,例如下面这张图片: 明显,无论是刻度标记...'k',lw=1,zorder=5) ax.plot([130,130],[10,50],[1000,1000],color='k',lw=1,zorder=5) ax.axis('off') 是不是具有三维视感而不是伪劣的窜入效果了呢...在三维图中实现栅格可视化 在之前的推文三维图形迁移中,我们已经介绍了如何使用收集collection的办法,来实现贴瓷砖式的数据可视化,这里我们仍然使用这种办法。...从上图不难看出,在东经90°剖面上,相对湿度在北纬25°以北突然降低到非常小的范围;在北纬50°剖面上,相对湿度在东经105°-115°有一个深厚的湿区。...=RH['r'].loc[1000:300,50:10,90:130] lev=RH['isobaricInhPa'].loc[1000:300] 随后,类似于二维平面的网格化,对三维坐标的lon,lat

    1.1K11

    腾讯位置服务开发应用-使用教程,案例分享,知识总结

    使用map地图组件开发,地图组件用于展示地图(使用的时腾讯地图),说一下它的属性有: longitude(类型为Number,没有默认值,表示中心经度) latitude(类型为Number,没有默认值...,比如说是 750rpx,不要设置百分比值,在uni-app中只支持gcj02坐标。...**介绍markers属性-类型为数组Array** 由之前描述可知,markers属性类型为Array,表示标记点用于在地图上显示标记的位置。...- 可识别换行符 anchor,经纬度在标注图标的锚点,默认底边中点,不是必填,{x, y},x表示横向(0-1),y表示竖向(0-1)。...在微信小程序中,当用户离开应用后,此接口无法调用,除非申请后台持续定位权限;当用户点击“显示在聊天顶部”时,此接口可继续调用。

    3K40

    腾讯位置服务开发应用-使用教程,案例分享,知识总结

    使用map地图组件开发,地图组件用于展示地图(使用的时腾讯地图),说一下它的属性有: longitude(类型为Number,没有默认值,表示中心经度) latitude(类型为Number,没有默认值...我们在写map组件时,组件的宽/高推荐写直接量,比如说是 750rpx,不要设置百分比值,在uni-app中只支持gcj02坐标。...介绍markers属性-类型为数组Array 由之前描述可知,markers属性类型为Array,表示标记点用于在地图上显示标记的位置。...,不是必填 - 可识别换行符 anchor,经纬度在标注图标的锚点,默认底边中点,不是必填,{x, y},x表示横向(0-1),y表示竖向(0-1)。...在微信小程序中,当用户离开应用后,此接口无法调用,除非申请后台持续定位权限;当用户点击“显示在聊天顶部”时,此接口可继续调用。

    6.4K51

    雷达系列 | 如何绘制极坐标下的雷达数据

    ,你必须知道雷达的坐标、方位角与库长 项目方法 azimuth_range_to_lat_lon 是 MetPy 库中的一个函数,用于将极坐标系统中的方位角和距离位置转换为经纬度坐标。...该函数的参数和返回值如下: 参数 azimuths (array-like):定义网格的一系列方位角。如果这不是一个 pint.Quantity 对象,则假定单位是度。...如果为 None,则使用默认的球形椭球体。 返回值 lon, lat (2D arrays):与原始位置相对应的经度和纬度二维数组。...在实际应用中,你可能需要安装 MetPy 和 pyproj 这两个库来使用这个功能。...如果你的数据不是以 pint.Quantity 的形式提供,确保它们是以正确的单位(例如,方位角为度,距离为米)给出的。

    9511

    关系运算符

    条件逻辑运算符是在类型logical和上定义的null。如果操作数值不是这些类型,"Expression.Error"则会引发带有原因代码的错误。 结果是一个逻辑值。...添加持续时间和某种类型的datetime值时,结果值具有相同的类型。 对于表中所列值以外的其他值组合,将"Expression.Error"引发带有原因代码的错误。以下各节介绍了每种组合。...下表列出了非零有限值、零、无穷大和 NaN 的所有可能组合的结果。在表中,x和y是非零有限值,z是 的结果x + y。如果x和y大小相同但符号相反,z则为正零。...如果x + y太大而无法在目标类型中表示,z则是与 具有相同符号的无穷大x + y。...Decimal Precision 中的总和是在不损失精度的情况下计算的。

    1K40

    基于腾讯地图定位组件实现周边公用厕所远近排序分布图

    在实际开发过程中,各有优劣。本次基于需求,使用腾讯位置服务作为一个公用厕所位置标注的H5页面开发。 本次使用版本: JavaScript API 2.0版本。...,默认值为10s; failTipFlag: 是否在定位失败时给出提示引导用户打开授权或打开定位开关。...2、WebService API地址解析(地址转坐标) 在项目完成测试后,如果遇到成千上百的地址时,一个一个的拾取,好像不是一个合格的开发者的所为。...v=2.exp&key=YOUR_KEY"> 在开发过程中,默认会这样引入到前端文件。测试环境和生成环境一致,或者更换环境也是一直,不会出现问题的。...在使用坐标拾取器时,一定要选择各个对应的工具,导航等牵涉到坐标的地方一定要注意。

    1.5K71

    用可视化地图讲照片的故事(Python+Leaflet)

    手机和数码相机拍的照片里除了我们能看到的RGB像元数据,还包含了拍摄时间、图像分辨率、感光值、GPS坐标等属性,记录在Exif(Exchangeable image file format)模块里。...随着手机像素越来越高,用手机记录身边的事(和自拍)已经变成很自然的动作,在一年里我们的手机肯定存了很多照片,照片和Exif数据块中的位置可以做哪些有趣的事情?...Exif块里,Exif信息以0xFFE1作为开头标记,采用TIFF格式,可以自己解析或直接用轮子exifread库,exifread是一个很方便使用的读取tiff和jpeg格式图片的Python库,在pypi...基于1中提取的坐标,保存为js文件,然后在浏览器打开html文件,就是上图中的效果了。另外需要说明的是,这些标记点(marker)点击之后都是能看到具体的文本的,展示的文本就是title里的内容。...',效果如下,是不是也很生动呢?

    2.3K30

    wrf-python 详解之如何使用

    wrf.to_np 函数按照以下流程执行: 如果没有缺省值或填充值,那么将直接调用 xarray.DataArray.values 属性返回值 如果有缺省值或填充值,那么会用 xarray.DataArray.attrs...属性 _FillValue 值替代 NaN 并返回 numpy.ma.MaskedArry # 获取3D对流有效位能(包含缺省值) cape_3d = getvar(ncfile, "cape_3d"...为了在输出数组中包含所有文件中的所有时间,设置 timeidx 参数为 wrf.ALL_TIMES(或设置为 None)。...当有多个文件并且每个文件具有多个时间时,如果最后一个文件的时间数少于之前文件的时间数,那么剩余的数组将用缺省值填充。...如果指定值的话,那么从每个文件中提取变量时,指定值将应用于每个文件。在具有多个时刻的多个文件中,这样做可能是没有意义的,因为每个文件的第 n 个索引可能表示不同的时刻。

    20.9K1012

    用可视化地图讲照片的故事(Python+Leaflet)

    本文转载自蛰虫始航 手机和数码相机拍的照片里除了我们能看到的RGB像元数据,还包含了拍摄时间、图像分辨率、感光值、GPS坐标等属性,记录在Exif(Exchangeable image file format...随着手机像素越来越高,用手机记录身边的事(和自拍)已经变成很自然的动作,在一年里我们的手机肯定存了很多照片,照片和Exif数据块中的位置可以做哪些有趣的事情?...基于1中提取的坐标,保存为js文件,然后在浏览器打开html文件,就是上图中的效果了。另外需要说明的是,这些标记点(marker)点击之后都是能看到具体的文本的,展示的文本就是title里的内容。...地图故事效果图 在html里可以根据自己的想法增加更多的内容,例如具体的地址文本,只需要调用百度/高德地图的Web服务 API中的逆地理编码服务就可以实现,逆地理编码就是指将经纬度转换为详细结构化的地址...',效果如下,是不是也很生动呢?

    2K20

    如何使用AngularJS和PHP为任何位置生成短而独特的数字地址

    此服务器应具有具有权限且配置了防火墙的具有sudo权限的非root用户。要进行此设置,没有服务器的同学可以在这里购买,不过我个人更推荐您使用免费的腾讯云开发者实验室进行试验,学会安装后再购买服务器。...第一列中,使用KEY命令来将digitaladdress编入索引。MySQL中的索引功能与它们在百科全书或其他参考工作中的工作方式类似。...您或您的应用程序发出包含的查询时WHERE声明,MySQL逐行读取每列中的每个条目,这可能成为一个资源密集程度极高的过程,因为您的表累积了越来越多的条目。...这样做的原因是为了加快启动应用程序运行的过程,而不是引导您完成创建每个文件并自行添加所有代码的过程。...如果您再次在浏览器中访问该应用程序,则不会看到其外观或行为的任何新变化。同样,如果您要输入地址并单击“ 生成”按钮,则应用程序仍然不会生成或显示地图代码。

    13.2K20

    google maps api_js调用谷歌浏览器接口

    注重得到的是地图类型GMapType对象,而不是对象名称 3.设置地图: setCenter(center, zoom?, type?)...() – southWest.lat();//维度差 //实例化一个在本张地图中的一个经纬度 var point = new GLatLng(southWest.lat() +...13.setImage(url) none 请求将 url 指定的图像设置为此标记的前景图。注:不调整打印图像和阴影图像。因此,此方法主要是为了达到高亮显示或变暗显示效果,而不是彻底改变标记的外观。...(自 2.46 开始) 2.getVertex(index) GLatLng 返回折线中给定索引值对应的顶点。...注重:在 v2 中,它不再是以地理坐标表示的地面上的一个点。现在,地理坐标可以用 GLatLng 表示。 在地图坐标系统中,x 坐标向右增大,y 坐标向下增大。

    5.7K10
    领券