#!/usr/bin/env python

# -*- coding: UTF-8 -*-

import tempfile

import shutil

import argparse

import sys

import os

from repipy import utilities

from lemon import methods

import repipy

import astropy.io.fits as fits

# Change to the directory where repipy is installed to load pyraf

with methods.tmp_chdir(repipy.__path__[0]):

import pyraf.iraf as iraf

from iraf import imcopy

def trim(args):

return args.output

# Create parser

parser = argparse.ArgumentParser(description='Trim a given section of an image, keeping WCS information true')

parser.add_argument("input", metavar='input', action='store', nargs="+", \

type=str, help='Images to be trimmed.')

parser.add_argument("--output", metavar='output', action='store', dest='output', \

default='', help='Name of output file.')

parser.add_argument("--suffix", metavar="suffix", dest='suffix', action='store', \

default='', help='suffix to be added at the end of the image' + \

'input list to generate the output.')

parser.add_argument("--region", metavar=('x0', 'x1', 'y0', 'y1'), action='store', nargs=4, type=int,

required=True, dest="region", \

help='Region of the image to be trimmed. x0,x1 is the range in the horizontal '

'axis as shown by ds9, while y0, y1 mean the same in the vertical axis.'

'Please, note that what ds9 shows as X (i.e. the horizontal axis) is the second index in '

'a numpy array, the Y axis as shown by ds9 being the first. '

'Usually x0,x1 correspond to Naxis1 and y0,y1 to Naxis2 in the header of fits images. ')

parser.add_argument("--overwrite", dest='overwrite', action='store_true',

help='Overwrite the input file if no suffix or output file is given. ')

parser.add_argument("--mask_key", metavar='mask_key', action='store',

dest='mask_key', default="",

help='Key where the name of the mask image is stored in the header. The mask image should be trimmed '

'to be the same size as the output image.')

def main(arguments=None):

return newfile

if __name__ == "__main__":

main()