diff options
author | Aaron LI <aaronly.me@outlook.com> | 2017-02-18 22:47:43 +0800 |
---|---|---|
committer | Aaron LI <aaronly.me@outlook.com> | 2017-02-18 22:47:43 +0800 |
commit | a411c50b460df7bd22e264dc149ee413a3ae6a6a (patch) | |
tree | 2bcbcb181b5485aa7cb3bf9c8233d0b16836e999 /bin | |
parent | 1fb2a0eefae77e242123100e7ac43aa157064bdc (diff) | |
download | atoolbox-a411c50b460df7bd22e264dc149ee413a3ae6a6a.tar.bz2 |
Remove bin/gen_points.py
Diffstat (limited to 'bin')
-rwxr-xr-x | bin/gen_points.py | 81 |
1 files changed, 0 insertions, 81 deletions
diff --git a/bin/gen_points.py b/bin/gen_points.py deleted file mode 100755 index 57733dc..0000000 --- a/bin/gen_points.py +++ /dev/null @@ -1,81 +0,0 @@ -#!/usr/bin/env python3 -# -*- coding: utf-8 -*- -# -# Aaron LI -# 2015/06/19 - -""" -Generate the required number of random points within the required region. -""" - -__version__ = "0.1.0" -__date__ = "2015/06/19" -DEBUG = True - -import sys -import argparse -import random -import time -import re - -from rand.sphere import sphere_point -from region.region import Region - -random.seed(time.time()) - - -def parse_region(regstring): - reg_par = re.sub(r"[(),]", " ", regstring).split() - regtype = reg_par[0].lower() - if regtype == "box": - xc = float(reg_par[1]) - yc = float(reg_par[2]) - width = parse_reg_value(reg_par[3]) - height = parse_reg_value(reg_par[4]) - rotation = float(reg_par[5]) - reg = Region(regtype, xc=xc, yc=yc, - width=width, height=height, rotation=rotation) - else: - raise ValueError("region type '%s' currently not implemented" % regtype) - return reg - - -def parse_reg_value(valstring): - if valstring[-1] == '"': - # arcsec -> deg - value = float(valstring.split('"')[0]) / 60.0 / 60.0 - elif valstring[-1] == "'": - # arcmin -> deg - value = float(valstring.split("'")[0]) / 60.0 - else: - value = float(valstring) - return value - - -def main(): - parser = argparse.ArgumentParser( - description="Generate random point within the given region.") - parser.add_argument("-V", "--version", action="version", - version="%(prog)s " + "%s (%s)" % (__version__, __date__)) - parser.add_argument("-n", "--number", dest="number", - type=int, default=1, - help="number of points to be generated") - parser.add_argument("-r", "--region", dest="region", required=True, - help="DS9 region") - args = parser.parse_args() - - reg = parse_region(args.region) - if DEBUG: - print("DEBUG: region: ", reg.dump(), file=sys.stderr) - - points = [] - while len(points) < args.number: - p = sphere_point(unit="deg") - if reg.is_inside(p): - points.append(p) - print("%s %s" % p) - - -if __name__ == "__main__": - main() - |