Designate class echo Find the maximum of aggregation values.
450 std::vector< std::vector< std::vector<HCA_Park> > > vol_Ai_filtered;
452 unsigned half_rg=0.5*(win_rg-1);
453 unsigned half_az=0.5*(win_az-1);
454 int pre_rg,post_rg,pre_az,post_az;
457 for(
unsigned el=0;el<
vol_Ai.size();el++)
465 for(
unsigned az=0;az<
vol_Ai[el].size();az++)
466 for(
unsigned rg=0;rg<
vol_Ai[el][az].size();rg++)
469 vol_Ai_filtered[el][az][rg].clearAi();
470 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][rg].Ai;
473 for(
unsigned j=1;j<half_az+1;j++)
478 if(pre_az<0)pre_az+=
vol_Ai[el].size();
479 if(post_az>=
vol_Ai[el].size())post_az-=
vol_Ai[el].size();
481 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][rg].Ai;
482 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][rg].Ai;
486 for(
unsigned i=1;i<half_rg+1;i++)
492 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][pre_rg].Ai;
495 if(post_rg<
vol_Ai[el][az].size())
497 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][post_rg].Ai;
501 for(
unsigned j=1;j<half_az+1;j++)
507 if(pre_az<0)pre_az+=
vol_Ai[el].size();
508 if(post_az>=
vol_Ai[el].size())post_az-=
vol_Ai[el].size();
511 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][pre_rg].Ai;
512 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][pre_rg].Ai;
515 if(post_rg<
vol_Ai[el][az].size())
517 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][post_rg].Ai;
518 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][post_rg].Ai;
524 vol_Ai_filtered[el][az][rg].Ai.array()/=(double)count;
526 vol_hca[el](az,rg)=vol_Ai_filtered[el][az][rg].echo(0.00001);
534 for(
unsigned el=0;el<
vol_z.size();el++)
538 for(
unsigned az=0;az<
vol_z.
scan(el).beam_count;az++)
539 for(
unsigned rg=0;rg<
vol_z.
scan(el).beam_size;rg++)
PolarScan< T > & scan(unsigned idx)
Access a polar scan.
Volume< EchoClass > vol_hca
PolarScan - structure to describe a polarScan containing a matrix of data and conversion factors...
std::vector< std::vector< std::vector< HCA_Park > > > vol_Ai