Commits
94 94 | //freqFrameValid_p=freqFrameValid; |
95 95 | //chanchunk may call the same object |
96 96 | if(initialized_p && nx_p==templateimage.shape()(0) && ny_p==templateimage.shape()(1)){ |
97 97 | CoordinateSystem cs=templateimage.coordinates(); |
98 98 | Double freq=cs.toWorld(IPosition(4,0,0,0,0))[3]; |
99 99 | if(freq==refFreq_p) |
100 100 | return; |
101 101 | } |
102 102 | //cerr << "in bgwt init " << endl; |
103 103 | //Need to save previous wieght scheme of vi |
104 - | VisImagingWeight visWgt_p=vi.getImagingWeightGenerator(); |
105 - | VisImagingWeight vWghtNat("natural"); |
104 + | visWgt_p=vi.getImagingWeightGenerator(); |
105 + | VisImagingWeight vWghtNat("natural"); |
106 106 | vi.useImagingWeight(vWghtNat); |
107 107 | vi::VisBuffer2 *vb=vi.getVisBuffer(); |
108 108 | Int nIndices=0; |
109 109 | for (vi.originChunks();vi.moreChunks();vi.nextChunk()) { |
110 110 | for (vi.origin(); vi.more(); vi.next()) { |
111 111 | String key=String::toString(vb->msId())+"_"+String::toString(vb->fieldId()(0)); |
112 112 | Int index=0; |
113 113 | if(multiField_p){ |
114 114 | //find how many indices will be needed |
115 115 | index=multiFieldMap_p.size(); |
292 292 | //ndrop++; |
293 293 | //} |
294 294 | } |
295 295 | else{ |
296 296 | imweight(chn,row)=0.0; |
297 297 | } |
298 298 | |
299 299 | } |
300 300 | } |
301 301 | |
302 + | |
302 303 | if(visWgt_p.doFilter()){ |
303 304 | visWgt_p.filter (imweight, flag, uvw, vb.getFrequencies(0), imweight); |
304 305 | |
305 306 | } |
306 307 | |
307 308 | } |
308 309 | |
309 310 | void BriggsCubeWeightor::initializeFTMachine(const uInt index, const ImageInterface<Complex>& templateimage, const RecordInterface& inRec){ |
310 311 | Int nchan=templateimage.shape()(3); |
311 312 | if(ft_p.nelements() <= index){ |