IDL将0值设置为nodata

最新推荐文章于 2023-06-20 03:56:19 发布

转载最新推荐文章于 2023-06-20 03:56:19 发布 · 7.2k 阅读

博客主要讲述了在IDL中把0值设置为nodata的相关内容，但具体操作等细节未给出。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

pro set_nodata

GPP=fltarr(720,360)

for i=1982L,2010L do begin 
 input_image='F:\P_SM_SPEI_GPP\Jung\GPP_JUNG\GPP_year\year\'+strtrim(string(i),02)+'.tif'
 GPP_image = read_tiff(input_image,Geotiff=Geotiff)
  for line=0,719 do begin
   for piex=0,359 do begin
   if (GPP_image[line,piex] le 0) then begin
    GPP[line,piex]=!VALUES.F_NAN
   endif else begin 
    GPP[line,piex]=GPP_image[line,piex]
   endelse
   endfor
  endfor
  output_file='F:\P_SM_SPEI_GPP\Jung\GPP_JUNG\GPP_year\year_nodata\'+strtrim(string(i),02)+'.tif'
  openw,lun,output_file,/get_lun
  write_tiff,output_file,GPP,Geotiff=Geotiff,/FLOAT
  free_lun,lun
  print,  i, '      OK!'
endfor

end