Python savestatsの例

コード例 #1

ファイル: dirty_selfcal_compare.py プロジェクト: nonyt/reduction

from spectral_cube import SpectralCube
from astropy.stats import mad_std
from astropy import log
import pylab as pl

from before_after_selfcal_quicklooks import make_comparison_image, get_selfcal_number

cwd = os.getcwd()
basepath = '/orange/adamginsburg/ALMA_IMF/2017.1.01355.L/imaging_results'
os.chdir(basepath)
prepostpath = '/orange/adamginsburg/web/secure/ALMA-IMF/October2020/comparisons/prepost/'
octoberpath = '/orange/adamginsburg/web/secure/ALMA-IMF/October2020/'

import imstats

tbl = imstats.savestats(basepath=octoberpath)

#tbl = Table.read('/orange/adamginsburg/web/secure/ALMA-IMF/October2020/metadata.ecsv')
tbl.add_column(
    Column(name='casaversion_pre', data=['             '] * len(tbl)))
tbl.add_column(
    Column(name='casaversion_post', data=['             '] * len(tbl)))
tbl.add_column(Column(name='scMaxDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMinDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMADDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMeanDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMedianDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_post', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='min_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='min_post', data=[np.nan] * len(tbl)))

コード例 #2

ファイル: 7m12m_comparison.py プロジェクト: xinglunju/reduction

from astropy import visualization
from astropy.table import Table, Column
from spectral_cube import SpectralCube
from astropy.stats import mad_std
from astropy import log
import pylab as pl

from compare_images import make_comparison_image

cwd = os.getcwd()
basepath = '/bio/web/secure/adamginsburg/ALMA-IMF/Feb2020'
os.chdir(basepath)

import imstats

tbl = imstats.savestats(basepath=basepath)

#tbl = Table.read('/bio/web/secure/adamginsburg/ALMA-IMF/Feb2020/metadata.ecsv')
tbl.add_column(
    Column(name='casaversion_7m12m', data=['             '] * len(tbl)))
tbl.add_column(
    Column(name='casaversion_cleanest', data=['             '] * len(tbl)))
tbl.add_column(Column(name='7m12mMaxDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='7m12mMinDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='7m12mMADDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='7m12mMeanDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='7m12mMedianDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='7m12mSumDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='shape', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='ppbeam', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_7m12m', data=[np.nan] * len(tbl)))

コード例 #3

from spectral_cube import SpectralCube
from astropy.stats import mad_std
from astropy import log
import pylab as pl

from before_after_selfcal_quicklooks import make_comparison_image, get_selfcal_number

cwd = os.getcwd()
basepath = '/orange/adamginsburg/ALMA_IMF/2017.1.01355.L/imaging_results'
os.chdir(basepath)
prepostpath = '/bio/web/secure/adamginsburg/ALMA-IMF/Feb2020/comparisons/prepost/'
febpath = '/bio/web/secure/adamginsburg/ALMA-IMF/Feb2020/'

import imstats

tbl = imstats.savestats(basepath=febpath)

#tbl = Table.read('/bio/web/secure/adamginsburg/ALMA-IMF/Feb2020/metadata.ecsv')
tbl.add_column(
    Column(name='casaversion_pre', data=['             '] * len(tbl)))
tbl.add_column(
    Column(name='casaversion_post', data=['             '] * len(tbl)))
tbl.add_column(Column(name='scMaxDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMinDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMADDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMeanDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMedianDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_post', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='min_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='min_post', data=[np.nan] * len(tbl)))

コード例 #4

from spectral_cube import SpectralCube
from astropy.stats import mad_std
from astropy import log
import pylab as pl

from before_after_selfcal_quicklooks import make_comparison_image, get_selfcal_number

cwd = os.getcwd()
basepath = '/orange/adamginsburg/ALMA_IMF/2017.1.01355.L/imaging_results'
os.chdir(basepath)
prepostpath = '/orange/adamginsburg/web/secure/ALMA-IMF/May2021/comparisons/prepost/'
Maypath = '/orange/adamginsburg/web/secure/ALMA-IMF/May2021/'

import imstats

tbl = imstats.savestats(basepath=Maypath)

tbl.add_column(Column(name='casaversion_pre', data=['             ']*len(tbl)))
tbl.add_column(Column(name='casaversion_post', data=['             ']*len(tbl)))
tbl.add_column(Column(name='scMaxDiff', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='scMinDiff', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='scMADDiff', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='scMeanDiff', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='scMedianDiff', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='dr_pre', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='dr_post', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='min_pre', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='min_post', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='max_pre', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='max_post', data=[np.nan]*len(tbl)))
tbl.add_column(Column(name='mad_pre', data=[np.nan]*len(tbl)))

コード例 #5

ファイル: before_after_selfcal_quicklooks_October31_2019_release.py プロジェクト: ctimea/reduction

from astropy.io import fits
from astropy import visualization
from astropy.table import Table, Column
from spectral_cube import SpectralCube
from astropy.stats import mad_std
from astropy import log
import pylab as pl

from before_after_selfcal_quicklooks import make_comparison_image, get_selfcal_number

cwd = os.getcwd()
os.chdir('/bio/web/secure/adamginsburg/ALMA-IMF/October31Release')

import imstats

tbl = imstats.savestats()

#tbl = Table.read('/bio/web/secure/adamginsburg/ALMA-IMF/October31Release/metadata.ecsv')
tbl.add_column(Column(name='scMaxDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMinDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMADDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMeanDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='scMedianDiff', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_post', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='max_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='max_post', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='mad_pre', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='mad_post', data=[np.nan] * len(tbl)))
tbl.add_column(Column(name='dr_improvement', data=[np.nan] * len(tbl)))