pro add_data_sst

  ;实现将SST增加1度，再将结果更新到SST中

  ;1. 获取SST索引

  ;2. 通过索引获取ID

  ;3. 获取SST数组

  ;4. 对数组进行操作

  ;5. 将数值进行更新

  cd, 'D:\McDelfino\dell-dev\IDL'

  file='D:\McDelfino\dell-dev\IDL\A20080200505006.L2_LAC'

  ;打开HDF文件，获取SD ID，注意使用关键字RDWR，可以读写HDF文件。

  sd_id=HDF_SD_START(file, /RDWR)

  ;获取SD中SST的索引

  sstIndex = HDF_SD_NAMETOINDEX(sd_id,'sst')

  ;获得SST所在索引的ID

  SSTSD_ID = HDF_SD_SELECT(sd_id, sstIndex)

  ;读取SST数据

  HDF_SD_GETDATA, SSTSD_ID, SSTData

  ;获取sst的实际值

  ;slope=0.005, intercept=0.0

  realSSTData = SSTData * 0.005 + 0.0

  ;处理SST数据

  tempSSTData = realSSTData+1

  ;从真实值返回

  newSSTData = (tempSSTData - 0.0)/0.005

  ;将新数据写出，更新原有数据

  HDF_SD_ADDDATA, SSTSD_ID, newSSTData

  ;关闭SD，关闭HDF文件

  HDF_SD_ENDACCESS, SSTSD_ID

  HDF_SD_END, sd_id

end

