function rot_optimize, x, y, number=number, range=range

widget_control, /hourglass

if n_elements(number) le 0 then number=10.

if n_elements(range) le 0 then range=10.

tot=fltarr(number)

index=-range/2.+findgen(number)/(number-1)*range

Szx=size(x)

Radius=(Szx(1) < Szx(2))*0.95/2

N=512

t=findgen(N)/(N-1)*!pi*2

xcir=Szx(1)/2.+Radius*cos(t)
ycir=Szx(2)/2.+Radius*sin(t)

disk_index=polyfillv(xcir, ycir, Szx(1), Szx(2))

for j=0, number-1 do tot(j)=total(x(disk_index)*(rot(y, index(j), /int, miss=0))(disk_index) )

Nsp=10

index1=-range/2.+findgen(number*Nsp)/(number*Nsp-1)*range

tot1=spline(index, tot-min(tot), index1)

amax=max(tot1, imax)

return, index1(imax)

	end