findlocmax.C

Go to the documentation of this file.

/* Gnome gxsm - Gnome X Scanning Microscopy
 * universal STM/AFM/SARLS/SPALEED/... controlling and
 * data analysis software
 * 
 * Gxsm Plugin Name: dummy.C
 * ========================================
 * 
 * Copyright (C) 1999 The Free Software Foundation
 *
 * Authors: Percy Zahl <zahl@fkp.uni-hannover.de>
 * additional features: Andreas Klust <klust@fkp.uni-hannover.de>
 *
 * 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.
 */


/* Please do not change the Begin/End lines of this comment section!
 * this is a LaTeX style section used for auto generation of the PlugIn Manual 
 * Chapter. Add a complete PlugIn documentation inbetween the Begin/End marks!
 * All "% PlugInXXX" commentary tags are mandatory
 * All "% OptPlugInXXX" tags are optional
 * --------------------------------------------------------------------------------
% BeginPlugInDocuSection
% PlugInDocuCaption: Find local maxima -- SARLS
% PlugInName: findlocmax
% PlugInAuthor: M. Langer?
% PlugInAuthorEmail: stefan_fkp@users.sf.net
% PlugInMenuPath: Math/Misc/Find Loc Max

% PlugInDescription
Finds local maximan in a scan. One (unknown to me, PZ) unknown guy
(M.Langer?) of the SARLS group in hannover wrote this PlugIn for
searching of local maxima. It looks somehow specialized to me.

% PlugInUsage
It is only available from \GxsmMenu{Math/Misc/Find Loc Max} if the
instrument type is set to SARLS or a reload of all PlugIns is forced.

% OptPlugInSources
The active channel is used as data source.

%% OptPlugInObjects
%A optional rectangle is used for data extraction...

%% OptPlugInDest
%The computation result is placed into an existing math channel, else
%into a new created math channel.

% OptPlugInFiles
Saves a list of local maxima in Ascii format (List of X Y Z) to a file.

%% OptPlugInRefs
%Any references?

%% OptPlugInKnownBugs
%Are there known bugs? List! How to work around if not fixed?

% OptPlugInNotes
I would appreciate if one of the current SARLS group members could
figure out more about this piece of code.

%% OptPlugInHints
%Any tips and tricks?

% EndPlugInDocuSection
 * -------------------------------------------------------------------------------- 
 */


#include <gtk/gtk.h>
#include "config.h"
#include "gxsm/plugin.h"
 
using namespace std;


static void findlocmax_init( void );
static void findlocmax_about( void );
static void findlocmax_configure( void );
static void findlocmax_cleanup( void );
static gboolean findlocmax_run( Scan *Src, Scan *Dest );

GxsmPlugin findlocmax_pi = {
  NULL,
  NULL,
  0,
  NULL,
  "Findlocmax-M1S-Misc",
  "+SARLS",
  NULL,
  "Nobody?, PZ",
  N_("_Math/_Misc/"),
  N_("Find Loc Max"),
  N_("Sorry, no help for findlocmax filter!"),
  "no more info",
  NULL,
  NULL,
  findlocmax_init,
  NULL,
  findlocmax_about,
  findlocmax_configure,
  NULL,
  findlocmax_cleanup
};

GxsmMathOneSrcPlugin findlocmax_m1s_pi = {
  findlocmax_run
};

static const char *about_text = N_("Gxsm Findlocmax Plugin\n\n"
                                   "Who has coded this one? -- PZ");

GxsmPlugin *get_gxsm_plugin_info ( void ){ 
  findlocmax_pi.description = g_strdup_printf(N_("Gxsm MathOneArg findlocmax plugin %s"), VERSION);
  return &findlocmax_pi; 
}

GxsmMathOneSrcPlugin *get_gxsm_math_one_src_plugin_info( void ) { 
  return &findlocmax_m1s_pi; 
}

static void findlocmax_init(void)
{
  PI_DEBUG (DBG_L2, "Findlocmax Plugin Init");
}

static void findlocmax_about(void)
{
  const gchar *authors[] = { findlocmax_pi.authors, NULL};
  gtk_widget_show(gnome_about_new ( findlocmax_pi.name,
                                    VERSION,
                                    N_("(C) 2000 the Free Software Foundation"),
                                    about_text,
                                    authors,
                                    NULL, NULL, NULL
                                    ));
}

static void findlocmax_configure(void)
{
  if(findlocmax_pi.app)
    findlocmax_pi.app->message("Findlocmax Plugin Configuration");
}

static void findlocmax_cleanup(void)
{
  PI_DEBUG (DBG_L2, "Findlocmax Plugin Cleanup");
}

