#include "hisview_view.h"
#include "hisview_misc.h"
#include "hisview_widgets.h"
#include <stdio.h>

static char dumstr[80];

int Hmin=1, Hmax=100, Smin=1, Smax=100, Vmin=1, Vmax=100;

void _2DfloatTo2Dbyte(startf,
		      xmin, xmax,
		      ymin, ymax,
		      startb, scalex, scaley,
		      clip)
float  *startf;      /* Addresse de depart des donnees 2Dfloat a */
		     /* extraire */
int
  xmin,              /* Borne inferieure sur la 1ere dimension */
  xmax,              /* Borne superieure sur la 1ere dimension */
  ymin,              /* Borne inferieure sur la 2eme dimension */
  ymax;              /* Borne superieure sur la 2eme dimension */

unsigned char *startb; /* Addresse de depart de la zone memoire */
                     /* contenant les valeurs de pixel */

int
  scalex,            /* Facteur d'echelle selon la 1ere dimension */
  scaley;            /* Facteur d'echelle selon la 2eme dimension */

float clip[2];
{
  int i, x, y, l;
  int width, height;
  float sc, tested;
  unsigned char *p;

  if (clip[1] == clip[0]) return;

  sc = ((float) (ncolors - 1)) / (clip[1] - clip[0]);
 
  width = (xmax - xmin + 1)*scalex;
  height = (ymax - ymin + 1)*scaley;

  startb = startb + (height-1)*width;

  for (y = ymin ; y <= ymax ; y++) {
    p = startb;
    for (x = xmin ; x <= xmax ; x++) {
      tested = *startf;

      l = (int) (sc*(tested - clip[0]));
      
      if (l < 0) {
	l = 0;
      }
      else {
	if (l > (ncolors-1)) l = ncolors - 1;
      }
      memset(startb, colmap[l], scalex);
      startb += scalex;
      startf++;
    }

    startb -=2*width;
    for (i = 1; i < scaley; i++) {
      (void) memcpy(startb, p, width);
      startb -= width;
    }
  }
}

int mode = 0;
extern int myowncmap;

void draw_wedge(float clip[2])
{
  int j, k;
  char c[12];
  Font fonte;
  float f;
  int offset;
  Pixel col;

  fonte = XLoadFont(display, "6x10");
  XSetFont(display, gc, fonte);

  offset = (cpixmapdim[1] - 2*ncolors)/2;

  XSetForeground(display, gc, BlackPixel(display, screen));
  XFillRectangle(display, cpixmap, gc, 0, 0, cpixmapdim[0], cpixmapdim[1]);

  for (k=0; k<ncolors; k++) {
    XSetForeground(display, gc, colmap[k]);
    XFillRectangle(display, cpixmap, gc,
		   cpixmapdim[0] - 30, 
		   cpixmapdim[1] - 2*k - 2 - offset,
		   30, 2);
  }

  XSetForeground(display, gc, WhitePixelOfScreen(ecran));

  switch(mode) {
  case 0: case 1:
    for (j = 0; j < nhue; j++) {
      XDrawLine(display, cpixmap, gc, 
		cpixmapdim[0] - 39, 
		cpixmapdim[1] - 2*nsat*nval*j - 1 - offset,
		cpixmapdim[0] - 31, 
		cpixmapdim[1] - 2*nsat*nval*j - 1 - offset);
      
      
      sprintf(c, "%3.0f", valhue[j]);
      XDrawString(display, cpixmap, gc, 
		  cpixmapdim[0] - 65, 
		  cpixmapdim[1] - 2*nsat*nval*j + 2 - offset,
		  c, strlen(c));
      
      for (k = 1; k < nval; k++)
	XDrawLine(display, cpixmap, gc, 
		  cpixmapdim[0] - 35, 
		  cpixmapdim[1] - 2*nsat*nval*j - 1 - 2*nsat*k - offset,
		  cpixmapdim[0] - 31, 
		  cpixmapdim[1] - 2*nsat*nval*j - 1 - 2*nsat*k - offset);
    }
    break;

  case 2: case 3:
    for (j = 0; j <= nhue; j++) {
      XDrawLine(display, cpixmap, gc, 
		cpixmapdim[0] - 39, 
		cpixmapdim[1] - 2*nsat*nval*j - offset,
		cpixmapdim[0] - 31, 
		cpixmapdim[1] - 2*nsat*nval*j - offset);
      
      
      sprintf(c, "%2.1f", (float) j/nhue);
      XDrawString(display, cpixmap, gc, 
		  cpixmapdim[0] - 65, 
		  cpixmapdim[1] - 2*nsat*nval*j + 2 - offset,
		  c, strlen(c));
    }
    break;
  }
}

