# Copyright (C) 2020 Cancer Care Associates and <NAME>
# Copyright (C) 2019 Cancer Care Associates

# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at

#     http://www.apache.org/licenses/LICENSE-2.0

# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

from pymedphys._imports import numpy as np

from pymedphys._experimental.vendor.pylinac_vendored import (
    winstonlutz as _pylinac_vendored_wlutz,
)
from pymedphys._experimental.vendor.pylinac_vendored._pylinac_installed import (
    pylinac as _pylinac_installed,
)

from . import transformation as _transformation


def _get_dx_dy_dpi(x, y):
    dx = _convert_grid_to_step_size(x)
    dy = _convert_grid_to_step_size(y)

    if dx == dy:
        dpmm = dx
        dpi = dpmm * 25.4
    else:
        dpi = None

    return dx, dy, dpi


def run_wlutz_raw(
    x, y, image, find_bb=True, pylinac_version=None, fill_errors_with_nan=False
):
    dx, dy, dpi = _get_dx_dy_dpi(x, y)

    WLImage = _get_class_for_version(pylinac_version)
    wl_image = WLImage(image, dpi=dpi)

    nan_coords = [np.nan, np.nan]

    try:
        field_centre = [
            wl_image.field_cax.x * dx + np.min(x),
            wl_image.field_cax.y * dy + np.min(y),
        ]
    except ValueError:
        if fill_errors_with_nan:
            field_centre = nan_coords

        else:
            raise

    if find_bb:
        try:
            bb_centre = [wl_image.bb.x * dx + np.min(x), wl_image.bb.y * dy + np.min(y)]
        except ValueError:
            if fill_errors_with_nan:
                bb_centre = nan_coords
            else:
                raise
    else:
        bb_centre = nan_coords

    return field_centre, bb_centre


def _get_class_for_version(pylinac_version=None):
    if pylinac_version is None:
        pylinac_version = _pylinac_installed.__version__

    VERSION_TO_CLASS_MAP = _pylinac_vendored_wlutz.get_version_to_class_map()
    WLImage = VERSION_TO_CLASS_MAP[pylinac_version]

    return WLImage


def find_bb_only_raw(x, y, image, padding):
    dx, dy, dpi = _get_dx_dy_dpi(x, y)

    WLImage = _pylinac_vendored_wlutz.WLImageCurrent
    wl_image = WLImage(image, dpi=dpi)
    wl_image.set_bounding_box_by_padding(padding)

    bb_centre = [wl_image.bb.x * dx + np.min(x), wl_image.bb.y * dy + np.min(y)]

    return bb_centre


def _convert_grid_to_step_size(x):
    diff_x = np.diff(x)

    dx_mean = np.mean(diff_x)
    dx_deviations = np.abs(diff_x - dx_mean)

    if np.any(dx_deviations > 0.00001):
        raise ValueError(
            "Exactly one grid step size required. Maximum deviation "
            f"from the mean was {np.max(dx_deviations)}."
        )

    return dx_mean


def find_bb_only(x, y, image, edge_lengths, penumbra, field_centre, field_rotation):
    extra_pixels_padding = 20
    out_of_field_padding_factor = 2
    in_field_padding_factor = 2
    bounding_box_padding_factor = out_of_field_padding_factor + in_field_padding_factor

    dx = _convert_grid_to_step_size(x)
    dy = _convert_grid_to_step_size(y)

    x_radius = (
        edge_lengths[0] / 2
        + penumbra * out_of_field_padding_factor
        + extra_pixels_padding * dx
    )
    y_radius = (
        edge_lengths[1] / 2
        + penumbra * out_of_field_padding_factor
        + extra_pixels_padding * dy
    )

    bounding_box_x_padding = (
        np.round(bounding_box_padding_factor * penumbra / dx) + extra_pixels_padding
    )
    bounding_box_y_padding = (
        np.round(bounding_box_padding_factor * penumbra / dy) + extra_pixels_padding
    )

    padding = [bounding_box_x_padding, bounding_box_y_padding]

    x_new = np.arange(-x_radius, x_radius + dx / 2, dx)
    y_new = np.arange(-y_radius, y_radius + dy / 2, dy)
    centralised_image = _transformation.create_centralised_image(
        x, y, image, field_centre, field_rotation, new_x=x_new, new_y=y_new
    )

    # try:
    raw_bb_centre = find_bb_only_raw(x_new, y_new, centralised_image, padding)
    # except Exception as e:
    #     plt.pcolormesh(x_new, y_new, centralised_image, shading="nearest")
    #     plt.axis("equal")

    #     print(e)
    #     plt.show()
    #     field_centre, bb_centre = run_wlutz_raw(x_new, y_new, centralised_image)
    #     print(field_centre)
    #     print(bb_centre)
    #     raise

    bb_centre = _transformation.transform_point(
        raw_bb_centre, field_centre, field_rotation
    )

    return bb_centre


