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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
|
package ij.process;
import ij.Prefs;
import java.awt.*;
/** This class does bit blitting of 32-bit floating-point images. */
public class FloatBlitter implements Blitter {
public static float divideByZeroValue;
private FloatProcessor ip;
private int width, height;
private float[] pixels;
static {
divideByZeroValue = (float)Prefs.getDouble(Prefs.DIV_BY_ZERO_VALUE, Float.POSITIVE_INFINITY);
if (divideByZeroValue==Float.MAX_VALUE)
divideByZeroValue = Float.POSITIVE_INFINITY;
}
/** Constructs a FloatBlitter from a FloatProcessor. */
public FloatBlitter(FloatProcessor ip) {
this.ip = ip;
width = ip.getWidth();
height = ip.getHeight();
pixels = (float[])ip.getPixels();
}
public void setTransparentColor(Color c) {
}
/** Copies the float image in 'ip' to (x,y) using the specified mode. */
public void copyBits(ImageProcessor ip, int xloc, int yloc, int mode) {
Rectangle r1, r2;
int srcIndex, dstIndex;
int xSrcBase, ySrcBase;
float[] srcPixels;
if (!(ip instanceof FloatProcessor))
ip = ip.convertToFloat();
int srcWidth = ip.getWidth();
int srcHeight = ip.getHeight();
r1 = new Rectangle(srcWidth, srcHeight);
r1.setLocation(xloc, yloc);
r2 = new Rectangle(width, height);
if (!r1.intersects(r2))
return;
srcPixels = (float [])ip.getPixels();
r1 = r1.intersection(r2);
xSrcBase = (xloc<0)?-xloc:0;
ySrcBase = (yloc<0)?-yloc:0;
boolean useDBZValue = !Float.isInfinite(divideByZeroValue);
float src, dst;
for (int y=r1.y; y<(r1.y+r1.height); y++) {
srcIndex = (y-yloc)*srcWidth + (r1.x-xloc);
dstIndex = y * width + r1.x;
switch (mode) {
case COPY: case COPY_INVERTED: case COPY_TRANSPARENT:
for (int i=r1.width; --i>=0;)
pixels[dstIndex++] = srcPixels[srcIndex++];
break;
case COPY_ZERO_TRANSPARENT:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++];
if (src==0f)
dst = pixels[dstIndex];
else
dst = src;
pixels[dstIndex++] = dst;
}
break;
case ADD:
for (int i=r1.width; --i>=0; srcIndex++, dstIndex++)
pixels[dstIndex] = srcPixels[srcIndex]+pixels[dstIndex];
break;
case AVERAGE:
for (int i=r1.width; --i>=0;) {
dst =(srcPixels[srcIndex++]+pixels[dstIndex])/2;
pixels[dstIndex++] = dst;
}
break;
case DIFFERENCE:
for (int i=r1.width; --i>=0; srcIndex++, dstIndex++) {
dst = pixels[dstIndex]-srcPixels[srcIndex];
pixels[dstIndex] = dst<0?-dst:dst;
}
break;
case SUBTRACT:
for (int i=r1.width; --i>=0; srcIndex++, dstIndex++)
pixels[dstIndex] = pixels[dstIndex]-srcPixels[srcIndex];
break;
case MULTIPLY:
for (int i=r1.width; --i>=0; srcIndex++, dstIndex++)
pixels[dstIndex] = srcPixels[srcIndex]*pixels[dstIndex];
break;
case DIVIDE:
for (int i=r1.width; --i>=0; srcIndex++, dstIndex++) {
src = srcPixels[srcIndex];
if (useDBZValue && src==0.0)
pixels[dstIndex] = divideByZeroValue;
else
pixels[dstIndex] = pixels[dstIndex]/src;
}
break;
case AND:
for (int i=r1.width; --i>=0;) {
dst = (int)srcPixels[srcIndex++]&(int)pixels[dstIndex];
pixels[dstIndex++] = dst;
}
break;
case OR:
for (int i=r1.width; --i>=0;) {
dst = (int)srcPixels[srcIndex++]|(int)pixels[dstIndex];
pixels[dstIndex++] = dst;
}
break;
case XOR:
for (int i=r1.width; --i>=0;) {
dst = (int)srcPixels[srcIndex++]^(int)pixels[dstIndex];
pixels[dstIndex++] = dst;
}
break;
case MIN:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++];
dst = pixels[dstIndex];
if (src<dst) dst = src;
pixels[dstIndex++] = dst;
}
break;
case MAX:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++];
dst = pixels[dstIndex];
if (src>dst) dst = src;
pixels[dstIndex++] = dst;
}
break;
}
if (y%20==0)
ip.showProgress((double)(y-r1.y)/r1.height);
}
ip.showProgress(1.0);
}
}
|
