Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/~kupry/SOFTWARE/polar.pro Дата изменения: Wed Aug 6 15:32:52 2003 Дата индексирования: Mon Oct 1 22:39:08 2012 Кодировка:

PRO gfunc, X, A, F, pder
bx = EXP(A[1] * X)
F = A[0] * bx + A[2]
IF N_PARAMS() GE 4 THEN $
pder= [[bx], [A[0] * X * bx], [replicate(1.0, N_ELEMENTS(X))]]
END

function intenzity,pole
vel=n_elements(pole(*,0))
intt=make_array(vel,2,/float)
intt(*,0)=pole(*,1)
atlas=make_array(4400,/float)
openr,unit,'c:\Rsi\Idl54\polarizace\liegout.dat',/get_lun
readf,unit,atlas
free_lun,unit

for i=0,vel-1 do begin
delka1=pole(i,0) & delka1=(round(delka1*100.0))/100.0
index=where(atlas eq delka1,count)
intt(i,1)=atlas(index+1)
endfor

return,intt
end

function vyber,pole,x1,x2
vel=n_elements(pole(*,0))
poleob=make_array(vel,/int)
intervaly=make_array(22,/float)
openr,unit,'c:\Rsi\Idl54\polarizace\Halfap.dat',/get_lun
readf,unit,intervaly
free_lun,unit

for j=x1,x2 do begin
poleob(j)=0
for i=0,10 do if (pole(j,0) gt intervaly(i*2)) and (pole(j,0) lt intervaly(i*2+1)) then poleob(j)=1
endfor

indexy=where(poleob eq 1)
cast=pole(indexy,*)

return,cast
end


pro list_event,event
COMMON spol,item,wave
COMMON spol1,pocet,delky
widget_control,event.id,get_uvalue=uvalue
case uvalue of
'LIST':begin
wave=float(item(event.index))
end
'OK':begin
widget_control,event.top,/destroy
delky(pocet)=wave
pocet=pocet+1
end
endcase
end

function najdi,pole,ind
i=ind
while (pole(195,i) lt 70) do i=i-1
return,i
end

function najdih,pole,ind
i=ind
while (pole(195,i) lt 70) do i=i+1
return,i
end

pro nakresli,zceho,kam,jedna,dva,diff,pozice,lab
COMMON disperze,dispx
if diff ne 0 then begin
if (((jedna-15) gt 0) and ((jedna+15) lt 619)) and (dva lt 569) and ((dva-diff) gt 0) then begin
!P.MULTI=[0,1,2]
widget_control,kam,get_value=nas & wset,nas

y1=make_array(31,value=0) & y2=make_array(31,value=0)

for i=-1,1 do begin
y1=y1+smooth(zceho((jedna-15):(jedna+15),dva+i),4)
y2=y2+smooth(zceho((jedna-15):(jedna+15),dva-diff+i),4)
endfor
y1=round(y1/3.0) & y2=round(y2/3.0)

min1=where(y1 eq min(y1),count1)
min2=where(y2 eq min(y2),count2)

if (count1 eq 1) and (count2 eq 1) then begin
dispx=0
dispx=min2(0)-min1(0)
val1=strcompress(string(dispx),/remove_all)
widget_control,lab,set_value=val1
endif else widget_control,lab,set_value=' '

plot,y1
if (count1 eq 1) then plots,min1(0),min(y1),psym=2,/data
plot,y2
if (count2 eq 1) then plots,min2(0),min(y2),psym=2,/data

device,set_graphics=6
arrow,pozice,0,pozice,399,hsize=0,color=!D.N_colors-1
device,set_graphics=3
endif
endif
end

pro nastav,co
widget_control,co,get_value=nas
wset,nas
end

