Документ взят из кэша поисковой машины. Адрес оригинального документа : http://lnfm1.sai.msu.ru/~kupry/SOFTWARE/foto2.pro Дата изменения: Thu Mar 18 17:19:50 2004 Дата индексирования: Mon Oct 1 22:39:05 2012 Кодировка: Поисковые слова: внешние планеты

pro prvni_event,event
COMMON profily,profil1,profil2
COMMON baze2,prvnilab,prvnidraw,prvnibut,prvni1,prvni
COMMON pozice,pozice1,pozice2
COMMON kolik,pocet2
COMMON data,kliny,obrazek
COMMON uloz,cesta
COMMON uloz2,obrazek2
COMMON novy,prvniposlab
;COMMON velikost,s1,s2
widget_control,event.id,get_uvalue=uvalue
case uvalue of
'CLICK1':begin
coord=round(convert_coord(event.x,event.y,/device,/to_data))
widget_control,prvniposlab,set_value=strcompress('('+string(coord(0))+','+string(coord(1))+')',/remove_all)
if (event.press eq 1) then begin
iks=coord(0)
index=where(pozice1 eq 0)
pozice1(index(0))=iks
print,iks
endif
end

'CLICK2':begin
coord=round(convert_coord(event.x,event.y,/device,/to_data))
widget_control,prvniposlab,set_value=strcompress('('+string(coord(0))+','+string(coord(1))+')',/remove_all)
if (event.press eq 1) then begin
pocet2=pocet2+1
index=where(pozice1 eq 0)
pocet1=index(0)

iks=coord(0)
index=where(pozice2 eq 0)
pozice2(index(0))=iks
print,iks

if pocet2 eq pocet1 then begin
widget_control,prvni1,/destroy
posun=round(float(total(pozice1-pozice2))/pocet1)
print,'shift = ',posun

if posun gt 0 then begin ;posun doprava (zvetseni obrazku)
s=size(obrazek) ;duvod - posledni cara v obrazku by se dostala posunem doleva
;obrazek2=bytarr(s(1)+posun,s(2))
obrazek2=lonarr(s(1)+posun,s(2))
obrazek2(posun:*,*)=obrazek
endif else begin ;posun doleva
obrazek2=shift(obrazek,posun)
endelse

s=size(obrazek2)
prvni1=widget_base(prvni,/column)
prvnilab=widget_label(prvni1,value='Shifted image')
prvnidraw=widget_draw(prvni1,xsize=s(1),ysize=s(2),x_scroll_size=400,y_scroll_size=400)
prvnibut=widget_button(prvni1,value='SAVE',uvalue='SAVE')
widget_control,prvnidraw,get_value=okno & wset,okno
tvscl,obrazek2

endif
endif
end

'NEXT':begin
widget_control,prvnilab,set_value='2nd profile'
widget_control,prvnibut,/destroy
widget_control,prvnidraw,set_uvalue='CLICK2'
widget_control,prvniposlab,set_value='position'
widget_control,prvnidraw,get_value=okno & wset,okno
plot,smooth(profil2,2),pos=[0,0,1,1],xstyle=4,ystyle=4
pocet2=0
print,'--------------'

end

'SAVE':begin
;widget_control,prvnidraw,get_value=okno & wset,okno
strput,cesta,'_s.ft',strlen(cesta)-4
cesta=cesta+'ts'
writefits,cesta,obrazek2 ;tvrd()
widget_control,event.top,/destroy
;****************************************************
end

endcase
end

