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