MODIS Ocean Level 3 Data读取程序.doc

一、MODIS Ocean Level 3 Binned Data读取程序 IDL语言程序 从HDF文件中读取SDS Function Get_SDS, file_name, SDS_Name If (N_Elements( file_name ) eq 0) then Message, Must specify filename. Id = HDF_SD_Start( file_name, /READ ) Index = HDF_SD_NameToIndex( Id, SDS_Name ) SDS = HDF_SD_Select( Id, Index ) HDF_SD_GetData, SDS, Data HDF_SD_End, Id Return, Data End 从一个byte中提取1个或多个bit的资料 ; The IDL function, Bit.pro, extracts bits from an array of byte values. ; ; i = array of byte values ; j = starting bit position ; k = number of bits to be extracted ; ; Kevin Turpie, SAIC GSO Function Bit, i, j, k If (N_Elements( k ) eq 0) then k = 1 k = Fix( k ) If (N_Elements( j ) eq 0) then j = 1 j = Fix( j ) Return, IShft ( IShft ( Long (i), 32-j-k+1 ), -(31-k+1) ) End 将二进制资料转换成经纬度坐标 ; This code converts bin number to (lat,lon) ; It is written in both C and IDL. ; The IDL routine, binloc.pro, calls a C handler, binloc.c, which in turn ; calls the bin2ll routine from a shareable library, created from bin_csub.c ; (originally from SeaDAS). Compilation instructions in the binloc.c file ; work for SGI. binloc.pro references a .so library file created using those ; instructions. The user must be careful to correctly specify the path and ; name of this library file for their system in binloc.pro. ;--------------------------- IDL Interface: ------------------------------- ;-------------------------------------------------------------------------- ; ; BinLoc ; ; This IDL procedure takes an array of L3 bin numbers (b) and returns the ; corresponding latitude and longitude values. ; ; ARGUMENTS: ; b = corresponding array of bin numbers ; Lat = lat/lon grid containing sum of values ; Lon = lat/lon grid containing number of values ; ; KEYWORDS: ; NROWS = number of rows in input L3 bin grid. ; default is 4320 for MODIS. for SeaWiFS, set to 2160. ; ; AUTHOR: Kevin Turpie, 08 May 2000 ; SAIC General Sciences Corp