app_profile.h

Go to the documentation of this file.

/* Gxsm - Gnome X Scanning Microscopy
 * universal STM/AFM/SARLS/SPALEED/... controlling and
 * data analysis software
 * 
 * Copyright (C) 1999,2000,2001,2002,2003 Percy Zahl
 *
 * Authors: Percy Zahl <zahl@users.sf.net>
 * additional features: Andreas Klust <klust@users.sf.net>
 * WWW Home: http://gxsm.sf.net
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307, USA.
 */

/* -*- Mode: C++; indent-tabs-mode: nil; c-basic-offset: 8 c-style: "K&R" -*- */


#ifndef __APP_PROFILE_H
#define __APP_PROFILE_H

typedef struct{
  gchar *mitem;
  int   msk;
  int   neg;
} menuoptionlist;

#define PROFILE_MODE_XGRID    (1<<0)
#define PROFILE_MODE_XLOG     (1<<1)
#define PROFILE_MODE_YGRID    (1<<2)
#define PROFILE_MODE_YLOG     (1<<3)
#define PROFILE_MODE_YLINREG  (1<<4)
#define PROFILE_MODE_POINTS   (1<<5)
#define PROFILE_MODE_CONNECT  (1<<6)
#define PROFILE_MODE_IMPULS   (1<<7)
#define PROFILE_MODE_NOTICS   (1<<8)
#define PROFILE_MODE_LEGEND   (1<<9)
#define PROFILE_MODE_SYMBOLS  (1<<10)
#define PROFILE_MODE_YPSD     (1<<11)

#define PROFILE_SCALE_XAUTO   (1<<0)
#define PROFILE_SCALE_YAUTO   (1<<1)
#define PROFILE_SCALE_XHOLD   (1<<3)
#define PROFILE_SCALE_YHOLD   (1<<4)
#define PROFILE_SCALE_XEXPAND (1<<5)
#define PROFILE_SCALE_YEXPAND (1<<6)

#define UTF8_DEGREE "\302\260"

class ProfileControl;

class ProfileElement{
#define PROFILE_ELEM_HIDE     (1<<0)
public:
  ProfileElement(Scan *sc, ProfileControl *pc, gchar *color=NULL, int Yy=0);
  ProfileElement(ProfileElement *pe, gchar *color=NULL, int Yy=0);
  ~ProfileElement();

  int GetNy(){ return scan->mem2d->GetNy(); };
  void SetY(int Yy=0){ yy=Yy; };
  void SetLastY(){ yy=scan->mem2d->GetNy()-1; };
  void SetMode(long Mode);
  void SetOptions(long Flg);
  
  void calc(int ymode=0);
  void draw(GtkWidget *canvas);

  void GetCurXYc(double *x, double *y, int i){
    if(pn)
      if(i<n && i>=0)
        *x=pn->coords[2*i], *y=pn->coords[2*i+1];
  };

  int GetCurX(double *x, double *y){
    if(pn){
      int i = (int)(n*(*x-pn->coords[0])/(pn->coords[2*n-2]-pn->coords[0]));
      //      while(*x > pn->coords[2*i] && i > 1) --i;
      //      while(*x < pn->coords[2*i] && i < n-2) ++i;
      if(i<n && i>=0)
        *x=pn->coords[2*i], *y=pn->coords[2*i+1];
      return i;
    }
    return 0;
  };

  gchar *GetInfo(int i, int ymode=0){
          double x,y;
          GetCurXYc (&x, &y, i);
          if(i<n && i>=0){
                  Scan *sc = ymode & PROFILE_MODE_YPSD ? psd_scan : scan;
                  x = sc->mem2d->data->GetXLookup(i);
                  return g_strconcat("Cursor (",
                                     ymode & PROFILE_MODE_YPSD ? 
                                     psd_scan->data.Xunit->UsrString (x) :
                                     scan->data.Xunit->UsrString (x),
                                     ", ",
                                     ymode & PROFILE_MODE_YPSD ? 
                                     scan->data.Zunit->UsrString (y) :
                                     ymode & PROFILE_MODE_YLINREG ?
                                     scan->data.Zunit->UsrString (scan->mem2d->GetDataPktLineReg(i,yy)*scan->data.s.dz) :
                                     scan->data.Zunit->UsrString (scan->mem2d->GetDataPkt(i,yy)*scan->data.s.dz),
                                     ")",
                                     NULL);
          }
    else
      return g_strdup("out of scan !");
  };

  double GetValue(int i){
    if(i<n && i>=0)
      return scan->mem2d->GetDataPktLineReg(i,yy)*scan->data.s.dz;
    return 0.;
  };

  int nextMax(int i, int dir, double eps=1.){
    if(i<n && i>=0){
      int k;
      double v=scan->mem2d->GetDataPktLineReg(i,yy);
      if(v > scan->mem2d->GetDataPktLineReg(i+dir,yy)+eps){
        i=nextMin(i+dir,dir,eps);
        v=scan->mem2d->GetDataPktLineReg(i,yy);
      }
      for(k=i+dir; k>=0 && k<n && v < scan->mem2d->GetDataPktLineReg(k,yy)+eps; k+=dir){
        v=scan->mem2d->GetDataPktLineReg(k,yy);
      }
      k-=dir;
      return k;
    }
    else return 0;
  };

