IDL画扇形图还是有些麻烦的，今天中午没午休，以红移和RA为例写了程序，希望以后能用上

pro sector

set_plot,'ps'

device,file='F:Aprilmap.ps'

READCOL, 'F:Aprilalfa.txt', FORMAT='(l,F,F,f)', index,ra, dec,redshift,/SILENT

radius=redshift*1000.0

angle=(ra-360.0)/360.0*!pi*2.0

Plot,radius,angle,PSym=3,symsize=0.8,xrange=[-20,80],yrange=[-15.0*!pi,15.0*!pi],/POLAR

a=[0.0,max(radius)]

b=[0.0,max(angle)]

OPLOT,a,b,/POLAR

c=[0.0,max(radius)]

d=[0.0,min(angle)]

OPLOT,c,d,/POLAR

e1=(fltarr(1000)+1.0)*max(radius)

e2=min(angle)+(max(angle)-min(angle))*randomu(10,1000)

oplot,e1,e2,PSym=3,/POLAR

e1=(fltarr(1000)+1.0)*max(radius)*1.0/3.0

e2=min(angle)+(max(angle)-min(angle))*randomu(10,1000)

oplot,e1,e2,PSym=3,/POLAR;画第一条等高线

e1=(fltarr(1000)+1.0)*max(radius)*2.0/3.0;画第二条等高线

e2=min(angle)+(max(angle)-min(angle))*randomu(10,1000)

oplot,e1,e2,PSym=3,/POLAR

print,max(radius)

print,max(radius)/3.0

print,max(radius)/3.0*2.0

xyouts,14,-12,'19.6',charsize=1,charthick=1,orientation=-20;画等高线19.6

xyouts,30.5,-21.5,'39.3',charsize=1,charthick=1,orientation=-20;画等高线39.3

;xyouts,48,-31,'58.9',charsize=1,charthick=1,orientation=-20

print,max(ra)-360

print,min(ra)-360

xyouts,50,-38,'-30',charsize=1,charthick=1,orientation=50;标记RA=-30

xyouts,38,48,'48',charsize=1,charthick=1,orientation=-40;标记RA=48

xyouts,62,2,'RA',charsize=1.5,charthick=2,orientation=90;标记RA

xyouts,20,-20,'z*1000.0',charsize=1,charthick=1,orientation=-20;标记z*1000.0

device,/close

set_plot,'win'

end

转载本文请联系原作者获取授权，同时请注明本文来自张国印科学网博客。

链接地址：http://blog.sciencenet.cn/blog-783377-679814.html

上一篇：Jy/beam到Jy/pixel的转化

下一篇：坐标轴是对数形式的时候均匀叠加上一些点