from pipeline import center_approx
import integration
import numpy as np
import peakfinding
from xicam import config

from pyFAI import calibrant
from functools import wraps
import msg
from xicam import threads
from scipy import signal


def calibrationAlgorithm(f):
    @wraps(f)
    def wrapped(*args, **kwargs):
        msg.showMessage('Calibrating...')
        msg.showBusy()

        runnable = threads.RunnableIterator(f, iterator_args=args, iterator_kwargs=kwargs,
                                            callback_slot=showProgress, finished_slot=calibrationFinish)
        threads.add_to_queue(runnable)

    return wrapped


def calibrationFinish():
    msg.showMessage('Calibration complete!', 4)
    msg.hideBusy()


def showProgress(value):
    msg.showProgress(value)


@calibrationAlgorithm
def fourierAutocorrelation(dimg, calibrantkey):
    yield 0
    if dimg.transformdata is None:
        return
    yield 20

    config.activeExperiment.center = center_approx.center_approx(dimg.transformdata,dimg.mask)

    yield 40

    radialprofile = integration.pixel_2Dintegrate(dimg, mask=dimg.mask)

    yield 60

    peaks = np.array(peakfinding.findpeaks(np.arange(len(radialprofile)), radialprofile)).T

    yield 80

    peaks = peaks[peaks[:, 1].argsort()[::-1]]

    for peak in peaks:
        if peak[0] > 15 and not np.isinf(peak[1]):  ####This thresholds the minimum sdd which is acceptable
            bestpeak = peak[0]
            # print peak
            break

    yield 85
    N = 1
    # guess which order is the first
    stds=[np.std((peaks[:,0]/(np.arange(len(peaks))+i))[:4]) for i in range(1,5)]
    if min(stds)<5: # threshold accepting peak as higher N
        N=np.argmin(stds)+1 # index of the first detected peak

    msg.logMessage(('Using peak #',N,' for calibration...'))
    yield 90

    calibrant1stpeak = calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing[N-1]

    # Calculate sample to detector distance for lowest q peak
    tth = 2 * np.arcsin(0.5 * config.activeExperiment.getvalue('Wavelength') / calibrant1stpeak / 1.e-10)
    tantth = np.tan(tth)
    sdd = bestpeak * config.activeExperiment.getvalue('Pixel Size X') / tantth

    config.activeExperiment.setvalue('Detector Distance', sdd)

    #        self.refinecenter()

    dimg.invalidatecache()

    yield 100


@calibrationAlgorithm
def rickerWavelets(dimg, calibrantkey):
    yield 0

    if dimg.transformdata is None:
        return

    yield 1

    radii = np.arange(30, 100)
    img = dimg.transformdata

    maxval = 0
    center = np.array([0, 0], dtype=np.int)
    for i in range(len(radii)):
        yield int(96. * i / len(radii)) + 1
        w = center_approx.tophat2(radii[i], scale=1000)
        im2 = signal.fftconvolve(img, w, 'same')
        if im2.max() > maxval:
            maxval = im2.max()
            center = np.array(np.unravel_index(im2.argmax(), img.shape))

    config.activeExperiment.center = center

    yield 97

    radialprofile = integration.pixel_2Dintegrate(dimg, mask=dimg.mask)

    yield 98

    peaks = np.array(peakfinding.findpeaks(np.arange(len(radialprofile)), radialprofile)).T

    yield 99

    peaks = peaks[peaks[:, 1].argsort()[::-1]]

    for peak in peaks:
        if peak[0] > 25 and not np.isinf(peak[1]):  ####This thresholds the minimum sdd which is acceptable
            bestpeak = peak[0]
            # print peak
            break

    calibrant1stpeak = calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing[0]

    # Calculate sample to detector distance for lowest q peak
    tth = 2 * np.arcsin(0.5 * config.activeExperiment.getvalue('Wavelength') / calibrant1stpeak / 1.e-10)
    tantth = np.tan(tth)
    sdd = bestpeak * config.activeExperiment.getvalue('Pixel Size X') / tantth

    config.activeExperiment.setvalue('Detector Distance', sdd)

    #        self.refinecenter()

    dimg.invalidatecache()

    yield 100


