IDL编程学习之读取nc/NetCDF文件

15195775117

------第一部分：基础知识------
nc文件简介：
1、nc文件=netCDF文件∈二进制文件；
2、一个NetCDF数据集包含维(dimensions)、变量(variables)、属性(attributes)，每种类型都会被分配一个名字和一个ID；
3、变量存储实际数据，维给出了变量维度信息，属性则给出了单位等信息；
------第二部分：读取netCDF案例------
nc文件见附件
;文件路径：
  file='C:\Users\Administrator\Desktop\BPLH.nc'
;获取文件ID：
  file_id = ncdf_open(file, /nowrite )
;下面需要获取结构体变量名，根据名字把数据读出来
;返回值file_info是个结构体，给出了维度数、变量数，根据变量数我们进行循环：
  file_info = ncdf_inquire( file_id )
  FOR varid=0, file_info.nvars-1 DO BEGIN
    var = ncdf_varinq( file_id, varid )
;看一下文件中的变量名字
    help,var.name
    lonid = ncdf_varid(file_ID,var.name)
    ncdf_varget, file_id, lonid,mid_var
;把文件中每个量，赋给相同名字的变量
    (scope_varfetch(var.name,/enter))=mid_var
;看看我转了之后的变量名，跟原始名字应该是一样的
    print,var.name
  ENDFOR
;在知道变量名后，help下，这样数据就读出来了，撒花！
  help,longitude,latitude,time,blh
;我们知道了：
LONGITUDE       FLOAT     = Array[253]
LATITUDE        FLOAT     = Array[205]
TIME            LONG      = Array[248]

BLH             INT       = Array[253, 205, 248]

15195775117

netCDF中，不能直接用“结构体名.成员名”取数据
数据和变量名、数据类型等放在一个结构体中
首先，我们用ncdf_inquire( file_id)了解有几个变量
这个函数的返回值是个结构体：
   NDIMS           LONG                 3
   NVARS           LONG                 4
   NGATTS          LONG                 2
   RECDIM          LONG                 2
它告诉我们，数据是三维的，有4个变量
(后面俩我也不知道有什么用，遇到了再说)
根据上楼的例子，NVars似乎是包含NDims的，
所以我们想获取所有变量，只要看nvars即可；
对于第i个变量，获取该变量的结构体（这个变量还是以结构体形式存的）

    var = ncdf_varinq( file_id,i)
返回值是个结构体，我们要的是变量名字的字符串
    var_id = ncdf_varid(file_ID,var.name)
我做个中间值mid_var用来装变量的数据
    ncdf_varget, file_id, var_id,mid_var
    (scope_varfetch(var.name,/enter))=mid_var

15195775117

根据一楼读取的数据绘图，叠加地图，方式有3种：
1、用icontour画图，另存为kmz格式；
2、用icontour画图，另存为kmz格式；
3、用icontour画图，用insert下拉菜单加载国界河流，
可惜没有省界，只有美国和加拿大可以加载州界和省界；

午后阳光

谢谢楼主分享，学习中

15195775117

午后阳光 发表于 2018-9-7 16:27
谢谢楼主分享，学习中

有什么资料啦、问题啦，都可以分享出来，很多进步都是在解决问题中获得的

15195775117

原来，IDL有nc文件的读取导航界面

SonGoku

请问楼主，可以自定义叠加地图吗，比如我只要青藏高原地区的背景地图，不要大洲大陆国家界线等

15195775117

SonGoku 发表于 2019-3-21 14:35
请问楼主，可以自定义叠加地图吗，比如我只要青藏高原地区的背景地图，不要大洲大陆国家界线等

可以直接铺在谷歌地球上，或者谷歌地图，或者百度地图，或者高德地图.........
参考我的这个帖子：
http://bbs.06climate.com/forum.p ... p;extra=&page=1

cgb

楼主，可否写下代码：如何把nc文件中许多年的数据输出为tiff格式？

15195775117

cgb 发表于 2019-10-8 16:38
楼主，可否写下代码：如何把nc文件中许多年的数据输出为tiff格式？

write_tiff函数即可输出，不知具体数据是什么样的