pro polar_event,event
COMMON profily,profzacx,profkonx,profzacy,profkony,lbcorx,rucorx,lbcory,rucory
COMMON box,box1,box2,calk1,calk2,st1,en1
COMMON kolik,rozhrani,pozicep
COMMON spol,item,wave
COMMON spol1,pocet,delky
COMMON pos,pozice
COMMON baze,base,basep,based,omenu,baset,basepet
COMMON kal,acoef1,bcoef1,ccoef1,acoef2,bcoef2,ccoef2
COMMON kal1,a1,b1,c1,a2,b2,c2
COMMON disperze,dispx
COMMON uloz,nazev,robraz
COMMON slid,line,rotat,rotval,f2,obr,prof,lineprof,roz2,o2,dis2 ;,line3
COMMON dalsi,pos2,pos3,pos4,pos5
COMMON inner,ypos,rotpos,xpos,prvni,druhy,dif;,xpos2,rozsah
COMMON inner2,n1,n2
COMMON inner3,cesta
COMMON posunobr,posun,zac1,kon1
COMMON oblasti,black1,black2,mensi
COMMON interpo,body,pocetint
COMMON knof,f1
widget_control,event.id,get_uvalue=uvalue
case uvalue of

'LOAD':begin
nastav,obr & erase & nastav,prof & erase
cesta=pickfile(filter='*.fts')
obraz=readfits(cesta) & robraz=obraz
widget_control,obr,sensitive=1
nastav,obr
tvscl,obraz
widget_control,line,sensitive=1
widget_control,rotat,sensitive=1
; widget_control,f2,sensitive=1
; widget_control,f22,sensitive=1
widget_control,f1,sensitive=0
widget_control,rotval,sensitive=1
widget_control,lineprof,sensitive=1
widget_control,o2,sensitive=1
widget_control,line,set_value=300 & ypos=300
widget_control,rotat,set_value=0 & rotpos=0.0
widget_control,rotval,set_value='0'
widget_control,lineprof,set_value=10
widget_control,roz2,set_value=' '
widget_control,dis2,set_value=' '
widget_control,pos3,set_value=' '

device,set_graphics=6

arrow,0,300,619,300,hsize=0,color=!D.N_colors-1
nastav,prof
arrow,10,0,10,399,hsize=0,color=!D.N_colors-1
device,set_graphics=3
xpos=10 & prvni=-1 & druhy=-1

dif=0 & n1=0 & n2=0
dispx=0 & posun=0 & pocetint=0
end

'SAVE':begin
nazev='p'+strmid(cesta,strlen(cesta)-10,10)
strput,cesta,nazev,strlen(cesta)-10
cesta=cesta+'s'
val='The result has been saved as a '+nazev
writefits,cesta,tvrd()
basesave=widget_base()
labelsave=widget_label(basesave,value=val)
widget_control,basesave,/realize & wait,2
widget_control,basesave,/destroy
widget_control,event.top,/destroy
polar
end

'QUIT':begin
!p.multi=0
widget_control,event.top,/destroy
if (!d.window ne -1) then wdelete,!d.window
end

'LINE':begin
nastav,obr
device,set_graphics=6
arrow,0,ypos,619,ypos,hsize=0,color=!D.N_colors-1
widget_control,line,get_value=hodnota
arrow,0,hodnota,619,hodnota,hsize=0,color=!D.N_colors-1
device,set_graphics=3
ypos=hodnota
end

'LINE2':begin
nastav,prof
device,set_graphics=6
arrow,xpos,0,xpos,399,hsize=0,color=!D.N_colors-1
widget_control,lineprof,get_value=hodn
arrow,hodn,0,hodn,399,hsize=0,color=!D.N_colors-1
device,set_graphics=3
xpos=hodn
end


'ROT':begin
widget_control,rotat,get_value=hodnp
hodnota=hodnp/10.0

if hodnota lt 0 then strhodnota=strmid(strcompress(string(hodnota),/remove_all),0,4) $
else strhodnota=strmid(strcompress(string(hodnota),/remove_all),0,3)
widget_control,rotval,set_value=strhodnota

