vig_optimize.cpp

Go to the documentation of this file.

// -*- c-basic-offset: 4 -*-

#include <hugin_config.h>
#include <fstream>
#include <sstream>

#include <algorithms/optimizer/PhotometricOptimizer.h>
#include "ExtractPoints.h"
#include <getopt.h>
#include <hugin_basic.h>
#include <nona/Stitcher.h>

#include <tiff.h>

static void usage(const char* name)
{
    std::cout << name << ": Try to determine the radial vignetting" << std::endl
         << "vig_optimize version " << hugin_utils::GetHuginVersion() << std::endl
         << std::endl
         << "Usage: " << name  << " [options] -o output.pto input.pto" << std::endl
         << "Valid options are:" << std::endl
         << "  -o|--output file   write results to output project" << std::endl
         << "  -v|--verbose       Verbose, print progress messages" << std::endl
         << "  -p n               Number of points to extract" << std::endl
         << "  -s level           Work on downscaled images, every step halves width and height" << std::endl
         << "  -h|--help          Display help (this text)" << std::endl
         << std::endl
         << " Expert and debugging options:" << std::endl
         << "  -i file   Read corresponding points from file" << std::endl
         << "  -w file   Dump corresponding points to file" << std::endl;
}


/*
template<class RIMG>
void importROIImage(ImageImportInfo & info, RIMG & img)
{
    vigra::Rect2D roi(Point2D(info.getPosition()), info.size());
    img.resize(roi);
    importImageAlpha(info, destImage(img.m_image), destImage(img.m_mask));
}
*/


void loadPointsC(FILE* f, std::vector<vigra_ext::PointPairRGB>& vec)
{
    double dummy1, dummy2;
    vigra_ext::PointPairRGB point;
    double i1, i2;

    int n=0;
    do
    {
        n = fscanf(f, " %lf %lf %lf %lf %lf %lf %f %f %f %f %f %f %lf %f %f %lf",
                   &i1, &i2, &(point.p1.x), &(point.p1.y),
                   &(point.p2.x), &(point.p2.y),
                   &(point.i1.red()), &(point.i1.green()), &(point.i1.blue()),
                   &(point.i2.red()), &(point.i2.green()), &(point.i2.blue()),
                   // matlab index.. set to zero
                   &dummy1,
                   &(point.r1), &(point.r2),
                   // edge score at this point, zero so far
                   &(dummy2) );


        point.imgNr1 = hugin_utils::roundi(i1);
        point.imgNr2 = hugin_utils::roundi(i2);
        if (n==16)
        {
            vec.push_back(point);
        }
    }
    while (n == 16);
}

void loadPoints(std::istream& i, std::vector<vigra_ext::PointPairRGB>& vec )
{
    while(!i.eof() && i.good())
    {
        double dummy1, dummy2;
        vigra_ext::PointPairRGB point;
        double i1, i2;
        i >> i1 >> i2
          >> point.p1.x >> point.p1.y
          >> point.p2.x  >> point.p2.y
          >> point.i1.red()  >> point.i1.green()  >> point.i1.blue()
          >> point.i2.red()  >> point.i2.green()  >> point.i2.blue()
          // matlab index.. set to zero
          >> dummy1
          >> point.r1  >> point.r2
          // edge score at this point, zero so far
          >> dummy2;
        point.imgNr1 = hugin_utils::roundi(i1);
        point.imgNr2 = hugin_utils::roundi(i2);
        if (i.good())
        {
            vec.push_back(point);
        }
    }
}

bool hasphotometricParams(HuginBase::Panorama& pano)
{
    HuginBase::OptimizeVector vars = pano.getOptimizeVector();

    for (HuginBase::OptimizeVector::const_iterator it = vars.begin(); it != vars.end(); ++it)
    {
        for (std::set<std::string>::const_iterator itv = (*it).begin();
                itv != (*it).end(); ++itv)
        {
            if ((*itv)[0] == 'E' || (*itv)[0] == 'R' || (*itv)[0] == 'V')
            {
                return true;
            }
        }
    }
    return false;
}


