示例#1
0
 def load_scene(self, fmt=None):
     if fmt is 'mitiff':
         self.scene = io.load_mitiff(self.file_path, self.config_path)
     elif fmt is 'netcdf':
         self.scene = io.load_netcdf(self.file_path, self.config_path)
     else:
         raise Exception('{0} reader not implemented'.format(fmt))
     self.bands = self.scene.bands
示例#2
0
def main():
    p = argparse.ArgumentParser()
    p.add_argument("-o", "--output-directory")
    p.add_argument("input_files", nargs="+", help="Input Mitiff Files")
    p.add_argument("-c", "--channels", type=int, nargs="+", help="Channels which should be saved into geotiff")
    args = p.parse_args()

    print "Channels: ", args.channels

    files_number = len(args.input_files)
    for i, input_file in enumerate(args.input_files):
        print "Converting file {} of {}".format(i + 1, files_number)
        output_filepath = make_output_filepath(input_file, args.output_directory)
        mitiff = io.load_mitiff(input_file)
        # mitiff.save_rgb_image('polar_lows.tif', args.channels)
        mitiff.save_geotiff(output_filepath, bands=args.channels)
示例#3
0
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)
示例#4
0
def step_impl(context, input_filepath):
    context.input_filepath = input_filepath
    context.scene = io.load_mitiff(input_filepath)
    assert context.scene.area_def is not None