otoceni=hodnota-rotpos
robraz=rot(robraz,otoceni,/INTERP)
nastav,obr
tvscl,robraz
device,set_graphics=6

arrow,0,ypos,619,ypos,hsize=0,color=!D.N_colors-1
device,set_graphics=3
rotpos=hodnota

if (prvni ne -1) then nakresli,robraz,prof,prvni,druhy,dif,xpos,dis2

if dif ne 0 then begin
widget_control,pos4,get_value=black1 & widget_control,pos5,get_value=black2
black1=fix(black1(0)) & black2=fix(black2(0))
n1=najdi(robraz,black1)
n2=najdi(robraz,black2)
dif=abs(n1-n2)
val1=strcompress(string(dif),/remove_all)
widget_control,roz2,set_value=val1
endif
end

'OKNO':begin
widget_control,pos2,set_value=STRCOMPRESS('('+string(event.x)+','+string(event.y)+')',/rem)
if event.press eq 1 then begin
nakresli,robraz,prof,event.x,event.y,dif,xpos,dis2
prvni=event.x & druhy=event.y
endif
end

'VER':begin
widget_control,pos4,get_value=black1 & widget_control,pos5,get_value=black2
black1=fix(black1(0)) & black2=fix(black2(0))
if (black1 ne 0) and (black2 ne 0) then begin
n1=najdi(robraz,black1)
n2=najdi(robraz,black2)
dif=abs(n1-n2)
val1=strcompress(string(dif),/remove_all)
widget_control,roz2,set_value=val1
endif
end


'NEXT1':begin
widget_control,profzacx,get_value=lbcorx & lbcorx=fix(lbcorx(0))
widget_control,profkonx,get_value=rucorx & rucorx=fix(rucorx(0))
widget_control,profzacy,get_value=lbcory & lbcory=fix(lbcory(0))
widget_control,profkony,get_value=rucory & rucory=fix(rucory(0))

widget_control,pos3,get_value=pozicep
widget_control,pos5,get_value=rozhrani
rozhrani=fix(rozhrani(0)) & pozicep=fix(pozicep(0))
;print,pozicep,rozhrani,dispx ;!!!!!!!!!!!!!!!!!!!!!!!!!!
if (dispx ne 0) and (dif ne 0) and (rozhrani ne 0) and (pozicep ne 0)then begin
part=shift(robraz(*,0:rozhrani),-dispx,0)
robraz(*,0:rozhrani)=part
widget_control,basep,/destroy
widget_control,omenu,/destroy
based=widget_base(base,/column)
based1=widget_base(based,/row)
drawd=widget_draw(based1,/button_events,/motion_events,uvalue='OBR2',xsize=620,ysize=570)
based11=widget_base(based1,/column)
plab=widget_label(based11,value='POSITION')
pos2=widget_label(based11,value=' ',/frame,xsize=70,/align_center)
pos4=cw_field(based11,title='start for calibration')
pos5=cw_field(based11,title='end for calibration')
box1=cw_field(based11,title='start of box')
box2=cw_field(based11,title='end of box')
based2=widget_base(based,xsize=100)
nbut=widget_button(based2,value='CALIBRATE',uvalue='CAL')
plot,robraz(*,pozicep),position=[0,0.1,1,0.9],xrange=[0,620],xstyle=1
;writefit,'h:\idl\polar\o1.fts',robraz
; writefit,'h:\idl\polar\1.fts',float(robraz(*,pozicep)) ;!!!!!!!!!!!!!!!!!!!!!!!!!111
widget_control,based,/realize
pocet=0
endif
end

