DQI",
                       region, &toa_brf_data, &toa_brf_map_info);
status = MtkReadData("MISR_AM1_AS_LAND_P037_O029058_F06_0017.hdf", "SubregParamsLnd", "LandBRF", region, &land_brf_data, &land_brf_map_info);
status = MtkDataBufferAllocate(land_brf_data.nline, land_brf_data.nsample, MTKe_float,&land_brf_sigma);
for (iline = 0; iline < land_brf_data.nline; iline++) {
  for (isample = 0; isample < land_brf_data.nsample; isample++) {
    land_brf_sigma.data.f[iline][isample] = 1.0;
  }
}                  