  int nextMin(int i, int dir, double eps=1.){
    if(i<n && i>=0){
      int k;
      double v=scan->mem2d->GetDataPktLineReg(i,yy);
      if(v < scan->mem2d->GetDataPktLineReg(i+dir,yy)-eps){
        i=nextMax(i+dir,dir,eps);
        v=scan->mem2d->GetDataPktLineReg(i,yy);
      }
      for(k=i+dir; k>=0 && k<n && v > scan->mem2d->GetDataPktLineReg(k,yy)-eps; k+=dir)
        v=scan->mem2d->GetDataPktLineReg(k,yy);
      k-=dir;
      return k;
    }
    else return 0;
  };

  gchar *GetDeltaInfo(int i, int j, int ymode=0){
          double x1,y1, x2,y2, dx,dz;
          GetCurXYc (&x1, &y1, i);
          GetCurXYc (&x2, &y2, j);
          if(i<n && i>=0 && j<n && j>=0){
                  Scan *sc = ymode & PROFILE_MODE_YPSD ? psd_scan : scan;
                  x1 = sc->mem2d->data->GetXLookup(i);
                  x2 = sc->mem2d->data->GetXLookup(j);
                  dx = x2-x1;
                  gchar *s1 = sc->data.Xunit->UsrString (dx);
                  gchar *s2 = ymode & PROFILE_MODE_YLINREG ?
                          sc->data.Zunit->UsrString (dz = (( sc->mem2d->GetDataPktLineReg(i,yy)
                                                              - sc->mem2d->GetDataPktLineReg(j,yy))
                                                             *sc->data.s.dz))
                          :
                          sc->data.Zunit->UsrString (dz = (( sc->mem2d->GetDataPkt(i,yy) 
                                                               - sc->mem2d->GetDataPkt(j,yy) )
                                                             *sc->data.s.dz));
                  gchar *s3 = g_strdup_printf ("%g"UTF8_DEGREE, 180./M_PI*atan(fabs(dz/dx)));
                  gchar *s = g_strconcat ("Cursor Delta (", s1, ", ", s2, ")",
                                          " Slope: ", s3,
                                          NULL);
                  g_free (s1);
                  g_free (s2);
                  g_free (s3);
                  return s;
          }
          else
                  return g_strdup("out of scan !");
  };

  double ylocmin, ylocmax;
  double xlocmin, xlocmax;
private:
  int n;
  int yy;

  ProfileControl    *pctrl;
  GnomeCanvasPoints *pn;
  GnomeCanvasItem   *pathitem;

  Scan *scan;
  Scan *psd_scan;
  gchar *color;
  long flg;
  long mode;
};

class ProfileControl : public AppBase, public LineProfile1D{

public:
  ProfileControl(gchar *titlestring=NULL, int ChNo=-1);
  ProfileControl(gchar *titlestring, int n, UnitObj *ux, UnitObj *uy, double xmin=0., double xmax=1., gchar *resid=NULL);
  ProfileControl(gchar *filename, gchar *resource_id_string);
  ~ProfileControl();

  void Init(gchar *titlestring, int ChNo, gchar *resid=NULL);

  static void file_open_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_save_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_save_as_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_print1_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_print2_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_print3_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_print4_callback (GtkWidget *widget, ProfileControl *pc);
// experimental 5th & 6th print for xmgrace
  static void file_print5_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_print6_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_activate_callback (GtkWidget *widget, ProfileControl *pc);
  static void file_close_callback (GtkWidget *widget, ProfileControl *pc);

  static void logy_callback (GtkWidget *widget, ProfileControl *pc);
  static void linreg_callback (GtkWidget *widget, ProfileControl *pc);
  static void psd_callback (GtkWidget *widget, ProfileControl *pc);
  static void tics_callback (GtkWidget *widget, ProfileControl *pc);
  static void symbols_callback (GtkWidget *widget, ProfileControl *pc);
  static void legend_callback (GtkWidget *widget, ProfileControl *pc);
  static void nogrid_callback (GtkWidget *widget, ProfileControl *pc);

  static void yhold_callback (GtkWidget *widget, ProfileControl *pc);
  static void yexpand_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yauto_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yupperup_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yupperdn_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Ylowerup_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Ylowerdn_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yzoomin_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yzoomout_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Yset_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Xauto_callback (GtkWidget *widget, ProfileControl *pc);
  static void skl_Xset_callback (GtkWidget *widget, ProfileControl *pc);

  static void sel_next_p1_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_p5_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_p10_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_p50_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_p100_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_p500_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m500_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m100_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m50_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m10_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m5_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_next_m1_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_all_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_first_callback(GtkWidget *widget, ProfileControl *pc);
  static void sel_last_callback(GtkWidget *widget, ProfileControl *pc);