'OBR2':begin
widget_control,pos2,set_value=STRCOMPRESS('('+string(event.x)+','+string(event.y)+')',/rem)
if (event.press eq 1) and (pocet lt 10) then begin
if pocet eq 0 then begin
pozice=make_array(11,/int,value=0) & delky=make_array(11,/float,value=0.0)
endif
pozice(pocet)=event.x
wave=6552.629
basecal=widget_base(/column)
item=['6552.629','6553.785','6554.238','6555.466','6556.077','6556.806','6557.171',$
'6558.149','6559.576','6560.555','6561.007','6562.808','6564.206','6565.545','6569.22']
list = widget_droplist(basecal,value=item,uvalue='LIST',ysize=100)
butt=widget_button(basecal,value='OK',uvalue='OK')
widget_control,basecal,/realize
xmanager,'list',basecal

endif
end

'CAL':begin
widget_control,pos4,get_value=st & st=fix(st(0))
widget_control,box1,get_value=st1 & st1=fix(st1(0))
widget_control,pos5,get_value=en & en=fix(en(0))
widget_control,box2,get_value=en1 & en1=fix(en1(0))

if (st ne 0) and (en ne 0) and (pocet gt 0) and $
(st1 ne 0) and (en1 ne 0) then begin
widget_control,based,/destroy

baset=widget_base(base,/column)

baset1=widget_base(baset,/row)
drawt1=widget_draw(baset1,xsize=500,ysize=500)
drawt2=widget_draw(baset1,xsize=500,ysize=500)
baset2=widget_base(baset,/row)
al1=widget_label(baset2,value='1st region',/frame)
al2=widget_label(baset2,value=' ',xsize=100)
bl2=widget_label(baset2,value=' ',xsize=100)
cl2=widget_label(baset2,value=' ',xsize=100)
al11=widget_label(baset2,value='2nd region',/frame)
al22=widget_label(baset2,value=' ',xsize=100)
bl22=widget_label(baset2,value=' ',xsize=100)
cl22=widget_label(baset2,value=' ',xsize=100)

basen=widget_base(baset,xsize=100)
nbut=widget_button(basen,value='NEXT',uvalue='NEXT2')

widget_control,baset,/realize
!p.multi=0

for oblast=1,2 do begin

pozice=pozice(0:pocet-1) & delky=delky(0:pocet-1)
coef=linfit(pozice,delky)
profil=make_array(620,2,/float)
if oblast eq 1 then profil(*,1)=robraz(*,pozicep) else profil(*,1)=robraz(*,pozicep-dif)
for i=0,619 do profil(i,0)=coef(0)+coef(1)*i
if oblast eq 1 then begin
nastav,drawt1
plot,profil(*,0),yrange=[6552,6572],ystyle=1
oplot,pozice,delky,psym=2
endif

profil2=vyber(profil,st,en)
inten=intenzity(profil2)
c=10
X=smooth(inten(*,0),c)
Y=smooth(inten(*,1),c)
W=1.0/Y
A=[20.0,0.0085,-30.0]
result=curvefit(X,Y,W,A,sigma,function_name='gfunc')
x=findgen(150)+80
nastav,drawt2
if oblast eq 1 then plot,x,a(0)*exp(x*a(1))+a(2) else oplot,x,a(0)*exp(x*a(1))+a(2)
oplot,smooth(inten(*,0),c),smooth(inten(*,1),c),psym=1
if oblast eq 1 then calk1=a else calk2=a
if oblast eq 1 then begin
widget_control,al2,set_value=string(a(0)) & widget_control,bl2,set_value=string(a(1))
widget_control,cl2,set_value=string(a(2))
endif else begin
widget_control,al22,set_value=string(a(0)) & widget_control,bl22,set_value=string(a(1))
widget_control,cl22,set_value=string(a(2))
endelse

endfor
endif
end

