版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

zzahkj

木蟲(chóng) (著名寫(xiě)手)

龍

應(yīng)助: 42 (小學(xué)生)
金幣: 1931.7
散金: 2554
紅花: 30
帖子: 1522
在線: 465.8小時(shí)
蟲(chóng)號(hào): 964751
注冊(cè): 2010-03-08
性別: GG
專(zhuān)業(yè): 地理信息系統(tǒng)

[求助] IDL批量求取某坐標(biāo)點(diǎn)DN值已有3人參與

我的數(shù)據(jù)如下：
TIFF格式——含有很多波段
我自己試著寫(xiě)的代碼如下

CODE:

;批量讀取點(diǎn)的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);打開(kāi)影像的個(gè)數(shù)

  FOR i=0,n-1 do begin;對(duì)每個(gè)影像，讀取指定坐標(biāo)下的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;讀取影像基本信息，獲取行列數(shù)

     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

運(yùn)行出現(xiàn)下列錯(cuò)誤
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$

代碼寫(xiě)得好像也不對(duì)，哪位大神幫我改改？非常感謝！

回復(fù)此樓

» 猜你喜歡

曹建勁團(tuán)隊(duì)在各種介質(zhì)中發(fā)現(xiàn)含金屬微粒，發(fā)明微粒找礦技術(shù)，獲得8項(xiàng)國(guó)家發(fā)明專(zhuān)利授權(quán) 已經(jīng)有11人回復(fù)
首次在隱伏礦體上方發(fā)現(xiàn)金、硫酸鉛、硝酸鉛、三氧化鎢等含金屬納米微粒已經(jīng)有8人回復(fù)
地球物理學(xué)和空間物理學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有94人回復(fù)
求助：有沒(méi)有大神有流體包裹體計(jì)算軟件已經(jīng)有0人回復(fù)
282求調(diào)劑已經(jīng)有0人回復(fù)

不要迷戀哥，哥只是一個(gè)過(guò)客！

1樓 2014-07-02 22:53:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Liwenkai

