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
145
146
147
148
|
package ij.process;
import java.awt.*;
/** This class does bit blitting of 16-bit images. */
public class ShortBlitter implements Blitter {
private ShortProcessor ip;
private int width, height;
private short[] pixels;
public void setTransparentColor(Color c) {
}
/** Constructs a ShortBlitter from a ShortProcessor. */
public ShortBlitter(ShortProcessor ip) {
this.ip = ip;
width = ip.getWidth();
height = ip.getHeight();
pixels = (short[])ip.getPixels();
}
/** Copies the 16-bit 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;
short[] srcPixels;
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 = (short [])ip.getPixels();
//new ij.ImagePlus("srcPixels", new ShortProcessor(srcWidth, srcHeight, srcPixels, null)).show();
r1 = r1.intersection(r2);
xSrcBase = (xloc<0)?-xloc:0;
ySrcBase = (yloc<0)?-yloc:0;
int 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++]&0xffff;
if (src==0)
dst = pixels[dstIndex];
else
dst = src;
pixels[dstIndex++] = (short)dst;
}
break;
case ADD:
for (int i=r1.width; --i>=0;) {
dst = (srcPixels[srcIndex++]&0xffff)+(pixels[dstIndex]&0xffff);
if (dst<0) dst = 0;
if (dst>65535) dst = 65535;
pixels[dstIndex++] = (short)dst;
}
break;
case AVERAGE:
for (int i=r1.width; --i>=0;) {
dst = ((srcPixels[srcIndex++]&0xffff)+(pixels[dstIndex]&0xffff))/2;
pixels[dstIndex++] = (short)dst;
}
break;
case DIFFERENCE:
for (int i=r1.width; --i>=0;) {
dst = (pixels[dstIndex]&0xffff)-(srcPixels[srcIndex++]&0xffff);
if (dst<0) dst = -dst;
if (dst>65535) dst = 65535;
pixels[dstIndex++] = (short)dst;
}
break;
case SUBTRACT:
for (int i=r1.width; --i>=0;) {
dst = (pixels[dstIndex]&0xffff)-(srcPixels[srcIndex++]&0xffff);
if (dst<0) dst = 0;
if (dst>65535) dst = 65535;
pixels[dstIndex++] = (short)dst;
}
break;
case MULTIPLY:
for (int i=r1.width; --i>=0;) {
dst = (srcPixels[srcIndex++]&0xffff)*(pixels[dstIndex]&0xffff);
if (dst<0) dst = 0;
if (dst>65535) dst = 65535;
pixels[dstIndex++] = (short)dst;
}
break;
case DIVIDE:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++]&0xffff;
if (src==0)
dst = 65535;
else
dst = (pixels[dstIndex]&0xffff)/src;
pixels[dstIndex++] = (short)dst;
}
break;
case AND:
for (int i=r1.width; --i>=0;) {
dst = srcPixels[srcIndex++]&pixels[dstIndex]&0xffff;
pixels[dstIndex++] = (short)dst;
}
break;
case OR:
for (int i=r1.width; --i>=0;) {
dst = srcPixels[srcIndex++]|pixels[dstIndex];
pixels[dstIndex++] = (short)dst;
}
break;
case XOR:
for (int i=r1.width; --i>=0;) {
dst = srcPixels[srcIndex++]^pixels[dstIndex];
pixels[dstIndex++] = (short)dst;
}
break;
case MIN:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++]&0xffff;
dst = pixels[dstIndex]&0xffff;
if (src<dst) dst = src;
pixels[dstIndex++] = (short)dst;
}
break;
case MAX:
for (int i=r1.width; --i>=0;) {
src = srcPixels[srcIndex++]&0xffff;
dst = pixels[dstIndex]&0xffff;
if (src>dst) dst = src;
pixels[dstIndex++] = (short)dst;
}
break;
}
if (y%20==0)
ip.showProgress((double)(y-r1.y)/r1.height);
}
ip.showProgress(1.0);
}
}
|
