projects / ctsim.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
r328: *** empty log message ***
[ctsim.git] / include / imagefile.h
diff --git a/include/imagefile.h b/include/imagefile.h
index c5a68e726b3073a950c3006d68e8e0b8a4666968..3587f1cde01e5cc90b905ea5bafd9874200ca4d5 100644 (file)
--- a/include/imagefile.h
+++ b/include/imagefile.h
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (C) 1983-2000 Kevin Rosenberg
 **
-**  $Id: imagefile.h,v 1.15 2000/06/26 21:15:24 kevin Exp $
+**  $Id: imagefile.h,v 1.27 2001/01/02 05:34:57 kevin Exp $
 **
 **  This program is free software; you can redistribute it and/or modify
 **  it under the terms of the GNU General Public License (version 2) as
@@ -28,17 +28,17 @@
 #ifndef IMAGEFILE_H
 #define IMAGEFILE_H
 
+#ifndef MSVC\r
+#include <unistd.h>\r
+#endif\r
 #include <string>
 #include <sys/types.h>
-#include <unistd.h>
 #include <fstream>
 #include <iostream>
 #include "ctsupport.h"
 #include "fnetorderstream.h"
 #include "array2dfile.h"
 
-using namespace std;
-
 #ifdef HAVE_MPI
 #include <mpi++.h>
 #endif