import saxs_calibration as sc
@calibrationAlgorithm
def dpdakRefine(dimg,calibrantkey):

    if dimg.transformdata is None:
        return

    # Refine calibration
    # d-spacing for Silver Behenate
    d_spacings = np.array(sorted(calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing,key=float,reverse=True))

    geometry = config.activeExperiment.getAI()

    data = dimg.rawdata

    print 'Start parameter:'
    print geometry.getFit2D()

    fit_param = ['distance', 'rotation', 'tilt', 'center_x', 'center_y']
    # calculate maxima for every d_spacing
    center = (geometry.getFit2D()['centerX'],
              geometry.getFit2D()['centerY'])
    radial_pos = sc.radial_array(center, data.shape)
    x_data, y_data, = [], []

    # calculate maxima for every d_spacing
    for i in range(len(d_spacings)):
        yield 100/len(d_spacings)*i
        maxima_x, maxima_y, radial_pos = sc.ring_maxima(geometry,
                                                     d_spacings[i],
                                                     data,
                                                     radial_pos,
                                                     10) # TODO: parameterize this
        x_data.extend(maxima_x)
        y_data.extend(maxima_y)



    # start fit
    sc.fit_geometry(geometry,
                        (np.array(x_data), np.array(y_data)),
                        d_spacings,
                        fit_param)


    print 'Final parameter:'
    print geometry.getFit2D()

    config.activeExperiment.setvalue('Detector Distance', geometry.get_dist())
    config.activeExperiment.center = geometry.getFit2D()['centerX'],geometry.getFit2D()['centerY']
    yield 100

from pipeline import center_approx
import integration
import numpy as np
import peakfinding
from xicam import config

from pyFAI import calibrant
from functools import wraps
import msg
from xicam import threads
from scipy import signal


def calibrationAlgorithm(f):
    @wraps(f)
    def wrapped(*args, **kwargs):
        msg.showMessage('Calibrating...')
        msg.showBusy()

        runnable = threads.RunnableIterator(f, iterator_args=args, iterator_kwargs=kwargs,
                                            callback_slot=showProgress, finished_slot=calibrationFinish)
        threads.add_to_queue(runnable)

    return wrapped


def calibrationFinish():
    msg.showMessage('Calibration complete!', 4)
    msg.hideBusy()


def showProgress(value):
    msg.showProgress(value)


@calibrationAlgorithm
def fourierAutocorrelation(dimg, calibrantkey):
    yield 0
    if dimg.transformdata is None:
        return
    yield 20

    config.activeExperiment.center = center_approx.center_approx(dimg.transformdata,dimg.mask)

    yield 40

    radialprofile = integration.pixel_2Dintegrate(dimg, mask=dimg.mask)

    yield 60

    peaks = np.array(peakfinding.findpeaks(np.arange(len(radialprofile)), radialprofile)).T

    yield 80

    peaks = peaks[peaks[:, 1].argsort()[::-1]]

    for peak in peaks:
        if peak[0] > 15 and not np.isinf(peak[1]):  ####This thresholds the minimum sdd which is acceptable
            bestpeak = peak[0]
            # print peak
            break

    yield 85
    N = 1
    # guess which order is the first
    stds=[np.std((peaks[:,0]/(np.arange(len(peaks))+i))[:4]) for i in range(1,5)]
    if min(stds)<5: # threshold accepting peak as higher N
        N=np.argmin(stds)+1 # index of the first detected peak

    msg.logMessage(('Using peak #',N,' for calibration...'))
    yield 90

    calibrant1stpeak = calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing[N-1]

    # Calculate sample to detector distance for lowest q peak
    tth = 2 * np.arcsin(0.5 * config.activeExperiment.getvalue('Wavelength') / calibrant1stpeak / 1.e-10)
    tantth = np.tan(tth)
    sdd = bestpeak * config.activeExperiment.getvalue('Pixel Size X') / tantth

    config.activeExperiment.setvalue('Detector Distance', sdd)

    #        self.refinecenter()

    dimg.invalidatecache()

    yield 100


