@@ -149,13 +149,13 @@ pic1 = phantom(x, y, [ellipse1])
149
149
150
150
marker = :star
151
151
p0 = jim (x, y, pic0, " Original ellipse" ;
152
- xlabel= " x" , ylabel= " y" , prompt= :false )
152
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 300 ), prompt= :false )
153
153
x0,y0 = 7 ,0
154
154
scatter! ([x0], [y0], color= :blue ; marker)
155
155
point1 = [cos (ϕ1) - sin (ϕ1); sin (ϕ1) cos (ϕ1)] * [x0; y0] # rotate point
156
156
x1,y1 = point1[1 ], point1[2 ]
157
157
p1 = jim (x, y, pic1, " Rotated by ϕ = $ϕ1s " ;
158
- xlabel= " x" , ylabel= " y" , prompt= :false )
158
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 300 ), prompt= :false )
159
159
scatter! ([x1], [y1], color= :blue ; marker)
160
160
jim (p0, p1)
161
161
@@ -221,15 +221,15 @@ pic1 = phantom(x, y, z, [ellipsoid1])
221
221
222
222
p0z = jim (x, y, pic0,
223
223
" Original ellipsoid:\n (x,y) slice" ;
224
- xlabel= " x" , ylabel= " y" , prompt= :false )
224
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 300 ), prompt= :false )
225
225
x0,y0 = 7 ,0
226
226
scatter! ([x0], [y0], color= :blue ; marker)
227
227
Rz (ϕ) = [cos (ϕ) - sin (ϕ) 0 ; sin (ϕ) cos (ϕ) 0 ; 0 0 1 ]
228
228
point1 = Rz (ϕ1) * [x0; y0; 0 ] # rotate point
229
229
x1,y1 = point1[1 ], point1[2 ]
230
230
p1z = jim (x, y, pic1,
231
231
" Rotated about z\n by ϕ = $ϕ1s \n (z out of 'board')" ;
232
- xlabel= " x" , ylabel= " y" , prompt= :false )
232
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 350 ), prompt= :false )
233
233
scatter! ([x1], [y1], color= :blue ; marker)
234
234
jim (p0z, p1z)
235
235
@@ -285,15 +285,15 @@ pic1 = phantom(x, y, z, [ellipsoid1]); pic1 = selectdim(pic1, 2, 1)
285
285
286
286
p0y = jim (x, z, pic0,
287
287
" Original ellipsoid:\n (x,z) slice" ;
288
- xlabel= " x" , ylabel= " z" , prompt= false )
288
+ xlabel= " x" , ylabel= " z" , size = ( 700 , 350 ), prompt= false )
289
289
x0,z0 = 7 ,0
290
290
scatter! ([x0], [z0], color= :green ; marker)
291
291
Ry (θ) = [cos (θ) 0 sin (θ); 0 1 0 ; - sin (θ) 0 cos (θ)]
292
292
point1 = Ry (θ1) * [x0; 0 ; z0] # rotate point
293
293
x1,z1 = point1[1 ], point1[3 ]
294
294
p1y = jim (x, z, pic1,
295
295
" Rotated about y\n by θ = $θ1s \n (y into 'board')" ;
296
- xlabel= " x" , ylabel= " z" , prompt= false )
296
+ xlabel= " x" , ylabel= " z" , size = ( 700 , 350 ), prompt= false )
297
297
scatter! ([x1], [z1], color= :green ; marker)
298
298
jim (p0y, p1y)
299
299
@@ -347,15 +347,15 @@ pic1 = phantom(x, y, z, [ellipsoid1]); pic1 = selectdim(pic1, 1, 1)
347
347
348
348
p0x = jim (y, z, pic0,
349
349
" Original ellipsoid:\n (y,z) slice)" ;
350
- xlabel= " y" , ylabel= " z" , prompt= false )
350
+ xlabel= " y" , ylabel= " z" , size = ( 700 , 350 ), prompt= false )
351
351
y0,z0 = 3 ,0
352
352
scatter! ([y0], [z0], color= :red ; marker)
353
353
Rx (ψ) = [1 0 0 ; 0 cos (ψ) - sin (ψ); 0 sin (ψ) cos (ψ)]
354
354
point1 = Rx (ψ1) * [0 ; y0; z0] # rotate point
355
355
y1,z1 = point1[2 ], point1[3 ]
356
356
p1x = jim (y, z, pic1,
357
357
" Rotated about x\n by ψ = $ψ1s \n (x out of 'board')" ;
358
- xlabel= " y" , ylabel= " z" , prompt= false )
358
+ xlabel= " y" , ylabel= " z" , size = ( 700 , 350 ), prompt= false )
359
359
scatter! ([y1], [z1], color= :red ; marker)
360
360
jim (p0x, p1x)
361
361
@@ -441,10 +441,10 @@ pic1 = phantom(x, y, z, [ellipsoid1])
441
441
442
442
p0a = jim (mid3 (pic0),
443
443
" Original ellipsoid\n (central slices)" ;
444
- xlabel= " x" , ylabel= " y" , prompt= :false )
444
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 320 ), prompt= :false )
445
445
p1a = jim (mid3 (pic1),
446
446
" Rotated\n ϕ = $ϕ1s , θ = $θ1s , ψ = $ψ1s " ;
447
- xlabel= " x" , ylabel= " y" , prompt= :false )
447
+ xlabel= " x" , ylabel= " y" , size = ( 700 , 320 ), prompt= :false )
448
448
jim (p0a, p1a)
449
449
450
450
0 commit comments