Esempio n. 1
0
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')
Esempio n. 2
0
 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')