paint_utils.c

#if defined(__GNUC__)
#ident "University of Edinburgh $Id$"
#else
static char _paint_utils_c[] = "University of Edinburgh $Id$";
#endif
/*!
* \file         paint_utils.c
* \author	Richard Baldock
* \date		April 2009
* \version      $Id$
* \par
* Address:
*               MRC Human Genetics Unit,
*               MRC Institute of Genetics and Molecular Medicine,
*               University of Edinburgh,
*               Western General Hospital,
*               Edinburgh, EH4 2XU, UK.
* \par
* Copyright (C), [2012],
* The University Court of the University of Edinburgh,
* Old College, Edinburgh, UK.
* 
* 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., 51 Franklin Street, Fifth Floor,
* Boston, MA  02110-1301, USA.
* \ingroup      MAPaint
* \brief        
*/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include <MAPaint.h>
#include <X11/Xproto.h>

void myHGU_XmHelpStandardCb(Widget	w,
			  XtPointer	client_data,
			  XtPointer	call_data)
{
  /* dummy routine to switch off help */
  HGU_XmUserMessage(w,
		    "Button and context help has been\n"
		    "disabled to switch to a more usable\n"
		    "system based on JavaHelp.\n\n"
		    "For now please see help pages on the\n"
		    "Mouse Atlas web site:\n"
		    "http://genex.hgu.mrc.ac.uk/Resources/"
		    "or email ma-tech@hgu.mrc.ac.uk.\n"
		    "Thank you",
		    XmDIALOG_FULL_APPLICATION_MODAL);

  return;
}



void popup_dialog_cb(
Widget	w,
XtPointer	client_data,
XtPointer	call_data)
{
    XtManageChild( (Widget) client_data );
    XtPopup( XtParent((Widget) client_data), XtGrabNone );

    /* kludge to update the file selections */
    if( XmIsFileSelectionBox((Widget) client_data) ){
      XtCallCallbacks((Widget) client_data, XmNmapCallback, call_data);
    }
    return;
}

void UnmanageChildCb(
  Widget	w,
  XtPointer	client_data,
  XtPointer	call_data)
{
  if( client_data != NULL ){
    XtUnmanageChild( (Widget) client_data );
  }
  else {
    XtUnmanageChild( w );
  }

  return;
}

void destroy_cb(
Widget          w,
XtPointer       client_data,
XtPointer       call_data)
{
    XtDestroyWidget( (Widget) client_data );
    return;
}
 
int InteractDispatchEvent(
Display			*dpy,
Window			win,
XtPointer		client_data,
XEvent			*event)
{
  XtDispatchEvent( event );

  return( globals.paint_action_quit_flag );
}

int InteractHelpCb(
Display			*dpy,
Window			win,
XtPointer		client_data,
String			help_str)
{
    HGU_XmHelpShowString( globals.topl, (String) client_data, help_str );
    return( 0 );
}