int main(int argc, char* argv[])
{
    AppBase::StreamProgressDisplay progressDisplay(std::cout);

    // parse arguments
    const char* optstring = "hi:o:p:s:vw:";
    static struct option longOptions[] =
    {
        { "output", required_argument, NULL, 'o' },
        { "verbose", no_argument, NULL, 'v' },
        { "help", no_argument, NULL, 'h' },
        0
    };
    int c;

    int pyrLevel=3;
    int verbose = 0;
    int nPoints = 200;
    std::string outputFile;
    std::string outputPointsFile;
    std::string inputPointsFile;
    while ((c = getopt_long(argc, argv, optstring, longOptions, nullptr)) != -1)
    {
        switch (c)
        {
            case 'i':
                inputPointsFile = optarg;
                break;
            case 'o':
                outputFile = optarg;
                break;
            case 'p':
                nPoints = atoi(optarg);
                break;
            case 's':
                pyrLevel = atoi(optarg);
                break;
            case 'v':
                verbose++;
                break;
            case 'h':
                usage(hugin_utils::stripPath(argv[0]).c_str());
                return 1;
            case 'w':
                outputPointsFile = optarg;
                break;
            case ':':
            case '?':
                // missing argument or invalid switch
                return 1;
                break;
            default:
                // this should not happen
                abort();
        }
    };

    unsigned nFiles = argc - optind;
    if (nFiles != 1)
    {
        if (nFiles< 1)
        {
            std::cerr << hugin_utils::stripPath(argv[0]) << ": No project file given." << std::endl;
        }
        else
        {
            std::cerr << hugin_utils::stripPath(argv[0]) << ": Only one project file expected." << std::endl;
        };
        return 1;
    };

    const char* scriptFile = argv[optind];
    HuginBase::Panorama pano;
    if (!pano.ReadPTOFile(scriptFile, hugin_utils::getPathPrefix(scriptFile)))
    {
        return 1;
    };

    // Ensure photometric parameters are selected for optimizaion
    if (!hasphotometricParams(pano))
    {
        std::cerr << "ERROR:no photometric parameters were selected for optimization" << std::endl;
        std::cerr << "please update 'v' line in .pto script and try again." << std::endl;
        return 1;
    }

    // suppress tiff warnings
    TIFFSetWarningHandler(0);

    // todo: handle images with alpha masks
    try
    {
        std::vector<vigra_ext::PointPairRGB> points;

        float imageStepSize = 1 / 255.0f;
        if (inputPointsFile != "" )
        {
            //ifstream ifs(inputPointsFile.c_str());
            //loadPoints(ifs, points);
            FILE* f = fopen(inputPointsFile.c_str(), "r");
            if (f == 0)
            {
                perror("Could not open input point file");
                return 1;
            }
            loadPointsC(f, points);
            fclose(f);
        }
        else
        {
            loadImgsAndExtractPoints(pano, nPoints, pyrLevel, true, progressDisplay, points, verbose, imageStepSize);
        }
        if (verbose)
        {
            std::cout << "\rSelected " << points.size() << " points" << std::endl;
        }
        if (points.empty())
        {
            std::cerr << "Error: no overlapping points found, exiting" << std::endl;
            return 1;
        }

        if (!outputPointsFile.empty())
        {
            std::ofstream of(outputPointsFile.c_str());
            for (std::vector<vigra_ext::PointPairRGB>::iterator it = points.begin(); it != points.end(); ++it)
            {
                of << (*it).imgNr1 << " " << (*it).imgNr2 << " "
                   << (*it).p1.x << " " << (*it).p1.y<< " "
                   << (*it).p2.x << " " << (*it).p2.y<< " "
                   << (*it).i1.red() << " " << (*it).i1.green() << " " << (*it).i1.blue() << " "
                   << (*it).i2.red() << " " << (*it).i2.green() << " " << (*it).i2.blue() << " "
                   // matlab index.. set to zero
                   << 0 << " "
                   << (*it).r1 << " " << (*it).r2 << " "
                   // edge score at this point, zero so far
                   << 0 << std::endl;
            }
        }


        progressDisplay.setMessage("Vignetting Optimization");
        HuginBase::PhotometricOptimizer photoopt(pano, &progressDisplay, pano.getOptimizeVector(), points, imageStepSize);
        photoopt.run();
        //              double error = photoopt.getResultError();

        progressDisplay.taskFinished();

        // save project
        outputFile = hugin_utils::GetOutputFilename(outputFile, scriptFile, "vig");
        pano.WritePTOFile(outputFile, hugin_utils::getPathPrefix(outputFile));

    }
    catch (std::exception& e)
    {
        std::cerr << "caught exception: " << e.what() << std::endl;
        return 1;
    }
    return 0;
}