특정성단에 대한 자료를 받아서 그 성단의 구성천체들에 대한 성도를 만들려고 합니다.
만드는 방법은 중심별을 정해놓고 이 별에 대한 상대위치로 다른 별들의 위치를 표시하고
각 별들의 밝기에 따라 별의 크기를 다르게 넣고 싶습니다.
현재까지 readcol을 이용해서 ascii파일을 불러와 별의 상대위치에따라 plot을 하도록 pro를 짰는데요
문제는 이 plot의 symbol size를 별의 밝기에 따라서 다르게 넣는 부분입니다.
pro starchart
readcol,’파일명’, A,B,C,D,E,F,G,H,I,J,K
DEC0=(26*60)+28(43.62/60)
Q0=((DEC0*0.016667)/360)*2*3.14
RA0=(20+(12.0/60)+(00.699/3600))*900*cos(Q0)
;중심별의 좌표
DECi=(E*60)+F(G/60)
Q=((DECi*0.016667)/360)*2*3.14
RAi=(B+(C/60)+(D/3600))*900*cos(Q)
;각 별의 적위 적경을 분단위로 변경
DEC=DECi-DEC0
RA=RAi-RA0
;중심별로부터 각 별의 상대 위치 지정
q=plot(RA,DEC)
현재 제가 짠 pro가 대충 이런 상태입니다.
여기서 별의 밝기 등급 SIZ=0.25*(17-K)의 값에 따라 plot된 점의 크기를 다르게 하고 싶은건데요
단순히 plot(RA,DEC,sym_size=SIZ)로 넣으면 RA,DEC변수값과 SIZ변수 사이의 관계를 IDL이 알수가 없는지
plot되지 않네요.
RA,DEC와 SIZ변수 사이의 관계를 지정할수 있는 방법이 있을런지요…
이걸 어떻게 해야 되면 좋을까요?
(참고로 ascii파일은 txt형태의 파일이고 BCD는 적경의 h:m:s값 EFG는 적위의 deg:’:”값이고 K가 밝기 등급관련 값입니다. 자료는 아래로 쭉 이어지는 열 형태 자료입니다)
-
이 게시글은 orion03에 의해 10 years, 7 months 전에 수정됐습니다.