我是靠谱客的博主 悦耳皮带,这篇文章主要介绍IDL:读栅格数据出现-Inf值,现在分享给大家,希望可以做个参考。

栅格数据中存有明确的负值,-3.0607,但是读数据出现错误,这个时候只需要按照数值设置范围提取就行

代码

pro cal_city
 ;计算城市指数
 ;读取数据
 cd,'G:城市发展指数计算1992'
 fns=file_search('*.tif',count=fnums)
 ;print,fnums
 for k=0,fnums-1 do begin
   img=read_tiff(fns[k])
   sz=size(img)
   nb=sz[1]
   ns=sz[2]
   w=where(img gt 0,count)
   ma=max(img)
   sum=0.0
   for i=0,nb-1 do begin
     for j=0,ns-1 do begin
       if img[i,j] gt -4 then begin
        sum=sum+img[i,j]
       endif
     endfor
   endfor
   I=sum/(count*ma)
   print,I
 endfor
end

其中设置数值大于-4即可

最后

以上就是悦耳皮带最近收集整理的关于IDL:读栅格数据出现-Inf值的全部内容,更多相关IDL内容请搜索靠谱客的其他文章。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(68)

评论列表共有 0 条评论

立即
投稿
返回
顶部