IDL读取数据 方法及程序实例 编程示例.docxVIP

自己看到的方法，不知道是谁写的，就没办法注明出处了。作者见谅！ 1========== 1 pro readhdf 2 ; Set some constants 3 FILE_NAME=E:/sssbak/MODSST.hdf 4 SDS_NAME=sst 5 X_LENGTH=1354 6 Y_LENGTH=4856 7 8 ; Open the file and initialize the SD interface 9 sd_id = HDF_SD_START( FILE_NAME, /read ) 10 ; Find the index of the sds to read using its name 11 sds_index = HDF_SD_NAMETOINDEX(sd_id,SDS_NAME) 12 ; Select it 13 sds_id = HDF_SD_SELECT( sd_id, sds_index ) 14 15 ; Set the data subset limits. Actually, it is read from the first element, so start is [0,0]. X_LENGTH elements will be read along the X axis and Y_LENGTH along Y. 16 start=INTARR(2) ; the start position of the data to be read 17 start[0] = 0 18 start[1] = 0 19 edges=INTARR(2) ; the number of elements to read in each direction 20 edges[0] = X_LENGTH 21 edges[1] = Y_LENGTH 22 ; Read the data : you can notice that here, it is not needed to allocate the data array yourself 23 HDF_SD_GETDATA, sds_id, data 24 dims = size(data) 25 ;HDF_SD_GETDATA, sds_id, data, start = start, count = edges 26 27 ; print them on the screen. You can notice that the IDL method HDF_SD_GETDATA swaps the HDF indexes convention [Z,Y,X] to [X,Y,Z]. This method is more efficient on IDL. If you prefer the usual HDF convention, you should better use the set the NOREVERSE keyword when calling HDF_SD_GETDATA 28 ;FOR i=0,(X_LENGTH-1),1 DO BEGIN ; crosses X axis 29 ; FOR j=0,(Y_LENGTH-1),1 DO BEGIN ; crosses Y axis 30 ; PRINT, FORMAT=(I, ,$), data[i,j] 31 ; ENDFOR 32 ; PRINT, 33 ;ENDFOR 34 35 ; end access to SDS 36 HDF_SD_ENDACCESS, sds_id 37 ; close the hdf file 38 HDF_SD_END, sd_id 39 print,dims[1];1354 也就是说，实际上X_LENGTH，X_LENGTH可以通过size读出。 40 print,dims[2];4856 41 print,dims[3] 42 print,dims[4] 43 filter_data = uintarr(dims[1], dims[2],dims[3]) 44 ;imagedata = BYTSCL(data) 45 imagedata = bytscl( congrid(data,800,600) );缩小图片。 46 47 ;print,imagedata 48 loadct, 15;调色板编号不一样。 49 window, 0, xsize=800, ysize=600, retain=2 50 tv, imagedata 2.============= Readin