zgyin的个人博客分享 http://blog.sciencenet.cn/u/zgyin

博文

画扇形图(idl程序)

已有 4824 次阅读 2013-4-13 14:23 |个人分类:idl使用方法|系统分类:科研笔记


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




https://wap.sciencenet.cn/blog-783377-679814.html

上一篇:Jy/beam到Jy/pixel的转化
下一篇:坐标轴是对数形式的时候均匀叠加上一些点
收藏 IP: 202.101.72.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-2-4 13:39

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部