makemovie_gaze.pro 1.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104
  1. pro makemovie_gaze
  2. default, filename, "$HOME/prog/objsim/data/movies/gazemovie.idlmov"
  3. retw=50
  4. reth=1
  5. modw=50
  6. modh=1
  7. a = obj_new('movieinput', /create, filename=filename)
  8. a->generatefile, owarray=[retw,modw], oharray=[reth,modh], nfilters=2
  9. pic = fltarr(retw,reth)
  10. pic2 = fltarr(modw,modh)
  11. nretneurons=retw*reth
  12. nmodneurons=modw*modh
  13. nkombi=nretneurons*nmodneurons
  14. g=indgen(7)-3
  15. b=fltarr(7)
  16. for i=0,6 do begin
  17. b(i)=exp(-g(i)*g(i)/5.)
  18. endfor
  19. h=indgen(11)-5
  20. c=fltarr(11)
  21. for i=0,10 do begin
  22. c(i)=exp(-h(i)*h(i)/10.)
  23. endfor
  24. objposition=intarr(3)
  25. objposition(0)=0
  26. objposition(1)=20
  27. objposition(2)=40
  28. for objpos=0,2 do begin
  29. for modpos=0,nmodneurons-9 do begin
  30. ; for retpos=0,nretneurons-9 do begin
  31. rettemppos=35
  32. modtemppos=40
  33. pic(*)=0
  34. pic2(*)=0
  35. ;if (modpos-retpos eq objposition(objpos)) then begin
  36. for j=0,6 do begin
  37. pic(objposition(objpos)+j)=1*b(j)
  38. endfor
  39. for j=0,6 do begin
  40. pic2(modpos+j)=1*b(j)
  41. endfor
  42. a->WriteFrame, pic
  43. a->WriteFrame, pic2
  44. ;endif else begin
  45. ; pic(retpos)=0
  46. ; pic2(modpos)=0
  47. ;endelse
  48. tvscl, rebin(pic, 10*retw, reth*10, /sample)
  49. ; endfor
  50. endfor
  51. wait, 0.15
  52. endfor
  53. ;for u=0,nmodneurons-5 do begin
  54. ; for i=0, nretneurons-5 do begin
  55. ; frame = pic
  56. ; frame2 = pic2
  57. ; for j=0,4 do begin
  58. ; frame(i+j)=1*b(j)
  59. ; frame2(u+j)=1*b(j)
  60. ; ;;frame2(i+j)=1+b(j)
  61. ; endfor
  62. ; ;frame2(u)=1
  63. ; ;wait, 0.15
  64. ; tvscl, rebin(frame, 10*retw, reth*10, /sample)
  65. ; ;tvscl, 0, 0.3, /norm, rebin(frame2, 10*modw, modh*10, /sample)
  66. ; a->WriteFrame, frame
  67. ; a->WriteFrame, frame2
  68. ; endfor
  69. ;endfor
  70. ;endfor
  71. obj_destroy, a
  72. end