Designate class echo Find the maximum of aggregation values.
444 std::vector< std::vector< std::vector<HCA_Park> > > vol_Ai_filtered;
446 unsigned half_rg=0.5*(win_rg-1);
447 unsigned half_az=0.5*(win_az-1);
448 int pre_rg,post_rg,pre_az,post_az;
451 for(
unsigned el=0;el<
vol_Ai.size();el++)
459 for(
unsigned az=0;az<
vol_Ai[el].size();az++)
460 for(
unsigned rg=0;rg<
vol_Ai[el][az].size();rg++)
463 vol_Ai_filtered[el][az][rg].clearAi();
464 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][rg].Ai;
467 for(
unsigned j=1;j<half_az+1;j++)
472 if(pre_az<0)pre_az+=
vol_Ai[el].size();
473 if(post_az>=
vol_Ai[el].size())post_az-=
vol_Ai[el].size();
475 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][rg].Ai;
476 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][rg].Ai;
480 for(
unsigned i=1;i<half_rg+1;i++)
486 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][pre_rg].Ai;
489 if(post_rg<
vol_Ai[el][az].size())
491 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][az][post_rg].Ai;
495 for(
unsigned j=1;j<half_az+1;j++)
501 if(pre_az<0)pre_az+=
vol_Ai[el].size();
502 if(post_az>=
vol_Ai[el].size())post_az-=
vol_Ai[el].size();
505 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][pre_rg].Ai;
506 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][pre_az][pre_rg].Ai;
509 if(post_rg<
vol_Ai[el][az].size())
511 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][post_rg].Ai;
512 vol_Ai_filtered[el][az][rg].Ai+=
vol_Ai[el][post_az][post_rg].Ai;
518 vol_Ai_filtered[el][az][rg].Ai.array()/=(double)count;
520 vol_hca[el](az,rg)=vol_Ai_filtered[el][az][rg].echo(0.00001);
528 for(
unsigned el=0;el<
vol_z.size();el++)
532 for(
unsigned az=0;az<
vol_z.
scan(el).beam_count;az++)
533 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