// Gmsh - Copyright (C) 1997-2011 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// bugs and problems to <gmsh@geuz.org>.

#include <string.h>
#include "ReadImg.h"
#include "GmshMessage.h"
#include "PView.h"
#include "PViewDataList.h"
#include "GmshConfig.h"

#if defined(HAVE_FLTK)

#include <FL/Fl_JPEG_Image.H>
#include <FL/Fl_PNM_Image.H>
#include <FL/Fl_PNG_Image.H>
#include <FL/Fl_BMP_Image.H>

// from an image, we create a post-procession view

static PViewDataList *Img2Data(Fl_RGB_Image &img_init, int quads=1,
                               int resizex=0, int resizey=0) 
{
  img_init.desaturate(); // convert to grayscale

  // resize if necessary
  Fl_RGB_Image *img;
  if(!resizex || !resizey)
    img = (Fl_RGB_Image*)img_init.copy();
  else
    img = (Fl_RGB_Image*)img_init.copy(resizex, resizey);

  const unsigned char *data = img->array;
  int height = img->h();
  int width = img->w();
  int dim = img->d();

  if(dim != 1) {
    Msg::Error("Unable to obtain one-channel image");
    return 0;
  }

  PViewDataList *d = new PViewDataList();

  double z = 0.;
  for(int i = 0; i < height - 1; i++) {
    const unsigned char *a = data + i * width * dim;
    const unsigned char *a1 = data + (i + 1) * width * dim;
    double y = height - i - 1;
    double y1 = height - i - 2;
    for(int j = 0; j < width - 1; j++) {
      double x = j;
      double x1 = j + 1;
      double val1 = (double)a[j]/255.;
      double val2 = (double)a1[j]/255.;
      double val3 = (double)a1[j + 1]/255.;
      double val4 = (double)a[j + 1]/255.;
      if(quads){ // generate quads
        d->SQ.push_back(x); d->SQ.push_back(x); 
        d->SQ.push_back(x1); d->SQ.push_back(x1);
        d->SQ.push_back(y); d->SQ.push_back(y1);
        d->SQ.push_back(y1); d->SQ.push_back(y);
        d->SQ.push_back(z); d->SQ.push_back(z);
        d->SQ.push_back(z); d->SQ.push_back(z);
        d->SQ.push_back(val1); d->SQ.push_back(val2);
        d->SQ.push_back(val3); d->SQ.push_back(val4);
        d->NbSQ++;
      }
      else{ // generate triangles
        d->ST.push_back(x); d->ST.push_back(x); d->ST.push_back(x1);
        d->ST.push_back(y); d->ST.push_back(y1); d->ST.push_back(y1);
        d->ST.push_back(z); d->ST.push_back(z); d->ST.push_back(z);
        d->ST.push_back(val1); d->ST.push_back(val2); d->ST.push_back(val3);
        d->NbST++;
        d->ST.push_back(x); d->ST.push_back(x1); d->ST.push_back(x1);
        d->ST.push_back(y); d->ST.push_back(y1); d->ST.push_back(y);
        d->ST.push_back(z); d->ST.push_back(z); d->ST.push_back(z);
        d->ST.push_back(val1); d->ST.push_back(val3); d->ST.push_back(val4);
        d->NbST++;
      }
    }
  }
  delete img;
  return d;
}

static int EndPos(const char *name, PViewData *d)
{
  if(!d) return 0;
  char name_pos[256], title[256];
  strcpy(name_pos, name);
  strcat(name_pos, ".pos");
  int i;
  for(i = strlen(name) - 1; i >= 0; i--){
    if(name[i] == '/' || name[i] == '\\')
      break;
  }
  if(i <= 0)
    strcpy(title, name);
  else
    strcpy(title, &name[i + 1]);
  d->setName(title);
  d->setFileName(name_pos);
  if(d->finalize()){
    new PView(d);
    return 1;
  }
  else{
    delete d;
    return 0;
  }
}

int read_pnm(std::string fileName) 
{
  Fl_PNM_Image img(fileName.c_str());
  return EndPos(fileName.c_str(), Img2Data(img));
}

int read_jpeg(std::string fileName) 
{
  Fl_JPEG_Image img(fileName.c_str());
  return EndPos(fileName.c_str(), Img2Data(img));
}

int read_png(std::string fileName) 
{
  Fl_PNG_Image img(fileName.c_str());
  return EndPos(fileName.c_str(), Img2Data(img));
}

int read_bmp(std::string fileName) 
{
  Fl_BMP_Image img(fileName.c_str());
  return EndPos(fileName.c_str(), Img2Data(img));
}

#else

int read_pnm(std::string fileName){ return 0; }
int read_jpeg(std::string fileName){ return 0; }
int read_png(std::string fileName){ return 0; }
int read_bmp(std::string fileName){ return 0; }

#endif
