IDL 读取葵花8（Himawari-8）HSD数据

可以实现直接输入.DAT文件，读取获得定标完成的影像、时间、太阳位置的三个数组
inputfile为.DAT文件，例如HS_H08_20170623_0250_B01_FLDK_R10_S0110.DAT
resolution为HS_H08_20170623_0250_B01_FLDK_R后的数字，代表了分辨率
代码会自动对数据进行对应的定标（可见光波段定标为表观反射率，红外定标为亮温）
时间提前转为double型，防止后面麻烦
会将500m的波段采样到一千米，然后将数据按分辨率划分，合并

;输入文件夹（存储全部DAT文件）并进行预处理
Function H8_Preprocess,inputfolder=inputfolder,e=eCOMPILE_OPT idl2R10_data=make_array(11000,11000,4,/FLOAT)R20_data=make_array(5500,5500,12,/FLOAT)bandnames=['01','02','03','04','05','06','07','08','09','10','11','12','13','14','15','16']for i=0,15 do beginband_files=file_search(inputfolder,'*B'+bandnames[i]+'_FLDK*.DAT',/test_regular)if N_ELEMENTS(band_files) EQ 10 THEN BEGINband_data=[]obstime=[]sunpos=[]for j=0,9 do begininfile=band_files[j]segment=read_hsd(inputfile=infile[0])imgdata=segment['pixels']imgtime=segment['time']sun_pos=segment['sun_pos']band_data=[[band_data],[imgdata]]obstime=[[obstime],[imgtime]]sunpos=[[sunpos],[sun_pos]];销毁哈希表OBJ_DESTROY,segmentendforif i lt 2 then beginR10_data[*,*,i]=band_dataR10_obstime=obstimeR10_sunpos=sunposendif else if i eq 2 then begin       band_raster=ENVIRaster(band_data,URI=e.GetTemporaryFilename(CLEANUP_ON_EXIT='True'))band_raster.SavesampleRaster=ENVIResampleRaster(band_raster, DIMENSIONS=[11000,11000])band_data=sampleRaster.GetData()sampleRaster.Closeband_raster.CloseR10_data[*,*,i]=band_dataendif else if i lt 4 then beginR10_data[*,*,i]=band_dataendif else beginR20_data[*,*,i-4]=band_dataR20_obstime=obstimeR20_sunpos=sunposendelseendifendforoutdata=hash('R10',R10_data,'R10_time',R10_obstime,'R10_sunpos',R10_sunpos,'R20',R20_data,'R20_time',R20_obstime,'R20_sunpos',R20_sunpos)return,outdataEND
;读取HSD数据并定标为albedo/BT
Function read_hsd,inputfile=inputfileCOMPILE_OPT idl2spli=strsplit(file_basename(inputfile),'_');获取分辨率resolution=fix(strmid(file_basename(inputfile),spli[6]+1,2))if resolution eq 5 then begincols=22000rows=2200endif else if resolution eq 10 then begincols=11000rows=1100endif else if resolution eq 20 then begincols=5500rows=550endifopenr,h8_lun,inputfile,/get_lun;获取Observation start timepoint_lun,h8_lun,46imgtime=dblarr(1);R8readu,h8_lun,imgtime;获取Total header lengthpoint_lun,h8_lun,70header_length=uintarr(1);I4readu,h8_lun,header_length;获取影像point_lun,h8_lun,header_lengthimgdata=uintarr(cols,rows)readu,h8_lun,imgdata ;获取Sun's positionpoint_lun,h8_lun,510sun_pos=dblarr(3);R8readu,h8_lun,sun_pos;获取Band numberpoint_lun,h8_lun,601band=uintarr(1);I2readu,h8_lun,band  ;获取Gainpoint_lun,h8_lun,617Gain=dblarr(1);R8readu,h8_lun,Gain;获取Offsetpoint_lun,h8_lun,625Offset=dblarr(1);R8readu,h8_lun,Offset;计算radianceradiance=imgdata*Gain[0]+Offset[0]if band le 6 then begin;计算反射率;获取radiance to albedopoint_lun,h8_lun,633cc=dblarr(1);R8readu,h8_lun,cc;计算反射率albedo=radiance*cc[0]   ;返回值outdataoutdata=albedoendif else begin;计算亮温;获取Central wave lengthpoint_lun,h8_lun,603wv=dblarr(1);R8readu,h8_lun,wv;获取radiance to brightness temperature(c0)point_lun,h8_lun,633c0=dblarr(1);R8readu,h8_lun,c0;获取radiance to brightness temperature(c1)point_lun,h8_lun,641c1=dblarr(1);R8readu,h8_lun,c1;获取radiance to brightness temperature(c2)point_lun,h8_lun,649c2=dblarr(1);R8readu,h8_lun,c2;获取Speed of light (c)point_lun,h8_lun,681c=dblarr(1);R8readu,h8_lun,c    ;获取Planck constant (h)point_lun,h8_lun,689h=dblarr(1);R8readu,h8_lun,h;获取Boltzmann constant(k)point_lun,h8_lun,697k=dblarr(1);R8readu,h8_lun,k;计算亮温wv=wv[0]*1e-6rad=radiance*1e6Te=h[0]*c[0]/k[0]/wv[0]/(ALOG(2*h[0]*c[0]^2/(wv[0]^5*rad)+1))BT=c0[0]+c1[0]*Te+c2[0]*Te^2;返回值 outdataoutdata=BTendelsefree_lun,h8_lun ;返回：albedo/BT,时间,太阳坐标obstime=make_array(cols,rows,value=imgtime[0]*1.0d)sunpos=make_array(cols,rows,3,/float)sunpos[*,*,0]=make_array(cols,rows,value=sun_pos[0])sunpos[*,*,1]=make_array(cols,rows,value=sun_pos[1])sunpos[*,*,2]=make_array(cols,rows,value=sun_pos[2])out=hash('pixels',outdata,'time',obstime,'sun_pos',sunpos)return,outEND

