source: branches/alma/src/MathUtils.h@ 1589

Last change on this file since 1589 was 1514, checked in by Kana Sugimoto, 16 years ago

New Development: No

JIRA Issue: Yes (CAS-1079)

Ready to Release: Yes

Interface Changes: Yes

What Interface Changed:
Added the new functions to get min/max position (channel) of spectra.

  1. std::vector<int> asap::python::stmath::_minmaxpos(ScantableWrapper& in,

std::vector<bool>& mask, std::string& which)

@python_STMath.cpp

  1. std::vector<int> STMathWrapper::minMaxPos(ScantableWrapper& in,

std::vector<bool>& mask, std::string& which)

@STMathWrapper.h

  1. std::vector<int> STMath::minMaxPos(CountedPtr< Scantable > & in,

std::vector< bool > & mask, std::string& which)

@STMath.h & .cpp

  1. IPosition mathutil::minMaxPos(String& which, MaskedArray<Float>& data) @MathUtils.h & .cpp

Test Programs:

Run scantable.stats() with the parameter stat='minpos' or 'maxpos',
and you'll get min/max value with its psition.

Put in Release Notes: No

Module(s): scantable.stats()

Description:

These modifications are to return min/max value with its
position (channel/frequency/velocity) by running scantable.stats().

Diagram:
scantable.stats() ->asap::python::stmath::_minmaxpos
-> STMathWrapper::minMaxPos -> STMath::minMaxPos() -> mathutil::minMaxPos()
-> casa::minMax (@casacore/casa/casa/Arrays/MaskArrMath.tcc)


  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.6 KB
Line 
1//#---------------------------------------------------------------------------
2//# MathUtilities.h: General math operations
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2004
5//# ATNF
6//#
7//# This program is free software; you can redistribute it and/or modify it
8//# under the terms of the GNU General Public License as published by the Free
9//# Software Foundation; either version 2 of the License, or (at your option)
10//# any later version.
11//#
12//# This program is distributed in the hope that it will be useful, but
13//# WITHOUT ANY WARRANTY; without even the implied warranty of
14//# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
15//# Public License for more details.
16//#
17//# You should have received a copy of the GNU General Public License along
18//# with this program; if not, write to the Free Software Foundation, Inc.,
19//# 675 Massachusetts Ave, Cambridge, MA 02139, USA.
20//#
21//# Correspondence concerning this software should be addressed as follows:
22//# Internet email: Malte.Marquarding@csiro.au
23//# Postal address: Malte Marquarding,
24//# Australia Telescope National Facility,
25//# P.O. Box 76,
26//# Epping, NSW, 2121,
27//# AUSTRALIA
28//#
29//# $Id: MathUtils.h 1514 2009-02-23 09:55:35Z KanaSugimoto $
30//#---------------------------------------------------------------------------
31#ifndef MATHUTILS_H
32#define MATHUTILS_H
33
34#include <string>
35#include <vector>
36#include <casa/aips.h>
37#include <casa/Arrays/Vector.h>
38#include <casa/BasicSL/String.h>
39#include <casa/Arrays/IPosition.h>
40
41namespace mathutil {
42
43// Hanning smoothing
44/**
45 * Hanning smooth a masked vector
46 * @param out the smoothed vector
47 * @param outmask the smoothed mask
48 * @param in the input vector
49 * @param mask the input mask
50 * @param relaxed a weighting scheme
51 * @param ignoreOther drop every second channel (NYI)
52 */
53void hanning(casa::Vector<casa::Float>& out,
54 casa::Vector<casa::Bool>& outmask,
55 const casa::Vector<casa::Float>& in,
56 const casa::Vector<casa::Bool>& mask,
57 casa::Bool relaxed=casa::False,
58 casa::Bool ignoreOther=casa::False);
59
60/**
61 * Apply a running median to a masked vector.
62 * Edge solution: The first and last hwidth channels will be replicated
63 * from the first/last value from a full window.
64 * @param out the smoothed vector
65 * @param outmask the smoothed mask
66 * @param in the input vector
67 * @param mask the input mask
68 * @param hwidth half-width of the smoothing window
69 */
70 void runningMedian(casa::Vector<casa::Float>& out,
71 casa::Vector<casa::Bool>& outmask,
72 const casa::Vector<casa::Float>& in,
73 const casa::Vector<casa::Bool>& mask,
74 float hwidth);
75
76// Generate specified statistic
77float statistics(const casa::String& which,
78 const casa::MaskedArray<casa::Float>& data);
79
80// Return a position of min or max value
81 casa::IPosition minMaxPos(const casa::String& which,
82 const casa::MaskedArray<casa::Float>& data);
83
84// Replace masked value by zero
85void replaceMaskByZero(casa::Vector<casa::Float>& data,
86 const casa::Vector<casa::Bool>& mask);
87
88/**
89 * Convert a std::vector of std::string
90 * to a casa::Vector casa::String
91 * @param in
92 * @return
93 */
94std::vector<std::string> tovectorstring(const casa::Vector<casa::String>& in);
95
96/**
97 * Convert a casa::Vector of casa::String
98 * to a stl std::vector of stl std::string
99 * @param in
100 * @return
101 */
102casa::Vector<casa::String> toVectorString(const std::vector<std::string>& in);
103
104}
105
106#endif
Note: See TracBrowser for help on using the repository browser.