'NEXT2':begin
widget_control,baset,/destroy
basepet=widget_base(base,/column)
intdraw=widget_draw(basepet,xsize=800,ysize=400,/button_events,uvalue='INTER')
basepet1=widget_base(basepet,xsize=100)
nbut=widget_button(basepet1,value='NEXT',uvalue='NEXT3')
widget_control,basepet,/realize
nastav,intdraw
zac1=najdi(robraz,black1)
kon1=najdih(robraz,black2)
mensi=robraz & mensi(*,0:(zac1-kon1))=robraz(*,kon1-dif:zac1-dif)
mensi(*,zac1-kon1+1:(2*zac1-2*kon1+1))=robraz(*,kon1:zac1)
mensi=mensi(*,0:(2*zac1-2*kon1+1))


boxx2=float(mensi(st1:en1,*))
mensi(*,0:zac1-kon1)=calk2(0)*exp(float(mensi(*,0:zac1-kon1))*calk2(1))+calk2(2)
mensi(*,zac1-kon1+1:(2*zac1-2*kon1+1))=calk1(0)*exp(float(mensi(*,zac1-kon1+1:(2*zac1-2*kon1+1)))*calk1(1))+calk1(2)

s=size(boxx2)
boxx1=make_array(s(2),/float)

for i=0,s(2)-1 do begin
boxx1(i)=total(boxx2(*,i))/float(n_elements(boxx2(*,i)))
;mensi(*,i)=(90*mensi(*,i)/boxx1(i))
endfor

plot,boxx1
body=make_array(100,2)
end

'NEXT3':begin
if pocetint gt 0 then begin
s=size(mensi)
U=findgen(s(2))
X=body(0:pocetint-1,0)
Y=body(0:pocetint-1,1)
R=interpol(Y,X,U)
oplot,U,R,color=100
wait,2
for i=0,s(2)-1 do mensi(*,i)=90.0*float(mensi(*,i))/float(R(i))

widget_control,basepet,/destroy
basec=widget_base(base)
drawc=widget_draw(basec,xsize=620,ysize=s(2)/2)
widget_control,basec,/realize

mensi1=mensi(*,0:zac1-kon1)
mensi2=mensi(*,zac1-kon1+1:(2*zac1-2*kon1+1))
mensivysl=make_array(s(1),zac1-kon1,/byte)
mensivysl=float(mensi1)-float(mensi2)
tvscl,mensivysl

;-------------------profily-------
if (lbcorx ne 0) and (lbcory ne 0) and (rucorx ne 0) and (rucory ne 0) then begin
horni=mensi2(lbcorx:rucorx,(lbcory-kon1):(rucory-kon1))
ips1=make_array(rucorx-lbcorx+1,value=0,/double)
;ips1=horni(*,0)
for i=0,(rucory-lbcory) do ips1=ips1+double(horni(*,i))
ips1=ips1/double(rucory-lbcory+1) & ips1=0.9*ips1/90.0

dolni=mensi1(lbcorx:rucorx,(lbcory-kon1):(rucory-kon1))
ips2=make_array(rucorx-lbcorx+1,value=0,/double)
;ips2=dolni(*,0)
for i=0,(rucory-lbcory) do ips2=ips2+double(dolni(*,i))
ips2=ips2/double(rucory-lbcory+1) & ips2=0.9*ips2/90.0
res=poly_fit(pozice,delky,1)
iks=float(findgen(620)) & iks=res(0)+iks*res(1) & iks=iks(lbcorx:rucorx)
profil1=dblarr(n_elements(iks),2) & profil2=dblarr(n_elements(iks),2)
profil1(*,0)=iks & profil2(*,0)=iks & profil1(*,1)=ips1 & profil2(*,1)=ips2
writefits,'profil1.fts',profil1 & writefits,'profil2.fts',profil2
basehlas=widget_base()
labhlas=widget_label(basehlas,value='Profiles have been saved.')
widget_control,basehlas,/realize & wait,2
widget_control,basehlas,/destroy
;tvscl,mensi2 & plots,lbcorx,lbcory-kon1,psym=2,/device
;plots,rucorx,rucory-kon1,psym=2,/device & wait,2
;tvscl,mensi1 & plots,rucorx,rucory-kon1,psym=2,/device
;plots,lbcorx,lbcory-kon1,psym=2,/device

