main()

raster2pgsql.main

(

)

Definition at line 955 of file raster2pgsql.py.

def main():
 
    (opts, args) = parse_command_line()
 
    global VERBOSE
    VERBOSE = opts.verbose
 
    global SUMMARY
    SUMMARY = []
    
    saved_out = sys.stdout
    if isinstance(opts.output, str):
        filename = opts.output
        opts.output = open(filename, "w")
 
    # BEGIN
    opts.output.write('BEGIN;\n')
 
    # If overviews requested, CREATE TABLE raster_overviews
    if opts.create_raster_overviews_table:
        sql = make_sql_create_raster_overviews(opts)
        opts.output.write(sql)
 
    # Base raster schema
    if opts.overview_level == 1:
        # DROP TABLE
        if opts.drop_table:
            sql = make_sql_drop_raster_table(opts.table)
            opts.output.write(sql)
 
        # CREATE TABLE
        if opts.create_table and opts.overview_level == 1:
            sql = make_sql_create_table(opts)
            opts.output.write(sql)
              
    # INSERT
    i = 0
 
    # Burn all specified input raster files into single WKTRaster table
    gt = None
    for infile in opts.raster:
        filelist = glob.glob(infile)
        assert len(filelist) > 0, "No input raster files found for '" + str(infile) + "'"
 
        for filename in filelist:
            logit("MSG: Dataset #%d: %s\n" % (i + 1, filename))
            
            # Write raster data to WKB and send it to opts.output
            gt = wkblify_raster(opts, filename.replace( '\\', '/') , i, gt)
            i += 1
 
    # INDEX
    if opts.index and SUMMARY is not None:
        sql = make_sql_create_gist(SUMMARY[0][0], opts.column)
        opts.output.write(sql)
    
    # COMMIT
    opts.output.write('END;\n')
 
    # VACUUM
    if opts.vacuum and SUMMARY is not None:
        sql = make_sql_vacuum(SUMMARY[0][0])
        opts.output.write(sql)
 
    # Cleanup
    if opts.output != sys.stdout:
        sys.stdout = saved_out
 
        print("------------------------------------------------------------")
        print(" Summary of GDAL to PostGIS Raster processing:")
        print("------------------------------------------------------------")
        if i == 1:
            m = '%d (%s)' % (i, infile)
        else:
            m = '%d' % i
        print("Number of processed raster files: " + m)
        print("List of generated tables (number of tiles):")
        i = 0
        for s in SUMMARY:
            i += 1
            print("%d\t%s (%d)" % (i, s[0], s[1]))