;;;;make overlay map,integral intensity
;;;;you have to change 1,data dir   (defaul;../data/)
;;;;;;;;;;;;;;;;;;;;;;;2,output dir (defaul;../fig/multi/png/)
;;;;;;;;;;;;;;;;;;;;;;;3,map xrange,yrange ;;must
;;;;;;;;;;;;;;;;;;;;;;;4,corelation between 17G-34G


;;;windowsize
xsize=768
ysize=768
window,xsize=xsize,ysize=ysize


;;;make output dir
outdir='../figs/'
sortdir=outdir+'multi/'
pngdir=sortdir+'png/'
;case !d.name of
;   'Z' : device,set_resolution=[xsize,ysize]
;   'PS' : device,filename=filebd,/enc,xs=15.,ys=5.4,/color
;   else : window,1,xs=xsize,ys=ysize
;endcase
if (not file_test(outdir,/directory)) then file_mkdir,outdir
if (not file_test(sortdir,/directory)) then file_mkdir,sortdir
if (not file_test(pngdir,/directory)) then file_mkdir,pngdir


;;;set env
!p.font=1
!p.thick=2
!p.charsize=1.5
!p.charthick=2
!x.thick=2
!y.thick=2
!p.multi=[0,2,2,0,1]


;;;data dir
file17=findfile('../data/ipa*')
file34=findfile('../data/ipz*')

;;;file number
nt17=n_elements(file17[*])
nt34=n_elements(file34[*])

if(nt17 gt nt34) then nt0=nt34
if(nt17 lt nt34) then nt0=nt17

;;;read file
norh_rd_img,file17,ind17,dat17
norh_rd_img,file34,ind34,dat34


;;;intensity max
intmax17=max(dat17[*,*,*]) 
intmax34=max(dat34[*,*,*]) 


;;;;;;;;;;;auto corelation;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;for 17G
nx=n_elements(dat17[*,0,0])
ny=n_elements(dat17[0,*,0])
data_temp=fltarr(nx,ny,2)
dx_tot=0.d0
dy_tot=0.d0

for mm=1,nt17-1 do begin
data_temp[*,*,0]=dat17[*,*,mm-1]
data_temp[*,*,1]=dat17[*,*,mm]
xy_temp=get_correl_offsets(data_temp)
dx_tot=xy_temp[0,1]
dy_tot=xy_temp[1,1]
dat17[*,*,mm]=shift_img(dat17[*,*,mm],[dx_tot,dy_tot])
endfor

;;;for 34G
nx=n_elements(dat34[*,0,0])
ny=n_elements(dat34[0,*,0])
data_temp=fltarr(nx,ny,2)
dx_tot=0.d0
dy_tot=0.d0


for mm=1,nt34-1 do begin
data_temp[*,*,0]=dat34[*,*,mm-1]
data_temp[*,*,1]=dat34[*,*,mm]
xy_temp=get_correl_offsets(data_temp)
dx_tot=xy_temp[0,1]
dy_tot=xy_temp[1,1]
dat34[*,*,mm]=shift_img(dat34[*,*,mm],[dx_tot,dy_tot])
endfor

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


;;;make map
norh_index2map,ind17,dat17,nmap
norh_index2map,ind34,dat34,zmap


;;;correlation between 17G_34G
dx_34=-20
dy_34=-18
zmap[*].xc=zmap[*].xc+dx_34
zmap[*].yc=zmap[*].yc+dy_34


;;;rotaion map
foot_angle=-13
rnmap=rot_map(nmap,foot_angle)
rzmap=rot_map(zmap,foot_angle)




;;;;;;;;;;;;;;;;;;;;;;;;;;;plot image;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
loadct,5
;;;set xrange,yrange by arcsec
xr_min=870
xr_max=1020
yr_min=-200
yr_max=-50

for mm=0,nt0-1 do begin
;;;map
plot_map,rnmap[mm],/logscale,dmin=8e3,$
		xrange=[xr_min,xr_max],yrange=[yr_min,yr_max]
mapmax34=max(zmap[mm].data[*,*])
plot_map,rzmap[mm],/logscale,dmin=8e3,/over,$
	leveles=[mapmax34*0.1,mapmax34*0.5,mapmax34*0.9]


;;;integrate
ntx17=n_elements(rnmap[mm].data[*,0])
nty17=n_elements(rnmap[mm].data[0,*])
ntx34=n_elements(rzmap[mm].data[*,0])
nty34=n_elements(rzmap[mm].data[0,*])
intx17=fltarr(ntx17)
inty17=fltarr(nty17)
intx34=fltarr(ntx34)
inty34=fltarr(nty34)
for i=0,ntx17-1 do begin
intx17(i)=total(rnmap[mm].data[i,*])
inty17(i)=total(rnmap[mm].data[*,i])
endfor
for i=0,ntx34-1 do begin
intx34(i)=total(rzmap[mm].data[i,*])
inty34(i)=total(rzmap[mm].data[*,i])
endfor



xx17=findgen(ntx17)*rnmap[mm].dx+rnmap[mm].xc-rnmap[mm].dx*ntx17*0.5
yy17=findgen(nty17)*rnmap[mm].dy+rnmap[mm].yc-rnmap[mm].dx*nty17*0.5
xx34=findgen(ntx34)*rzmap[mm].dx+rzmap[mm].xc-rzmap[mm].dx*ntx34*0.5
yy34=findgen(nty34)*rzmap[mm].dy+rzmap[mm].yc-rzmap[mm].dx*nty34*0.5

plot,xx17,intx17,xrange=[xr_min,xr_max];,xtitle='17G'
oplot,xx34,intx34,linestyle=2
plot,inty17,yy17,yrange=[yr_min,yr_max]
oplot,inty34,yy34,linestyle=2

;;;fixed colorbar
plot_map,rnmap[mm],/logscale,dmin=1e5,dmax=intmax17,$
	xrange=[xr_min,xr_max],yrange=[yr_min,yr_max]
plot_map,rzmap[mm],/logscale,dmin=1e5,dmax=intmax34,/over,$
	leveles=[intmax34*0.1,intmax34*0.5,intmax34*0.9]
;	xrange=[xr_min,xr_max],yrange=[yr_min,yr_max]


;;;output png
filepng=pngdir+'idl'+string(mm,format='(i4.4)')+'.png'
img=tvrd(0,0,xsize,ysize,1,true=1)
tvlct,red,green,blue,/get
write_png,filepng,img,red,green,blue

endfor


end