pro readandcreate

  ;实现数据的计算以及新建数据

  ;1. 读取经纬度数据

  ;2. 读取SST数据，并存到数组中

  ;3. 读取Rrs数据，并进行计算获取SD、TSM、CHL结果，存到数组中

  ;4. 新建HDF，将经纬度以及结果写入进去

  file='D:\McDelfino\dell-dev\IDL\A20080200505006.L2_LAC'

  hdfid=HDF_SD_START(file, /rdwr)

  ;读取数据中的lat数据

  ;通过latitude来获取相应的索引值

  index=HDF_SD_NAMETOINDEX(hdfid, 'latitude')

  ;通过索引值获取ID值

  varid=HDF_SD_SELECT(hdfid, index)

  ;通过ID值获取数组值

  HDF_SD_GETDATA, varid, latdata

  help,latdata

  ;读取数据中的lon数据

  index=hdf_sd_nametoindex(hdfid, 'longitude')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, londata

  ;读取数据中的sst数据，有数据的进行计算，mask的部分赋值为NaN

  sstdata1=dindgen(1354,2040)

  index=hdf_sd_nametoindex(hdfid, 'sst')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, sstdata

  for i=0, 1354-1 do begin

    for j=0, 2040-1 do begin

      if(sstdata[i,j] gt -32767) then begin

        sstdata1[i,j]=sstdata[i,j]*0.005

      endif else begin

        sstdata1[i,j]='nan'

      endelse

    endfor

  endfor

  ;读取数据中的Rrs_667数据，有数据的进行SD的计算，mask的部分赋值为NaN

  rrs667data1=dindgen(1354,2040)

  sddata1=dindgen(1354,2040)

  index=hdf_sd_nametoindex(hdfid, 'Rrs_667')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, rrs667data

  for i=0, 1354-1 do begin

    for j=0, 2040-1 do begin

      if(rrs667data[i,j] gt -32767) then begin

        rrs667data1[i,j]=rrs667data[i,j]*(2.0E-6) + 0.05

        ;只对正数进行计算

        if(rrs667data1[i,j] gt 0) then begin

          sddata1[i,j]=-85.666-108.807*alog10(rrs667data1[i,j])-45.591*(alog10(rrs667data1[i,j]))^2.0-6.426*(alog10(rrs667data1[i,j]))^3.0

        endif else begin

          sddata1[i,j]='nan'

        endelse

      endif else begin

        sddata1[i,j]='nan'

      endelse

    endfor

  endfor

  ;读取数据中的Rrs_645数据，有数据的进行TSM的计算，mask的部分赋值为NaN

  rrs645data1=dindgen(1354,2040)

  tsmdata1=dindgen(1354,2040)

  index=hdf_sd_nametoindex(hdfid, 'Rrs_645')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, rrs645data

  for i=0, 1354-1 do begin

    for j=0, 2040-1 do begin

      if(rrs645data[i,j] gt -32767) then begin

        rrs645data1[i,j]=rrs645data[i,j]*(2.0E-6) + 0.05

        ;只对正数进行计算

        if(rrs645data1[i,j] gt 0) then begin

          tsmdata1[i,j]=1.572+2172.239*rrs645data1[i,j]

        endif else begin

          tsmdata1[i,j]='nan'

        endelse

      endif else begin

        tsmdata1[i,j]='nan'

      endelse

    endfor

  endfor

  ;读取数据中的Rrs_488/Rrs_547数据，有数据的进行CHL的计算，mask的部分赋值为NaN

  rrs488data1=dindgen(1354,2040)

  rrs547data1=dindgen(1354,2040)

  chldata1=dindgen(1354,2040)

  index=hdf_sd_nametoindex(hdfid, 'Rrs_488')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, rrs488data

  index=hdf_sd_nametoindex(hdfid, 'Rrs_547')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, rrs547data

  for i=0, 1354-1 do begin

    for j=0, 2040-1 do begin

      if(rrs488data[i,j] gt -32767) then begin

        rrs488data1[i,j]=rrs488data[i,j]*(2.0E-6) + 0.05

        rrs547data1[i,j]=rrs547data[i,j]*(2.0E-6) + 0.05

        ;只对正数进行计算

        if(rrs488data1[i,j] gt 0 && rrs547data1[i,j] gt 0) then begin

          R=alog10(rrs488data1[i,j]/rrs547data1[i,j])

          chldata1[i,j]=10^(4.608-14.742*R+22.025*R^2-13.509*R^3)

        endif else begin

          chldata1[i,j]='nan'

        endelse

      endif else begin

        chldata1[i,j]='nan'

      endelse

    endfor

  endfor

  ;读取数据中的chlor_a数据，有数据的进行计算，mask的部分赋值为NaN

  chlordata1=dindgen(1354,2040)

  index=hdf_sd_nametoindex(hdfid, 'chlor_a')

  varid=hdf_sd_select(hdfid, index)

  hdf_sd_getdata, varid, chlordata

  for i=0, 1354-1 do begin

    for j=0, 2040-1 do begin

      if(chlordata[i,j] gt -32767) then begin

        chlordata1[i,j]=chlordata[i,j]*1.0

      endif else begin

        chlordata1[i,j]='nan'

      endelse

    endfor

  endfor

  hdf_sd_endaccess, varid

  hdf_sd_end, hdfid

  ;新建HDF的文件路径

  file='D:\McDelfino\dell-dev\IDL\sst23.hdf'

  ;以create的形式打开文件

  hdfid=HDF_SD_START(file, /create)

  ;新建latitude变量，在HDF上显示1354列，2040行，默认是浮点型数组

  sds_id=HDF_SD_CREATE(hdfid, 'latitude', [1354, 2040], /double)

  ;由于从HDF中读取的数据与实际数据存在一个中心对称的关系，因此通过两个reverse实现

  ;latdata就是从其他HDF中读取的数组，将数组的结果添加到变量latitude中去

  HDF_SD_ADDDATA, sds_id, reverse(reverse(latdata), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'longitude', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(londata), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'sst', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(sstdata1), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'sd', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(sddata1), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'tsm', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(tsmdata1), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'chl', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(chldata1), 2)

  sds_id=HDF_SD_CREATE(hdfid, 'chlor_a', [1354, 2040], /double)

  HDF_SD_ADDDATA, sds_id, reverse(reverse(chlordata1), 2)

  HDF_SD_ENDACCESS, sds_id

  HDF_SD_END, hdfid

end