function rh_ronron,ima,n,ma,xc,yc,larg
;
;calcul des i,j sur lesquels interpoler les points pour tv
; ces pts sont ceux pour lesquels les x,y sont reguliers entre
; -larg/2. et +larg/2.

; version acceleree (une seule boucle de remplissage des x,y)
; le 7/1/+98 par ak



xi=fltarr(n,n)
yi=fltarr(n,n)
imama=fltarr(n,n)
;larg=4. ; champ affichˇ

det=ma(0)*ma(3)-ma(2)*ma(1)
 
nnorm=larg/float(n); champ de 4 rayons

; a refaire: ils sont tous pareils par ligne et colonnes, les xi yi
; puis ss doute appeler interpolate avec l'option grid 

;for i=0 , n-1 do begin
;	ii=-larg/2.+float(i)*nnorm-xc
;	for j=0 , n-1 do begin
;		jj=-larg/2.+float(j)*nnorm-yc
;		xi(i,j)=ma(3)*ii-ma(2)*jj
;		yi(i,j)=ma(0)*jj-ma(1)*ii
;		print,i,j,yi(i,j),x(i,j),y(i,j)
;	endfor
;endfor

indic=findgen(n)
ii=-larg/2.-xc+indic*nnorm
jj=ii+xc-yc
for ik = 0,n-1 do begin
	xi(0:n-1,ik)=ma(3)*ii-ma(2)*jj(ik)
	yi(0:n-1,ik)=ma(0)*jj(ik)-ma(1)*ii
endfor



nn2=float(n)/2.
; on calcule avec 0,0 a nn2,nn2, plus le mvt du centre du soleil xc,yc
;print,'maille xc,yc',xc,yc

xi=xi/det+nn2
yi=yi/det+nn2

; interpolation avec raz des pts en dehors du parallelogramme xi,yi

imama=interpolate(ima,xi,yi,miss=0)

return,imama
end