;;; generate mosaic of filament.  Use Ha wing data

GOTO, weiter

fh1 = findfile('Filament-ha-08*')
fh2 = findfile('Filament-ha-00*') & fh2=fh2(1:*)

fg = 'Filament-gb-0000'+strmid(fh1, 16)
fb = 'Filament-bc-0000'+strmid(fh1, 16)
fr = 'Filament-rc-0000'+strmid(fh1, 16)
fc = 'Filament-ca-0000'+strmid(fh1, 16)

par_bc = [-0.4102, 1.00254]
par_rc = [-0.8612, 1.00474]
par_ca = [-0.2845, 1.00385]
par_ha = [-0.7507, 1.00389]


im = 0
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh1 = [[0, 0], [4, 7], [11, 49], [-18, -15], [0, 7]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))

x = 1206 & y = 899
part1=[[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
       [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
       [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
       [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
       [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

im = 1
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh2 = [[0, 0], [3, 9], [0, 58], [-21, -15], [-2, 21]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1212 & y=917
part2 = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
         [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
         [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
         [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
         [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

im = 2
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh3 = [[0, 0], [3, 6], [9, 43], [-20, -16], [8, 6]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1206 & y=904
part3 = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
         [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
         [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
         [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
         [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

im = 3
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh4 = [[0, 0], [4, 8], [11, 44], [-19, -16], [10, 7]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1205 & y=903
part4 = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
         [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
         [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
         [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
         [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

im = 4
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh5 = [[0, 0], [-3, 6], [8, 43], [-20, -16], [6, 6]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1207 & y=904
part5 = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
         [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
         [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
         [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
         [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

im = 5
rdfits, p, fg(im) & g=reform(p(*, *, 0))
rdfits, p, fb(im) & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, fr(im) & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, fc(im) & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
rdfits, p, fh1(im) & h=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

sh = intarr(2, 5)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh(*, 4) = shc(g, h, /filt)
sh6 = [[0, 0], [3, 7], [7, 56], [-21, -15], [6, 18]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1208 & y=919
part6 = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
         [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
         [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
         [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]], $
         [[(shift(h, sh(0, 4), sh(1, 4)))(0:x, 0:y)]]]

;;; parts 2 and 3 don't overlap!  guiding error. create one (without h-alpha
;;; wing) from the other timestep

im = '-0000-20041006-083518.fits'
rdfits, p, 'F*gb'+im & g=reform(p(*, *, 0))
rdfits, p, 'F*bc'+im & b=rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
rdfits, p, 'F*rc'+im & r=rot(reform(p(*, *, 0)), par_rc(0), par_rc(1), cub=-.5)
rdfits, p, 'F*ca'+im & c=rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
sh = intarr(2, 4)
sh(*, 1) = shc(g, b, /filt, /N2)
sh(*, 2) = shc(g, r, /filt, /N2)
sh(*, 3) = shc(g-smooth(g, 7, /edge), c-smooth(c, 7, /edge))
sh3a = [[0, 0], [4, 8], [11, 58], [-19, -16]]
sh(0, *) = sh(0, *)-max(sh(0, *))
sh(1, *) = sh(1, *)-max(sh(1, *))
x = 1205 & y=916
part3a = [[[(shift(g, sh(0, 0), sh(1, 0)))(0:x, 0:y)]], $
          [[(shift(b, sh(0, 1), sh(1, 1)))(0:x, 0:y)]], $
          [[(shift(r, sh(0, 2), sh(1, 2)))(0:x, 0:y)]], $
          [[(shift(c, sh(0, 3), sh(1, 3)))(0:x, 0:y)]]]

print, shc(part6(*, 756:*, 0), part5(1:*, *, 0), /filt)
print, shc(part5(8:*, 702:*, 0), part4(*, *, 0), /filt)
print, shc(part4(9:*, 706:*, 0), part3a(*, *, 0), /filt)
print, shc(part3a(*, 795:*, 0), part2(1:*, *, 0), /filt)
print, shc(part2(20:*, 690:*, 0), part1(*, *, 0), /filt)

Weiter:

x = [10, 9, 17, 26, 25, 45]
y = [10, 766, 1468, 2174, 2969, 3659]

mos = fltarr(869, 4549, 6)
FOR im=0, 4 DO BEGIN
    big = fltarr(1300, 4560)
    w = big
    paste, big, part6(*, *, im), w, x(0), y(0)
    paste, big, part5(*, *, im), w, x(1), y(1)
    paste, big, part4(*, *, im), w, x(2), y(2)
    IF im EQ 4 THEN $
      paste, big, part3(*, *, im), w, x(3)-14, y(3)-272 $
    ELSE $
      paste, big, part3a(*, *, im), w, x(3), y(3)
    paste, big, part2(*, *, im), w, x(4), y(4)
    paste, big, part1(*, *, im), w, x(5), y(5)
    ix = where(w NE 0)
    sum = big
    sum(ix) = big(ix)/w(ix)
    sum = sum(45:913, 10:4558)
    show, sum, im+1, sc=.2
    IF im EQ 4 THEN sum(where(sum eq 0)) = 1.
    mos(*, *, im) = sum
ENDFOR

big = fltarr(1300, 4603)
w = big

rdfits, p, fh2(5) & h1=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)
rdfits, p, fh2(4) & h2=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)
rdfits, p, fh2(3) & h3=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)
rdfits, p, fh2(2) & h4=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)
rdfits, p, fh2(1) & h5=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)
rdfits, p, fh2(0) & h6=rot(reform(p(*, *, 0)), par_ha(0), par_ha(1), cub=-.5)

print, shc(h1(*, 751:*), h2(3:*, *))
print, shc(h2(5:*, 700:*), h3(*, *))
print, shc(h3(19:*, 698:*), h4(*, *))
print, shc(h4(2:*, 820:*), h5(*, *))
print, shc(h5(17:*, 657:*), h6(*, *))

paste, big, h1, w, 10, 0
paste, big, h2, w, 7, 751
paste, big, h3, w, 12, 1451
paste, big, h4, w, 31, 2149
paste, big, h5, w, 33, 2969
paste, big, h6, w, 50, 3626
ix = where(w NE 0)
sum = big
sum(ix) = big(ix)/w(ix)
sum = sum(50:982, *)

mos(1:*, *, 5) = sum(0:867, 27:4575)
mos = mos(1:*, *, *)

END