def main():
from argparse import ArgumentParser
description = """
Plot binary files using matplotlib.
"""
parser = ArgumentParser(description=description)
parser.add_argument("scene_file", help="the data file to be sampled from")
parser.add_argument("product_name", help="product in scene to inspect")
parser.add_argument("lonVal", type=float, help="Specify the longitude to find a the point nearest")
parser.add_argument("latVal", type=float, help="Specify the latitude to find a the point nearest")
args = parser.parse_args()
# check what the user gave us on the command line
if len(args.data_file) != 1 :
print("Data file required for point selection. Skipping.")
elif args.latVal is None or args.lonVal is None:
print("Longitude and latitude to find must be selected. Skipping.")
else:
# if we got to here, we have the minimum we needed to try
scene = SwathScene.load(args.scene_file)
data = scene[args.product_name].get_data_array()
swath_def = scene[args.product_name]["swath_definition"]
lon_data = swath_def.get_longitude_array()
lat_data = swath_def.get_longitude_array()
value, lonPt, latPt = _find_nearest_data_pt(args.lonVal, args.latVal, lon_data, lat_data, data)
print("Found data value %f at longitude %f and latitude %f" % (value, lonPt, latPt))
def main():
from polar2grid.core.script_utils import create_basic_parser, create_exc_handler, setup_logging
from polar2grid.core.containers import SwathScene
parser = create_basic_parser(description="Remap a SwathScene to the provided grids")
subgroup_titles = add_remap_argument_groups(parser)
parser.add_argument("--scene", required=True,
help="JSON SwathScene filename to be remapped")
parser.add_argument('-o', dest="output_filename", default="gridded_scene_{grid_name}.json",
help="Output filename for JSON scene (default is to 'gridded_scene_{grid_name}.json')")
global_keywords = ("keep_intermediate", "overwrite_existing", "exit_on_error")
args = parser.parse_args(subgroup_titles=subgroup_titles, global_keywords=global_keywords)
levels = [logging.ERROR, logging.WARN, logging.INFO, logging.DEBUG]
setup_logging(console_level=levels[min(3, args.verbosity)], log_filename=args.log_fn)
sys.excepthook = create_exc_handler(LOG.name)
LOG.debug("Starting script with arguments: %s", " ".join(sys.argv))
if args.output_filename and args.output_filename != "-" and os.path.isfile(args.output_filename):
LOG.error("JSON file '%s' already exists, will not overwrite." % (args.output_filename,))
raise RuntimeError("JSON file '%s' already exists, will not overwrite." % (args.output_filename,))
scene = SwathScene.load(args.scene)
remapper = Remapper(**args.subgroup_args["Remapping Initialization"])
remap_kwargs = args.subgroup_args["Remapping"]
for grid_name in remap_kwargs.pop("forced_grids", ["wgs84_fit"]):
gridded_scene = remapper.remap_scene(scene, grid_name, **remap_kwargs)
if args.output_filename is None or args.output_filename == "-":
print(gridded_scene.dumps(persist=True))
else:
fn = args.output_filename.format(grid_name=grid_name)
LOG.info("Saving gridded scene to file: %s", fn)
gridded_scene.save(fn)
def main():
from polar2grid.core.script_utils import create_basic_parser, create_exc_handler, setup_logging
from polar2grid.core.containers import SwathScene
parser = create_basic_parser(description="Remap a SwathScene to the provided grids")
subgroup_titles = add_remap_argument_groups(parser)
parser.add_argument("--scene", required=True,
help="JSON SwathScene filename to be remapped")
parser.add_argument('-o', dest="output_filename", default="gridded_scene_{grid_name}.json",
help="Output filename for JSON scene (default is to 'gridded_scene_{grid_name}.json')")
global_keywords = ("keep_intermediate", "overwrite_existing", "exit_on_error")
args = parser.parse_args(subgroup_titles=subgroup_titles, global_keywords=global_keywords)
levels = [logging.ERROR, logging.WARN, logging.INFO, logging.DEBUG]
setup_logging(console_level=levels[min(3, args.verbosity)], log_filename=args.log_fn)
sys.excepthook = create_exc_handler(LOG.name)
LOG.debug("Starting script with arguments: %s", " ".join(sys.argv))
if args.output_filename and args.output_filename != "-" and os.path.isfile(args.output_filename):
LOG.error("JSON file '%s' already exists, will not overwrite." % (args.output_filename,))
raise RuntimeError("JSON file '%s' already exists, will not overwrite." % (args.output_filename,))
scene = SwathScene.load(args.scene)
remapper = Remapper(**args.subgroup_args["Remapping Initialization"])
remap_kwargs = args.subgroup_args["Remapping"]
for grid_name in remap_kwargs.pop("forced_grids", ["wgs84_fit"]):
gridded_scene = remapper.remap_scene(scene, grid_name, **remap_kwargs)
if args.output_filename is None or args.output_filename == "-":
print(gridded_scene.dumps(persist=True))
else:
fn = args.output_filename.format(grid_name=grid_name)
LOG.info("Saving gridded scene to file: %s", fn)
gridded_scene.save(fn)
def main():
from argparse import ArgumentParser
description = """
Plot binary files using matplotlib.
"""
parser = ArgumentParser(description=description)
parser.add_argument("scene_file", help="the data file to be sampled from")
parser.add_argument("product_name", help="product in scene to inspect")
parser.add_argument(
"lonVal",
type=float,
help="Specify the longitude to find a the point nearest")
parser.add_argument(
"latVal",
type=float,
help="Specify the latitude to find a the point nearest")
args = parser.parse_args()
# check what the user gave us on the command line
if len(args.data_file) != 1:
print("Data file required for point selection. Skipping.")
elif args.latVal is None or args.lonVal is None:
print("Longitude and latitude to find must be selected. Skipping.")
else:
# if we got to here, we have the minimum we needed to try
scene = SwathScene.load(args.scene_file)
data = scene[args.product_name].get_data_array()
swath_def = scene[args.product_name]["swath_definition"]
lon_data = swath_def.get_longitude_array()
lat_data = swath_def.get_longitude_array()
value, lonPt, latPt = _find_nearest_data_pt(args.lonVal, args.latVal,
lon_data, lat_data, data)
print("Found data value %f at longitude %f and latitude %f" %
(value, lonPt, latPt))