IDL 读取葵花8（Himawari-8）HSD数据相关推荐

xml文件c语言读取函数,IDL读取XML文件
使用IDL读取RADARSAT-2的数据,需要用到lutSigma.xml文件中的定标常数来计算相关参量.本文需要提取lutSigma.xml中的offset和gains参数.使用IDL来读取xml文 ...
IDL读取TXT文件出现 READF: End of file encountered. Unit: l, File: H: 错误
IDL读取TXT文件出现 READF: End of file encountered. Unit: l, File: H: 错误 1.出现原因出现错误图在读取txt文件时出现这个错误,研究了一天 ...
php读取大文件某行内容,php读取大文件最后几行数据的实现代码
/* * 读取大文件最后几行数据 * by www.jbxue.com */ $file = $handledir.'/venocap.1'; $fp = fopen($file, "r&q ...
如何用matlab读取npz文件,Python Numpy中数据的常用的保存与读取方法
在经常性读取大量的数值文件时(比如深度学习训练数据),可以考虑现将数据存储为Numpy格式,然后直接使用Numpy去读取,速度相比为转化前快很多. 下面就常用的保存数据到二进制文件和保存数据到文本文件 ...
微信新版支持读取iPhone M7/M8协处理器运动数据与好友PK一下运动量吧
iPhone的创新是有目共睹的,Healthkit的推出预示着苹果进军健康领域,iPhone M7/M8协处理器可以收集和分析用户的健康数据,那么好的硬件自然不会被势在打造完整生态圈的微信给错过,这不 ...
如何用C语言编写wav读取函数,C++读取WAV音频文件的头部数据的实现方法
C++读取WAV音频文件的头部数据的实现方法前言: 在这里分享一下自己的心得,希望和大家一起分享技术,如果有什么不足,还请大家指正.写出这篇目的,就是希望大家一起成长,我也相信技术之间没有高低,只有 ...
python读取文件某一行-使用python读取.text文件特定行的数据方法
如何用python循环读取下面.txt文件中,用红括号标出来的数据呢? 首先,观察数据可知,不同行的第一个数据元素不一样,所以考虑直接用正则表达式. 再加上,对读和写文件的操作,就行了注:我用的是p ...
python读取文件多行内容-使用python读取.text文件特定行的数据方法
如何用python循环读取下面.txt文件中,用红括号标出来的数据呢? 首先,观察数据可知,不同行的第一个数据元素不一样,所以考虑直接用正则表达式. 再加上,对读和写文件的操作,就行了注:我用的是p ...
stm32 iic接口进入busy_STM32通过IIC接口读取JY61模块MPU6050陀螺仪芯片数据核心程序...
1 简述最近,想学角度融合算法在网上买一个JY61的模块.他们家的模块用起来还不错.模块分为串口通讯和IIC通讯的.串口读取数据他们家有例程,我就不说了.想分享给大家这个模块的IIC是怎么去读取MP ...
java通过poi读取excel中的日期类型数据或自定义类型日期
java通过poi读取excel中的日期类型数据或自定义类型日期 Java 读取Excel表格日期类型数据的时候,读出来的是这样的 12-十月-2019,而Excel中输入的是 2019/10/12 ...

IDL 读取葵花8（Himawari-8）HSD数据

IDL 读取葵花8（Himawari-8）HSD数据相关推荐

最新文章

热门文章