def run_wlutz(
    x,
    y,
    image,
    edge_lengths,
    field_rotation,
    find_bb=True,
    interpolated_pixel_size=0.25,
    pylinac_versions=None,
    fill_errors_with_nan=False,
):
    offset_iter = 10
    current_pylinac_version = _pylinac_installed.__version__

    # By defining a search offset and radius, artefacts that can cause
    # offsets in the pylinac algorithm can be cropped out. See the
    # following issue for more details:

    #    <https://github.com/jrkerns/pylinac/issues/333>

    # By defining the search radius to equal the maximum side length
    # the interpolation region being search over by PyLinac is twice
    # that of the maximum field edge.
    pylinac_offset_calculation = run_wlutz_with_manual_search_definition(
        x,
        y,
        image,
        field_rotation,
        search_radius=None,
        find_bb=False,
        pylinac_versions=[current_pylinac_version],
    )
    previous_search_offset = pylinac_offset_calculation[current_pylinac_version][
        "field_centre"
    ]
    for _ in range(offset_iter):
        pylinac_offset_calculation = run_wlutz_with_manual_search_definition(
            x,
            y,
            image,
            field_rotation,
            search_radius=np.max(edge_lengths),
            search_offset=previous_search_offset,
            find_bb=False,
            pylinac_versions=[current_pylinac_version],
        )
        search_offset = pylinac_offset_calculation[current_pylinac_version][
            "field_centre"
        ]
        if np.allclose(search_offset, previous_search_offset, atol=0.2):
            break

        previous_search_offset = search_offset

    pylinac_calculation_with_offset = run_wlutz_with_manual_search_definition(
        x,
        y,
        image,
        field_rotation,
        search_radius=np.max(edge_lengths),
        search_offset=search_offset,
        find_bb=find_bb,
        interpolated_pixel_size=interpolated_pixel_size,
        pylinac_versions=pylinac_versions,
        fill_errors_with_nan=fill_errors_with_nan,
    )

    return pylinac_calculation_with_offset


def run_wlutz_with_manual_search_definition(
    x,
    y,
    image,
    field_rotation,
    search_radius=None,
    search_offset=None,
    find_bb=True,
    interpolated_pixel_size=0.25,
    pylinac_versions=None,
    fill_errors_with_nan=False,
):
    if search_offset is None:
        search_offset = [np.mean(x), np.mean(y)]

    if search_radius is None:
        new_x = x - search_offset[0]
        new_y = y - search_offset[1]
    else:
        new_x = np.arange(
            -search_radius,  # pylint: disable = invalid-unary-operand-type
            search_radius + interpolated_pixel_size / 2,
            interpolated_pixel_size,
        )
        new_y = new_x

    centralised_image = _transformation.create_centralised_image(
        x, y, image, search_offset, field_rotation, new_x=new_x, new_y=new_y
    )

    if pylinac_versions is None:
        VERSION_TO_CLASS_MAP = _pylinac_vendored_wlutz.get_version_to_class_map()
        pylinac_versions = VERSION_TO_CLASS_MAP.keys()

    results = {}
    for pylinac_version in pylinac_versions:
        raw_field_centre, raw_bb_centre = run_wlutz_raw(
            new_x,
            new_y,
            centralised_image,
            find_bb=find_bb,
            pylinac_version=pylinac_version,
            fill_errors_with_nan=fill_errors_with_nan,
        )

        bb_centre = _transformation.transform_point(
            raw_bb_centre, search_offset, field_rotation
        )
        field_centre = _transformation.transform_point(
            raw_field_centre, search_offset, field_rotation
        )

        results[pylinac_version] = {
            "field_centre": field_centre,
            "bb_centre": bb_centre,
        }

    return results

# Copyright (C) 2020 Cancer Care Associates and <NAME>
# Copyright (C) 2019 Cancer Care Associates

# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at

#     http://www.apache.org/licenses/LICENSE-2.0

# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

from pymedphys._imports import numpy as np

from pymedphys._experimental.vendor.pylinac_vendored import (
    winstonlutz as _pylinac_vendored_wlutz,
)
from pymedphys._experimental.vendor.pylinac_vendored._pylinac_installed import (
    pylinac as _pylinac_installed,
)

from . import transformation as _transformation


