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')
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')