summaryrefslogtreecommitdiff
path: root/src/Stitch.cxx
blob: 535995f63db58fb4b0f57435b84d2b3f1a6fbb96 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
// 
// Stitch routines.
//
// Copyright 2006 by Johannes Hofmann
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Library General Public
// License as published by the Free Software Foundation; either
// version 2 of the License, or (at your option) any later version.
//
// This library 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
// Library General Public License for more details.
//
// You should have received a copy of the GNU Library General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA.
//

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

#include <Fl/Fl.H>

#include "Stitch.H"

static double pi_d = asin(1.0) * 2.0;

Stitch::Stitch() {
	for (int i=0; i<MAX_PICS; i++) {
		gipf[i] = NULL;
	}
}

Stitch::~Stitch() {
	for (int i=0; i<MAX_PICS; i++) {
		if (gipf[i]) {
			delete(gipf[i]);
		} else {
			break;
		}
	}
}


int
Stitch::load_image(char *file) {
	for (int i=0; i<MAX_PICS; i++) {
		if (gipf[i] == NULL) {
			gipf[i] = new GipfelWidget(0, 0, 800, 600);
			if (gipf[i]->load_image(file) != 0) {
				delete gipf[i];
				gipf[i] = NULL;
			}
			break;
		}
	}


}

int
Stitch::resample(DataImage *img,
            double view_start, double view_end) {
	double step_view = (view_end - view_start) / img->w();
	char r, g, b;
	int y_off = img->h() / 2;
 	double radius = (double) img->w() / (view_end -view_start);

	for (int y=0; y<img->h(); y++) {
		double a_nick = atan((double)(y_off - y)/radius);

		for (int x=0; x<img->w(); x++) {
			double a_view;
			a_view = view_start + x * step_view;

			for (int i=0; i<MAX_PICS; i++) {
				if (gipf[i] == NULL) {
					break;
				} else if (gipf[i]->get_pixel(a_view, a_nick, &r, &g, &b)==0) {
					img->set_pixel(x, y, r, g, b);
					break;
				}
			}
		}

		if (y % (img->h() / 100 + 1) == 0) {
			img->redraw();
			Fl::check();
		}
	}
}