endif
;----------------------------
widget_control,f2,sensitive=1
;widget_control,f22,sensitive=1
endif
end

'INTER':begin
if event.press eq 1 then begin
plots,event.x,event.y,psym=2,/device
coord=convert_coord(event.x,event.y,/device,/to_data)
body(pocetint,0)=coord(0) & body(pocetint,1)=coord(1)
pocetint=pocetint+1
endif
end
;'SAVE2':begin




; end
endcase
end

pro polar
COMMON profily,profzacx,profkonx,profzacy,profkony,lbcorx,rucorx,lbcory,rucory
COMMON slid,line,rotat,rotval,f2,obr,prof,lineprof,roz2,o2,dis2;line3
COMMON dalsi,pos2,pos3,pos4,pos5
COMMON baze,base,basep,based,omenu,basepet
COMMON knof,f1
device,set_graphics=3

base=widget_base(/row,mbar=bar,xsize=1020,ysize=700)
fmenu=widget_button(bar,value='FILE',/menu)
f1=widget_button(fmenu,value='LOAD',uvalue='LOAD')
f2=widget_button(fmenu,value='SAVE',uvalue='SAVE')
; f22=widget_button(fmenu,value='SAVE PROFILES',uvalue='SAVE2')
f3=widget_button(fmenu,value='QUIT',uvalue='QUIT')
omenu=widget_button(bar,value='OPTIONS',/menu)
o2=widget_button(omenu,value='VERTICAL DIFFERENCE',uvalue='VER')

basep=widget_base(base,/row)
base1=widget_base(basep,/column)
obr=widget_draw(base1,xsize=620,ysize=570,/button_events,/motion_events,uvalue='OKNO')
base3=widget_base(base1,/row)
roz1=widget_label(base3,value='VERTICAL DIFFERENCE OF REGIONS')
roz2=widget_label(base3,value=' ',/frame,xsize=70,/align_center)
dis1=widget_label(base3,value='DISPERSION (px)')
dis2=widget_label(base3,value=' ',/frame,xsize=70,/align_center)
pos1=widget_label(base3,value='POSITION')
pos2=widget_label(base3,value=' ',/frame,xsize=70,/align_center)
base4=widget_base(base1,/row,xsize=500)
pos3=cw_field(base4,title='y profile position')
pos4=cw_field(base4,title='1st region')
pos5=cw_field(base4,title='2nd region')


base5=widget_base(base1,xsize=100)
nbut1=widget_button(base5,value='NEXT',uvalue='NEXT1')
base2=widget_base(basep,/column)
prof=widget_draw(base2,xsize=370,ysize=400)
lineprof=widget_slider(base2,title='Line motion (profiles)',uvalue='LINE2', $
value=10,min=0,max=369)

line=widget_slider(base2,title='Line motion (image)',uvalue='LINE', $
value=300,min=0,max=569)

rotval=widget_label(base2,value='0.0',xsize=200,/align_center)
rotat=widget_slider(base2,title='Image rotation',uvalue='ROT', $
value=0,min=-30,max=30,/suppress_value)
base6=widget_base(base2,/row)
profzacx=cw_field(base6,title='1st corner (x)')
profzacy=cw_field(base6,title='1nd corner (y)')
base7=widget_base(base2,/row)
profkonx=cw_field(base7,title='2nd corner (x)')
profkony=cw_field(base7,title='2nd corner (y)')

widget_control,line,sensitive=0
widget_control,rotat,sensitive=0
widget_control,f2,sensitive=0
;widget_control,f22,sensitive=0
widget_control,rotval,sensitive=0
widget_control,lineprof,sensitive=0
widget_control,base,/realize
widget_control,obr,sensitive=0
widget_control,o2,sensitive=0
xmanager,'polar',base
end