def _get_dx_dy_dpi(x, y):
    dx = _convert_grid_to_step_size(x)
    dy = _convert_grid_to_step_size(y)

    if dx == dy:
        dpmm = dx
        dpi = dpmm * 25.4
    else:
        dpi = None

    return dx, dy, dpi


def run_wlutz_raw(
    x, y, image, find_bb=True, pylinac_version=None, fill_errors_with_nan=False
):
    dx, dy, dpi = _get_dx_dy_dpi(x, y)

    WLImage = _get_class_for_version(pylinac_version)
    wl_image = WLImage(image, dpi=dpi)

    nan_coords = [np.nan, np.nan]

    try:
        field_centre = [
            wl_image.field_cax.x * dx + np.min(x),
            wl_image.field_cax.y * dy + np.min(y),
        ]
    except ValueError:
        if fill_errors_with_nan:
            field_centre = nan_coords

        else:
            raise

    if find_bb:
        try:
            bb_centre = [wl_image.bb.x * dx + np.min(x), wl_image.bb.y * dy + np.min(y)]
        except ValueError:
            if fill_errors_with_nan:
                bb_centre = nan_coords
            else:
                raise
    else:
        bb_centre = nan_coords

    return field_centre, bb_centre


def _get_class_for_version(pylinac_version=None):
    if pylinac_version is None:
        pylinac_version = _pylinac_installed.__version__

    VERSION_TO_CLASS_MAP = _pylinac_vendored_wlutz.get_version_to_class_map()
    WLImage = VERSION_TO_CLASS_MAP[pylinac_version]

    return WLImage


def find_bb_only_raw(x, y, image, padding):
    dx, dy, dpi = _get_dx_dy_dpi(x, y)

    WLImage = _pylinac_vendored_wlutz.WLImageCurrent
    wl_image = WLImage(image, dpi=dpi)
    wl_image.set_bounding_box_by_padding(padding)

    bb_centre = [wl_image.bb.x * dx + np.min(x), wl_image.bb.y * dy + np.min(y)]

    return bb_centre


def _convert_grid_to_step_size(x):
    diff_x = np.diff(x)

    dx_mean = np.mean(diff_x)
    dx_deviations = np.abs(diff_x - dx_mean)

    if np.any(dx_deviations > 0.00001):
        raise ValueError(
            "Exactly one grid step size required. Maximum deviation "
            f"from the mean was {np.max(dx_deviations)}."
        )

    return dx_mean


def find_bb_only(x, y, image, edge_lengths, penumbra, field_centre, field_rotation):
    extra_pixels_padding = 20
    out_of_field_padding_factor = 2
    in_field_padding_factor = 2
    bounding_box_padding_factor = out_of_field_padding_factor + in_field_padding_factor

    dx = _convert_grid_to_step_size(x)
    dy = _convert_grid_to_step_size(y)

    x_radius = (
        edge_lengths[0] / 2
        + penumbra * out_of_field_padding_factor
        + extra_pixels_padding * dx
    )
    y_radius = (
        edge_lengths[1] / 2
        + penumbra * out_of_field_padding_factor
        + extra_pixels_padding * dy
    )

    bounding_box_x_padding = (
        np.round(bounding_box_padding_factor * penumbra / dx) + extra_pixels_padding
    )
    bounding_box_y_padding = (
        np.round(bounding_box_padding_factor * penumbra / dy) + extra_pixels_padding
    )

    padding = [bounding_box_x_padding, bounding_box_y_padding]

    x_new = np.arange(-x_radius, x_radius + dx / 2, dx)
    y_new = np.arange(-y_radius, y_radius + dy / 2, dy)
    centralised_image = _transformation.create_centralised_image(
        x, y, image, field_centre, field_rotation, new_x=x_new, new_y=y_new
    )

    # try:
    raw_bb_centre = find_bb_only_raw(x_new, y_new, centralised_image, padding)
    # except Exception as e:
    #     plt.pcolormesh(x_new, y_new, centralised_image, shading="nearest")
    #     plt.axis("equal")

    #     print(e)
    #     plt.show()
    #     field_centre, bb_centre = run_wlutz_raw(x_new, y_new, centralised_image)
    #     print(field_centre)
    #     print(bb_centre)
    #     raise

    bb_centre = _transformation.transform_point(
        raw_bb_centre, field_centre, field_rotation
    )

    return bb_centre


