菜鸟速成记-IDL读取显示dat文件

DAT文件的读取

问题描述

之前，博主没有理解dat文件导致对文件获取产生问题和误解，还可笑的认为没办法获取demission的话，是没有办法读取dat文件的。

由于dat文件是遥感影像的常用格式之一。所以，首先认识一番，dat文件。

1. 先找个文件夹看看它在哪
   

2.看看它们到底长什么样子

• hdr 头文件
  

• dat 文件
  

3.思考怎么去读取dat文件

读取dat文件的思路:
- 首先获取到dat文件的头文件信息。

• 调用函数进行显示READ_BINARY()

代码实现过程

由于代码实现过于简单，所以没有写相应的注释，有什么不会的直接留言~

pro read_dat

hdr_path=DIALOG_PICKFILE(/READ)
hdr_path=string(hdr_path)

head_info=strarr(6)
openr,lun,hdr_path,/get_lun
readf,lun,head_info
;print,head_info
free_lun,lun
 I=5;
 while (i gt 0) do begin
   temp=string(head_info[i])
   IF(temp.contains('samples'))THEN BEGIN
     xs=temp
   ENDIF
   IF(temp.contains('bands'))THEN BEGIN
     band=temp
   ENDIF
   IF(temp.contains('lines'))THEN BEGIN
     ys=temp
   ENDIF
   IF i GT 5 THEN BEGIN
     BREAK
   ENDIF
   i--
 endwhile
  ;print,head_info[j]

;xs=string(head_info[3])
xs=xs.split("=")
xs=xs[1]
;ys=string(head_info[4])
ys=ys.split("=")
ys=ys[1]
;band=string(head_info[4])
band=band.split("=")
band=band[1]

img_path=dialog_pickfile(/read)
img_path=string(img_path)
imageSize=[uint(xs),uint(ys)]
mbells = READ_BINARY(img_path, DATA_DIMS = imagesize)

; TrueColor display
DEVICE, DECOMPOSED = 0
LOADCT, 0
WINDOW, 0, XSIZE = imageSize[0], YSIZE = imageSize[1], $
   TITLE = 'A Binary Image, Scaled'
TVSCL, mbells
end

代码实现的效果

灰度图实现效果

送你们个惊喜

国外优秀网站

在GIS的沙场，传播有价值的东西!