请选择 进入手机版 | 继续访问电脑版

遥感集市论坛

 
登录  立即注册
查看: 1124|回复: 3

[数据预处理] netcdf在IDL中的处理

[复制链接]

8

主题

15

帖子

39

积分

RS学酥

Rank: 2

积分
39
发表于 2015-10-12 18:59:26 | 显示全部楼层 |阅读模式
PRO read_asc
   COMPILE_OPT idl2
   ENVI, /RESTORE_BASE_SAVE_FILES
   ENVI_BATCH_INIT, LOG_FILE = 'batch_log.txt'

   finame = 'M:\asc\cru_ts3.21.1991.2000.tmp.dat.nc'
   NCid = NCDF_OPEN(finame)
   NCinfo = NCDF_INQUIRE(NCid)
   print, NCinfo.Nvars
   ; 查询NC文件中的变量,以及变量的 att 信息
   FOR iVar = 0, NCinfo.Nvars-1 DO BEGIN
     Varinfo = NCDF_VARINQ(NCid, iVar)
     print, "Var Name: ", Varinfo.Name
     print, "Att Number: ", Varinfo.Natts
     FOR iAtt = 0, Varinfo.Natts-1 DO BEGIN
       AttName = NCDF_ATTNAME(NCid, iVar, iAtt)
       NCDF_ATTGET, NCid, iVar, AttName, Att
       print, AttName, ":  ", STRING(Att)
     ENDFOR
     print, ""
   ENDFOR
   ; 读取变量名为tmp的数据
   Dataid = NCDF_VARID(NCid,'tmp')
   NCDF_VARGET, NCid, Dataid, Data
   ;help, Data
   ;   data=(data_t gt 100)*(-999)+(data_t lt 100)*data
   tags={GEOGANGULARUNITSGEOKEY:9102,$
     GEOGCITATIONGEOKEY:'GCS_WGS_1984',$
     GEOGINVFLATTENINGGEOKEY:298.25722356300003,$
     GEOGRAPHICTYPEGEOKEY:4326,$
     GEOGSEMIMAJORAXISGEOKEY:6378137.0000000000,$
     GTMODELTYPEGEOKEY:2,$
     GTRASTERTYPEGEOKEY:1,$
     MODELPIXELSCALETAG:[0.50000000,0.50000000,0.00000000],$
     MODELTIEPOINTTAG:[0.00000000,0.00000000,0.00000000,-180.00000,90.000000,0.00000000]}

   date=strarr(120);文件是一个含有120波段的数组,建立一个string array,用date接收返回值
   index=0
    for i=0,9,1 do begin
     for j=0,11,1 do begin     
       date[index]=string(1991+i)+'.'+string(j+1)      
       out_name='M:\temp\'+'tmp'+date[index]+'.tif'
       out_name=strcompress(out_name,/REMOVE_ALL)
       data_1=data[*,*,index]
       data_2=(data[*,*,index] gt 100)*(-999)+(data[*,*,index] lt 100)*(data[*,*,index])
       data_3=data_2*float(data_2 ne (-999))/(data_2 ne (-999))
       write_tiff, out_name, rotate(data_3,7), /float, geotiff=tags;rotate()
       index=index+1
     endfor
   endfor
end

回复

使用道具 举报

0

主题

6

帖子

6

积分

RS学渣

Rank: 1

积分
6
发表于 2015-10-14 19:54:03 | 显示全部楼层
给你个顶
回复

使用道具 举报

0

主题

4

帖子

4

积分

RS学渣

Rank: 1

积分
4
发表于 2017-5-3 23:22:12 | 显示全部楼层
mark
回复

使用道具 举报

0

主题

35

帖子

35

积分

RS学酥

Rank: 2

积分
35
发表于 2017-12-29 17:15:05 | 显示全部楼层
谢谢楼主!
回复

使用道具 举报

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

1 人签到
将意见发给我们
  • 通过email将您的想法和建议发给我们

    编辑部:bbs@rscloudmart.com

联系我们
  • 论坛官方QQ群:128651389

    论坛客服QQ:2739549450

遥感集市服务热线
  • 400-674-7498

QQ|小黑屋|遥感集市论坛 ( 粤ICP备12003192号-2

GMT+8, 2018-1-22 06:39 , Processed in 0.103443 second(s), 30 queries, Gzip On.

Powered by Discuz! X3.1 Copyright
© 2001-2013 Comsenz Inc.    All Rights Reserved.

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