천문학적 데이터 분석 입문 2 (Montage의 사용 방법)

3354 단어 AstronomyPython
이 글은 Montage의 사용 방법을 소개한다.

설치 모듈


http://montage.ipac.caltech.edu/docs/download.html
에서 기술한 장면은 다음과 같은 절차를 이용하여 명세표를 작성하여 개념 디자인에서 체량의 부피를 분석하도록 한다.
그리고 경로를 통해 터미널에서 어디서든 "mImgtbl"등의 명령을 사용할 수 있는 상태를 확인합니다.

데이터 준비


이번에는 Spitzer24 µm의 아카이브 데이터를 사용합니다.
https://irsa.ipac.caltech.edu/data/SPITZER/MIPSGAL/
W43 근처의 데이터를 다운로드합니다.
리본에서 를 클릭합니다.
https://irsa.ipac.caltech.edu/data/SPITZER/MIPSGAL/images/mosaics24/
에서 설명한 대로 해당 매개변수의 값을 수정합니다.
MG0290n005_024.fits
MG0290p005_024.fits
MG0300n005_024.fits
MG0300p005_024.fits
MG0310n005_024.fits
MG0310p005_024.fits
의 총 6개의fits를 연결합니다.

모듈 사용(2D)


우선,'imagedir'라는 디렉터리를 만들어서 모든fits를 넣습니다.
터미널에서fits가 있는 디렉터리로 이동
mkdir imagedir
mkdir resultdir
mv MG*_024.fits imagedir/
에서 설명한 대로 해당 매개변수의 값을 수정합니다.
그리고
mImgtbl imagedir images.tbl
에서 설명한 대로 해당 매개변수의 값을 수정합니다.images.tbl이라는 파일을 생성합니다.파일 목록.
리본에서 를 클릭합니다.
mMakeHdr images.tbl template.hdr GAL
에서 설명한 대로 해당 매개변수의 값을 수정합니다.template.hdr라는 파일을 생성합니다.새 머리글.뒤에 GAL을 붙이면 좌표계는 은경위가 될 수 있다.
리본에서 를 클릭합니다.
mProjExec -p imagedir images.tbl template.hdr resultdir stats.tbl
에서 설명한 대로 해당 매개변수의 값을 수정합니다.resultdir에서reprojectfits를 생성합니다.나는 계산할 시간이 좀 필요하다.
리본에서 를 클릭합니다.
mImgtbl resultdir/ resultimages.tbl
에서 설명한 대로 해당 매개변수의 값을 수정합니다.resultimages.tbl이라는 파일을 생성합니다.
리본에서 를 클릭합니다.
mAdd -p resultdir/ resultimages.tbl template.hdr W43_24um.fits
에서 설명한 대로 해당 매개변수의 값을 수정합니다.resultdir의fits를 추가합니다.W43_24um.fits는 저장할 이름입니다.나는 계산할 시간이 좀 필요하다.W43_24um_area.fits도 생성되지만 중간 파일과 같아서 삭제할 수 있습니다.
생성된fits 값은float64로 쓰입니다.python의 astropy.io.피츠 등으로float32를 열면 반쯤 가벼워집니다.
마지막으로 ds9에서 생성된fits를 검사합니다.

Montage 사용(3D)


내 생각에도 3D의fits(cube)를 서로 덧붙이고 싶은 일이 있다.
Montage도 cube를 지원하지만 속도축이 때때로 이상한 행동을 할 수 있기 때문에 가능한 casa의'imregrid'등 속도축을 맞추고 덧붙이는 것을 권장합니다.
카사를 어떻게 사용하는지에 대한 정보는 다른 문장을 참고하십시오.
mImgtbl imagedir/ images.tbl 

mMakeHdr images.tbl template.hdr

for infile in imagedir/*.fits; do
mProjectCube $infile resultdir/${infile##*/}_proj.fits template.hdr ;
done

mImgtbl resultdir/ resultimages.tbl

mAddCube -p resultdir resultimages.tbl template.hdr result.fits
여기서 마치겠습니다.
링크
카탈로그

좋은 웹페이지 즐겨찾기