銅蟲(chóng) (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 193.5
散金: 60
帖子: 32
在線: 12.9小時(shí)
蟲(chóng)號(hào): 2011460
注冊(cè): 2012-09-18
性別: GG
專(zhuān)業(yè): 海洋遙感

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

非得IDL嗎，arcgis也可以啊

贊一下

回復(fù)此樓

2樓2014-07-03 19:00:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Liwenkai

銅蟲(chóng) (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 193.5
散金: 60
帖子: 32
在線: 12.9小時(shí)
蟲(chóng)號(hào): 2011460
注冊(cè): 2012-09-18
性別: GG
專(zhuān)業(yè): 海洋遙感

【答案】應(yīng)助回帖

求助一些seadas的資料，網(wǎng)上好難找

贊一下

回復(fù)此樓

3樓2014-07-03 19:02:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zzahkj

木蟲(chóng) (著名寫(xiě)手)

龍

應(yīng)助: 42 (小學(xué)生)
金幣: 1931.7
散金: 2554
紅花: 30
帖子: 1522
在線: 465.8小時(shí)
蟲(chóng)號(hào): 964751
注冊(cè): 2010-03-08
性別: GG
專(zhuān)業(yè): 地理信息系統(tǒng)

引用回帖:

2樓: Originally posted by Liwenkai at 2014-07-03 19:00:08
非得IDL嗎，arcgis也可以啊

ArcGIS如何批量一個(gè)TIFF文件有365個(gè)波段，多個(gè)TIFF文件

贊一下

回復(fù)此樓

不要迷戀哥，哥只是一個(gè)過(guò)客！

4樓2014-07-03 20:55:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zzahkj

木蟲(chóng) (著名寫(xiě)手)

龍

應(yīng)助: 42 (小學(xué)生)
金幣: 1931.7
散金: 2554
紅花: 30
帖子: 1522
在線: 465.8小時(shí)
蟲(chóng)號(hào): 964751
注冊(cè): 2010-03-08
性別: GG
專(zhuān)業(yè): 地理信息系統(tǒng)

引用回帖:

3樓: Originally posted by Liwenkai at 2014-07-03 19:02:32
求助一些seadas的資料，網(wǎng)上好難找

海洋遙感軟件？你做海洋遙感？

贊一下(1人)

回復(fù)此樓

不要迷戀哥，哥只是一個(gè)過(guò)客！

5樓2014-07-03 20:56:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

瀟湘es

新蟲(chóng) (小有名氣)

應(yīng)助: 21 (小學(xué)生)
金幣: 1102.6
散金: 50
帖子: 151
在線: 57小時(shí)
蟲(chóng)號(hào): 2001190
注冊(cè): 2012-09-14
性別: GG
專(zhuān)業(yè): 遙感機(jī)理與方法

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

;不會(huì)，我這里有一個(gè)批處理模板，你看著修改下
;---------------------------------------------------------------
;+
;

escription:
; ENVI二次開(kāi)發(fā)的批處理模版
; 默認(rèn)為數(shù)據(jù)格式轉(zhuǎn)換為tiff格式
;
; Author: DYQ
;-
;析構(gòu)函數(shù)
PRO ENVI_BATCH_TEMPLATE_CLEANUP,tlb
  WIDGET_CONTROL,tlb,get_UValue = pState
  PTR_FREE,pState
END
;事件響應(yīng)函數(shù)
PRO ENVI_BATCH_TEMPLATE_EVENT,event
  COMPILE_OPT idl2
  WIDGET_CONTROL,event.TOP, get_UValue = pState

  ;關(guān)閉事件
  IF TAG_NAMES(event, /Structure_Name) EQ 'WIDGET_KILL_REQUEST' THEN BEGIN
;
status = DIALOG_MESSAGE('關(guān)閉?',/Question)
IF status EQ 'No' THEN RETURN
;銷(xiāo)毀指針
; PTR_FREE, pState
WIDGET_CONTROL, event.TOP,/Destroy
RETURN;
  ENDIF
  ;根據(jù)系統(tǒng)的uname進(jìn)行判斷點(diǎn)擊的組件
  uName = WIDGET_INFO(event.ID,/uName)
  ;
  CASE uname OF
;打開(kāi)文件
'open': BEGIN
   files = DIALOG_PICKFILE(/MULTIPLE_FILES, $
      title = !SYS_Title+' 打開(kāi)文件', $
      path = (*pState).ORIROOT)
   IF N_ELEMENTS(files) EQ 0 THEN RETURN
   ;設(shè)置顯示文件
   WIDGET_CONTROL, (*pState).WLIST, set_value = files
   (*pState).INPUTFILES = PTR_NEW(files)
   (*pState).ORIROOT = FILE_DIRNAME(files[0])
   ;重置進(jìn)度條進(jìn)度
   IDLITWDPROGRESSBAR_SETVALUE,(*pState).PRSBAR,0

END
;退出
'exit': BEGIN
   status = DIALOG_MESSAGE('關(guān)閉?',$
      title = !SYS_Title, $
      /Question)
   IF status EQ 'No' THEN RETURN
   WIDGET_CONTROL, event.TOP,/Destroy
END
;關(guān)于
'about': BEGIN
   void = DIALOG_MESSAGE(!SYS_Title+' V1.0'+STRING(13b)+'歡迎使用，問(wèn)題討論請(qǐng)去bbs.esrichina-bj.cn！' ,/information)
END
;
;路徑選擇按鈕
'filepathsele': BEGIN
   WIDGET_CONTROL, event.ID,get_value = value
   WIDGET_CONTROL,(*pState).WSELE, Sensitive= value
   WIDGET_CONTROL,(*pState).OUTPATH, Sensitive= value
END
;選擇輸出路徑
'selePath' : BEGIN
   outroot = DIALOG_PICKFILE(/dire,title = !SYS_Title)
   WIDGET_CONTROL,(*pState).OUTPATH,set_value = outRoot
END

;功能執(zhí)行
'execute': BEGIN
   ;獲取選擇的方法
   WIDGET_CONTROL,(*pState).BGROUP, get_Value = mValue
   IF PTR_VALID((*pState).INPUTFILES) EQ 0 THEN RETURN
   ;初始化ENVI
   ENVI, /restore_base_save_files
   ENVI_BATCH_INIT,/NO_Status_Window

   ;獲取文件名
   files = *((*pState).INPUTFILES)
   per = 100./N_ELEMENTS(files)
   ;判斷是否需要選擇路徑
   IF mValue NE 0 THEN BEGIN
      ;構(gòu)建輸出文件名
      WIDGET_CONTROL, (*pState).OUTPATH,get_value= outfiledir
      IF (outfiledir[0] EQ ' ') THEN  outfiledir = DIALOG_PICKFILE(/dire, title =!SYS_Title+' 輸出路徑')
   ENDIF  ELSE outfiledir = FILE_DIRNAME(files[0])

   FOR i=0,N_ELEMENTS(files)-1 DO BEGIN
      ;構(gòu)建輸出文件名
      fileName = FILE_BASENAME(files)
      pointPos = STRPOS(fileName,'.')
      ;查找文件名中點(diǎn)的位置
      IF pointPos[0] NE -1 THEN BEGIN
      fileName= STRMID(fileName,0,pointPos)
      ENDIF
      out_name = outfiledir+PATH_SEP()+fileName+'.tiff'

      ENVI_OPEN_FILE, files, r_fid=fid
      IF (fid EQ -1) THEN BEGIN
      tmp = DIALOG_MESSAGE(files+'文件讀取錯(cuò)誤',$
         title = !sys_title, /error)
      CONTINUE
      ENDIF
      ;文件信息
      ENVI_FILE_QUERY, fid, dims=dims, nb=nb,bnames = bnames
      ;設(shè)置tiff文件輸出參數(shù)
      ;如果波段小于3個(gè)
      IF nb LE 3 THEN bandList = INDGEN(nb)ELSE $
      bandList = [3,2,1]
      ;調(diào)用ENVI功能函數(shù)另存數(shù)據(jù)
      ENVI_OUTPUT_TO_EXTERNAL_FORMAT,fid = fid,dims = dims, out_name=out_name,pos = bandList, $
      out_bname=bnames[bandlist],/TIFF
      ;輸出完成
      ENVI_FILE_MNG, id=fid, /remove
      ;設(shè)置進(jìn)度條
      IDLITWDPROGRESSBAR_SETVALUE,(*pState).PRSBAR,(i+1)*per
   ENDFOR
   void = DIALOG_MESSAGE('處理完成 ',title = !sys_title,/infor)
   ;關(guān)閉ENVI二次開(kāi)發(fā)模式
   ENVI_BATCH_EXIT
END
ELSE:
  ENDCASE
END
;
;--------------------------
;ENVI二次開(kāi)發(fā)批處理模版
PRO ENVI_BATCH_TEMPLATE
  ;
  COMPILE_OPT idl2
  ;初始化組件大小
  sz = [600,400]
  ;設(shè)置系統(tǒng)變量，可方便修改系統(tǒng)標(biāo)題
  DEFSYSV,'!SYS_Title','ENVI批處理模版'
  ;創(chuàng)建界面的代碼
  tlb = WIDGET_BASE(MBAR= mBar, $
/COLUMN , $
title = !SYS_Title, $
/Tlb_Kill_Request_Events, $
tlb_frame_attr = 1, $
Map = 0)
  ;創(chuàng)建菜單
  fMenu = WIDGET_BUTTON(mBar, value ='文件',/Menu)
  wButton = WIDGET_BUTTON(fMenu,value ='打開(kāi)數(shù)據(jù)文件', $
uName = 'open')
  fExit = WIDGET_BUTTON(fMenu, value = '退出', $
uName = 'exit',/Sep)
  eMenu = WIDGET_BUTTON(mBar,value ='功能',/Menu)
  wButton = WIDGET_BUTTON(eMenu,$
value ='運(yùn)行批處理', $
uName = 'execute')
  hMenu =  WIDGET_BUTTON(mBar, value ='幫助',/Menu)
  hHelp = WIDGET_BUTTON(hmenu, value = '關(guān)于', $
uName = 'about',/Sep)
  ;上面的輸入base
  wInputBase = WIDGET_BASE(tlb, $
xSize =sz[0], $
/Frame, $
/Align_Center,$
/Column)

  wLabel= WIDGET_LABEL(wInputBase, $
value ='文件列表')
  wList = WIDGET_LIST(wInputBase, $
YSize = sz[1]/(2*15),$
XSize = sz[0]/8)

  ;輸出路徑設(shè)置
  wLabel= WIDGET_LABEL(tlb, $
value ='輸出參數(shù)設(shè)置')

  ;輸出參數(shù)控制界面
  wSetBase = WIDGET_BASE(tlb, $
xSize =sz[0], $
/Row)
  values = ['源文件路徑', $
'另選擇路徑']
  bgroup = CW_BGROUP(wSetBase, values, $
/ROW, /EXCLUSIVE, $
/No_Release, $
SET_VALUE=1, $
uName = 'filepathsele', $
/FRAME)
  outPath = WIDGET_TEXT(wSetBase, $
value =' ', $
xSize =30, $
/Editable, $
uName = 'outroot')
  wSele = WIDGET_BUTTON(wSetBase, $
value ='選擇路徑', $
uName ='selePath')
  ;
  ;執(zhí)行按鈕base
  wExecuteBase = WIDGET_BASE(tlb,$
/align_center,$
/row)
  wButton = WIDGET_BUTTON(wExecuteBase, $
ysize =40,$
value ='打開(kāi)數(shù)據(jù)文件', $
uName = 'open')
  wButton = WIDGET_BUTTON(wExecuteBase,$
value ='運(yùn)行批處理', $
uName = 'execute')
  ;狀態(tài)欄，僅顯示進(jìn)度條
  wStatus = WIDGET_BASE(tlb,/align_right)
  prsbar = IDLITWDPROGRESSBAR( wExecuteBase ,$
title ='進(jìn)度', $
CANCEL =0)
  ;結(jié)構(gòu)體傳遞參數(shù)
  state = {wButton:wButton, $
tlb : tlb, $
oriRoot: '', $
outPath: outPath, $
wSele : wSele, $
bgroup : bgroup , $
inputFiles : PTR_NEW(), $
prsbar : prsbar , $
wList : WLIST }

  pState = PTR_NEW(state,/no_copy)
  ;操作界面居中
  CENTERTLB,tlb
  ;
  WIDGET_CONTROL, tlb,/Realize,/map,set_uValue = pState
  XMANAGER,'ENVI_Batch_Template',tlb,/No_Block,$
cleanup ='ENVI_Batch_Template_Cleanup'
END

贊一下(1人)

回復(fù)此樓

6樓2014-07-04 11:36:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hyPang

銅蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1309.5
帖子: 136
在線: 58.4小時(shí)
蟲(chóng)號(hào): 4206052
注冊(cè): 2015-11-08
性別: GG
專(zhuān)業(yè): 遙感機(jī)理與方法

【答案】應(yīng)助回帖

不知道你在
envi,/restore_base_save_files 之后加
envi_batch_init
會(huì)不會(huì)管事

贊一下

回復(fù)此樓

路雖遠(yuǎn)我用心永恒

7樓2015-11-21 14:03:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 zzahkj 的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求材料，環(huán)境專(zhuān)業(yè)調(diào)劑 +3	18567500178 2026-03-18	3/150	2026-03-23 23:50 by 熱情沙漠
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	11/550	2026-03-23 23:29 by 呆呆師姐
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 269專(zhuān)碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4	石瑞0426 2026-03-19	4/200	2026-03-22 08:50 by Wei_ren
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 材料工程專(zhuān)碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 材料工程（專(zhuān)）一志愿985 初試335求調(diào)劑 +3	hiloiy 2026-03-17	4/200	2026-03-21 03:04 by JourneyLucky
[考研] 二本跨考鄭大材料306英一數(shù)二 +3	z1z2z3879 2026-03-17	3/150	2026-03-21 02:29 by JourneyLucky
[考研] 085700資源與環(huán)境308求調(diào)劑 +12	墨墨漠 2026-03-18	13/650	2026-03-21 01:42 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿西南交通專(zhuān)碩材料355 本科雙非求調(diào)劑 +5	西南交通專(zhuān)材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 295復(fù)試調(diào)劑 +8	簡(jiǎn)木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing

24小時(shí)熱門(mén)版塊排行榜

[求助] IDL批量求取某坐標(biāo)點(diǎn)DN值 已有3人參與

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] IDL批量求取某坐標(biāo)點(diǎn)DN值已有3人參與