请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4990|回复: 9

[其他] (焦急求助)如何提取nc文件的经纬度?

[复制链接]

新浪微博达人勋

发表于 2022-11-16 16:01:13 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 顺利毕业! 于 2022-11-17 16:47 编辑

一个nc数据文件(全球格点)中,lat(维度)有96个值,lon(经度)有144个值,现在想利用ncl将经纬度提取成一个txt文件(或者excel文件),共两列,第一列为纬度数据,第二列为经度数据,并且一个纬度下分别对应144个经度,即输出文件为96*144=13824行,2列的数据。请问各位大佬,这应该如何实现呀?如果能提供相应脚本就更感谢啦!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-11-17 08:02:33 | 显示全部楼层
用 conform 函数 将lon和 lat 扩展成两个 array(96,144)
之后 用 ndtoond  将  array(96,144) 变为1维
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-17 09:13:33 | 显示全部楼层
不会ncl,提供一个python的,很简单
  1. import xarray as xr
  2. # --------------------------------------------------------

  3. nc_filein = "era5.temperature.20200901.nc"

  4. ds = xr.open_dataset(nc_filein)
  5. lngs = ds.longitude.data
  6. lats = ds.latitude.data
  7. with open("lnglat_output.txt", "w") as f:
  8.     for lat in lats:
  9.         for lng in lngs:
  10.             f.write(f"{lat},{lng}\n")
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-17 12:12:38 | 显示全部楼层
betasy 发表于 2022-11-17 09:13
不会ncl,提供一个python的,很简单

谢谢,目前还不会python 呜呜呜
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-17 12:15:31 | 显示全部楼层
houtian 发表于 2022-11-17 08:02
用 conform 函数 将lon和 lat 扩展成两个 array(96,144)
之后 用 ndtoond  将  array(96,144) 变为1 ...

这两个array用ndtoond会分别变成13824个经度或者纬度吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-17 18:16:59 | 显示全部楼层
betasy 发表于 2022-11-17 09:13
不会ncl,提供一个python的,很简单

我在conda(miniconda)下安装了python3.9,运行该脚本后,报错如下:
   line 1, in <module>
    import xarray as xr
ModuleNotFoundError: No module named 'xarray'
请问这是怎么回事?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-17 18:28:11 来自手机 | 显示全部楼层
顺利毕业! 发表于 2022-11-17 18:16
我在conda(miniconda)下安装了python3.9,运行该脚本后,报错如下:
&nbsp; &nbsp;line 1, in &lt;module&gt;
&nbsp; &nbsp; import xarray as xr
ModuleNotFoundError: No module named 'xarray'
请问这是怎么回事?

没有库,按网上的教程装一下库就行
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-17 19:54:28 | 显示全部楼层
小余 发表于 2022-11-17 18:28
没有库,按网上的教程装一下库就行

好的,谢谢啦~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-11-21 04:32:15 | 显示全部楼层
顺利毕业! 发表于 2022-11-17 12:15
这两个array用ndtoond会分别变成13824个经度或者纬度吗?

lz 注意多看NCL官网,对于所用的函数,在NCL官网上查查相应说明:http://ncl.ucar.edu/Document/Functions/list_alpha.shtml
另外有个建议:学习程序语言,不理解的用法,可以在ncl命令行下,输入输出看看结果。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-11-21 10:33:10 | 显示全部楼层
houtian 发表于 2022-11-21 04:32
lz 注意多看NCL官网,对于所用的函数,在NCL官网上查查相应说明:http://ncl.ucar.edu/Document/Function ...

好的 谢谢!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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