2021-10-29

IDL的学习总结

spdsw包的使用

1. load数据

直接用IDL下载mms数据(时间较长，建议在网站直接下载，但是文件夹一直要按照mms的文件顺序建立)
mms_load_mec, probes=probe, trange=trange, local_data_dir = ‘G:\mms’
mms_load_fgm, trange = trange, probes = probe, data_rate = [‘srvy’,‘brst’], local_data_dir = ‘G:\mms’
mms_load_fpi, trange = trange, data_rate=‘fast’, level=‘l2’, probes=probe, datatype=[‘des-moms’,‘dis-moms’], local_data_dir = ‘G:\mms’
其中mec:卫星数据
fgm:磁场数据
fpi:粒子数据
scm:波动的数据，扰动数据，burst-mode 扰动磁场数据 分析哨声波的

2. 加载数据
   file = filepath(‘mms1_fgm_srvy_l2_20150901_v4.18.0.cdf’,ROOT_DIR=‘F:\mms_data\data\mms\mms1\fgm\srvy\l2\2015\09’)
   cdf2tplot,file,varformat=’*’,/all,/get_support_data

3. 卫星时间转换
   t1 = time_double(‘2015-09-01/12:15:28’) ;将世界时（string）转化为卫星上面的时间（double）
   t2 = time_double(‘2015-09-01/12:24:00’)
   time_string(a.x[1095263]) ;将卫星中的时间变为世界时，即地球上的时间

4. 数据的查找、储存
   pos_t = where(a.x ge t1 and a.x le t2) ;where函数 用于查找指定范围的数据
   store_data,‘shi’,data={x:a.x[pos_t],y:a.y[pos_t,*]}；将数据选定范围存储进选定的变量
   get_data,‘mms1_fgm_b_dmpa_srvy_l2_bvec’,data=a ;提取文件中的数据，并将数据存储到变量a中
   tplot_names ;用于查看文件的名称

5. 画图
   tplot,[‘mms1_fgm_b_bcs_srvy_l2_btot’,‘mms1_fgm_b_bcs_srvy_l2_bvec’] ;同时画出不同的数据
   tlimit,[‘2015-09-01/12:15:28’,‘2015-09-01/12:24:00’] ;限制画图区间的函数并绘画出来
   限制画图区间
   options, ‘mms_scm_fac_powspec’, zlog=1 ，￥
   （z轴设置为对数轴）, ylog=1（y轴设置为对数轴）, yrange=yrange（y轴范围）, ytitle=‘f (Hz)’ （y轴标题）, ztitle=‘B!ds!n (nT!u2!n/Hz)’ （z轴标题）, zrange=[1e-8,1e-4]

6. 将磁场数据转换为场向坐标系之后，将频谱画出来
   twavpol,‘shi’ $
   ,nopfft=128,steplength=128

第一步是load数据，第二步将数据转到场向坐标系（fac_matrix_make, tvector_rotate），第三步做频谱分析（twavpol）
背景磁场可以先不用滤出来，滤波是 output_data=time_domain_filter(input_data,f1,f2)，带通滤波 滤除的信号是在[f1, f2]范围

7. 带感叹号的一般都是IDL系统的变量，所以其实这个只要知道应该就够了。比如你想把某个值去掉，那就改成这个值就行了，比如我想去除的部分，那我就赋值。
   a = !VALUES.F_NAN ；NAN
   IDL> b = -!VALUES.F_NAN ;-NAN
   IDL> c = !VALUES.F_INFINITY ;INF
   IDL> d = -!VALUES.F_INFINITY ;-INF