dim=[2707  4059 ; %  500 m resolution [pixels_across by scanlines] bands 1-2
     1354  2030 ; % 1000 m resolution                              bands 3-7
      677  1015]; % 2000 m resolution                              bands 17-36
      
fid = fopen('./MODIS_band31.be.bin','r','ieee-be');
  band31 = fread(fid,[677 1015],'float32');
fclose(fid);

fid = fopen('./lat_2000m.be.bin','r','ieee-be');
  lat2000m = fread(fid,[677 1015],'float32');
fclose(fid);      

fid = fopen('./long_2000m.le.bin','r','ieee-le');
  long2000m = fread(fid,[677 1015],'float32');
fclose(fid);  