@@ -46,24 +46,21 @@ using namespace std;
 class F32Image : public Array2dFile
 {
 public:
-  F32Image (int nx, int ny)
-      : Array2dFile::Array2dFile (nx, ny, sizeof(kfloat32), Array2dFile::PIXEL_FLOAT32)
-  {
-  }
-
-  F32Image (void)
-      : Array2dFile::Array2dFile()
-  {
-      setPixelFormat (Array2dFile::PIXEL_FLOAT32);
-      setPixelSize (sizeof(kfloat32));
-  }
-
-  kfloat32** getArray (void)
-      { return (kfloat32**) (m_arrayData); }
-
-  const kfloat32** getArray (void) const
-       { return (const kfloat32**) (m_arrayData); }
-
+  F32Image (int nx, int ny, int dataType = Array2dFile::DATA_TYPE_REAL);     
+  F32Image (void);\r
+
+  kfloat32** getArray (void)\r
+      { return (kfloat32**) (m_arrayData); }\r
+\r
+  kfloat32** const getArray (void) const\r
+       { return (kfloat32** const) (m_arrayData); }\r
+\r
+  kfloat32** getImaginaryArray (void)\r
+      { return (kfloat32**) (m_imaginaryArrayData); }\r
+\r
+  kfloat32** const getImaginaryArray (void) const\r
+       { return (kfloat32** const) (m_imaginaryArrayData); }\r
+\r
 #ifdef HAVE_MPI
   MPI::Datatype getMPIDataType (void) const
       { return MPI::FLOAT; }
@@ -79,24 +76,21 @@ class F64Image : public Array2dFile
 {
  public:
 
-  F64Image (int nx, int ny)
-      : Array2dFile::Array2dFile (nx, ny, sizeof(kfloat64), Array2dFile::PIXEL_FLOAT64)
-  {
-  }
-
-  F64Image (void)
-      : Array2dFile::Array2dFile ()
-  {
-      setPixelFormat (PIXEL_FLOAT64);
-      setPixelSize (sizeof(kfloat64));
-  }
-
-  kfloat64** getArray (void)
-      { return (kfloat64**) (m_arrayData); }
-
-  const kfloat64** getArray (void) const
-      { return (const kfloat64**) (m_arrayData); }
-
+   F64Image (int nx, int ny, int dataType = Array2dFile::DATA_TYPE_REAL);
+  F64Image (void);
+
+  kfloat64** getArray (void)\r
+      { return (kfloat64**) (m_arrayData); }\r
+\r
+  kfloat64** const getArray (void) const\r
+      { return (kfloat64** const) (m_arrayData); }\r
+\r
+  kfloat64** getImaginaryArray (void)\r
+      { return (kfloat64**) (m_imaginaryArrayData); }\r
+\r
+  kfloat64** const getImaginaryArray (void) const\r
+      { return (kfloat64** const) (m_imaginaryArrayData); }\r
+\r
 #ifdef HAVE_MPI
   MPI::Datatype getMPIDataType (void) const
       { return MPI::DOUBLE; }
@@ -112,17 +106,43 @@ typedef F64Image   ImageFileBase;
 typedef kfloat64   ImageFileValue;
 typedef kfloat64*  ImageFileColumn;
 typedef kfloat64** ImageFileArray;
+typedef kfloat64** const ImageFileArrayConst;\r
+typedef const kfloat64* ImageFileColumnConst;
 #else
 typedef F32Image   ImageFileBase;
 typedef kfloat32   ImageFileValue;
 typedef kfloat32*  ImageFileColumn;
 typedef kfloat32** ImageFileArray;
+typedef kfloat32** const ImageFileArrayConst;
+typedef const kfloat32* ImageFileColumnConst;\r
 #endif
 
 
 class ImageFile : public ImageFileBase
-{
- public:
+{\r
+private:\r
+\r
+  static const char* s_aszFormatName[];\r
+  static const char* s_aszFormatTitle[];\r
+  static const int s_iFormatCount;\r
+
+public:\r
+\r
+  static const int FORMAT_INVALID;\r
+  static const int FORMAT_PGM;\r
+  static const int FORMAT_PGMASCII;\r
+#if HAVE_PNG\r
+  static const int FORMAT_PNG;\r
+  static const int FORMAT_PNG16;\r
+#endif\r
+\r
+  static const int getFormatCount() {return s_iFormatCount;}\r
+  static const char** getFormatNameArray() {return s_aszFormatName;}\r
+  static const char** getFormatTitleArray() {return s_aszFormatTitle;}\r
+  static int convertFormatNameToID (const char* const formatName);\r
+  static const char* convertFormatIDToName (const int formatID);\r
+  static const char* convertFormatIDToTitle (const int formatID);\r
+
   ImageFile (int nx, int ny)
       : ImageFileBase (nx, ny)
   {}
@@ -130,21 +150,53 @@ class ImageFile : public ImageFileBase
   ImageFile (void)
       : ImageFileBase ()
   {}
+\r
+  bool convertRealToComplex ();\r
+  bool convertComplexToReal ();\r
 
-  void filterResponse (const char* const domainName, double bw, const char* const filterName, double filt_param);
-
-  void statistics (double& min, double& max, double& mean, double& mode, double& median, double& stddev) const;
-
-  void printStatistics (ostream& os) const;
+  void filterResponse (const char* const domainName, double bw, const char* const filterName, double filt_param, double dInputScale = 1., double dOutputScale = 1.);
 
+  void statistics (double& min, double& max, double& mean, double& mode, double& median, double& stddev) const;\r
+  void statistics (ImageFileArrayConst v, double& min, double& max, double& mean, double& mode, double& median, double& stddev) const;\r
+  void getMinMax (double& min, double& max) const;
+  void printStatistics (std::ostream& os) const;
   bool comparativeStatistics (const ImageFile& imComp, double& d, double& r, double& e) const;
-
-  bool printComparativeStatistics (const ImageFile& imComp, ostream& os) const;
-
-  int display (void);
-
-  int displayScaling (const int scaleFactor, ImageFileValue pmin, ImageFileValue pmax);
-
+  bool printComparativeStatistics (const ImageFile& imComp, std::ostream& os) const;
+\r
+  bool subtractImages (const ImageFile& rRHS, ImageFile& result) const;
+  bool addImages (const ImageFile& rRHS, ImageFile& result) const;\r
+  bool multiplyImages (const ImageFile& rRHS, ImageFile& result) const;\r
+  bool divideImages (const ImageFile& rRHS, ImageFile& result) const;\r
+\r
+  bool scaleImage (ImageFile& result) const;\r
+\r
+  bool invertPixelValues (ImageFile& result) const;\r
+  bool sqrt (ImageFile& result) const;\r
+  bool square (ImageFile& result) const;\r
+  bool log (ImageFile& result) const;\r
+  bool exp (ImageFile& result) const;\r
+  bool fourier (ImageFile& result) const;\r
+  bool inverseFourier (ImageFile& result) const;\r
+#ifdef HAVE_FFTW\r
+  bool fft (ImageFile& result) const;\r
+  bool ifft (ImageFile& result) const;\r
+#endif\r
+  bool magnitude (ImageFile& result) const;\r
+  bool phase (ImageFile& result) const;\r
+\r
+  int display (void) const;
+  int displayScaling (const int scaleFactor, ImageFileValue pmin, ImageFileValue pmax) const;
+\r
+  bool exportImage (const char* const pszFormat, const char* const pszFilename, int nxcell, int nycell, double densmin, double densmax);\r
+
+#if HAVE_PNG
+  bool writeImagePNG (const char* const outfile, int bitdepth, int nxcell, int nycell, double densmin, double densmax);
+#endif
+#if HAVE_GD
+  bool writeImageGIF (const char* const outfile, int nxcell, int nycell, double densmin, double densmax);
+#endif
+  bool writeImagePGM (const char* const outfile, int nxcell, int nycell, double densmin, double densmax);
+  bool writeImagePGMASCII (const char* const outfile, int nxcell, int nycell, double densmin, double densmax);
 };
 
 
Computed Tomography simulator and reconstructor