IDL_原函数分析3

pro test
PRINT, '----------------------------Read programme start!-------------------------'


;--------------------------------批量读取部分---------------------------------------------
;MyRootDir='H:\Q';
;filearr = file_search(MyRootDir,'*.dat',count=num);
;FOR fileindex=0,num-num,1 DO BEGIN
;  OPENR , LUN , filearr[fileindex] , /COMPRESS,/GET_LUN
;  filename=FILE_BASENAME(filearr[fileindex]);  ;获取文件名
;  PRINT,'filename is '+filename
;  FleArr = BYTARR(90,90)                     ;创建1440列720行3个字节的随机数组
;  READU , LUN , FleArr                            ;读LUN标识符的文件到IDL变量库
;  TV, FleArr
  ;用Contour函数可以绘制数据的等高线。在构造等高线之前要先确定每个轮廓线的高程值、标注的内容及宽度，根据首曲线和计曲线的概念，
;--------------------------------批量读取部分---------------------------------------------


;----------------------------参数计算部分-------------------------------------
;  for h=ceil(min(grid)),ceil(max(grid)),d do begin; 
;  ;grid为高程值，d为高程间距
;  hd=[hd,h]
;  ;每个轮廓线高程值
;    if h mod 5*d eq 0 then begin & hs=[hs,h] & hc=[hc,3]
;  ;计曲线的标注内容及宽度
;    endif else begin & hs=[hs,''] & hc=[hc,1]
;  ;首曲线的标注内容及宽度
;    endelse 
;  endfor
;  contour,grid,levels=hd,c_labels=hs,C_THICK=hc,c_color=['00ffff'xl];
;----------------------------参数计算部分-------------------------------------




;; Open the file for reading as file unit 1:
OPENR, 1, FILEPATH('SLA19930101.dat')
;; The image is a x by y byte array, so make B that size:
B = BYTARR(480, 480, /NOZERO)  ;这里填入数据的长和宽
;; Read the data into B:
READU, 1, B
;; Close the file:
CLOSE, 1
; Display the image:
TV, B


end