应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 地学 » 3S技术 » IDL批量求取某坐标点DN值
51/1返回列表
查看: 2303  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zzahkj

木虫 (著名写手)

龙

[求助] IDL批量求取某坐标点DN值 已有3人参与

我的数据如下：
TIFF格式——含有很多波段
我自己试着写的代码如下
CODE:
;批量读取点的DN值
pro txtpoints
  envi,/restore_base_save_files
  imges=dialog_pickfile(/read,/multiple_files,filter='*.tif',title='please choose your images',path='C:\Users\BDQS\Desktop\zz');读取影像文件路径
  n= N_ELEMENTS(imges);打开影像的个数
  FOR i=0,n-1 do begin;对每个影像，读取指定坐标下的DN值
    envi_open_data_file,imges(i),/IMAGINE,r_fid=fid;
    ENVI_FILE_QUERY, fid, dims=dims, nb=nb ,nl=nl,ns=ns,FNAME=fname;读取影像基本信息，获取行列数
     for curBand =0,nb-1 do begin
      xmap=84.4
      ymap=43.266667   
      ENVI_CONVERT_FILE_COORDINATES,fid,xf,yf,xmap,ymap
      yf=nl-round(yf)
      xf=round(xf)-1
      dims=[-1,xf,xf,yf,yf]
      pointArray(curBand).fname=fname
      pointArray(curBand).data = ENVI_GET_DATA(fid=fid, dims=dims, pos=curBand)
     endfor     
  ENDFOR
  outfile='E:\zz.txt'
  openw,lun,outfile,/get_lun,width=9999999
  printf,lun,pointArray,format='(17a)'
  free_lun,lun
  print,'completed

end

运行出现下列错误
txtpoints

      pointArray(curBand).data = ENVI_GET_DATA(fid=fid, dims=dims, pos=curBand)
                                                   ^
% Syntax error.
  At: C:\Users\BDQS\Desktop\txtpoints.pro, Line 17
% 1 Compilation error(s) in module TXTPOINTS.
% Attempt to call undefined procedure/function: 'TXTPOINTS'.
% Execution halted at: $MAIN$      

代码写得好像也不对，哪位大神帮我改改？非常感谢！
回复此楼

» 猜你喜欢
不要迷恋哥，哥只是一个过客！
1楼 2014-07-02 22:53:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Liwenkai

铜虫 (初入文坛)

【答案】应助回帖

求助一些seadas的资料，网上好难找
赞一下 回复此楼
3楼2014-07-03 19:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

Liwenkai

铜虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
非得IDL吗，arcgis也可以啊
赞一下 回复此楼
2楼2014-07-03 19:00:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzahkj

木虫 (著名写手)

龙
引用回帖:
2楼: Originally posted by Liwenkai at 2014-07-03 19:00:08
非得IDL吗，arcgis也可以啊

ArcGIS如何批量 一个TIFF文件有365个波段，多个TIFF文件
赞一下 回复此楼
不要迷恋哥，哥只是一个过客！
4楼2014-07-03 20:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzahkj

木虫 (著名写手)

龙
引用回帖:
3楼: Originally posted by Liwenkai at 2014-07-03 19:02:32
求助一些seadas的资料，网上好难找

海洋遥感软件？你做海洋遥感？
赞一下(1人) 回复此楼
不要迷恋哥，哥只是一个过客！
5楼2014-07-03 20:56:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