【円弧-22】緯線・経線による球の表現

【円弧-20】と【円弧-21】を足せば地球儀のような球が表現できます。

参考：数式によるパソコン・グラッフィクスのテクニック／著者：依田谹（ヨダコダマ）／昭和61年８月１日初版発行／発行所：一橋出版株式会社

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

%!PS-Adobe-3.0 EPSF-3.0
%%BoundingBox: 0 0 720 720
%%HiResBoundingBox: 0 0 720 720
%%CropBox: 0 0 720 720
%%Title:緯線・経線による球の表現
%%Creator:Studio Fruit Jam / Toyokazu Nishi
%%Copyright:Studio Fruit Jam / Toyokazu Nishi
%%CreationDate:2014年11月30日日曜日11:35:39
%%ModificationDate:2018年10月4日木曜日21:24:04

% 数式によるパソコン・グラッフィクスのテクニック
% 著者：依田谹（ヨダコダマ）
% 昭和61年８月１日初版発行
% 発行所：一橋出版株式会社

% ================ 座標変換 ================
360 360 translate
1 dup scale
% ================ 背景 =================
0 0 0 setrgbcolor
-360 360 720 -720 rectfill
% ================ 定義 ================
1 1 1 setrgbcolor
1 setlinewidth
/p 180 def % 円の角度
/n 20 def % 緯度（横）・経度(縦) の数
/k p n div def % th2の角度
/xa -45 def % 球のX軸角度 a:angle
/ya -45 def % 球のY軸角度
/za 0 def % 球のZ軸角度
/r 300 def % 球のサイズ

% ================ プロシージャ ================
% 図形の回転
/rotation {
% x軸角度
/x1 x def
/y1 y xa cos mul z xa sin mul sub def
/z1 y xa sin mul z xa cos mul add def
% y軸角度
/x2 x1 ya cos mul z1 ya sin mul add def
/y2 y1 def
/z2 x1 neg ya sin mul z1 ya cos mul add def
% Z軸角度
/x3 x2 za cos mul y2 za sin mul sub def
/y3 x2 za sin mul y2 za cos mul add def
/z3 z2 def
} def
% 描画
/line {
/px x3 def
/py y3 neg def
sw 1 eq { px py lineto }{ px py moveto } ifelse
/sw 1 def
} def

% ================ メイン================
newpath % パスの初期化
0 0 r 0 360 arc % 輪郭線

% ================ 球の関数 % 緯度（横）
0 k p { % 緯線の角度
/th2 exch def
/sw 0 def % 余計な線を出さない
0 1 2 p mul { % 緯線の円
/th exch def
/x r th sin mul th2 sin mul def
/y r th2 cos mul def
/z r th cos mul th2 sin mul def
rotation % 図形の回転
% 陰線処理 z3の負の値は回転角度によって変える
z3 -10 lt { /sw 0 def } { line } ifelse
} for
} for
% ================ 球の関数 % %経度（縦）
0 k p 1 sub { % %経線の角度
/th2 exch def
/sw 0 def % 余計な線を出さない
0 1 2 p mul { % 緯線の円
/th exch def
/x r th sin mul th2 cos mul def
/y r th cos mul def
/z r th sin mul th2 sin mul def
rotation % 図形の回転
% 陰線処理 z3の負の値は回転角度によって変える
z3 -10 lt { /sw 0 def } { line } ifelse
} for
} for
stroke