Changeset 504 for trunk/src/SDMemTable.cc
- Timestamp:
- 02/23/05 14:02:08 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SDMemTable.cc
r499 r504 43 43 #include <casa/Arrays/Vector.h> 44 44 #include <casa/BasicMath/Math.h> 45 #include <casa/BasicSL/Constants.h> 45 46 #include <casa/Quanta/MVAngle.h> 46 47 … … 377 378 378 379 379 std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol, 380 Float paOffset) const 380 std::vector<float> SDMemTable::getStokesSpectrum(Int whichRow, Bool doPol) const 381 381 // 382 382 // Gets one STokes parameter depending on cursor polSel location … … 394 394 if (doPol && (polSel_==1 || polSel_==2)) { // Q,U --> P, P.A. 395 395 396 // Set current cursor location 396 // Set current cursor location 397 397 398 const IPosition& shape = arr.shape(); 398 399 IPosition start, end; 399 400 getCursorSlice(start, end, shape); 400 401 401 402 // Get Q and U slices 402 403 403 404 Array<Float> Q = SDPolUtil::getStokesSlice(arr,start,end,"Q"); 404 405 Array<Float> U = SDPolUtil::getStokesSlice(arr,start,end,"U"); 405 406 406 407 // Compute output 407 408 408 409 Array<Float> out; … … 410 411 out = SDPolUtil::polarizedIntensity(Q,U); 411 412 } else if (polSel_==2) { // P.A. 412 out = SDPolUtil::positionAngle(Q,U) + paOffset;413 out = SDPolUtil::positionAngle(Q,U); 413 414 } 414 415
Note: See TracChangeset
for help on using the changeset viewer.