input moments.sh
show.pol
input moments.sh
show.pol
print e1 e2
input moments.sh
show.pol
print theta e1 e2
set fr = f0
rotate fr 190
tv -n tv fr -0.01 0.5
rotate fr 90
set fr = f0
rotate fr 90
tv -n tv fr -0.01 0.5
set fr = f0
rotate fr 95
tv -n tv fr -0.01 0.5
rotate fr 85
tv -n tv fr -0.01 0.5
set fr = f0
rotate fr 85
tv -n tv fr -0.01 0.5
star fr {fr[][0]/2} {fr[0][]/2} 71
echo star fr {fr[][0]/2} {fr[0][]/2} 71
buff
input moments.sh; test.rot
tv fr -0.01 0.2
set fr = f0
rotate fr 5
tv fr
input moments.sh; test.rot
tv fr -0.01 0.2
tv f0
set fr = f0
rotate fr 2
tv fr
input moments.sh; test.rot; tv f0
input moments.sh; test.rot; tv f0
tv fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
resize 800 800
set fr = f0
rotate fr 10
tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
buffers
echo {fr[][0]/2} {fr[0][]/2}
$x = 360 - 720/2
$y = 360 - 720/2
$xo = 360 - 720/2
$yo = 360 - 720/2
$x = $xo*dcos(2) + $yo*dsin(2)
$y = $yo*dcos(2) - $xo*dsin(2)
echo $x
echo $y
echo $xo
echo quit
echo {dcos(2)} {dsin(2)}
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
echo {fr[][0]/2}  {fr[0][]/2}
input moments.sh; test.rot; tv -ch 1 f0; tv -ch 2 fr
header
head f0
head fr
show.pol
lim -n 2 theta e1; clear; box; plot theta e1
lim -n 2 theta e2; clear; box; plot theta e2
input moments.sh; show.pol
lim -n 2 theta e2; clear; box; plot theta e2
plot -c blue theta e1
vstat e1
tv -n 0 -0.01 0.2 f0
tv -n 0 f0 -0.01 0.2
input moments.sh; show.pol
tv -n 0 f0 -0.01 0.2
input moments.sh; show.pol
tv -n 0 f0 -0.01 0.2
vstat e1
vstat e2
imsmooth.2d
imsmooth.generic
imsmooth
fft2d
input moments.sh
input moments.sh
opihi verbose on
input moments.sh
pwd
!cat moments.sh
input moments.sh
opihi verbose off
show.smear
show.smear 1.0
buffers
input moments.sh
show.smear 1.0
input moments.sh
show.smear 1.0
buffers
tv -n 0 -ch 1 f0 -0.01 0.5
tv -n 0 -ch 1 f0 -0.01 1.0
tv -n 0 -ch 2 psf -0.01 1.0
tv -n 0 -ch 2 psf -0.01 0.5
tv -n 0 -ch 2 fs -0.01 1.0
stat fs
tv -n 0 -ch 2 fs 0 2e-9
tv im_R -0.01 0.1
stat im_R
tv im_R -0.0005 0.001
tv im_I
stat im_I
tv kr_R -0.0005 0.001
tv ff_R
stat ff_R
stat kr_R
set ff_R = ff_R * 512^2
set ff_I = ff_I * 512^2
stat ff_R
tv ff_R 0 0.05
tv ff_I 0 0.05
stat ff_I
fft2d -inverse ff_R ff_I to ot_R ot_I
stat ot_R
set out_R = out_R * 512^2
set out_R = ot_R * 512^2
set out_I = ot_I * 512^2
stat out_R
stat out_I
tv out_R
set fs = fs * 512^4
tv fs
tv fs 0 152
tv fs 0 512
input moments.sh
show.smear 1.0
tv -n 0 -ch 2 fs 0 1
tv -n 0 -ch 2 psf 0 1
tv -n 0 -ch 2 psf 0 0.5
tv -n 0 -ch 1 f0 0 0.5
tv -n 0 -ch 3 fs 0 0.5
star
star f0 256 256 128
star -q f0 256 256 128; $Mxx = ($SXg/2.355)^2; $Myy = ($SYg/2.355)^2;
??
star -q f0 256 256 128; $Mxx = ($SXg/2.355)^2; $Myy = ($SYg/2.355)^2; echo $Zcg $Mxx $SXYg $Myy
star -q psf 256 256 128; $Mxx = ($SXg/2.355)^2; $Myy = ($SYg/2.355)^2; echo $Zcg $Mxx $SXYg $Myy
star -q fs 256 256 128; $Mxx = ($SXg/2.355)^2; $Myy = ($SYg/2.355)^2; echo $Zcg $Mxx $SXYg $Myy
echo {1256.63706807*25.1327415272}
input moments.sh
show.smear
show.smear 1
input moments.sh
show.smear
show.smear 1
input moments.sh
show.smear 1
lim -n 2 theta e1_o; clear; box; plot theta e1_o -c blue; plot theta e1_s -c red
vectors
create t 0 720 5
vectors
create t 0 360 5
vectors
lim -n 2 t e1_o; clear; box; plot t e1_o -c blue; plot t e1_s -c red
lim -n 2 t e1_o; clear; box; plot t e1_o -c blue; plot t e1_s -c red; plot t e2_o -c blue -pt box; plot t e2_s -c red -pt box
lim -n 2 t e1_o; clear; box; plot t e1_o -c blue; plot t e1_s -c red; plot t e2_o -c blue -pt obox; plot t e2_s -c red -pt obox
input moments.sh
show.smear 1
star psf {psf[][0]/2} {psf[0][]/2} 128
$Mxx_p = ($SXg/2.355)^2; $Myy_p = ($SYg/2.355)^2; $Mxy_p = $SXYg
lim -n 2 theta Mxx_s; clear; box; plot t Mxx_o -c blue; plot t Mxx_s -c red
set Mxx_S = Mxx_o + $Mxx_p
lim -n 2 theta Myy_s; clear; box; plot theta Myy_o -c blue; plot theta Myy_s -c red
lim -n 2 theta 0 200; clear; box; plot theta Myy_o -c blue; plot theta Myy_s -c red
lim -n 2 theta 0 150; clear; box; plot theta Myy_o -c blue; plot theta Myy_s -c red
set dMxx = Mxx_s - Mxx_o
set dMxy = Mxy_s - Mxy_o
set dMyy = Myy_s - Myy_o
lim -n 2 theta dMxx; clear; box; plot theta dMxx
line -c red $Mxx_p 0 to $Mxx_p 360
line -c red 0 $Mxx_p to 360 $Mxx_p
lim -n 2 theta dMxx; clear; box; plot theta dMxx; line -c red 0 $Mxx_p to 360 $Mxx_p
lim -n 2 theta dMyy; clear; box; plot theta dMyy; line -c red 0 $Myy_p to 360 $Myy_p
vstat dMyy
lim -n 2 theta 35 37; clear; box; plot theta dMyy; line -c red 0 $Myy_p to 360 $Myy_p
lim -n 2 theta dMxx; clear; box; plot theta dMxx; line -c red 0 $Mxx_p to 360 $Mxx_p
lim -n 2 theta dMxy; clear; box; plot theta dMxy; line -c red 0 $Mxy_p to 360 $Mxy_p
lim -n 1 theta Mxy_s; clear; box; plot theta Mxy_s -c red -pt cir -sz 2; plot theta Mxy_o -c blue -pt cir -sz 2; line -c red 0 $Mxy_p to 360 $Mxy_p
input moments.sh
show.smear 1
lim -n 1 theta Mxy_s; clear; box; plot theta Mxy_s -c red -pt cir -sz 2; plot theta Mxy_o -c blue -pt cir -sz 2; line -c red 0 $Mxy_p to 360 $Mxy_p
set dMxy = Mxy_s - Mxy_o
lim -n 2 theta dMxy; clear; box; plot theta dMxy; line -c red 0 $Mxy_p to 360 $Mxy_p
input moments.sh
test.convolve
test.convolve 1
lim theta Mxx_s; clear; box; plot theta Mxx_s
lim theta Mxx_s; clear; box; plot theta Mxx_s; plot -c red theta Mxx_o
set dMxx = Mxx_s - Mxx_o
lim -n 1 theta dMxx; clear; box; plot theta dMxx
line 0 $Mxx_p to 360 $Mxx_p -c red
set dMyy = Myy_s - Myy_o
lim theta Myy_s; clear; box; plot theta Myy_s; plot -c red theta Myy_o
lim -n 0 theta 0 110; clear; box; plot theta Myy_s; plot -c red theta Myy_o
lim -n 0 theta 0 160; clear; box; plot theta Myy_s; plot -c red theta Myy_o
lim -n 0 theta 0 150; clear; box; plot theta Myy_s; plot -c red theta Myy_o
lim -n 1 theta dMyy; clear; box; plot theta dMyy
line 0 $Myy_p to 360 $Myy_p -c red
lim -n 1 theta 30 40; clear; box; plot theta dMyy; line 0 $Myy_p to 360 $Myy_p -c red
lim -n 1 theta 35 37; clear; box; plot theta dMyy; line 0 $Myy_p to 360 $Myy_p -c red
lim -n 1 theta 35.9 36.1; clear; box; plot theta dMyy; line 0 $Myy_p to 360 $Myy_p -c red
lim -n 1 theta 35.99 36.01; clear; box; plot theta dMyy; line 0 $Myy_p to 360 $Myy_p -c red
lim -n 0 theta -10 10; clear; box; plot theta Mxy_s; plot -c red theta Mxy_o
lim -n 0 theta Mxy_s; clear; box; plot theta Mxy_s; plot -c red theta Mxy_o
lim -n 0 theta Mxy_o; clear; box; plot theta Mxy_s; plot -c red theta Mxy_o
lim -n 0 theta Mxy_o; clear; box; plot theta Mxy_s; plot -c red theta Mxy_o -pt ocir -sz 2
input moments.sh
test.convolve
test.convolve 0
lim -n 0 theta 0 150; clear; box; plot theta Myy_s; plot -c red theta Myy_o
set dMyy = Myy_s - Myy_o
set dMxy = Mxy_s - Mxy_o
set dMxx = Mxx_s - Mxx_o
lim -n 1 theta dMxx; clear; box; plot theta dMxx
lim -n 1 theta dMxx; clear; box; plot theta dMxx; line 0 $Mxx_p to 360 $Mxx_p -c red
lim -n 0 theta 0 150; clear; box; plot theta Mxx_s; plot -c red theta Mxx_o
lim -n 0 theta 0 120; clear; box; plot theta Mxx_s; plot -c red theta Mxx_o
lim -n 0 theta 0 150; clear; box; plot theta Myy_s; plot -c red theta Myy_o
lim -n 1 theta dMyy; clear; box; plot theta dMyy; line 0 $Myy_p to 360 $Myy_p -c red
echo $Mxx_p $Myy_p
lim -n 0 theta Mxy_s; clear; box; plot theta Mxy_s; plot -c red theta Mxy_o
lim -n 1 theta dMxy; clear; box; plot theta dMxy; line 0 $Mxy_p to 360 $Mxy_p -c red
pwd
