# This file is part of meas_base.
#
# Developed for the LSST Data Management System.
# This product includes software developed by the LSST Project
# (https://www.lsst.org).
# See the COPYRIGHT file at the top-level directory of this distribution
# for details of code ownership.
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see <https://www.gnu.org/licenses/>.

import unittest

import lsst.geom
import lsst.daf.base as dafBase
from lsst.afw.geom import makeSkyWcs
import lsst.utils.tests

from lsst.meas.base.tests import AlgorithmTestCase
from lsst.meas.base.sfm import SingleFrameMeasurementConfig


class JacobianTestCase(AlgorithmTestCase, lsst.utils.tests.TestCase):

    def setUp(self):
        # Pick arbitrary numbers to create a detector object, and a synthetic
        # dataset. The particular numbers have no special meaning to the test
        # and be anything as long as they are self consistent (i.e. the
        # fake source is inside the bounding box)
        self.center = lsst.geom.Point2D(50.1, 49.8)
        self.bbox = lsst.geom.Box2I(lsst.geom.Point2I(-20, -30),
                                    lsst.geom.Extent2I(140, 160))
        self.dataset = lsst.meas.base.tests.TestDataset(self.bbox)
        self.dataset.addSource(100000.0, self.center)

        md = dafBase.PropertyList()

        for k, v in (
            ("EQUINOX", 2000.0),
            ("CRPIX1", 5353.0),
            ("CRPIX2", -35.0),
            ("CD1_1", 0.0),
            ("CD1_2", -5.611E-05),
            ("CD2_1", -5.611E-05),
            ("CD2_2", -0.0),
            ("CRVAL1", 4.5789875),
            ("CRVAL2", 16.30004444),
            ("CUNIT1", 'deg'),
            ("CUNIT2", 'deg'),
            ("CTYPE1", 'RA---TAN'),
            ("CTYPE2", 'DEC--TAN'),
        ):
            md.set(k, v)

        self.wcs = makeSkyWcs(md)

    def tearDown(self):
        del self.center
        del self.bbox
        del self.dataset
        del self.wcs

    def testJacobianPlugin(self):
        config = SingleFrameMeasurementConfig()
        config.plugins.names |= ['base_Jacobian']
        # Pixel scale chosen to approximately match the scale defined in the
        # above WCS
        config.plugins['base_Jacobian'].pixelScale = 0.2
        task = self.makeSingleFrameMeasurementTask(config=config)
        exposure, catalog = self.dataset.realize(10.0, task.schema, randomSeed=0)
        exposure.setWcs(self.wcs)
        task.run(catalog, exposure)
        record = catalog[0]
        self.assertFalse(record.get("base_Jacobian_flag"))
        self.assertAlmostEqual(record.get("base_Jacobian_value"), 1.0200183929088285, 4)


class TestMemory(lsst.utils.tests.MemoryTestCase):
    pass


def setup_module(module):
    lsst.utils.tests.init()


if __name__ == "__main__":
    lsst.utils.tests.init()
    unittest.main()