imagemap.cs

Image Fusion Techniues

// Waveblend - complex dualtree based image fusion
// (C) Copyright 2004 -- Sebastian Nowozin <[e-mail]>
//
// This file is part of Waveblend.
//
// Waveblend 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; version 2 of the License.
//
// Waveblend 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.
//
// The license is included with the distribution in the file 'LICENSE'.
//

/* ImageMap.cs
 *
 * Abstract floating point image map processing functionality.
 */

using System;

public class ImageMap : ICloneable
{
	private ImageMap ()
	{
	}

	public object Clone ()
	{
		ImageMap cl = new ImageMap (xDim, yDim);

		for (int y = 0 ; y < yDim ; ++y)
			for (int x = 0 ; x < xDim ; ++x)
				cl[x, y] = this[x, y];

		return (cl);
	}

	int yDim, xDim;
	double[,] valArr;
	public double[,] ValueArray {
		get {
			return (valArr);
		}
	}

	public int YDim {
		get {
			return (yDim);
		}
	}
	public int XDim {
		get {
			return (xDim);
		}
	}

	public double MaxElem {
		get {
			double max = 0.0;

			for (int y = 0 ; y < yDim ; ++y) {
				for (int x = 0 ; x < xDim ; ++x) {
					if (this[x, y] > max)
						max = this[x, y];
				}
			}
			return (max);
		}
	}

	public double this[int x, int y] {
		get {
			return (valArr[y, x]);
		}
		set {
			valArr[y, x] = value;
		}
	}

	public ImageMap (int xDim, int yDim)
	{
		this.xDim = xDim;
		this.yDim = yDim;

		valArr = new double[yDim, xDim];
	}

	public ImageMap EnlargeCanvas (int newX, int newY)
	{
		if (newX < xDim || newY < yDim)
			throw (new ArgumentException ("New canvas dimensions too small."));

		ImageMap nmap = new ImageMap (newX, newY);
		for (int y = 0 ; y < yDim ; ++y)
			for (int x = 0 ; x < xDim ; ++x)
				nmap[x, y] = this[x, y];

		return (nmap);
	}

	public ImageMap ShrinkCanvas (int newX, int newY)
	{
		if (newX > xDim || newY > yDim)
			throw (new ArgumentException ("New canvas dimensions too large."));

		ImageMap nmap = new ImageMap (newX, newY);
		for (int y = 0 ; y < newY ; ++y)
			for (int x = 0 ; x < newX ; ++x)
				nmap[x, y] = this[x, y];

		return (nmap);
	}

	public ImageMap (int[,] values)
	{
		valArr = new double[values.GetLength (0), values.GetLength (1)];
		xDim = values.GetLength (1);
		yDim = values.GetLength (0);

		for (int y = 0 ; y < yDim ; ++y)
			for (int x = 0 ; x < xDim ; ++x)
				valArr[y, x] = values[y, x];
	}

	public ImageMap (double[,] values)
	{
		valArr = new double[values.GetLength (0), values.GetLength (1)];
		xDim = values.GetLength (1);
		yDim = values.GetLength (0);

		for (int y = 0 ; y < yDim ; ++y)
			for (int x = 0 ; x < xDim ; ++x)
				valArr[y, x] = values[y, x];
	}

	public ImageMap ScaleHalf ()
	{
		if ((xDim / 2) == 0 || (yDim / 2) == 0)
			return (null);

		ImageMap res = new ImageMap (xDim / 2, yDim / 2);

		for (int y = 0 ; y < res.yDim ; ++y) {
			for (int x = 0 ; x < res.xDim ; ++x) {
				res[x, y] = this[2 * x, 2 * y];
			}
		}

		return (res);
	}

	// Double the size of an imagemap using linear interpolation.
	// It is not a real doubling as the last line is omitted (the image size
	// would always be odd otherwise and we have no second line for
	// interpolation).
	public ImageMap ScaleDouble ()
	{
		// Doubling an image with x/y dimensions less or equal than 2 will
		// result in an image with just (2, 2) dims, so its useless.
		if (xDim <= 2 || yDim <= 2)
			return (null);

		ImageMap res = new ImageMap (xDim * 2 - 2, yDim * 2 - 2);

		// fill four pixels per step, except for the last line/col, which will
		// be omitted
		for (int y = 0 ; y < (yDim - 1) ; ++y) {
			for (int x = 0 ; x < (xDim - 1) ; ++x) {
				// pixel layout:
				// A B
				// C D

				// A
				res[2 * x + 0, 2 * y + 0] = this[x, y];
				// B
				res[2 * x + 1, 2 * y + 0] =
					(this[x, y] + this[x + 1, y]) / 2.0;
				// C
				res[2 * x + 0, 2 * y + 1] =
					(this[x, y] + this[x, y + 1]) / 2.0;
				// D
				res[2 * x + 1, 2 * y + 1] = (this[x, y] + this[x + 1, y] +
					this[x, y + 1] + this[x + 1, y + 1]) / 4.0;
			}
		}

		return (res);
	}

	static public ImageMap operator* (ImageMap f1, ImageMap f2)
	{
		if (f1.xDim != f2.xDim || f1.yDim != f2.yDim) {
			throw (new ArgumentException ("Mismatching dimensions"));
		}

		ImageMap resultMap = new ImageMap (f1.xDim, f1.yDim);

		for (int y = 0 ; y < f1.yDim ; ++y) {
			for (int x = 0 ; x < f1.xDim ; ++x) {
				resultMap[x, y] = f1[x, y] * f2[x, y];
			}
		}

		return (resultMap);
	}

	static public ImageMap operator+ (ImageMap f1, ImageMap f2)
	{
		if (f1.xDim != f2.xDim || f1.yDim != f2.yDim) {
			throw (new ArgumentException ("Mismatching dimensions"));
		}

		ImageMap resultMap = new ImageMap (f1.xDim, f1.yDim);

		for (int y = 0 ; y < f1.yDim ; ++y) {
			for (int x = 0 ; x < f1.xDim ; ++x) {
				resultMap[x, y] = f1[x, y] + f2[x, y];
			}
		}

		return (resultMap);
	}

	static public ImageMap operator- (ImageMap f1, ImageMap f2)
	{
		if (f1.xDim != f2.xDim || f1.yDim != f2.yDim) {
			throw (new ArgumentException ("Mismatching dimensions"));
		}

		ImageMap resultMap = new ImageMap (f1.xDim, f1.yDim);

		for (int y = 0 ; y < f1.yDim ; ++y) {
			for (int x = 0 ; x < f1.xDim ; ++x) {
				resultMap[x, y] = f1[x, y] - f2[x, y];
			}
		}

		return (resultMap);
	}

	// Normalize: Find the minimum to maximum range, then stretch and limit
	// those to exactly 0.0 to 1.0. If both the minimum and maximum values are
	// equal, no normalization takes place.
	public void Normalize ()
	{
		double min = 0.0;
		double max = 0.0;

		for (int y = 0 ; y < yDim ; ++y) {
			for (int x = 0 ; x < xDim ; ++x) {
				if (min > this[x, y])
					min = this[x, y];

				if (max < this[x, y])
					max = this[x, y];
			}
		}
		if (min == max)
			return;

		double diff = max - min;

		for (int y = 0 ; y < yDim ; ++y) {
			for (int x = 0 ; x < xDim ; ++x) {
				this[x, y] = (this[x, y] - min) / diff;
			}
		}
	}
}