def run_wlutz(
    x,
    y,
    image,
    edge_lengths,
    field_rotation,
    find_bb=True,
    interpolated_pixel_size=0.25,
    pylinac_versions=None,
    fill_errors_with_nan=False,
):
    offset_iter = 10
    current_pylinac_version = _pylinac_installed.__version__

    # By defining a search offset and radius, artefacts that can cause
    # offsets in the pylinac algorithm can be cropped out. See the
    # following issue for more details:

    #    <https://github.com/jrkerns/pylinac/issues/333>

    # By defining the search radius to equal the maximum side length
    # the interpolation region being search over by PyLinac is twice
    # that of the maximum field edge.
    pylinac_offset_calculation = run_wlutz_with_manual_search_definition(
        x,
        y,
        image,
        field_rotation,
        search_radius=None,
        find_bb=False,
        pylinac_versions=[current_pylinac_version],
    )
    previous_search_offset = pylinac_offset_calculation[current_pylinac_version][
        "field_centre"
    ]
    for _ in range(offset_iter):
        pylinac_offset_calculation = run_wlutz_with_manual_search_definition(
            x,
            y,
            image,
            field_rotation,
            search_radius=np.max(edge_lengths),
            search_offset=previous_search_offset,
            find_bb=False,
            pylinac_versions=[current_pylinac_version],
        )
        search_offset = pylinac_offset_calculation[current_pylinac_version][
            "field_centre"
        ]
        if np.allclose(search_offset, previous_search_offset, atol=0.2):
            break

        previous_search_offset = search_offset

    pylinac_calculation_with_offset = run_wlutz_with_manual_search_definition(
        x,
        y,
        image,
        field_rotation,
        search_radius=np.max(edge_lengths),
        search_offset=search_offset,
        find_bb=find_bb,
        interpolated_pixel_size=interpolated_pixel_size,
        pylinac_versions=pylinac_versions,
        fill_errors_with_nan=fill_errors_with_nan,
    )

    return pylinac_calculation_with_offset


def run_wlutz_with_manual_search_definition(
    x,
    y,
    image,
    field_rotation,
    search_radius=None,
    search_offset=None,
    find_bb=True,
    interpolated_pixel_size=0.25,
    pylinac_versions=None,
    fill_errors_with_nan=False,
):
    if search_offset is None:
        search_offset = [np.mean(x), np.mean(y)]

    if search_radius is None:
        new_x = x - search_offset[0]
        new_y = y - search_offset[1]
    else:
        new_x = np.arange(
            -search_radius,  # pylint: disable = invalid-unary-operand-type
            search_radius + interpolated_pixel_size / 2,
            interpolated_pixel_size,
        )
        new_y = new_x

    centralised_image = _transformation.create_centralised_image(
        x, y, image, search_offset, field_rotation, new_x=new_x, new_y=new_y
    )

    if pylinac_versions is None:
        VERSION_TO_CLASS_MAP = _pylinac_vendored_wlutz.get_version_to_class_map()
        pylinac_versions = VERSION_TO_CLASS_MAP.keys()

    results = {}
    for pylinac_version in pylinac_versions:
        raw_field_centre, raw_bb_centre = run_wlutz_raw(
            new_x,
            new_y,
            centralised_image,
            find_bb=find_bb,
            pylinac_version=pylinac_version,
            fill_errors_with_nan=fill_errors_with_nan,
        )

        bb_centre = _transformation.transform_point(
            raw_bb_centre, search_offset, field_rotation
        )
        field_centre = _transformation.transform_point(
            raw_field_centre, search_offset, field_rotation
        )

        results[pylinac_version] = {
            "field_centre": field_centre,
            "bb_centre": bb_centre,
        }

    return results

# Copyright (C) 2020 Cancer Care Associates and <NAME> # Copyright (C) 2019 Cancer Care Associates # Licensed under the Apache License, Version 2.0 (the "License"); # you may not use this file except in compliance with the License. # You may obtain a copy of the License at #     http://www.apache.org/licenses/LICENSE-2.0 # Unless required by applicable law or agreed to in writing, software # distributed under the License is distributed on an "AS IS" BASIS, # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. # See the License for the specific language governing permissions and # limitations under the License. # try: # except Exception as e: #     plt.pcolormesh(x_new, y_new, centralised_image, shading="nearest") #     plt.axis("equal") #     print(e) #     plt.show() #     field_centre, bb_centre = run_wlutz_raw(x_new, y_new, centralised_image) #     print(field_centre) #     print(bb_centre) #     raise # By defining a search offset and radius, artefacts that can cause # offsets in the pylinac algorithm can be cropped out. See the # following issue for more details: #    <https://github.com/jrkerns/pylinac/issues/333> # By defining the search radius to equal the maximum side length # the interpolation region being search over by PyLinac is twice # that of the maximum field edge. # pylint: disable = invalid-unary-operand-type