pro foto2_event,event
COMMON inner1,base2,kliknuti1,kliknuti2
COMMON data,kliny,obrazek
COMMON profily,profil1,profil2
COMMON baze1,ibase,base,f2
COMMON baze2,prvnilab,prvnidraw,prvnibut,prvni1,prvni
COMMON pozice,pozice1,pozice2
COMMON uloz,cesta
COMMON pozice2,pos1,pos2
COMMON novy,prvniposlab
;COMMON velikost,s1,s2
widget_control,event.id,get_uvalue=uvalue
case uvalue of
'LOAD1':begin
widget_control,ibase,/destroy
cestaklin=pickfile()
pripona=strmid(cestaklin,strlen(cestaklin)-3,3)
if pripona eq 'fts' then kliny=readfits(cestaklin)
if pripona eq 'tif' then kliny=read_tiff(cestaklin,/unsigned)
;kliny=read_tiff(pickfile(filter='*.tif'),/unsigned)
s1=size(kliny)
ibase=widget_base(base,/row)
base1=widget_base(ibase,/column)
lab1=widget_label(base1,value='WEDGES')
draw1=widget_draw(base1,xsize=s1(1),ysize=s1(2),x_scroll_size=400,$
y_scroll_size=400,/button_events,uvalue='OBR1',/motion_events)
pos1=widget_label(base1,value=' ',xsize=150)
base2=widget_base(ibase)
widget_control,draw1,get_value=okno & wset,okno
tvscl,kliny
widget_control,f2,sensitive=1
end
'LOAD2':begin
widget_control,base2,/destroy

cesta=pickfile()
pripona=strmid(cesta,strlen(cesta)-3,3)
if pripona eq 'fts' then obrazek=readfits(cesta)
if pripona eq 'tif' then obrazek=read_tiff(cesta,/unsigned)
s2=size(obrazek)
base2=widget_base(ibase,/column)
lab2=widget_label(base2,value='IMAGE')
draw2=widget_draw(base2,xsize=s2(1),ysize=s2(2),x_scroll_size=400,$
y_scroll_size=400,/button_events,uvalue='OBR2',/motion_events)
pos2=widget_label(base2,value=' ',xsize=150)
shiftbut=widget_button(base2,value='COMPUTE SHIFT',uvalue='SHIFT')
widget_control,draw2,get_value=okno & wset,okno
tvscl,obrazek
kliknuti1=0 & kliknuti2=0
end

'OBR1':begin
position=strcompress('('+string(event.x)+','+string(event.y)+')',/remove_all)
widget_control,pos1,set_value=position

if (event.press eq 1) then begin
s=size(kliny)
profil1=make_array(s(1),value=0)
for i=-5,5 do profil1=profil1+float(kliny(*,event.y+i))
profil1=profil1/11
kliknuti1=1
endif
end

'OBR2':begin
position=strcompress('('+string(event.x)+','+string(event.y)+')',/remove_all)
widget_control,pos2,set_value=position

if (event.press eq 1) then begin
s=size(obrazek)
profil2=make_array(s(1),value=0)
for i=-5,5 do profil2=profil2+float(obrazek(*,event.y+i))
profil2=profil2/11
kliknuti2=1
endif
end


'SHIFT':begin
if (kliknuti1 ne 0) and (kliknuti2 ne 0) then begin
prvni=widget_base()
prvni1=widget_base(prvni,/column)
prvnilab=widget_label(prvni1,value='1st profile')
prvnidraw=widget_draw(prvni1,xsize=600,ysize=400,/button_event,uvalue='CLICK1',/motion_events)
prvniposlab=widget_label(prvni1,value='position',xsize=150)
prvnibut=widget_button(prvni1,value='Next profile',uvalue='NEXT')
widget_control,prvni,/realize
widget_control,prvnidraw,get_value=okno & wset,okno
plot,smooth(profil1,2),pos=[0,0,1,1],xstyle=4,ystyle=4
pozice1=make_array(10,value=0)
pozice2=make_array(10,value=0)
xmanager,'prvni',prvni
endif else begin
chyba=widget_base()
chlab=widget_label(chyba,value='You are supposed to take profiles by clicking in images')
widget_control,chyba,/realize
endelse
end
'QUIT':begin
widget_control,event.top,/destroy
end
endcase
end


pro foto2
COMMON baze1,ibase,base,f2
base=widget_base(mbar=bar,xsize=850,ysize=500)
ibase=widget_base(base)
fmenu=widget_button(bar,value='FILE',/menu)
f1=widget_button(fmenu,value='LOAD WEDGES ',uvalue='LOAD1')
f2=widget_button(fmenu,value='LOAD IMAGE',uvalue='LOAD2')
f3=widget_button(fmenu,value='QUIT',uvalue='QUIT')
widget_control,base,/realize
widget_control,f2,sensitive=0
xmanager,'foto2',base
end