||
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
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-19 22:55
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社