  static void cur_Ashow_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Bshow_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Aleft_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Almax_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Almin_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Armax_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Armin_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Aright_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Bleft_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Blmax_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Blmin_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Brmax_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Brmin_callback (GtkWidget *widget, ProfileControl *pc);
  static void cur_Bright_callback (GtkWidget *widget, ProfileControl *pc);

  static gint cursor_event(GnomeCanvasItem *item, GdkEvent *event, ProfileControl *pc);

  void file_print_callback (GtkWidget *widget, ProfileControl *pc, int index);

  void showCur(int id, int show);
  void moveCur(int id, int dir, int search=0, double ix=0.);

  void nextLine(int n);

  void drawScans(int flg=0);
  void UpdateArea();

  static void destroy_item(GnomeCanvasItem *gci, ProfileControl *pc){ gtk_object_destroy(GTK_OBJECT(gci)); };
  static void kill_elem(ProfileElement *pe, ProfileControl *pc){ delete pe; };
  static void draw_elem(ProfileElement *pe, ProfileControl *pc);

  gint NewData(Scan* sc, int line=-1, int cpyscan=TRUE);
  void AddScan(Scan* scan, int line=0, gchar *col=NULL);
  void AddLine(int line=0, gchar *col=NULL);
  void RemoveScans();

  void SetActive(int flg);

  void SetYrange(double y1, double y2);
  void SetXrange(double x1, double x2);

  void SetXlabel(gchar *xlab = NULL);
  void SetYlabel(gchar *ylab = NULL);
  void SetTitle(gchar *tit);

  void SetMode(long flg) { mode=flg; };
  void SetScaling(long flg) { scaleing = flg; };

  void scan2canvas(double sx, double sy, double &cx, double &cy){
                  cx = (sx-xmin)*cxwidth/xrange;
                  XSM_DEBUG (DBG_L2, "s2cxy"<< (mode& PROFILE_MODE_YLOG));
                  if (mode & PROFILE_MODE_YLOG){
                                  if (fabs (sy)>0.)
                                                  cy = cywidth*(1.-(log (sy)-lmin)/lmaxmin);
                                  else
                                                  cy = 1.;
                  }
                  else
                                  cy = cywidth*(1.-(sy-ymin)/yrange);
  };

  void scan2canvas(double &scx, double &scy){
                  scx = (scx-xmin)*cxwidth/xrange;
                  if(mode & PROFILE_MODE_YLOG){
                                  if (fabs(scy) > 0.)
                                                  scy = cywidth*(1.-(log(scy)-lmin)/lmaxmin);
                                  else
                                                  scy = 1.;
                  }
                  else
                                  scy = cywidth*(1.-(scy-ymin)/yrange);
  };

  double scan2canvasX(double sx){
    return (sx-xmin)*cxwidth/xrange;
  };

  double scan2canvasY(double sy){
                  if (mode & PROFILE_MODE_YLOG){
                                  if (fabs(sy)>0.)
                                                  return cywidth*(1.-(log (sy)-lmin)/lmaxmin);
                                  else
                                                  return 1.;
                  }
                  else
                                  return cywidth*(1.-(sy-ymin)/yrange);
  };
  void ref (){ ++ref_count; };
  void unref (){ --ref_count; };

private:

  void SetSize();
  void drawFrame();
  gint drawTics();
  void addTic(GnomeCanvasItem **tic, GnomeCanvasItem **lab,
              double val, double len, int pos, gchar *fmt=NULL, double tval=0.);

  int    chno;
  int    working;
  int    SklOnly;

  int    statusheight;
  double border;
  double pasize;
  double cxwidth, cywidth;
  double aspect;
  int    papixel;

  double xmin, xmax, xrange;
  double ymin, ymax, yrange;
  double lmaxmin, lmin;

  int CursorsIdx[2];

  long mode, scaleing;

  gchar *xlabel, *xlabel0;
  gchar *ylabel, *ylabel0;
  gchar *title;

  gchar *profile_res_id;

  double tic_x1, tic_x2, tic_y1, tic_y2;
  int tic_ym;
  gchar *xticfmt;
  gchar *yticfmt;

  GtkWidget *canvas;
  GtkWidget *statusbar;
  GtkWidget *linecounter;
  gint statusid;
  int new_data_nl, new_data_lastl; 

  GSList *ScanList;
  ProfileElement *last_pe;

  GnomeCanvasItem *frame, *background, *xaxislabel, *yaxislabel;

#define PC_XTN (2*Xticn+100)
#define PC_YTN (2*Xticn+100)
#define PC_XLN (Xticn+1+100)
#define PC_YLN (Yticn+1+100)
  int Yticn, Xticn;
  int ixt, iyt, ixl, iyl;
  GnomeCanvasItem **Xtics, **Ytics;
  GnomeCanvasItem **Xlabels, **Ylabels;

  GnomeCanvasItem *Cursor[2][2];
  int scount;

  int ref_count;
};

#endif

---