int set_grey_values_from_ref_object(
  WlzObject	*obj,
  WlzObject	*ref_obj)
{
  WlzObject		*objs[3], *new_obj;
  WlzPlaneDomain	*pdom, *pdom1, *pdom2;
  WlzVoxelValues	*voxvals1, *voxvals2;
  WlzIntervalWSpace	iwsp1, iwsp2;
  WlzGreyWSpace	gwsp1, gwsp2;
  int			p;
  WlzPixelV		min, max, Min, Max;
  WlzErrorNum		errNum=WLZ_ERR_NONE;

  /* check the object */
  if((obj == NULL) || (obj->type != ref_obj->type) ||
     ((obj->type != WLZ_2D_DOMAINOBJ) && (obj->type != WLZ_3D_DOMAINOBJ))){
    return( -1 );
  }

  /* set the required grey ranges */
  min.type = WLZ_GREY_UBYTE;
  max.type = WLZ_GREY_UBYTE;
  Min.type = WLZ_GREY_UBYTE;
  Max.type = WLZ_GREY_UBYTE;
  min.v.ubv = 0;
  max.v.ubv = 255;
  Min.v.ubv = globals.cmapstruct->gmin;
  Max.v.ubv = globals.cmapstruct->gmax;

  /* scan through the new domain */
  objs[0] = WlzAssignObject(obj, NULL);
  objs[1] = WlzAssignObject(globals.obj, NULL);
  objs[2] = WlzAssignObject(ref_obj, NULL);
  new_obj = WlzIntersectN(3, objs, 0, &errNum);
  WlzFreeObj(objs[0]);
  WlzFreeObj(objs[1]);
  WlzFreeObj(objs[2]);

  if( errNum == WLZ_ERR_NONE ){
    switch( obj->type ){
    default:
      errNum = WLZ_ERR_OBJECT_TYPE;
      break;

    case WLZ_2D_DOMAINOBJ:
      if((objs[0] = WlzMakeMain(WLZ_2D_DOMAINOBJ,
			       new_obj->domain,
			       globals.obj->values,
				NULL, NULL, &errNum))){

	if((objs[1] = WlzMakeMain(WLZ_2D_DOMAINOBJ,
				 new_obj->domain,
				 ref_obj->values,
				  NULL, NULL, &errNum))){
	  errNum = WlzInitGreyScan( objs[0], &iwsp1, &gwsp1 );
	  errNum |= WlzInitGreyScan( objs[1], &iwsp2, &gwsp2 );
	  while((errNum == WLZ_ERR_NONE) &&
		((errNum = WlzNextGreyInterval(&iwsp1)) == WLZ_ERR_NONE) &&
		(WlzNextGreyInterval(&iwsp2) == WLZ_ERR_NONE) )
	  {
	    memcpy((void *) gwsp1.u_grintptr.ubp,
		   (const void *) gwsp2.u_grintptr.ubp,
		   (size_t) (iwsp1.rgtpos - iwsp1.lftpos + 1));
	  }
	  if( errNum == WLZ_ERR_EOO ){
	    errNum = WLZ_ERR_NONE;
	  }
	  WlzGreySetRange(objs[0], min, max, Min, Max, 0);
	  WlzFreeObj( objs[1] );
	}
	WlzFreeObj( objs[0] );
      }
      break;

    case WLZ_3D_DOMAINOBJ:
      pdom = new_obj->domain.p;
      pdom1 = globals.obj->domain.p;
      pdom2 = ref_obj->domain.p;
      voxvals1 = globals.obj->values.vox;
      voxvals2 = ref_obj->values.vox;

      for(p=pdom->plane1; p <= pdom->lastpl; p++){

	if( pdom->domains[p - pdom->plane1].core == NULL ){
	  continue;
	}

	if((objs[0] = WlzMakeMain(WLZ_2D_DOMAINOBJ,
				  pdom->domains[p - pdom->plane1],
				  voxvals1->values[p - pdom1->plane1],
				  NULL, NULL, &errNum))){

	  if((objs[1] = WlzMakeMain(WLZ_2D_DOMAINOBJ,
				    pdom->domains[p - pdom->plane1],
				    voxvals2->values[p - pdom2->plane1],
				    NULL, NULL, &errNum))){

	    errNum = WlzInitGreyScan( objs[0], &iwsp1, &gwsp1 );
	    errNum |= WlzInitGreyScan( objs[1], &iwsp2, &gwsp2 );
	    while((errNum == WLZ_ERR_NONE) &&
		  ((errNum = WlzNextGreyInterval(&iwsp1)) == WLZ_ERR_NONE) &&
		  (WlzNextGreyInterval(&iwsp2) == WLZ_ERR_NONE) )
	    {
	      memcpy((void *) gwsp1.u_grintptr.ubp,
		     (const void *) gwsp2.u_grintptr.ubp,
		     (size_t) (iwsp1.rgtpos - iwsp1.lftpos + 1));
	    }
	    if( errNum == WLZ_ERR_EOO ){
	      errNum = WLZ_ERR_NONE;
	    }
	    WlzGreySetRange(objs[0], min, max, Min, Max, 0);
	    WlzFreeObj( objs[1] );
	  }
	  WlzFreeObj( objs[0] );
	}
      }
      break;
    }
  }

  if( new_obj ){
    WlzFreeObj( new_obj );
  }

  if( errNum != WLZ_ERR_NONE ){
    MAPaintReportWlzError(globals.topl, "set_grey_values_from_ref_object", errNum);
  }
  return( 0 );
}

