PRO plotaz_el,file

rdkimg,file,hb,img
col=REBIN(LINDGEN(256),256,256)
row=TRANSPOSE(col)

c0=130.302
r0=126.375
r=199.895
k=0.638308
ca=0.

l=256.

dc=col-REPLICATE(c0,l,l)
dr=row-REPLICATE(r0,l,l)
rho=SQRT(dc^2+dr^2)
phi=ATAN(dc,-dr)
xp=COS(phi)
yp=SIN(phi)
cc=COS(ca)
sc=SIN(ca)
x=xp*cc+yp*sc
y=-xp*sc+yp*cc
el=90.-!RADEG*(ASIN(rho/r)/k)
az=!RADEG*ATAN(-y,x)
wn=WHERE(az LT 0.)
az(wn)=az(wn)+360.
angx=!RADEG*ATAN(x,-y)

p=where(img GT max(img)/2.0)
one2two,p,img,cp,rp
si=img(min(cp,max=cm):cm,min(rp,max=rm):rm)
TRIANGULATE,az(p),el(p),triangles
ti=TRIGRID(az(p),el(p),img(p),triangles,[.25,.25],[180.,10.,240.,60.])
MAP_SET,35.,215.,0.,/ORTHOGRAPHIC,LIMIT=[15.,185.,55.,245.]
mi=MAP_IMAGE(ti,sx,sy,LATMIN=15.,LATMAX=55.,LONMIN=185.,LONMAX=245. $
	    ,COMPRESS=4)
bmi=BYTSCL(mi,MIN=0,TOP=254)
IF !D.NAME EQ 'PS' THEN BEGIN
    bmi(WHERE(bmi EQ 0B))=255B
    TVLCT,255B,255B,255B,!D.N_COLORS,/RGB
ENDIF
TVSCL,mi,sx,sy
MAP_GRID,LATDEL=5,LONDEL=5,/LABEL
stop
RETURN
END