24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zzahkj

木虫 (著名写手)

龙

应助: 42 (小学生)
金币: 1923.7
散金: 2554
红花: 30
帖子: 1517
在线: 465.1小时
虫号: 964751
注册: 2010-03-08
性别: GG
专业: 地理信息系统

[求助] 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$

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

回复此楼