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的轉化

下一篇：坐標軸是對數形式的時候均勻疊加上一些點

總結

以上是生活随笔為你收集整理的hdc mfc 画扇形图_科学网—画扇形图（idl程序） - 张国印的博文的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。