爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10366|回复: 13

如何用matlab或者arcgis提取netcdf格式的经纬度为双精度的区域数据

[复制链接]

新浪微博达人勋

发表于 2015-1-22 11:30:48 | 显示全部楼层 |阅读模式
GrADS
系统平台:
问题截图:
问题概况: 如何用matlab或者arcgis提取netcdf格式的经纬度为双精度的区域数据
我看过提问的智慧: 看过
自己思考时长(天): 7

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

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

x
本帖最后由 musudsay 于 2015-1-22 11:48 编辑

我下载的RegCM的数据他本身具有地图投影,是兰伯特投影,然后这是他网站上的图。。。
94cad1c8a786c917d1546d91ca3d70cf3ac757c6.jpg
我下载下来的数据用matlab读出来发现lat,lon是双精度的,不知道该如何提取区域。。。。用matlab与ArcGIS都试过了。。。都不行
用matlab读取了一下发现,lat:241*195   lon:241*195  不能直接提取区域。。。
d0c8a786c9177f3e67ad70b673cf3bc79e3d56eb.jpg
用ArcGIS打开也是,他的经纬度数据都显示为1,2,3,。。。到241,lon显示为1,2,3,。。。。195。感觉好像是格网的编号。。。
QQ图片20150122111511.jpg
4ec2d5628535e5dd2b5a296475c6a7efce1b6213.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-23 13:39:17 | 显示全部楼层
musudsay 发表于 2015-1-23 12:49
matlab有地图重投影的功能吗?

MATLAB有没有不知道,不过MeteoInfo有哇
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2015-1-22 14:18:51 | 显示全部楼层
又是那隻貓 发表于 2015-1-22 12:19
先看下数据里面有没有经纬度信息,没有的话一般在说明文档中会给范围,根据给的范围自己写经纬度也不难吧

他的数据经纬度lat与lon我用matlab读出来了,但是都是241*195。没有办法getvar......提取区域。。。说明文档里面说它是模式的中心格点为Central Lat. and Lon. : 22.04°N, 118.96°E  。模式格网为197 (Lat.) × 243 (Lon.)。。。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-1-22 12:19:45 | 显示全部楼层

回帖奖励 +5 金钱

先看下数据里面有没有经纬度信息,没有的话一般在说明文档中会给范围,根据给的范围自己写经纬度也不难吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-22 14:31:39 | 显示全部楼层
musudsay 发表于 2015-1-22 14:18
他的数据经纬度lat与lon我用matlab读出来了,但是都是241*195。没有办法getvar......提取区域。。。说明 ...

那您得学习如何用matlab读取netcdf了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-22 14:42:17 | 显示全部楼层
又是那隻貓 发表于 2015-1-22 14:31
那您得学习如何用matlab读取netcdf了

不好意思我是刚刚开始学习netcdf,不是很懂。而且也查阅了网站的资料,发现没有办法解决我的问题。所以请求帮助。我知道当lat,与lon为单精度的时候可以直接用sfcwind=netcdf.getVar (ncid,vid,[,,,],[,,,]);但是我现在的lat与lon是双精度的, 是不是要写一个循环语句??
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2015-1-22 15:24:28 | 显示全部楼层
musudsay 发表于 2015-1-22 14:42
不好意思我是刚刚开始学习netcdf,不是很懂。而且也查阅了网站的资料,发现没有办法解决我的问题。所以请 ...

没有问题,你再试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-22 16:05:24 | 显示全部楼层
如果是想提取等经纬度的区域数据,需要向想原始投影(Lambert)的数据投影至等经纬度投影,然后按经纬度范围提取。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-22 16:30:52 | 显示全部楼层
显示的错误应该是因为电脑内存不足,你直接meshgrid的话,出来的X,Y都是(241*195)*(241*195 )的数据,电脑没有相应的内存是放不下的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-22 16:35:03 | 显示全部楼层
很有助于编程学习!多谢楼主!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-23 12:47:45 | 显示全部楼层
asmon 发表于 2015-1-22 16:30
显示的错误应该是因为电脑内存不足,你直接meshgrid的话,出来的X,Y都是(241*195)*(241*195 )的数据, ...

不是的,我是直接打开lat,lon他是241*195,我截图放错了,这个不需要再定义网格
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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