Changeset 773 for branches/Release12


Ignore:
Timestamp:
12/06/05 14:03:55 (19 years ago)
Author:
mar637
Message:

Request: added recalc_azel

Location:
branches/Release12
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/Release12/python/__init__.py

    r737 r773  
    268268
    269269__date__ = '$Date$'.split()[1]
    270 __version__  = '1.2'
     270__version__  = '1.2.1'
    271271
    272272if rcParams['verbose']:
     
    340340            poly_baseline   - fit a polynomial baseline to all Beams/IFs/Pols
    341341            auto_poly_baseline - automatically fit a polynomial baseline
     342            recalc_azel     - recalculate azimuth and elevation based on
     343                              the pointing
    342344            gain_el         - apply gain-elevation correction
    343345            opacity         - apply opacity correction
  • branches/Release12/python/scantable.py

    r760 r773  
    13721372            else: raise
    13731373
     1374    def recalc_azel(self):
     1375        varlist = vars()
     1376        """
     1377        Recalculate the azimuth and elevation for each position.
     1378        Parameters:
     1379            none
     1380        Example:
     1381        """
     1382        self._recalc_azel()
     1383        self._add_history("recalc_azel", varlist)
     1384        return
     1385
    13741386    def __add__(self, other):
    13751387        varlist = vars()
     
    13861398        print_log()
    13871399        return s
    1388 
    13891400
    13901401    def __sub__(self, other):
  • branches/Release12/src/SDMemTable.cc

    r754 r773  
    16311631    mdr = MDirection::J2000;
    16321632    pushLog("WARNING: Unknown equinox using J2000");
    1633 
    16341633  }
    16351634
     
    18371836}
    18381837
     1838void SDMemTable::calculateAZEL()
     1839{
     1840  MPosition mp = getAntennaPosition();
     1841
     1842  for (uInt i=0; i<nRow();++i) {
     1843    MEpoch me = getEpoch(i);
     1844    MDirection md = getDirection(i,False);
     1845    MeasFrame frame(mp, me);
     1846    Vector<Double> azel =
     1847        MDirection::Convert(md, MDirection::Ref(MDirection::AZEL,
     1848                                                frame)
     1849                            )().getAngle("rad").getValue();
     1850    azCol_.put(i,azel[0]);
     1851    elCol_.put(i,azel[1]);
     1852  }
     1853}
  • branches/Release12/src/SDMemTable.h

    r745 r773  
    265265              const casa::Vector<casa::Int>& c);
    266266
     267
     268  void calculateAZEL();
    267269
    268270private:
  • branches/Release12/src/SDMemTableWrapper.h

    r745 r773  
    217217  }
    218218
     219  void calculateAZEL() { table_->calculateAZEL(); };
    219220
    220221private:
  • branches/Release12/src/python_SDMemTable.cc

    r745 r773  
    108108    .def("_addfit", &SDMemTableWrapper::addFit)
    109109    .def("_getfit", &SDMemTableWrapper::getSDFitTable)
     110    .def("_recalc_azel", &SDMemTableWrapper::calculateAZEL)
    110111  ;
    111112};
Note: See TracChangeset for help on using the changeset viewer.