[IDL] 畫假星星程式

又是IDL文XD 沒辦法，最近在整理程式嘛
希望能把雜亂的程式檔們減肥，讓電腦真正空一些~

嚴格說來這個程式是傅老寫的(驚)，是在meeting(還是上課?)時口述要我們寫下的
總之就是write_bmp這個碼沒什麼特別的
但很實用，有它立刻就做出三維平面了~





pro test

n=200
xa=100.
ya=100.
sa=15.
ia=250.

xb=80.
yb=80.
sb=15.
ib=150.

i12=fltarr(n,n)

for i=0,n-1 do begin
  for j=0,n-1 do begin
    d1=(i-xa)^2+(j-ya)^2
    d2=(i-xb)^2+(j-yb)^2
    i1=ia*exp(-(d1/sa^2))
    i2=ib*exp(-(d2/sb^2))
    i12(i,j)=i1+i2
  endfor
endfor

surface,i12
write_bmp,'d:\write_bmp\ixy.bmp',i12
ixy=read_bmp('d:\write_bmp\ixy.bmp')

i12_fft=abs(1+shift(fft(ixy),100,100))*10000
write_bmp,'d:\write_bmp\ifft.bmp',i12_fft

end