@calibrationAlgorithm
def rickerWavelets(dimg, calibrantkey):
    yield 0

    if dimg.transformdata is None:
        return

    yield 1

    radii = np.arange(30, 100)
    img = dimg.transformdata

    maxval = 0
    center = np.array([0, 0], dtype=np.int)
    for i in range(len(radii)):
        yield int(96. * i / len(radii)) + 1
        w = center_approx.tophat2(radii[i], scale=1000)
        im2 = signal.fftconvolve(img, w, 'same')
        if im2.max() > maxval:
            maxval = im2.max()
            center = np.array(np.unravel_index(im2.argmax(), img.shape))

    config.activeExperiment.center = center

    yield 97

    radialprofile = integration.pixel_2Dintegrate(dimg, mask=dimg.mask)

    yield 98

    peaks = np.array(peakfinding.findpeaks(np.arange(len(radialprofile)), radialprofile)).T

    yield 99

    peaks = peaks[peaks[:, 1].argsort()[::-1]]

    for peak in peaks:
        if peak[0] > 25 and not np.isinf(peak[1]):  ####This thresholds the minimum sdd which is acceptable
            bestpeak = peak[0]
            # print peak
            break

    calibrant1stpeak = calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing[0]

    # Calculate sample to detector distance for lowest q peak
    tth = 2 * np.arcsin(0.5 * config.activeExperiment.getvalue('Wavelength') / calibrant1stpeak / 1.e-10)
    tantth = np.tan(tth)
    sdd = bestpeak * config.activeExperiment.getvalue('Pixel Size X') / tantth

    config.activeExperiment.setvalue('Detector Distance', sdd)

    #        self.refinecenter()

    dimg.invalidatecache()

    yield 100


import saxs_calibration as sc
@calibrationAlgorithm
def dpdakRefine(dimg,calibrantkey):

    if dimg.transformdata is None:
        return

    # Refine calibration
    # d-spacing for Silver Behenate
    d_spacings = np.array(sorted(calibrant.ALL_CALIBRANTS[calibrantkey].dSpacing,key=float,reverse=True))

    geometry = config.activeExperiment.getAI()

    data = dimg.rawdata

    print 'Start parameter:'
    print geometry.getFit2D()

    fit_param = ['distance', 'rotation', 'tilt', 'center_x', 'center_y']
    # calculate maxima for every d_spacing
    center = (geometry.getFit2D()['centerX'],
              geometry.getFit2D()['centerY'])
    radial_pos = sc.radial_array(center, data.shape)
    x_data, y_data, = [], []

    # calculate maxima for every d_spacing
    for i in range(len(d_spacings)):
        yield 100/len(d_spacings)*i
        maxima_x, maxima_y, radial_pos = sc.ring_maxima(geometry,
                                                     d_spacings[i],
                                                     data,
                                                     radial_pos,
                                                     10) # TODO: parameterize this
        x_data.extend(maxima_x)
        y_data.extend(maxima_y)



    # start fit
    sc.fit_geometry(geometry,
                        (np.array(x_data), np.array(y_data)),
                        d_spacings,
                        fit_param)


    print 'Final parameter:'
    print geometry.getFit2D()

    config.activeExperiment.setvalue('Detector Distance', geometry.get_dist())
    config.activeExperiment.center = geometry.getFit2D()['centerX'],geometry.getFit2D()['centerY']
    yield 100

####This thresholds the minimum sdd which is acceptable # print peak # guess which order is the first # threshold accepting peak as higher N # index of the first detected peak #',N,' for calibration...')) # Calculate sample to detector distance for lowest q peak #        self.refinecenter() ####This thresholds the minimum sdd which is acceptable # print peak # Calculate sample to detector distance for lowest q peak #        self.refinecenter() # Refine calibration # d-spacing for Silver Behenate # calculate maxima for every d_spacing # calculate maxima for every d_spacing # TODO: parameterize this # start fit