ima=read_tiff('PIA06175_ring_panorama.tif')
image=reform(ima(2,5789-2800:*,*))
yscale=1./sin(4./!radeg)

image2=congrid(image,2800,1000*yscale)

image2=image2(*,0:13999)
nwin,xs=1400,ys=1050


;IMAGE2           BYTE      = Array[2800, 14000]


before_ansa=rebin(image2(0:799,13999-1999:13999),400,1000)
dx=150
after_ansa=rebin(image2(dx+0:dx+799,0:1999),400,1000)

dy2=-600
dx2=2000
ansa=rebin(image2(dx2+0:dx2+799,dy2+7000:dy2+7999),400,1000)

ansa(0,0)=0
ansa(1,0)=255

before_ansa(0,0)=0
before_ansa(1,0)=255
after_ansa(0,0)=0
after_ansa(1,0)=255


mini=100
maxi=180
tvscl,rebin(image,200,1000)<maxi>mini,0,0
tvscl,before_ansa>mini<maxi,201,0
tvscl,ansa>mini<maxi,601,0
tvscl,after_ansa>mini<maxi,1001,0

xyouts,0.5,0.97,/nor,"PIA06175.tif   12.12.2004  B=-4, B'=-23.2  alpha=42  (theta_sun=-40) scale=10.5km/pixel",chars=2.5,ali=.5


;------------------------------------------------------
;this is another Cassini image
;"RINGS FROM AFAR" = 572_1170_3.tif  

;just demonstrating how different the resonances can look depending
;on geometry!

imageb=read_tiff('572_1170_3.tif')
imageb=imageb(0:549,360:826)
yscale=1./sin(13.5/!radeg)

imageb2=congrid(imageb,550,yscale*467,/cubic)


nwin,xs=1000,ys=1000
tvscl,rebin(imageb2,275,1000)<140>100

imageb3=imageb2(*,500:1499)
tvscl,imageb3,400,0


ots," 572_1170_3.tif   1.11.2004  B=-13.5, B'=-23.5  alpha=112 !c!c  (theta_sun=-120) scale=129km/pixel" 




rmax=140.
rmin=60.
show_iss,[rmin,rmax],xs=600,ys=50,lit=lit,/noplot
!p.position=[0.5,0.03,.97,0.1]
!p.position=[0.472,0.50,.959,0.53]
!p.position=[0.472,0.03,.959,0.1]

xx=rmin+findgen(600)/600.*(rmax-rmin)
tvplot2,sigrange(lit),xx=xx,/nonew,/noerase
!p.position=0


end