scale_hsv()
{
  XImage *image;
  char *pixels;

  pixels = (char*) malloc(size[0]*size[1]*sizeof(char));
  
  hsvimage(Himage, Simage, Vimage,
	   1, npix[0],
	   1, npix[1],
	   pixels, zoom, zoom,
	   Hclip, Sclip, Vclip);
  
  image = XCreateImage(display, visual, depth,
		       ZPixmap, 0, pixels, size[0], size[1], 8, 0);
  
  free(pixels);

  XPutImage(display, HSVpixmap, gc, image, 0, 0, 0, 0, size[0], size[1]);

}

scale_hue(Boolean bool)
{
  XImage *image;
  char *pixels;


  pixels = (char*) malloc(size[0]*size[1]*sizeof(char));

  _2DfloatTo2Dbyte(Himage,
		   1, npix[0],
		   1, npix[1],
		   pixels, zoom, zoom,
		   Hclip); 

  image = XCreateImage(display, visual, depth,
		       ZPixmap, 0, pixels, size[0], size[1], 8, 0);

  free(pixels);
  
  XPutImage(display, Hpixmap, gc, image, 0, 0, 0, 0, size[0], size[1]);

  scale_hsv();

  if (bool) 
    XCopyArea(display, Hpixmap, pixmap, gc,
	      0, 0, pixmapdim[0], pixmapdim[1], 0, 0);
  else 
    XCopyArea(display, HSVpixmap, pixmap, gc,
	      0, 0, pixmapdim[0], pixmapdim[1], 0, 0);

  XCopyArea(display, pixmap, XtWindow(radec), gc,
	    0, 0, pixmapdim[0], pixmapdim[1], 0, 0);


}

scale_sat()
{
  XImage *image;
  char *pixels;

  pixels = (char*) malloc(size[0]*size[1]*sizeof(char));

  _2DfloatTo2Dbyte(Simage,
		   1, npix[0],
		   1, npix[1],
		   pixels, zoom, zoom,
		   Sclip); 

  image = XCreateImage(display, visual, depth,
		       ZPixmap, 0, pixels, size[0], size[1], 8, 0);
  
  free(pixels);

  XPutImage(display, Spixmap, gc, image, 0, 0, 0, 0, size[0], size[1]);

  scale_hsv();
  
  XCopyArea(display, Spixmap, pixmap, gc,
	    0, 0, pixmapdim[0], pixmapdim[1], 0, 0);

  XCopyArea(display, pixmap, XtWindow(radec), gc,
	    0, 0, pixmapdim[0], pixmapdim[1], 0, 0);

}

scale_val()
{
  XImage *image;
  char *pixels;

  pixels = (char*) malloc(size[0]*size[1]*sizeof(char));

  _2DfloatTo2Dbyte(Vimage,
		   1, npix[0],
		   1, npix[1],
		   pixels, zoom, zoom,
		   Vclip); 

  image = XCreateImage(display, visual, depth,
		       ZPixmap, 0, pixels, size[0], size[1], 8, 0);
  
  free(pixels);

  XPutImage(display, Vpixmap, gc, image, 0, 0, 0, 0, size[0], size[1]);

  scale_hsv();
  
  XCopyArea(display, Vpixmap, pixmap, gc,
	    0, 0, pixmapdim[0], pixmapdim[1], 0, 0);

  XCopyArea(display, pixmap, XtWindow(radec), gc,
	    0, 0, pixmapdim[0], pixmapdim[1], 0, 0);
}

static int nosense;
 
