| 2 |
- |
1 |
% Copyright (C) 1995 Aladdin Enterprises. All rights reserved.
|
|
|
2 |
%
|
|
|
3 |
% This software is provided AS-IS with no warranty, either express or
|
|
|
4 |
% implied.
|
|
|
5 |
%
|
|
|
6 |
% This software is distributed under license and may not be copied,
|
|
|
7 |
% modified or distributed except as expressly authorized under the terms
|
|
|
8 |
% of the license contained in the file LICENSE in this distribution.
|
|
|
9 |
%
|
|
|
10 |
% For more information about licensing, please refer to
|
|
|
11 |
% http://www.ghostscript.com/licensing/. For information on
|
|
|
12 |
% commercial licensing, go to http://www.artifex.com/licensing/ or
|
|
|
13 |
% contact Artifex Software, Inc., 101 Lucas Valley Road #110,
|
|
|
14 |
% San Rafael, CA 94903, U.S.A., +1(415)492-9861.
|
|
|
15 |
|
|
|
16 |
% $Id: docie.ps,v 1.5 2003/04/12 18:08:18 ray Exp $
|
|
|
17 |
% docie.ps
|
|
|
18 |
% Emulate CIE algorithms in PostScript.
|
|
|
19 |
|
|
|
20 |
% ---------------- Auxiliary procedures ---------------- %
|
|
|
21 |
|
|
|
22 |
/r1default [0 1] def
|
|
|
23 |
/r3default [0 1 0 1 0 1] def
|
|
|
24 |
|
|
|
25 |
/apply3 % <u> <v> <w> [<pu> <pv> <pw>] apply3 <u'> <v'> <w'>
|
|
|
26 |
{ { 4 -1 roll exch exec } forall
|
|
|
27 |
} bind def
|
|
|
28 |
|
|
|
29 |
/restrict % <u> <min> <max> restrict <u'>
|
|
|
30 |
{ 3 1 roll .max .min
|
|
|
31 |
} bind def
|
|
|
32 |
|
|
|
33 |
/restrict3 % <u> <v> <w> [<minu> ... <maxw>] restrict3 <u'> <v'> <w'>
|
|
|
34 |
{ aload pop
|
|
|
35 |
7 -1 roll 3 1 roll restrict 7 1 roll
|
|
|
36 |
5 -1 roll 3 1 roll restrict 5 1 roll
|
|
|
37 |
restrict 3 1 roll
|
|
|
38 |
} bind def
|
|
|
39 |
|
|
|
40 |
/rescale % <u> <min> <max> rescale <u'>
|
|
|
41 |
{ 1 index sub 3 1 roll sub exch div 0 .max 1 .min
|
|
|
42 |
} bind def
|
|
|
43 |
|
|
|
44 |
/rescale3 % <u> <v> <w> [<minu> ... <maxw>] rescale3 <u'> <v'> <w'>
|
|
|
45 |
{ aload pop
|
|
|
46 |
7 -1 roll 3 1 roll rescale 7 1 roll
|
|
|
47 |
5 -1 roll 3 1 roll rescale 5 1 roll
|
|
|
48 |
rescale 3 1 roll
|
|
|
49 |
} bind def
|
|
|
50 |
|
|
|
51 |
/mmult3 % <u> <v> <w> [<uu> <uv> ... <wv> <ww>] mmult3
|
|
|
52 |
% <u'> <v'> <w'>
|
|
|
53 |
{ 4 -1 roll dup dup 6 -1 roll dup dup 8 -1 roll dup dup
|
|
|
54 |
10 -1 roll { 10 -1 roll mul } forall
|
|
|
55 |
% Stack: u1 v1 w1 u2 v2 w2 u3 v3 w3
|
|
|
56 |
4 -1 roll add 6 -1 roll add
|
|
|
57 |
% Stack: u1 v1 u2 v2 u3 v3 w'
|
|
|
58 |
7 1 roll 3 -1 roll add 4 -1 roll add
|
|
|
59 |
% Stack: w' u1 u2 u3 v'
|
|
|
60 |
5 1 roll add add 3 1 roll
|
|
|
61 |
} bind def
|
|
|
62 |
|
|
|
63 |
/minvert3 % [<uu> <uv> ... <wv> <ww>] minvert3
|
|
|
64 |
% [<uu'> <uv'> ... <wv'> <ww'>]
|
|
|
65 |
{ 16 dict begin
|
|
|
66 |
aload pop { I H G F E D C B A } { exch def } forall
|
|
|
67 |
/coa E I mul F H mul sub def
|
|
|
68 |
/cob F G mul D I mul sub def
|
|
|
69 |
/coc D H mul E G mul sub def
|
|
|
70 |
/det A coa mul B cob mul add C coc mul add def
|
|
|
71 |
[ coa det div
|
|
|
72 |
C H mul B I mul sub det div
|
|
|
73 |
B F mul C E mul sub det div
|
|
|
74 |
cob det div
|
|
|
75 |
A I mul C G mul sub det div
|
|
|
76 |
C D mul A F mul sub det div
|
|
|
77 |
coc det div
|
|
|
78 |
B G mul A H mul sub det div
|
|
|
79 |
A E mul B D mul sub det div
|
|
|
80 |
]
|
|
|
81 |
end
|
|
|
82 |
} bind def
|
|
|
83 |
|
|
|
84 |
/print1
|
|
|
85 |
{ print dup ==
|
|
|
86 |
} bind def
|
|
|
87 |
|
|
|
88 |
/print3
|
|
|
89 |
{ print 3 array astore dup == aload pop
|
|
|
90 |
} bind def
|
|
|
91 |
|
|
|
92 |
% ---------------- Mapping to XYZ ---------------- %
|
|
|
93 |
|
|
|
94 |
/csmap % <csdict> <l> <m> <n> csmap <csdict> <x> <y> <z>
|
|
|
95 |
{ 3 index /RangeLMN .knownget not { r3default } if restrict3
|
|
|
96 |
DOCIEDEBUG { (After RangeLMN Decode: ) print3 } if
|
|
|
97 |
3 index /DecodeLMN .knownget { apply3 } if
|
|
|
98 |
DOCIEDEBUG { (After DecodeLMN Decode: ) print3 } if
|
|
|
99 |
3 index /MatrixLMN .knownget { mmult3 } if
|
|
|
100 |
DOCIEDEBUG { (After MatrixLMN Decode: ) print3 } if
|
|
|
101 |
} bind def
|
|
|
102 |
|
|
|
103 |
/csciea % <csdict> <a> csciea <csdict> <x> <y> <z>
|
|
|
104 |
{ 1 index /RangeA .knownget not { r1default } if restrict
|
|
|
105 |
DOCIEDEBUG { (After RangeA Decode: ) print1 } if
|
|
|
106 |
1 index /DecodeA .knownget { exec } if
|
|
|
107 |
DOCIEDEBUG { (After DecodeA Decode: ) print1 } if
|
|
|
108 |
1 index /MatrixA .knownget
|
|
|
109 |
{ { 1 index mul exch } forall pop }
|
|
|
110 |
{ dup dup }
|
|
|
111 |
ifelse
|
|
|
112 |
DOCIEDEBUG { (After MatrixA Decode: ) print3 } if
|
|
|
113 |
csmap
|
|
|
114 |
} bind def
|
|
|
115 |
|
|
|
116 |
/cscieabc % <csdict> <a> <b> <c> cscieabc <csdict> <x> <y> <z>
|
|
|
117 |
{ 3 index /RangeABC .knownget not { r3default } if restrict3
|
|
|
118 |
DOCIEDEBUG { (After RangeABC Decode: ) print3 } if
|
|
|
119 |
3 index /DecodeABC .knownget { apply3 } if
|
|
|
120 |
DOCIEDEBUG { (After DecodeABC Decode: ) print3 } if
|
|
|
121 |
3 index /MatrixABC .knownget { mmult3 } if
|
|
|
122 |
DOCIEDEBUG { (After MatrixABC Decode: ) print3 } if
|
|
|
123 |
csmap
|
|
|
124 |
} bind def
|
|
|
125 |
|
|
|
126 |
% ---------------- Rendering from XYZ ---------------- %
|
|
|
127 |
|
|
|
128 |
/lookup3 % <rtable> <a[0..1]> <b[0..1]> <c[0..1]> lookup3
|
|
|
129 |
% <rtable> <bytes>
|
|
|
130 |
{ 3 -1 roll 3 index 0 get 1 sub mul
|
|
|
131 |
3 -1 roll 3 index 1 get 1 sub mul
|
|
|
132 |
3 -1 roll 3 index 2 get 1 sub mul
|
|
|
133 |
% Stack: rtable ia ib ic
|
|
|
134 |
DOCIEDEBUG { (RenderTable indices: ) print3 mark 5 1 roll } if
|
|
|
135 |
3 -1 roll round cvi 3 index 3 get exch get
|
|
|
136 |
% Stack: rtable ib ic string
|
|
|
137 |
3 -1 roll round cvi 3 index 2 get mul
|
|
|
138 |
% Stack: rtable ic string ib*nc
|
|
|
139 |
3 -1 roll round cvi add 2 index 4 get mul
|
|
|
140 |
% Stack: rtable string index
|
|
|
141 |
2 index 4 get getinterval
|
|
|
142 |
% Stack: rtable bytes
|
|
|
143 |
DOCIEDEBUG { (RenderTable values: ) print (<) print (%stdout) (w) file 1 index writehexstring (>) = } if
|
|
|
144 |
} bind def
|
|
|
145 |
|
|
|
146 |
/bpdefault [0 0 0] def
|
|
|
147 |
|
|
|
148 |
/crmap % <csdict> <crdict> <x> <y> <z> crmap <v1> ...
|
|
|
149 |
{
|
|
|
150 |
DOCIEDEBUG { (CIE XYZ = ) print3 } if
|
|
|
151 |
3 index /MatrixPQR .knownget { mmult3 } if
|
|
|
152 |
DOCIEDEBUG { (After MatrixPQR: ) print3 } if
|
|
|
153 |
4 index /WhitePoint get
|
|
|
154 |
5 index /BlackPoint .knownget not { bpdefault } if
|
|
|
155 |
5 index /WhitePoint get
|
|
|
156 |
6 index /BlackPoint .knownget not { bpdefault } if
|
|
|
157 |
4
|
|
|
158 |
{ 4 -1 roll aload pop
|
|
|
159 |
% Stack: csdict crdict x y z pt pt pt px py pz
|
|
|
160 |
3 copy 12 index /MatrixPQR .knownget { mmult3 } if 6 array astore
|
|
|
161 |
}
|
|
|
162 |
repeat
|
|
|
163 |
% Stack: csdict crdict x y z wps+ bps+ wpd+ bpd+
|
|
|
164 |
9 -1 roll pop % get rid of csdict
|
|
|
165 |
7 4 roll
|
|
|
166 |
7 index /TransformPQR get
|
|
|
167 |
{ % Stack: crdict wps+ bps+ wpd+ bpd+ u v w proc
|
|
|
168 |
8 copy exch pop exch pop
|
|
|
169 |
exec exch pop 4 -1 roll pop
|
|
|
170 |
}
|
|
|
171 |
forall
|
|
|
172 |
7 3 roll pop pop pop pop % get rid of White/BlackPoints
|
|
|
173 |
DOCIEDEBUG { (After TransformPQR: ) print3 } if
|
|
|
174 |
3 index /MatrixPQR .knownget { minvert3 mmult3 } if
|
|
|
175 |
DOCIEDEBUG { (After MatrixPQR': ) print3 } if
|
|
|
176 |
3 index /MatrixLMN .knownget { mmult3 } if
|
|
|
177 |
DOCIEDEBUG { (After MatrixLMN Encode: ) print3 } if
|
|
|
178 |
3 index /EncodeLMN .knownget { apply3 } if
|
|
|
179 |
DOCIEDEBUG { (After EncodeLMN Encode: ) print3 } if
|
|
|
180 |
3 index /RangeLMN .knownget not { r3default } if restrict3
|
|
|
181 |
DOCIEDEBUG { (After RangeLMN Encode: ) print3 } if
|
|
|
182 |
3 index /MatrixABC .knownget { mmult3 } if
|
|
|
183 |
DOCIEDEBUG { (After MatrixABC Encode: ) print3 } if
|
|
|
184 |
3 index /EncodeABC .knownget { apply3 } if
|
|
|
185 |
DOCIEDEBUG { (After EncodeABC Encode: ) print3 } if
|
|
|
186 |
3 index /RangeABC .knownget not { r3default } if
|
|
|
187 |
5 -1 roll /RenderTable .knownget
|
|
|
188 |
{ % Stack: u v w ranges rtable
|
|
|
189 |
5 1 roll rescale3
|
|
|
190 |
DOCIEDEBUG { (Rescaled ABC: ) print3 } if
|
|
|
191 |
% Stack: rtable a b c
|
|
|
192 |
lookup3
|
|
|
193 |
% Stack: rtable bytes
|
|
|
194 |
|
|
|
195 |
{ % Stack: values rtable bytes c
|
|
|
196 |
2 copy get 255 div
|
|
|
197 |
% Stack: values rtable bytes c v
|
|
|
198 |
3 index 3 -1 roll 5 add get exec 3 1 roll
|
|
|
199 |
}
|
|
|
200 |
for pop pop
|
|
|
201 |
DOCIEDEBUG { (After RenderTableT: ) print ] dup == aload pop } if
|
|
|
202 |
}
|
|
|
203 |
{ restrict3
|
|
|
204 |
DOCIEDEBUG { (After RangeABC Encode: ) print3 } if
|
|
|
205 |
}
|
|
|
206 |
ifelse
|
|
|
207 |
} bind def
|
|
|
208 |
|
|
|
209 |
% ---------------- Top level control ---------------- %
|
|
|
210 |
|
|
|
211 |
/mapdict mark
|
|
|
212 |
/CIEBasedA { 1 get exch csciea currentcolorrendering 4 1 roll crmap } bind
|
|
|
213 |
/CIEBasedABC { 1 get 4 1 roll cscieabc currentcolorrendering 4 1 roll crmap } bind
|
|
|
214 |
.dicttomark def
|
|
|
215 |
/mapcie % <a> mapcie <v1> ...
|
|
|
216 |
% <a> <b> <c> mapcie <v1> ...
|
|
|
217 |
{ currentcolorspace dup 0 get //mapdict exch get exec
|
|
|
218 |
} bind def
|