AddCleanerQuantities.cpp

#include <iostream>
#include <cstring>
#include <radarelab/volume.h>
#include <radarelab/odim.h>
#include <radarlib/radar.hpp>
#include <sstream>
#include <radarelab/image.h>
#include <radarelab/algo/azimuth_resample.h>
#include <radarelab/algo/cleaner.h>

#include <radarelab/algo/elabora_volume.h>

using namespace radarelab;
using namespace std;

using namespace volume;
namespace odim = OdimH5v21;

int main(int argc,char* argv[])
{
        if (argc < 1)
        {
                fprintf(stderr, "Usage: %s <h5-volume-input> [--Use_undetect] \n", argv[0]);
                exit(1);
        }
        
        std::string pathname = argv[1]; 

        printf("il nome del mio file è %s\n", pathname.c_str());

        volume::ODIMLoader loader_all;

        volume::Scans<double> full_volume_z;
        volume::Scans<double> full_volume_zdr;
        volume::Scans<double> full_volume_vrad;
        volume::Scans<double> full_volume_wrad;
        volume::Scans<double> full_volume_sqi;
        volume::Scans<unsigned char> full_volume_cleanID;

        volume::Scans<double> Z_SD2D;
        Z_SD2D.quantity="DBZH_SD2D";
        volume::Scans<double> Z_SDRay_9;
        Z_SDRay_9.quantity="DBZH_SDRay_9";
        volume::Scans<double> Z_SDRay_21;
        Z_SDRay_21.quantity="DBZH_SDRay_21";
        volume::Scans<double> Z_SDAz;
        Z_SDAz.quantity="DBZH_SDAz";
        volume::Scans<double> ZDR_SD2D;
        ZDR_SD2D.quantity="ZDR_SD2D";

        std::string task;

        loader_all.request_quantity(odim::PRODUCT_QUANTITY_DBZH,&full_volume_z);
        loader_all.request_quantity(odim::PRODUCT_QUANTITY_ZDR,&full_volume_zdr);
        loader_all.request_quantity(odim::PRODUCT_QUANTITY_VRAD,&full_volume_vrad);
        loader_all.request_quantity(odim::PRODUCT_QUANTITY_WRAD,&full_volume_wrad);
        loader_all.request_quantity(odim::PRODUCT_QUANTITY_SQI,&full_volume_sqi);



        loader_all.load(argv[1]);

    
        for (unsigned i = 0; i < full_volume_z.size(); ++i){
            volume::Scans<double> Input, Input_zdr, Texture;
            Input.push_back(full_volume_z.at(i));
// Texture Z 2D 
            radarelab::volume::textureSD( Input, Texture, 1000. , 3.,false);
            Texture.at(0).nodata=65535.;
            Texture.at(0).undetect=0.;
            Z_SD2D.push_back(Texture.at(0));
// Texture Z 1D Ray 9 bins
            radarelab::volume::textureSD( Input,Texture, Input.at(0).cell_size*9 , 360./Input.at(0).beam_count,true);
            Texture.at(0).nodata=65535.;
            Texture.at(0).undetect=0.;
            Z_SDRay_9.push_back(Texture.at(0));
// Texture Z 1D Ray 21 bins
            radarelab::volume::textureSD( Input,Texture, Input.at(0).cell_size*21 , 360./Input.at(0).beam_count,true);
            Texture.at(0).nodata=65535.;
            Texture.at(0).undetect=0.;
            Z_SDRay_21.push_back(Texture.at(0));
// Texture Z 1D Az 5 rays
            radarelab::volume::textureSD( Input,Texture, Input.at(0).cell_size , 5*360./Input.at(0).beam_count,true);
            Texture.at(0).nodata=65535.;
            Texture.at(0).undetect=0.;
            Z_SDAz.push_back(Texture.at(0));
// Texture ZDR 2D 
            Input_zdr.push_back(full_volume_zdr.at(i));
            radarelab::volume::textureSD( Input_zdr, Texture, 1000. , 3.,false);
            Texture.at(0).nodata=65535.;
            Texture.at(0).undetect=0.;
            ZDR_SD2D.push_back(Texture.at(0));
        }
      std::cout<<"Finito Cleaner, salvo risultati"<<std::endl;
        volume::ODIMStorer storer;
        storer.store_quantity_fp((Volume<double>*)(&Z_SD2D));
        storer.store_quantity_fp((Volume<double>*)(&Z_SDRay_9));
        storer.store_quantity_fp((Volume<double>*)(&Z_SDRay_21));
        storer.store_quantity_fp((Volume<double>*)(&Z_SDAz));
        storer.store_quantity_fp((Volume<double>*)(&ZDR_SD2D));
        storer.store(argv[1]);
        cout<<endl<<"Fine"<<endl;
}


//
//    radarelab::volume::Scans<double>   Z_S,  SD2D,SD_Ray,SD_Az;
//    Z_S.push_back(scan_z);
//    radarelab::volume::textureSD( Z_S,SD2D, 1000. , 3,false);

//    radarelab::volume::textureSD( Z_S,SD_Ray, scan_z.cell_size*9 , 360./scan_z.beam_count,true);
//    radarelab::volume::textureSD( Z_S,SD_Az, scan_z.cell_size , 5*360./scan_z.beam_count,true);
//