MATLAB: How to use MULTIBAND read to read a BIL file of elevation data in Mapping Toolbox 3.6 (R2012b)

banddatafileinterleavedMapping Toolbox

I have downloaded data from USGS and have the a .BIL file and the following header file (.HDR). From this information, how do I determine the arguments to be passed to MULTIBANDREAD, to read and display the image ?
 BYTEORDER      I
 LAYOUT         BIL
 NROWS          1383
 NCOLS          1621
 NBANDS         1
 NBITS          16
 BANDROWBYTES   3242
 TOTALROWBYTES  3242
 PIXELTYPE      SIGNEDINT
 ULXMAP         -121.386990740857
 ULYMAP         43.7869907403301
 XDIM           9.25925925999961e-005
 YDIM           9.25925925999957e-005
 NODATA         32767

Best Answer

The raster file in .BIL format can be read using the MULTIBAND function in the Mapping Toolbox. The arguments to the function are all obtained from the corresponding fields of the header file.
filename ='MyFiles\Bil\45765067\45765067.bil';
% From the .hdr file, I obtained the following arguments for the MUTLIBANDREAD function.
size =      [ 1383 1621 1];      % [ NROWS NCOLS NBANDS]
precision = 'int16';                  % from NBITS and PIXEL TYPE = int
offset = 0;                              % since the header is not included in this file
interleave = 'bil';                     % LAYOUT
byteorder = 'ieee-le';              % BYTEORDER = I (refers to Intel order or little endian format)
X_bil = multibandread(filename, size, precision, offset, interleave, byteorder);
figure, imagesc(X_bil)               % Display the image file using IMAGESC
demcmap(X_bil)                        % Changing the colormap
Related Question