def main(): p = argparse.ArgumentParser() p.add_argument("-s", "--satellite_name") p.add_argument("-o", "--output_dir", default='.', nargs=1) p.add_argument('-b', '--channels', nargs='+', help='Input channels', type=str) p.add_argument('-a', '--area-name', nargs=1, help='Name of the area definition', type=str) p.add_argument("-i", "--input-files", nargs='+', help="Input Mitiff Files") args = p.parse_args() scene_list = [] for input_file in args.input_files: print "Loading %s" % (input_file) amsr2_scene = io.load_osisaf_amsr2_netcdf(input_file, bands=args.channels) scene_list.append(amsr2_scene) mosaic = MosaicScene() mosaic.get_area_def(args.area_name[0]) mosaic.compose_mosaic(scene_list, resample_method='gaussian') output_filename = "{}_12hr-mosaic_{}.tiff".format(args.satellite_name, mosaic.end_timestamp_string) output_filepath = os.path.join(args.output_dir[0], output_filename) mosaic.save_geotiff(output_filepath, bands=args.channels, cmap='istjenesten')
class TestMosaicScene(unittest.TestCase): def setUp(self): self.bands_list = ['lat_h'] self.scene_1 = load_osisaf_amsr2_netcdf('test_data/amsr2.nc', bands=self.bands_list) self.scene_2 = load_osisaf_amsr2_netcdf('test_data/amsr2.nc', bands=self.bands_list) self.mosaic = MosaicScene() self.mosaic.get_area_def(area_name='nsidc_stere_north_300k') def test_MosaicScene_ComposeMosaic(self): self.mosaic.compose_mosaic([self.scene_1, self.scene_2], resample_method='nn') self.mosaic.save_geotiff('out.tif', bands=self.bands_list) self.assertTrue(os.path.exists('out.tif'))
def main(): p = argparse.ArgumentParser() p.add_argument("-s", "--satellite_name") p.add_argument("-o", "--output_dir", default='.', nargs=1) p.add_argument('-c', '--channels', nargs='+', help='Input channels', type=int) p.add_argument('-a', '--area-name', nargs=1, help='Name of the area definition', type=str) p.add_argument("-i", "--input-files", nargs='+', help="Input Mitiff Files") args = p.parse_args() scene_list = [] for input_file in args.input_files: print "Loading %s" % (input_file) mitiff = io.load_mitiff(input_file, bands=args.channels) scene_list.append(mitiff) mosaic = MosaicScene() mosaic.get_area_def(args.area_name[0]) mosaic.compose_mosaic(scene_list, resample_method='nn') # Add '1' to channel id's for readability channels_string = "ch"+'-'.join(map(str, numpy.array(args.channels)+1)) output_filename = "{}_mosaic_{}-{}_{}.tiff".format(args.satellite_name, mosaic.start_timestamp_string, mosaic.end_timestamp_string, channels_string) output_filepath = os.path.join(args.output_dir[0], output_filename) mosaic.save_geotiff(output_filepath, bands=args.channels)
def setUp(self): self.bands_list = ['lat_h'] self.scene_1 = load_osisaf_amsr2_netcdf('test_data/amsr2.nc', bands=self.bands_list) self.scene_2 = load_osisaf_amsr2_netcdf('test_data/amsr2.nc', bands=self.bands_list) self.mosaic = MosaicScene() self.mosaic.get_area_def(area_name='nsidc_stere_north_300k')