网上好像没有PIOMAS数据处理的文章,最近项目需要进行了处理,因此给大家提供以下参考,不过自己摸索一下应该很快就会了
PIOMAS数据下载
http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model_grid
在表格中的Sea ice thickness (Volume per unit Area) , daily mean即为我下载的海冰厚度数据
点击右侧的 hiday.H<yyyy>
跳转到https://pscfiles.apl.washington.edu/zhang/PIOMAS/data/v2.1/hiday/下载数据
我这里是下载了hiday.H2021
这个数据明显比其他数据小,因为其实只有1和2月的数据
PIOMAS数据处理
其实官网提供了读取数据的IDL程序,就是叫做read_draw_hi.pro这个程序,但是这个程序包含了制图,因此我觉得对我来说稍微有点复杂了。
我按照它程序提供的信息和思路,重新用IDL写了代码,因此以下代码就是能够将我上述下载的数据读取出来,存成每天一个波段的envi格式的数据(不带地理信息),我也将mask之外的数据存储成nan值了
除了上述的数据外,还需要下载grid.dat数据和io.dat_360_120.output数据,http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model_grid 这个网址都有。
以下程序应该还算是比较简明,值得注意的是一些数据的地址要修改,另外就是如果读取的数据改变导致波段改变,59要改成别的值。因为我处理的是北极的信息,因此nx和ny是这个值,但是如果不是的话,可以从read_draw_hi.pro这个程序再找出进行修改
之后我使用envi的Georeference from IGM建立地理信息。这也是以下代码需要保存lon和lat的原因。
之后可能会更新批量处理的代码,代码也会加上建立地理信息的步骤。
pro read_piomas_2021
nx=360 & ny=120;120 ;Arctic Parallel Ocean and sea Ice Model (POIM) grid
max_levels=30
desired_level =1
mask_value=99999.0
barwidth=200.0
lon=fltarr(nx,ny) & lat=fltarr(nx,ny)
data4=fltarr(nx,ny,59, /NOZERO);注意这里的59其实是指数据为1和2月数据一起59天,最后生成59各波段的数据
kmt=intarr(nx,ny)
;读取经纬度
openr, lunin, 'H:\mission\PIOMAS\grid.dat',/Get_lun ;打开数据文件
readf, lunin, lon, format='(10f8.2)'
readf, lunin, lat, format='(10f8.2)'
close,lunin
free_lun,lunin
;read grid mask (ocean levels)
openr,lunin,'H:\mission\PIOMAS\io.dat_360_120.output',/Get_lun ;打开数据文件
readf, lunin, kmt,format='(360i2)'
close,lunin
free_lun,lunin
;读取数据
openr,lunin, 'H:\mission\PIOMAS\hiday.H2021', /swap_if_big_endian,/Get_lun ;打开数据文件'H:\mission\LNG\PIOMAS\monthly\heff.H2021'
readu,lunin, data4
close,lunin
free_lun,lunin
location=where(kmt lt 1,count)
if count gt 0 then begin
numFiles=size(data4)
numFiles=numFiles[3]
FOR fidx=0, numFiles-1 DO BEGIN
temp=data4[*,*,fidx]
temp[location]=0.0/0
data4[*,*,fidx]=temp[*,*]
ENDFOR
endif
print,'count',count
;储存数据
openw,lunou,'H:\mission\PIOMAS\lon',/get_lun
writeu,lunou,lon;
close,lunou
free_lun,lunou
openw,lunou,'H:\mission\PIOMAS\lat',/get_lun
writeu,lunou,lat;
close,lunou
free_lun,lunou
openw,lunou,'H:\mission\PIOMAS\piomas2021',/get_lun
writeu,lunou,data4;
close,lunou
free_lun,lunou
print,'finished'
end
PIOMAS的一些详细信息,包括格网可以参考:http://psc.apl.uw.edu/research/projects/projections-of-an-ice-diminished-arctic-ocean/model/