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=!PI/2.-(ASIN(rho/r)/k);<(89.5*!DTOR)
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),!RADEG*el(p),triangles
ti=TRIGRID(az(p),!RADEG*el(p),img(p),triangles,[.25,.25],[180.,0.,270.,70.])
MAP_SET,30.,200.,0.,/ORTHOGRAPHIC,LIMIT=[15.,175.,45.,225.]
mi=MAP_IMAGE(ti,sx,sy,LATMIN=15.,LATMAX=45.,LONMIN=175.,LONMAX=225. $
	    ,COMPRESS=4)
TVSCL,mi,sx,sy
MAP_GRID,LATDEL=5,LONDEL=5,/LABEL
stop
RETURN
END