casa6
Public
Commits
| 2085 2085 | //csys.obsInfo().telescope() << endl; |
| 2086 2086 | std::tie(cell,convnx)=getBeamCellSize(bandname); |
| 2087 2087 | convnx = Int(ceil(Float(convnx) / 8.0)) * 8; |
| 2088 2088 | //////////////////// |
| 2089 2089 | //cerr << "@@@cell " << cell << " pbnpix " << convnx << " imnpix "<< npix << " imcell"<< csys_p.increment()<< endl; |
| 2090 2090 | csys_p=csys; |
| 2091 2091 | Vector<String> units=csys_p.worldAxisUnits(); |
| 2092 2092 | Vector<Double> incr=csys_p.increment(); |
| 2093 2093 | Double inpFov=fabs(incr[0]*npix); |
| 2094 2094 | Double pbFov= fabs(cell.get(units[0]).getValue()*Double(convnx)); |
| 2095 - | cerr << "@@inpfov " << inpFov << " pbFov " << pbFov << endl; |
| 2096 2095 | if (inpFov > 0.125 * pbFov) { |
| 2097 2096 | incr[0] = cell.get(units[0]).getValue(); |
| 2098 2097 | incr[1] = cell.get(units[1]).getValue(); |
| 2099 2098 | //For small fov post fft resampling is not good enough...create alarge fov |
| 2100 2099 | // here to resample finer at calculation |
| 2101 2100 | if(inpFov < pbFov || isSingleField){ |
| 2102 2101 | incr[0]*=2.0; |
| 2103 2102 | incr[1]*=2.0; |
| 2104 2103 | convnx *=2.0; |
| 2105 2104 | pbFov *=4.0; |
| 2189 2188 | cerr << "@@@neArr shape "<< newArr.shape() << endl; |
| 2190 2189 | m.putMiddle(arr, newArr); |
| 2191 2190 | newArr.resize(); |
| 2192 2191 | newArr=m.resampleViaFFT(wtArr, fac, fac); |
| 2193 2192 | wtArr.set(0.0); |
| 2194 2193 | m.putMiddle(wtArr, newArr); |
| 2195 2194 | |
| 2196 2195 | |
| 2197 2196 | }*/ |
| 2198 2197 | //cerr << "Post FT MAX arr "<< max(wtArr) << " min "<< min(wtArr) << endl; |
| 2199 - | cerr << "inpFov/pbFov " << (inpFov / pbFov) << endl; |
| 2200 2198 | supportAndNormalizeAFunc(support, arr, wtArr, (inpFov/pbFov >1) ); |
| 2201 2199 | //cerr << "Post Norm MAX arr "<< max(arr) << " min "<< min(arr) << endl; |
| 2202 2200 | if(k==0){ |
| 2203 2201 | |
| 2204 2202 | IPosition tmpShape=shp; |
| 2205 2203 | tmpShape[3]=freqlist.nelements(); |
| 2206 2204 | //This is going to be returned to allow resacling etc later |
| 2207 2205 | //csys=refim::SynthesisUtils::makeUVCoords(csysA, tmpShape); |
| 2208 2206 | IPosition convShp(4, 2*support, 2*support, 4, freqlist.nelements()); |
| 2209 2207 | refpix[0]=refpix[1]=support; |