def main()

in utils/reproject_data.py [0:0]

• 41 lines of code
• 5 McCabe index (conditional complexity)

def main():
    parser = argparse.ArgumentParser(description="Reprojection script")

    parser.add_argument("-v", "--verbose", action="store_true", help="Enable verbose debugging", default=False)
    parser.add_argument("--input_fn", action="store", dest="input_fn", type=str, help="Path to input raster (this is the larger file that we want to chunk up)", required=True)
    parser.add_argument("--target_fn", action="store", dest="target_fn", type=str, help="Path to target raster (this is the shape that we want to copy)", required=True)
    parser.add_argument("--output_fn", action="store", dest="output_fn", type=str, help="Path to output", required=True)
   
    args = parser.parse_args(sys.argv[1:])

    start_time = float(time.time())

    if args.verbose:
        print("Starting reprojection")
        print("--input_fn: %s" % (args.input_fn))
        print("--target_fn: %s" % (args.target_fn))
        print("--output_fn: %s" % (args.output_fn))
        print("")

        print("Reading metadata from target_fn")
    f = rasterio.open(args.target_fn,"r")
    left, bottom, right, top = f.bounds
    crs = f.crs.to_string()
    height, width = f.height, f.width
    f.close()

    assert crs.startswith("EPSG")

    command = [
        "gdalwarp",
        "" if args.verbose else "-q",
        "-overwrite",
        "-ot", "Byte",
        "-t_srs", crs,
        "-r", "near",
        "-of", "GTiff",
        "-te", str(left), str(bottom), str(right), str(top),
        "-ts", str(width), str(height),
        "-co", "COMPRESS=LZW",
        "-co", "BIGTIFF=YES",
        args.input_fn,
        args.output_fn
    ]
    if args.verbose:
        print("Reprojecting input_fn to shape of target_fn")
    subprocess.call(command)

    if args.verbose:
        print("Finished in %0.4f seconds" % (time.time() - start_time))