게시판 답변이 생성됐습니다
-
글
-
답변 감사드립니다. 추가로 질문이 있는데요. for 구문 위에 window를 사용해서 하나의 창에서 작업하는 것은 수정하였으며, for 구문의 마지막에 w.erase를 넣어서 다음 i번째에서 그래프가 중복 전시되는 걸 막는 것까지는 했지만, 마우스 버튼을 눌러야 다음 i번째로 진행되게 할 수는 없을까요? 현재는 그냥 계속 순서대로 데이터를 처리하고 있는데 이를 마우스 클릭이 있어야만 다음으로 넘어가게 만들고 싶습니다. 또한 i가 대략 200여개 정도 되는데 for 구문동안 출력되는 모든 그래프를 개별적으로 저장하고 싶습니다. 말씀하신 것처럼 w1.save, ‘파일이름’으로 코딩하면 for에서 가장 마지막 파일만 저장이 됩니다.
우선 두번째에 대해 개인적으로 작성하여 봤는데, 정확하게 코딩이 된건지는 잘 모르겠습니다.
다음과 같이 하면 될까요?STT_RA = range_data.left_RA[0] ; 적경 시작점. 66.00000
STT_DEC = range_data.left_DEC[0] ; 적위 시작점. -55.410000
for i=0,101 do begin
xbin = 0.1d
ybin = 0.04d
w2=where(data.ra_des gt STT_RA and $
data.ra_Des lt STT_RA+xbin*(i+1))
data_RA = data[w2]
for j=0,101 do begin
w3=where(data_RA.DEC_des gt STT_DEC and $
data_RA.DEC_Des lt STT_DEC+ybin*(j+1))
data_RADEC = data_RA[w3]
density = n_elements(data_RADEC.sz)
print, density
endfor
endfor답변 정말 감사합니다. 먼저 첫번째로 드렸던 질문은 모든 경우에 대한(16912) 콘솔창에 출력되는 p 값을 하나의 텍스트 파일에 저장한 뒤, 텍스트 파일에 기록되어 있는 값에서 특정 조건을 만족하는 값을 찾는 것에 대하여 질문을 드린것입니다. 여기서 특정 조건은 p값 중 표준편차 시그마에 해당하는 값의 범위를 0.1에서 0.2 사이로 두는 작업을 하려고 한 것이었습니다.
또한 두번째는 현재 RA와 DEC에 따른 천체의 분포를 scatterplot으로 나타내었습니다. 문제는 ra,dec에 따른 천체의 위치만 알고 있지 단위면적당 밀도값이 주어지지 않았기 때문에 단위면적을 적경,적위를 기준으로 생성하고 (가로축 100, 세로축 100개로 총 10,000개의 격자를 생성하려고 하고 있습니다.) 각각의 격자의 RA,DEC 범위 안에 포함되는 데이터수를 where문으로 찾아(RA gt 74 and RA lt 75와 같이) 밀도, ra, dec의 contour map를 만드려고 하고 있습니다. 그런데 격자를 만드는게 쉽지 않아서 어떻게 하면 만들 수 있을지 질문을 드렸었습니다. 학생이어서 부족한 점이 너무 많음에도 불구하고 상세하게 답변해주셔서 정말 감사드립니다.