void change_clip(w, tag, cbs)
Widget w;
int *tag;
XmScaleCallbackStruct *cbs;
{

  Pixel col;
  static int firstpass = 1;
  static float newclip, Hclipsav[2], Sclipsav[2], Vclipsav[2];

  if (firstpass) {
    Hclipsav[0] = Hclip[0];
    Hclipsav[1] = Hclip[1];
    Sclipsav[0] = Sclip[0];
    Sclipsav[1] = Sclip[1];
    Vclipsav[0] = Vclip[0];
    Vclipsav[1] = Vclip[1];
    firstpass = 0;
  }

  nosense=0;

  switch (*tag) {
  case k_lowcut:
    switch(mode) {
    case 0: case 1:

      newclip = Hclipsav[0] + 
	((Hclipsav[1] - Hclipsav[0])/99.)*(cbs->value - 1);
      Hclip[0] = newclip;
      
      Hmin  = cbs->value;
      
      if (Hmin == Hmax) {
	XmTextSetString(lowcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "lowcut = %f", Hclip[0]);
      XmTextSetString(lowcut, dumstr);
      sprintf(dumstr, "highcut = %f", Hclip[1]);
      XmTextSetString(highcut, dumstr);
      break;
    case 2:
      newclip = Sclipsav[0] + 
	((Sclipsav[1] - Sclipsav[0])/99.)*(cbs->value - 1);
      Sclip[0] = newclip;

      Smin  = cbs->value;
      if (Smin == Smax) {
	XmTextSetString(lowcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "lowcut = %f", Sclip[0]);
      XmTextSetString(lowcut, dumstr);
      sprintf(dumstr, "highcut = %f", Sclip[1]);
      XmTextSetString(highcut, dumstr);

      break;
    case 3:
      newclip = Vclipsav[0] + 
	((Vclipsav[1] - Vclipsav[0])/99.)*(cbs->value - 1);
      Vclip[0] = newclip;

      Vmin  = cbs->value;
      if (Vmin == Vmax) {
	XmTextSetString(lowcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "lowcut = %f", Vclip[0]);
      XmTextSetString(lowcut, dumstr);
      sprintf(dumstr, "highcut = %f", Vclip[1]);
      XmTextSetString(highcut, dumstr);

      break;
    }      
    break;

  case k_highcut:
    switch(mode) {
    case 0: case 1:
      newclip = Hclipsav[0] + 
	((Hclipsav[1] - Hclipsav[0])/99.)*(cbs->value - 1);
      Hclip[1] = newclip;

      Hmax  = cbs->value;
      if (Hmin == Hmax) {
	XmTextSetString(highcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "highcut = %f", Hclip[1]);
      XmTextSetString(highcut, dumstr);
      sprintf(dumstr, "lowcut = %f", Hclip[0]);
      XmTextSetString(lowcut, dumstr);
      break;
    case 2:
      newclip = Sclipsav[0] + 
	((Sclipsav[1] - Sclipsav[0])/99.)*(cbs->value - 1);
      Sclip[1] = newclip;

      Smax  = cbs->value;
      if (Smin == Smax) {
	XmTextSetString(highcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "highcut = %f", Sclip[1]);
      XmTextSetString(highcut, dumstr);
      sprintf(dumstr, "lowcut = %f", Sclip[0]);
      XmTextSetString(lowcut, dumstr);

      break;
    case 3:
      newclip = Vclipsav[0] + 
	((Vclipsav[1] - Vclipsav[0])/99.)*(cbs->value - 1);
      Vclip[1] = newclip;

      Vmax  = cbs->value;
      if (Vmin == Vmax) {
	XmTextSetString(highcut, "Doesn't make sense ...");
	nosense = 1;
	return;
      }

      sprintf(dumstr, "highcut = %f", Vclip[1]);
      XmTextSetString(highcut, dumstr);
      sprintf(dumstr, "lowcut = %f", Vclip[0]);
      XmTextSetString(lowcut, dumstr);

      break;
    }      
    break;
  }
}

void scale_image(w, tag, cbs)
Widget w;
int *tag;
XmScaleCallbackStruct *cbs;
{

  Pixel col;

  if (nosense) return;

  XtVaSetValues(lowscale,
		XmNsensitive, 0,
		NULL);

  XtVaSetValues(highscale,
		XmNsensitive, 0,
		NULL);

  switch(mode) {
  case 0: case 1:
    scale_hue(mode);
    break;
  case 2:
    scale_sat();
    break;
  case 3:
    scale_val();
    break;
  }      

  XtVaSetValues(lowscale,
		XmNsensitive, 1,
		NULL);

  XtVaSetValues(highscale,
		XmNsensitive, 1,
		NULL);

}