xgray_tcl.cpp

Go to the documentation of this file.

/* Additional functions  for the Tcl interface */
#include "xgray_tcl.hpp"
#include <iostream> 

void getimage (const bild16& bmp, Tk_PhotoHandle vhandle) {
  Tk_PhotoImageBlock pb;
  pb.pixelPtr=new unsigned char [512*512];
  pb.width=512;
  pb.height=512;
  pb.pixelSize=1;
  pb.pitch=512;
  pb.offset[0]=0; // red, green and blue all point to the same value
  pb.offset[1]=0; 
  pb.offset[2]=0;
  pb.offset[3]=0; // Alpha should not been considered
  MAP (bmp) {
    int addr=(y*512+x);
    pb.pixelPtr[addr]=bmp(y,x)/256; // simple linear transform
 }   

 Tk_PhotoPutBlock(vhandle, &pb, 0, 0, 512, 512, TK_PHOTO_COMPOSITE_SET);
 delete[] pb.pixelPtr;
}

void getimage_contrast (const bild16& bmp, Tk_PhotoHandle vhandle, int low, int high) THROWSPEC
{
  if (low==high) STHROW("Error: Zero range (getimage_contrast), low==high=="<<low);
  
  Tk_PhotoImageBlock pb;
  pb.pixelPtr=new unsigned char [512*512];
  pb.width=512;
  pb.height=512;
  pb.pixelSize=1;
  pb.pitch=512;
  pb.offset[0]=0; // red, green and blue all point to the same value
  pb.offset[1]=0; 
  pb.offset[2]=0;
  pb.offset[3]=0; // Alpha should not been considered
  
  MAP (bmp) {
    int addr=(y*512+x);
    int val=int(255*double(bmp(y,x)-low)/double(high-low));
    if (val<0) val=0;
    if (val>255) val=255;
    pb.pixelPtr[addr]=val; // simple linear transform
 }   

 Tk_PhotoPutBlock(vhandle, &pb, 0, 0, 512, 512, TK_PHOTO_COMPOSITE_SET);
 delete[] pb.pixelPtr;
}

int /*low*/ getimage_range(const bild16& bmp, Tk_PhotoHandle vhandle, int x1, int y1, int x2, int y2, int *high) THROWSPEC
{
  int low=getminmax(bmp, x1,y1,x2,y2,high);

  if (low==*high) {
    getimage_contrast(bmp, vhandle, 0,65535);
  } else {
    getimage_contrast(bmp, vhandle, low, *high);
  }
  return low;
}  

int /*low*/ getminmax(const bild16& bmp, int x1, int y1, int x2, int y2, int *high) THROWSPEC
{
  if ((x1<0) || (y1<0) || (x2<0) || (y2<0) ||
      (x1>=bmp.SIZEX) || (x2>=bmp.SIZEX) ||
      (y1>=bmp.SIZEX) || (y2>=bmp.SIZEX)) {
         STHROW("Index error: getminmax x1="<<x1<<" x2="<<x2<<" y1="<<y1<<" y2="<<y2);
  }      
  if (x1>x2) std::swap(x1,x2);
  if (y1>y2) std::swap(y1,y2);

  int low=bmp(y1,x1);
  *high=low;

  for (int y=y1; y<=y2; y++) 
    for (int x=x1; x<=x2; x++) {
      int val=bmp(y,x);
      if (val<low) low=val;
      if (val>*high) *high=val;
  }
  return low;
}  



/* Variante für bildfloats. Geht nur mit Angabe des Bereichs */
        
void getimage_contrast (const bildfloat& bmp, Tk_PhotoHandle vhandle, float low, float high) THROWSPEC
{
  if (low==high) STHROW("Error: Zero range (getimage_contrast), low==high=="<<low);
  Tk_PhotoImageBlock pb;
  pb.pixelPtr=new unsigned char [512*512];
  pb.width=512;
  pb.height=512;
  pb.pixelSize=1;
  pb.pitch=512;
  pb.offset[0]=0; // red, green and blue all point to the same value
  pb.offset[1]=0; 
  pb.offset[2]=0;
  pb.offset[3]=0; // Alpha should not been considered
  
  MAP (bmp) {
    int addr=(y*512+x);
    int val=int(255*double(bmp(y,x)-low)/double(high-low));
    if (val<0) val=0;
    if (val>255) val=255;
    pb.pixelPtr[addr]=val; // simple linear transform
 }   

 Tk_PhotoPutBlock(vhandle, &pb, 0, 0, 512, 512, TK_PHOTO_COMPOSITE_SET);
 delete[] pb.pixelPtr;
}

void getimage (const bildfloat& bmp, Tk_PhotoHandle vhandle) {
   getimage_contrast(bmp, vhandle, 0.0, 1.0);
}   


float /*low*/ getimage_range (const bildfloat& bmp, Tk_PhotoHandle vhandle, int x1, int y1, int x2, int y2, float *high) THROWSPEC
{
 float low=getminmax(bmp,x1,y1,x2,y2,high);
 
 if (low==*high) {
    getimage_contrast(bmp, vhandle, 0.0,1.0);
  } else {
    getimage_contrast(bmp, vhandle, low, *high);
  }
  return low;
}  

float /*low*/ getminmax (const bildfloat& bmp, int x1, int y1, int x2, int y2, float *high) THROWSPEC
{
  if ((x1<0) || (y1<0) || (x2<0) || (y2<0) ||
      (x1>=bmp.SIZEX) || (x2>=bmp.SIZEX) ||
      (y1>=bmp.SIZEX) || (y2>=bmp.SIZEX)) {
         STHROW("Index error: getminmax x1="<<x1<<" x2="<<x2<<" y1="<<y1<<" y2="<<y2);
  }      
  if (x1>x2) std::swap(x1,x2);
  if (y1>y2) std::swap(y1,y2);

  float low=bmp(y1,x1);
  *high=low;

  for (int y=y1; y<=y2; y++) 
    for (int x=x1; x<=x2; x++) {
      float val=bmp(y,x);
      if (val<low) low=val;
      if (val>*high) *high=val;
  }
  return low;
}  

       
bild16 *cimg=NULL;
bildfloat *cimg_f=NULL;
bildcomplex *cimg_c=NULL;

---