static gboolean findlocmax_run(Scan *Src, Scan *Dest)
{
  findlocmax_pi.app->message("Sorry Findlocmax Plugin is buggy, fix it! PZ");
  return MATH_OK;

/* Suche lokaler Maxima im Fokusbereich (Kasten mit 10µm*10µm) */
/* Version 0.2 ML  (23/11/98)                                  */

  int line,col,fokline,fokcol,fokxpoints,fokypoints;
  int peakx,peaky,newpeakx,newpeaky;
  double referenzpeak;
  ofstream f;
  int   i;
  const char *fname;
  PI_DEBUG (DBG_L2, "Suche lokaler Maxima und Ausgabe in ASCII-Datei"); 
  fokxpoints=(int)(50000/Src->data.s.dx);
  if (fokxpoints==0) fokxpoints=1;
  fokypoints=(int)(50000/Src->data.s.dy);
  if (fokypoints==0) fokypoints=1;

  /* Empty-Scan erzeugen, alle Datenpunkte Null setzen */

  for(line=0;line<Dest->data.s.ny;line++)
   for(col=0;col<Dest->data.s.nx;col++)
      Dest->mem2d->PutDataPkt(0,col, line);     
  
  /* Peaks aufspüren, mit einem Mindestabstand = Fokus = 10µm <- Variablen: fokxpoints, fokypoints */
  /* doppelte Peaks werden leider noch nicht rausgeworfen */

  for(line=0;line<(Src->mem2d->GetNy()-fokypoints);line=line+(fokypoints/2))
    for(col=0;col<(Src->mem2d->GetNx()-fokxpoints);col=col+(fokxpoints/2)){      
      peakx=col;
      peaky=line;  
      referenzpeak=Src->mem2d->GetDataPkt(col,line);            
      do{
        newpeakx=peakx; newpeaky=peaky;
          for(fokcol=(int)(newpeakx-fokxpoints); fokcol<=(newpeakx+fokxpoints); fokcol++)
        for(fokline=(int)(newpeaky-fokypoints); fokline<=(newpeaky+fokypoints); fokline++)
            if((fokcol>=0) && (fokcol<Src->mem2d->GetNx()) && (fokline>=0) && (fokline<Src->mem2d->GetNy())){
              if(referenzpeak<=(double)Src->mem2d->GetDataPkt(fokcol,fokline)){
                referenzpeak=Src->mem2d->GetDataPkt(fokcol,fokline);
                peakx=fokcol;
                peaky=fokline;
                }                        
              }   
        }while(((newpeakx-peakx)!=0) && ((newpeaky-peaky)!=0));
      for(fokline=(newpeaky-fokypoints); fokline<=(newpeaky+fokypoints); fokline++)
        for(fokcol=(newpeakx-fokxpoints); fokcol<=(newpeakx+fokxpoints); fokcol++)
          if((fokcol>=0) && (fokcol<Dest->mem2d->GetNx()) && (fokline>=0) && (fokline<Dest->mem2d->GetNy()))
            if(0<Dest->mem2d->GetDataPkt(fokcol,fokline))
              referenzpeak=0;
      Dest->mem2d->PutDataPkt((SHT)referenzpeak,newpeakx,newpeaky);     
    }    
        
  /* File Selector */

  do{
    fname = findlocmax_pi.app->file_dialog("asc File to Save LOCAL MAXIMA","","*.asc","");
    if (fname == NULL) return 1;
    if (strncmp(fname+strlen(fname)-4,".asc",4))
      strcat((char*)fname, ".asc");
    PI_DEBUG (DBG_L2, "Save File: " << fname);
    i=2;
    f.open(fname, ios::in);
    if(f.good()){
      f.close();
      if((i=XSM_SHOW_CHOICE(WRN_WARNING, fname, WRN_FILEEXISTS, 3, L_CANCEL, L_OVERWRITE, L_RETRY, 1)) == 1)
        return 1;
    }
  }while(i!=2);
  
  time_t t;
  time(&t); 
  f.open(fname, ios::out | ios::trunc);
  if(!f.good()){
    XSM_SHOW_ALERT(ERR_SORRY, ERR_FILEWRITE,fname,1);
    return 1;
  }
  f << "# LOKALE MAXIMA ASCII-Data\n" << "# Date: " << ctime(&t) << "# " << fname << endl;
  f << "# Fokuspoints x = " << fokxpoints << endl;
  f << "# Fokuspoints y = " << fokypoints << endl;
  f << "# Datenpunktabstand [A]  dx = " << Src->data.s.dx << endl;
  f << "# Datenpunktabstand [A]  dy = " << Src->data.s.dy << endl;
  f << "# x   y   Intensität" << endl ;

  for(line=0;line<Dest->data.s.ny;line++)
   for(col=0;col<Dest->data.s.nx;col++)
      if (Dest->mem2d->GetDataPkt(col, line)!=0)     
         f << col << " " << line << " " << Dest->mem2d->GetDataPkt(col, line) << endl;     
  f.close();
  return MATH_OK;
}

---