/******************************************************************************
 * $Id: test_alg.cpp 5118c4d047129091c940f65a05582111aabf9a2b 2018-06-02 23:04:15 +0200 Even Rouault $
 *
 * Project:  GDAL algorithms
 * Purpose:  Test alg
 * Author:   Even Rouault, even.rouault at spatialys.com
 *
 ******************************************************************************
 * Copyright (c) 2016, Even Rouault <even.rouault at spatialys.com>
 *
 * Permission is hereby granted, free of charge, to any person obtaining a
 * copy of this software and associated documentation files (the "Software"),
 * to deal in the Software without restriction, including without limitation
 * the rights to use, copy, modify, merge, publish, distribute, sublicense,
 * and/or sell copies of the Software, and to permit persons to whom the
 * Software is furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included
 * in all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
 * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
 * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
 * DEALINGS IN THE SOFTWARE.
 ****************************************************************************/

#include "gdal_unit_test.h"

#include "cpl_conv.h"

#include "gdal_alg.h"
#include "gdalwarper.h"
#include "gdal_priv.h"

namespace tut
{
    // Common fixture with test data
    struct test_alg_data
    {
    };

    // Register test group
    typedef test_group<test_alg_data> group;
    typedef group::object object;
    group test_alg_group("ALG");

    typedef struct
    {
      double dfLevel;
      int nPoints;
      double x;
      double y;
    } writeCbkData;

    static CPLErr writeCbk( double dfLevel, int nPoints,
                            double *padfX, double *padfY, void * userData )
    {
        writeCbkData* data = (writeCbkData*)userData;
        data->dfLevel = dfLevel;
        data->nPoints = nPoints;
        if( nPoints == 1 )
        {
            data->x = padfX[0];
            data->y = padfY[0];
        }
        return CE_None;
    }

    // Dummy test
    template<>
    template<>
    void object::test<1>()
    {
        writeCbkData data;
        memset(&data, 0, sizeof(data));
        GDALContourGeneratorH hCG = GDAL_CG_Create(1,1,FALSE,0,1,0, writeCbk, &data);
        double scanline[] = { 0 };
        ensure_equals(GDAL_CG_FeedLine(hCG, scanline), CE_None);
        ensure_equals(data.dfLevel, 0);
        ensure_equals(data.nPoints, 0);
        ensure_approx_equals(data.x, 0.0);
        ensure_approx_equals(data.y, 0.0);
        GDAL_CG_Destroy(hCG);
    }

    // GDALWarpResolveWorkingDataType: default type
    template<>
    template<>
    void object::test<2>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );
        GDALDestroyWarpOptions(psOptions);
    }

    // GDALWarpResolveWorkingDataType: do not change user specified type
    template<>
    template<>
    void object::test<3>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->eWorkingDataType = GDT_CFloat64;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_CFloat64 );
        GDALDestroyWarpOptions(psOptions);
    }

    // GDALWarpResolveWorkingDataType: effect of padfSrcNoDataReal
    template<>
    template<>
    void object::test<4>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->nBandCount = 1;
        psOptions->padfSrcNoDataReal =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfSrcNoDataReal[0] = 0.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = -1.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Int16 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 2.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 256.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_UInt16 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 2.5;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Float32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 2.12345678;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Float64 );

        GDALDestroyWarpOptions(psOptions);
    }

    // GDALWarpResolveWorkingDataType: effect of padfSrcNoDataImag
    template<>
    template<>
    void object::test<5>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->nBandCount = 1;
        psOptions->padfSrcNoDataReal =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfSrcNoDataImag =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfSrcNoDataReal[0] = 0.0;
        psOptions->padfSrcNoDataImag[0] = 0.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 0.0;
        psOptions->padfSrcNoDataImag[0] = 1.0;
        GDALWarpResolveWorkingDataType(psOptions);
        // Could probably be CInt16
        ensure_equals( psOptions->eWorkingDataType, GDT_CInt32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 0.0;
        psOptions->padfSrcNoDataImag[0] = 1.5;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_CFloat32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfSrcNoDataReal[0] = 0.0;
        psOptions->padfSrcNoDataImag[0] = 2.12345678;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_CFloat64 );

        GDALDestroyWarpOptions(psOptions);
    }

    // GDALWarpResolveWorkingDataType: effect of padfDstNoDataReal
    template<>
    template<>
    void object::test<6>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->nBandCount = 1;
        psOptions->padfDstNoDataReal =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfDstNoDataReal[0] = 0.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = -1.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Int16 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = 2.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = 256.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_UInt16 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = 2.5;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Float32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = 2.12345678;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Float64 );

        GDALDestroyWarpOptions(psOptions);
    }

    // GDALWarpResolveWorkingDataType: effect of padfDstNoDataImag
    template<>
    template<>
    void object::test<7>()
    {
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->nBandCount = 1;
        psOptions->padfDstNoDataReal =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfDstNoDataImag =
            static_cast<double*>(CPLMalloc(sizeof(double)));
        psOptions->padfDstNoDataReal[0] = 0.0;
        psOptions->padfDstNoDataImag[0] = 0.0;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_Byte );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataReal[0] = 0.0;
        psOptions->padfDstNoDataImag[0] = 1.0;
        GDALWarpResolveWorkingDataType(psOptions);
        // Could probably be CInt16
        ensure_equals( psOptions->eWorkingDataType, GDT_CInt32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataImag[0] = 0.0;
        psOptions->padfDstNoDataImag[0] = 1.5;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_CFloat32 );

        psOptions->eWorkingDataType = GDT_Unknown;
        psOptions->padfDstNoDataImag[0] = 0.0;
        psOptions->padfDstNoDataImag[0] = 2.12345678;
        GDALWarpResolveWorkingDataType(psOptions);
        ensure_equals( psOptions->eWorkingDataType, GDT_CFloat64 );

        GDALDestroyWarpOptions(psOptions);
    }

    // Test GDALAutoCreateWarpedVRT() with creatino of an alpha band
    template<> template<> void object::test<8>()
    {
        GDALDatasetUniquePtr poDS(
            GDALDriver::FromHandle(
                GDALGetDriverByName("MEM"))->Create("", 1, 1, 1, GDT_Byte, nullptr));
        poDS->SetProjection( SRS_WKT_WGS84);
        double adfGeoTransform[6] = { 10, 1, 0, 20, 0, -1 };
        poDS->SetGeoTransform(adfGeoTransform);
        GDALWarpOptions* psOptions = GDALCreateWarpOptions();
        psOptions->nDstAlphaBand = 2;
        GDALDatasetH hWarpedVRT =
            GDALAutoCreateWarpedVRT(GDALDataset::ToHandle(poDS.get()), nullptr, nullptr,
                                GRA_NearestNeighbour, 0.0,
                                psOptions);
        ensure( hWarpedVRT != nullptr );
        ensure_equals( GDALGetRasterCount(hWarpedVRT), 2 );
        ensure_equals( GDALGetRasterColorInterpretation(
            GDALGetRasterBand(hWarpedVRT, 2)), GCI_AlphaBand );
        GDALDestroyWarpOptions(psOptions);
        GDALClose(hWarpedVRT);
    }


} // namespace tut