static int set_grey_values_from_domain3D(
  WlzObject		*obj,
  WlzObject		*destObj,
  unsigned int		col,
  unsigned int		planes)
{
  WlzObject		*obj1, *new_obj;
  WlzPlaneDomain	*pdom, *pdom1;
  WlzVoxelValues	*voxvaltb;
  WlzIntervalWSpace	iwsp;
  WlzGreyWSpace		gwsp;
  int			k, p;
  WlzErrorNum		errNum=WLZ_ERR_NONE;

  /* check the object */
  if( obj == NULL || obj->type != WLZ_3D_DOMAINOBJ )
    return( -1 );

  /* scan through the new domain */
  if(((new_obj = WlzIntersect2(obj, destObj, &errNum)) == NULL) ||
     (new_obj->type == WLZ_EMPTY_OBJ) ){
    HGU_XmUserMessage(globals.topl,
		      "The domain object just read has no\n"
		      "intersection with the reference object",
		      XmDIALOG_FULL_APPLICATION_MODAL);
    if( new_obj ){
      WlzFreeObj( new_obj );
    }
    if( errNum != WLZ_ERR_NONE ){
      MAPaintReportWlzError(globals.topl, "set_grey_values_from_domain3D",
			    errNum);
    }
    return( 0 );
  }
  pdom = new_obj->domain.p;
  pdom1 = destObj->domain.p;
  voxvaltb = destObj->values.vox;

  for(p=pdom->plane1; p <= pdom->lastpl; p++){

    if( pdom->domains[p - pdom->plane1].core == NULL )
      continue;

    if((obj1 = WlzMakeMain(WLZ_2D_DOMAINOBJ,
			   pdom->domains[p - pdom->plane1],
			   voxvaltb->values[p - pdom1->plane1],
			   NULL, NULL, &errNum))){

      errNum = WlzInitGreyScan( obj1, &iwsp, &gwsp );
      while((errNum == WLZ_ERR_NONE) &&
	    ((errNum = WlzNextGreyInterval(&iwsp)) == WLZ_ERR_NONE) ){
	for( k=0; k <= (iwsp.rgtpos - iwsp.lftpos); k++ ){
	  gwsp.u_grintptr.ubp[k] &= ~planes;
	  gwsp.u_grintptr.ubp[k] |= col;
	}
      }
      if( errNum == WLZ_ERR_EOO ){
	errNum = WLZ_ERR_NONE;
      }
      WlzFreeObj( obj1 );
    }
  }
  WlzFreeObj( new_obj );

  if( errNum != WLZ_ERR_NONE ){
    MAPaintReportWlzError(globals.topl, "set_grey_values_from_domain3D", errNum);
  }
  return( 0 );
}

int set_grey_values_from_domain(
  WlzObject		*obj,
  WlzObject		*destObj,
  unsigned int		col,
  unsigned int		planes)
{
  WlzObject		*new_obj;
  WlzIntervalWSpace	iwsp;
  WlzGreyWSpace		gwsp;
  int			k;
  WlzErrorNum		errNum=WLZ_ERR_NONE;

  /* check the object */
  if( obj == NULL ){
    return( -1 );
  }

  /* check the object type */
  if( obj->type != destObj->type ){
    return( -1 );
  }

  if( obj->type == WLZ_3D_DOMAINOBJ ){
    return set_grey_values_from_domain3D(obj, destObj, col, planes);
  }

  /* scan through the new domain */
  if((new_obj = WlzIntersect2(obj, destObj, &errNum))){

    new_obj = WlzAssignObject(new_obj, NULL);
    if( WlzArea(new_obj, &errNum) > 0 ){
      new_obj->values = WlzAssignValues(destObj->values, NULL);
      errNum = WlzInitGreyScan( new_obj, &iwsp, &gwsp );
      while((errNum == WLZ_ERR_NONE) &&
	    ((errNum = WlzNextGreyInterval(&iwsp)) == WLZ_ERR_NONE) ){
	for( k=0; k <= (iwsp.rgtpos - iwsp.lftpos); k++ ){
	  gwsp.u_grintptr.ubp[k] &= ~planes;
	  gwsp.u_grintptr.ubp[k] |= col;
	}
      }
      if( errNum == WLZ_ERR_EOO ){
	errNum = WLZ_ERR_NONE;
      }
    }
    WlzFreeObj(new_obj);
  }

  if( errNum != WLZ_ERR_NONE ){
    MAPaintReportWlzError(globals.topl, "set_grey_values_from_domain", errNum);
    return 1;
  }
  return 0;
}

void MAPaintReportWlzError(
  Widget	w,
  String	srcStr,
  WlzErrorNum	wlzErr)
{
  char		msgBuf[512];
  const char	*errStr, *msgStr;

  /* get the error and message strings */
  errStr = WlzStringFromErrorNum(wlzErr, &msgStr);

  /* create the message string to display */
  sprintf(msgBuf,
	  "A Woolz error has occurred in %s:\n\n"
	  "\t%s: %s\n\n"
	  "Safest action is to try and save domains\n"
	  "and restart paint. If that fails you can\n"
	  "recover data from the autosave file.\n\n"
	  "Please report this error to:\n"
	  "ma-tech@hgu.mrc.ac.uk",
	  srcStr, errStr, msgStr);
  HGU_XmUserError(w, msgBuf, XmDIALOG_FULL_APPLICATION_MODAL);
  
  return;
}