Ejemplo n.º 1
0
        print flt_directory_out

        test_flag = options.test_flag

        #        shutil.rmtree(flt_directory_out)
        #        os.mkdir(flt_directory_out)

        logger.log_status("Convert input data from ", str(input_year))

    lc = LandCoverDataset(in_storage=StorageFactory().get_storage(
        'flt_storage', storage_location=flt_directory_in),
                          out_storage=StorageFactory().get_storage(
                              'flt_storage',
                              storage_location=flt_directory_out))

    lc.get_header()  # added 23 june 2009 by mm
    mask = lc.get_mask()
    idx = where(mask == 0)[0]
    lcsubset = DatasetSubset(lc, idx)
    print "Converting:"
    lcsubset.write_dataset(attributes=["relative_x"],
                           out_table_name="land_covers")
    #lcsubset.write_dataset(attributes=["relative_x"], out_table_name="land_covers",
    #                            valuetypes=valuetypes)
    lc.delete_one_attribute("relative_x")
    lcsubset.write_dataset(attributes=["relative_y"],
                           out_table_name="land_covers")
    #lcsubset.write_dataset(attributes=["relative_y"], out_table_name="land_covers",
    #                            valuetypes=valuetypes)
    lc.delete_one_attribute("relative_y")
    #    srcdir = os.path.join(flt_directory_out, "land_covers", "computed")
Ejemplo n.º 2
0
        input_year = sys.argv[1]   
    
        flt_directory_in = options.input
        flt_directory_out = options.output
        
        print flt_directory_out
        
        test_flag = options.test_flag
        
        shutil.rmtree(flt_directory_out)
        os.mkdir(flt_directory_out)
        
        logger.log_status("Convert input data from ", str(input_year))
    
    lc =  LandCoverDataset(in_storage = StorageFactory().get_storage('flt_storage', storage_location = flt_directory_in), 
        out_storage = StorageFactory().get_storage('flt_storage', storage_location = flt_directory_out), debuglevel=4)
    
    lc.get_header()
    
#    mask = lc.get_mask()
#    idx = where(mask==0)[0]
#    lcsubset = DatasetSubset(lc, idx)
    print "Creating and writing relative_x and relative_y:"
    lc.write_dataset(attributes=["relative_x"], out_table_name="land_covers",
                                valuetypes=valuetypes)
    lc.delete_one_attribute("relative_x")
    lc.write_dataset(attributes=["relative_y"], out_table_name="land_covers",
                                valuetypes=valuetypes)
    lc.delete_one_attribute("relative_y")
    
    print "done."
Ejemplo n.º 3
0
        flt_directory_in = options.input
        flt_directory_out = options.output
        
        print flt_directory_out
        
        test_flag = options.test_flag
        
#        shutil.rmtree(flt_directory_out)
#        os.mkdir(flt_directory_out)
        
        logger.log_status("Convert input data from ", str(input_year))
    
    lc =  LandCoverDataset(in_storage = StorageFactory().get_storage('flt_storage', storage_location = flt_directory_in), 
        out_storage = StorageFactory().get_storage('flt_storage', storage_location = flt_directory_out))
    
    lc.get_header() # added 23 june 2009 by mm
    mask = lc.get_mask()
    idx = where(mask==0)[0]
    lcsubset = DatasetSubset(lc, idx)
    print "Converting:"
    lcsubset.write_dataset(attributes=["relative_x"], out_table_name="land_covers")
    #lcsubset.write_dataset(attributes=["relative_x"], out_table_name="land_covers",
    #                            valuetypes=valuetypes)
    lc.delete_one_attribute("relative_x")
    lcsubset.write_dataset(attributes=["relative_y"], out_table_name="land_covers")
    #lcsubset.write_dataset(attributes=["relative_y"], out_table_name="land_covers",
    #                            valuetypes=valuetypes)
    lc.delete_one_attribute("relative_y")
#    srcdir = os.path.join(flt_directory_out, "land_covers", "computed")
#    shutil.move(os.path.join(srcdir,"relative_x.li4"), os.path.join(flt_directory_out, "land_covers"))
#    shutil.move(os.path.join(srcdir,"relative_y.li4"), os.path.join(flt_directory_out, "land_covers"))