modularization/src/detection/vision_TensorRT_yolov3/main.cpp

#include <opencv2/opencv.hpp>
//txs
#include "opencv2/video/tracking.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "cameraobjectarray.pb.h"
//txs
#include "TrtNet.h"
#include "argsParser.h"
#include "configs.h"
#include <chrono>
#include "YoloLayer.h"
#include "dataReader.h"
#include "eval.h"

#include "imageBuffer.h"

#include "modulecomm.h"
#include "xmlparam.h"
#include "rawpic.pb.h"

#include <vector>
#include <QFile>
#include <QTextStream>
#include "localization.h"
using namespace std;
using namespace argsParser;
using namespace Tn;
using namespace Yolo;

//txs
using namespace cv;
iv::vision::cameraobjectarray gobj;
int gntemp = 0;
void * gpa;
void * gpcamout;
void * gdetec;
QTime gTime;
//txs

ConsumerProducerQueue<cv::Mat> *imageBuffer = new ConsumerProducerQueue<cv::Mat>(1, true);

static vector<string> GetNamesFromFile(string filePathAndName)
{
    vector<string> retVec;
    QFile qf(QString::fromStdString(filePathAndName));
    if(qf.exists())
    {
        if(qf.open(QIODevice::ReadOnly))
        {
            QTextStream in(&qf);
            while(!in.atEnd())
            {
                retVec.push_back(in.readLine(1024).toStdString());
            }
            qf.close();
        }
    }

    return retVec;
}


vector<float> prepareImage(cv::Mat& img)
{
    using namespace cv;

    int c = parser::getIntValue("C");
    int h = parser::getIntValue("H");   //net h
    int w = parser::getIntValue("W");   //net w

    float scale = min(float(w)/img.cols,float(h)/img.rows);
    auto scaleSize = cv::Size(img.cols * scale,img.rows * scale);

    cv::Mat rgb ;
    cv::cvtColor(img, rgb, CV_BGR2RGB);
    cv::Mat resized;
    cv::resize(rgb, resized,scaleSize,0,0,INTER_CUBIC);

    cv::Mat cropped(h, w,CV_8UC3, 127);
    Rect rect((w- scaleSize.width)/2, (h-scaleSize.height)/2, scaleSize.width,scaleSize.height);
    resized.copyTo(cropped(rect));

    cv::Mat img_float;
    if (c == 3)
        cropped.convertTo(img_float, CV_32FC3, 1/255.0);
    else
        cropped.convertTo(img_float, CV_32FC1 ,1/255.0);

    //HWC TO CHW
    vector<Mat> input_channels(c);
    cv::split(img_float, input_channels);

    vector<float> result(h*w*c);
    auto data = result.data();
    int channelLength = h * w;
    for (int i = 0; i < c; ++i) {
        memcpy(data,input_channels[i].data,channelLength*sizeof(float));
        data += channelLength;
    }

    return result;
}

void DoNms(vector<Detection>& detections,int classes ,float nmsThresh)
{
    auto t_start = chrono::high_resolution_clock::now();

    vector<vector<Detection>> resClass;
    resClass.resize(classes);

    for (const auto& item : detections)
        resClass[item.classId].push_back(item);

    auto iouCompute = [](float * lbox, float* rbox)
    {
        float interBox[] = {
            max(lbox[0] - lbox[2]/2.f , rbox[0] - rbox[2]/2.f), //left
            min(lbox[0] + lbox[2]/2.f , rbox[0] + rbox[2]/2.f), //right
            max(lbox[1] - lbox[3]/2.f , rbox[1] - rbox[3]/2.f), //top
            min(lbox[1] + lbox[3]/2.f , rbox[1] + rbox[3]/2.f), //bottom
        };
        
        if(interBox[2] > interBox[3] || interBox[0] > interBox[1])
            return 0.0f;

        float interBoxS =(interBox[1]-interBox[0])*(interBox[3]-interBox[2]);
        return interBoxS/(lbox[2]*lbox[3] + rbox[2]*rbox[3] -interBoxS);
    };

    vector<Detection> result;
    for (int i = 0;i<classes;++i)
    {
        auto& dets =resClass[i];
        if(dets.size() == 0)
            continue;

        sort(dets.begin(),dets.end(),[=](const Detection& left,const Detection& right){
            return left.prob > right.prob;
        });

        for (unsigned int m = 0;m < dets.size() ; ++m)
        {
            auto& item = dets[m];
            result.push_back(item);
            for(unsigned int n = m + 1;n < dets.size() ; ++n)
            {
                if (iouCompute(item.bbox,dets[n].bbox) > nmsThresh)
                {
                    dets.erase(dets.begin()+n);
                    --n;
                }
            }
        }
    }

    //swap(detections,result);
    detections = move(result);

    auto t_end = chrono::high_resolution_clock::now();
    float total = chrono::duration<float, milli>(t_end - t_start).count();
    mTensorIvlog->info("Time taken for nms is %f ms.", total);
//    cout << "Time taken for nms is " << total << " ms." << endl;
}


vector<Bbox> postProcessImg(cv::Mat& img,vector<Detection>& detections,int classes)
{
    using namespace cv;

    int h = parser::getIntValue("H");   //net h
    int w = parser::getIntValue("W");   //net w

    //scale bbox to img
    int width = img.cols;
    int height = img.rows;
    float scale = min(float(w)/width,float(h)/height);
    float scaleSize[] = {width * scale,height * scale};

    //correct box
    for (auto& item : detections)
    {
        auto& bbox = item.bbox;
        bbox[0] = (bbox[0] * w - (w - scaleSize[0])/2.f) / scaleSize[0];
        bbox[1] = (bbox[1] * h - (h - scaleSize[1])/2.f) / scaleSize[1];
        bbox[2] /= scaleSize[0];
        bbox[3] /= scaleSize[1];
    }
    
    //nms
    float nmsThresh = parser::getFloatValue("nms");
    if(nmsThresh > 0)
        DoNms(detections,classes,nmsThresh);

    vector<Bbox> boxes;
    for(const auto& item : detections)
    {
        auto& b = item.bbox;
        Bbox bbox =
        {
            item.classId,   //classId
            max(int((b[0]-b[2]/2.)*width),0), //left
            min(int((b[0]+b[2]/2.)*width),width), //right
            max(int((b[1]-b[3]/2.)*height),0), //top
            min(int((b[1]+b[3]/2.)*height),height), //bot
            item.prob       //score
        };
        boxes.push_back(bbox);
    }

    return boxes;
}

vector<string> split(const string& str, char delim)
{
    stringstream ss(str);
    string token;
    vector<string> container;
    while (getline(ss, token, delim)) {
        container.push_back(token);
    }

    return container;
}

qint64 gpictime;
qint64 gpicindex;
//cv::Mat gmatimage;
void Listenpic(const char * strdata,const unsigned int nSize,const unsigned int index,const QDateTime * dt,const char * strmemname)
{
    if(nSize<1000)return;
    iv::vision::rawpic pic;

    if(false == pic.ParseFromArray(strdata,nSize))
    {
//        std::cout<<"picview Listenpic fail."<<std::endl;
        mTensorIvlog->info("picview Listenpic fail.");
        return;
    }

    cv::Mat mat(pic.height(),pic.width(),pic.mattype());
    if(pic.type() == 1)
        memcpy(mat.data,pic.picdata().data(),mat.rows*mat.cols*mat.elemSize());
    else
    {
       mat.release();
       std::vector<unsigned char> buff(pic.picdata().data(),pic.picdata().data() + pic.picdata().size());
       mat = cv::imdecode(buff,CV_LOAD_IMAGE_COLOR);
    }
    gpictime = pic.time();
    gpicindex = pic.index();
    imageBuffer->add(mat);
    mat.release();

}
//std::vector<Obstacle> objvec,iv::lidar::objectarray & lidarobjvec
//txs
void GetCameraObj(std::vector<Bbox> objvec,iv::vision::cameraobjectarray &cameraobjvec)
{
    mTensorIvlog->info("Bbox count: %d", objvec.size());
//    std::cout<<"Bbox count: "<<objvec.size()<<std::endl;
    int i;
    for(i=0;i<objvec.size();i++){
        Bbox x;
        //iv::vision::cameraobject cameraobj;
        x = objvec.at(i);
        mTensorIvlog->info("left: %d\tright: %d", x.left, x.right);
//        std::cout<<"left: " <<x.left<<"\tright: "<<x.right<<std::endl;

        iv::vision::cameraobject obj;
        const char id = x.classId;
        obj.set_type(&id);
        obj.set_outpointx((x.left+x.right)/2);
        //2019-10-10 09-01AM
        //obj.set_outpointy((x.bot+x.top)/2);
        obj.set_outpointy(x.bot);
        obj.set_con(x.score);

        mTensorIvlog->info("raw output x: %f\traw output y: %f", obj.outpointx(), obj.outpointy());
//        std::cout<<"raw output x: "<<obj.outpointx()<<"\traw output y:" <<obj.outpointy()<<std::endl;

        iv::vision::cameraobject * po = cameraobjvec.add_obj();
        po->CopyFrom(obj);
    }

}

//txs


int main( int argc, char* argv[] )
{
    parser::ADD_ARG_STRING("prototxt",Desc("input yolov3 deploy"),DefaultValue(INPUT_PROTOTXT),ValueDesc("file"));
    parser::ADD_ARG_STRING("caffemodel",Desc("input yolov3 caffemodel"),DefaultValue(INPUT_CAFFEMODEL),ValueDesc("file"));
    parser::ADD_ARG_INT("C",Desc("channel"),DefaultValue(to_string(INPUT_CHANNEL)));
    parser::ADD_ARG_INT("H",Desc("height"),DefaultValue(to_string(INPUT_HEIGHT)));
    parser::ADD_ARG_INT("W",Desc("width"),DefaultValue(to_string(INPUT_WIDTH)));
    parser::ADD_ARG_STRING("calib",Desc("calibration image List"),DefaultValue(CALIBRATION_LIST),ValueDesc("file"));
    parser::ADD_ARG_STRING("mode",Desc("runtime mode"),DefaultValue(MODE), ValueDesc("fp32/fp16/int8"));
    parser::ADD_ARG_STRING("outputs",Desc("output nodes name"),DefaultValue(OUTPUTS));
    parser::ADD_ARG_INT("class",Desc("num of classes"),DefaultValue(to_string(DETECT_CLASSES)));
    parser::ADD_ARG_FLOAT("nms",Desc("non-maximum suppression value"),DefaultValue(to_string(NMS_THRESH)));
    parser::ADD_ARG_INT("batchsize",Desc("batch size for input"),DefaultValue("1"));
    parser::ADD_ARG_STRING("enginefile",Desc("load from engine"),DefaultValue(""));

    //input
    parser::ADD_ARG_STRING("input",Desc("input image file"),DefaultValue(INPUT_IMAGE),ValueDesc("file"));
    parser::ADD_ARG_STRING("evallist",Desc("eval gt list"),DefaultValue(EVAL_LIST),ValueDesc("file"));

    gTensorFault = new iv::Ivfault("TensorRT_yolov3");
    mTensorIvlog = new iv::Ivlog("TensorRT_youlov3");

//    //txs
//    KalmanFilter KF(2,1,0);
//    Mat state();
//    //txs

//    if(argc < 2){
//        parser::printDesc();
//        exit(-1);
//    }

//    parser::parseArgs(argc,argv);

    QString    strpath = "./vision_TensorRT_yolov3.xml";

    mTensorIvlog->info("%s", strpath.data());
//    std::cout<<strpath.toStdString()<<std::endl;
    iv::xmlparam::Xmlparam xp(strpath.toStdString());

    std::string xmlmsgname = xp.GetParam("msgname","image00");

    //string engineName = xp.GetParam("engineName","/home/nvidia/TensorRT-Yolov3/yolov3_fp16_608.engine");
    string engineName = "/home/nvidia/code/sd/src/detection/vision_TensorRT_yolov3/model/yolov3.engine";
//    string strwidth = xp.GetParam("width",608);
//    string strheight = xp.GetParam("height",608);
    mTensorIvlog->info("msgname is %s",xmlmsgname.data());
//    std::cout<<"msgname is "<<xmlmsgname<<std::endl;
//    void * mpa = iv::modulecomm::RegisterRecv(xmlmsgname.data(),Listenpic);
    void *pa;
    pa = iv::modulecomm::RegisterRecv("camera",Listenpic);

    gpcamout = iv::modulecomm::RegisterSend("cameraoutput",1000000,1);

    //txs
    gdetec = iv::modulecomm::RegisterSend("camera_detect",10000000,1);
    //txs

    std::unique_ptr<trtNet> net;
//    int batchSize = parser::getIntValue("batchsize");
//    string engineName =  parser::getStringValue("enginefile");
//    if(engineName.length() > 0)
//    {
//        net.reset(new trtNet(engineName));
//        assert(net->getBatchSize() == batchSize);
//    }
//    else
//    {
//    string deployFile = parser::getStringValue("prototxt");
//    string caffemodelFile = parser::getStringValue("caffemodel");

//    vector<vector<float>> calibData;
//    string calibFileList = parser::getStringValue("calib");
//    string mode = parser::getStringValue("mode");
//    if(calibFileList.length() > 0 && mode == "int8")
//    {
//        cout << "find calibration file,loading ..." << endl;

//        ifstream file(calibFileList);
//        if(!file.is_open())
//        {
//            cout << "read file list error,please check file :" << calibFileList << endl;
//            exit(-1);
//        }

//        string strLine;
//        while( getline(file,strLine) )
//        {
//            cv::Mat img = cv::imread(strLine);
//            auto data = prepareImage(img);
//            calibData.emplace_back(data);
//        }
//        file.close();
//    }

    int batchSize = 1;
    net.reset(new trtNet(engineName));
    assert(net->getBatchSize() == batchSize);

//    RUN_MODE run_mode = RUN_MODE::FLOAT32;
//    if(mode == "int8")
//    {
//        if(calibFileList.length() == 0)
//            cout << "run int8 please input calibration file, will run in fp32" << endl;
//        else
//            run_mode = RUN_MODE::INT8;
//    }
//    else if(mode == "fp16")
//    {
//        run_mode = RUN_MODE::FLOAT16;
//    }
    
//    string outputNodes = parser::getStringValue("outputs");
//    auto outputNames = split(outputNodes,',');
    
//        //save Engine name
//    string saveName = "yolov3_" + mode + ".engine";
//        net.reset(new trtNet(deployFile,caffemodelFile,outputNames,calibData,run_mode,batchSize));
//    cout << "save Engine..." << saveName <<endl;
//        net->saveEngine(saveName);
//    }

    RUN_MODE run_mode = RUN_MODE::FLOAT16;

    int outputCount = net->getOutputSize()/sizeof(float);
    unique_ptr<float[]> outputData(new float[outputCount]);

//    string listFile = parser::getStringValue("evallist");

    string listFile;

    list<string> fileNames;
    list<vector<Bbox>> groundTruth;

//    if(listFile.length() > 0)
//    {
//        std::cout << "loading from eval list " << listFile << std::endl;
//        tie(fileNames,groundTruth) = readObjectLabelFileList(listFile);
//    }
//    else
//    {
//        string inputFileName = parser::getStringValue("input");
//        fileNames.push_back(inputFileName);
//    }

    string inputFileName = "/home/nvidia/TensorRT-Yolov3/test.jpg";
    fileNames.push_back(inputFileName);

    list<vector<Bbox>> outputs;
//    int classNum = parser::getIntValue("class");
//    int c = parser::getIntValue("C");
//    int h = parser::getIntValue("H");
//    int w = parser::getIntValue("W");
    int classNum = 80;
    int c = 608;
    int h = 608;
    int w = 3;
    int batchCount = 0;
    vector<float> inputData;
    inputData.reserve(h*w*c*batchSize);
    vector<cv::Mat> inputImgs;

    auto iter = fileNames.begin();
    //txs
    cv::namedWindow("result",WINDOW_NORMAL);

    cv::VideoWriter ou;
    ou.open("test.avi",CV_FOURCC('X', 'V', 'I', 'D'),20.0,Size(1280, 720), true);
    //txs
    cv::Mat img2;
    int count = 0;
    while(1)
    {
        imageBuffer->consume(img2);
        if(!img2.rows)
        {
            gTensorFault->SetFaultState(1, 0, "img2.rows state Error");
            break;
        }

        vector<float> curInput = prepareImage(img2);

        net->doInference(curInput.data(), outputData.get(),1);

        mTensorIvlog->info("run");

//        std::cout<<"run"<<std::endl;

        auto output = outputData.get();

        int detCount = output[0];
    //later detect result
        vector<Detection> result;
        result.resize(detCount);
        memcpy(result.data(), &output[1], detCount*sizeof(Detection));

        outputs.clear();
        auto boxes = postProcessImg(img2,result,classNum);
        outputs.emplace_back(boxes);



        vector<string> cocoNames = GetNamesFromFile("/home/nvidia/code/sd/src/detection/vision_TensorRT_yolov3/model/coco.names");
//        cout<<"cocoNames: "<<endl<<endl;
        mTensorIvlog->info("cocoName:");
        for(int i=0; i<cocoNames.size(); i++)
        {
//            cout<<cocoNames[i]<<endl;
            mTensorIvlog->info("%s",cocoNames[i].data());
        }
//        cout<<"cocoNames end"<<endl<<endl;
        mTensorIvlog->info("cocoNames End");


        //txs
        iv::vision::cameraobjectarray cameraobjvec;
        GetCameraObj(boxes,cameraobjvec);

        int ntlen;
        std::string out = cameraobjvec.SerializeAsString();

        iv::modulecomm::ModuleSendMsg(gpcamout,out.data(),out.length());
        mTensorIvlog->info("--------------");
//        std::cout<<"------------------"<<std::endl;
        mTensorIvlog->info("lenth is %d",out.length());

        //txs


        auto bbox = *outputs.begin();
        for(const auto& item : bbox)
        {
            cv::rectangle(img2,cv::Point(item.left,item.top),cv::Point(item.right,item.bot),cv::Scalar(0,0,255),3,8,0);
            std::cout << "class=" << item.classId << " prob=" << item.score*100 << std::endl;
            cout << "left=" << item.left << " right=" << item.right << " top=" << item.top << " bot=" << item.bot << endl;

            cout<<cocoNames[item.classId]<<endl;
            //txs
            MatrixXd local(2,1);
            local<<(item.left +item.right)/2,(item.bot +item.bot)/2;
 //           local<<(item.bot +item.bot)/2,(item.left +item.right)/2;
            MatrixXd location=TranslateToStandardLocation::translateCameraToStandardLocation(local);
 //           string str = cocoNames[item.classId] +":" + std::to_string(item.score*100)+":"+std::to_string(location(0,0))+std::to_string(location(1,0));
            string str = std::to_string((int)(location(0,0)/100.0))+" " + std::to_string((int)(location(1,0)/100.0));
            string str1=cocoNames[item.classId];

            cv::putText(img2,str, cv::Point(item.left, item.top), cv::FONT_HERSHEY_PLAIN, 1, cv::Scalar(255, 0, 0));
            cv::putText(img2,str1, cv::Point(item.left, item.top-15), cv::FONT_HERSHEY_PLAIN, 1, cv::Scalar(0, 255, 0));
            //txs
        }

//        ou.write(img2);
        //txs
//        cv::imwrite("img/"+ std::to_string(count)+".jpg", img2);
        //txs

        iv::vision::rawpic pic;
        pic.set_time(gpictime);
        pic.set_index(gpicindex);
        pic.set_elemsize(img2.elemSize());
        pic.set_width(img2.cols);
        pic.set_height(img2.rows);
        pic.set_mattype(img2.type());

        std::vector<int> param = std::vector<int>(2);
        param[0] = CV_IMWRITE_JPEG_QUALITY;
        param[1] = 95; // default(95) 0-100
        std::vector<unsigned char> buff;
        cv::imencode(".jpg",img2,buff,param);
        pic.set_picdata(buff.data(),buff.size());
        buff.clear();
        pic.set_type(2);

        unsigned int nsize = pic.ByteSize();
        char * strser = new char[nsize];
        bool bser = pic.SerializeToArray(strser,nsize);
        if(bser)iv::modulecomm::ModuleSendMsg(gdetec,strser,nsize);

        delete[] strser;


        count += 1;
        cv::imshow("result",img2);


        cv::waitKey(1);
        curInput.clear();
//        free(output_cpu);
//        free(output2_cpu);
//        img2.release();
//        ou.release();
    }
        img2.release();
        ou.release();
    for (unsigned int i = 0;i < fileNames.size(); ++i ,++iter)
    {
        const string& filename  = *iter;

        std::cout << "process: " << filename << std::endl;

        cv::Mat img = cv::imread(filename);
        vector<float> curInput = prepareImage(img);
        if (!curInput.data())
            continue;
        inputImgs.emplace_back(img);

        inputData.insert(inputData.end(), curInput.begin(), curInput.end());
        batchCount++;

        if(batchCount < batchSize && i + 1 <  fileNames.size())
            continue;

        net->doInference(inputData.data(), outputData.get(),batchCount);

        //Get Output
        auto output = outputData.get();
        auto outputSize = net->getOutputSize()/ sizeof(float) / batchCount;
        for(int i = 0;i< batchCount ; ++i)
        {
        //first detect count
            int detCount = output[0];
        //later detect result
        vector<Detection> result;
            result.resize(detCount);
            memcpy(result.data(), &output[1], detCount*sizeof(Detection));

            auto boxes = postProcessImg(inputImgs[i],result,classNum);
        outputs.emplace_back(boxes);

            output += outputSize;
        }
        inputImgs.clear();
        inputData.clear();

        batchCount = 0;
    }
    
    
    net->printTime();

    if(groundTruth.size() > 0)
    {
        //eval map
        evalMAPResult(outputs,groundTruth,classNum,0.5f);
        evalMAPResult(outputs,groundTruth,classNum,0.75f);
    }

    if(fileNames.size() == 1)
    {
        //draw on image
        cv::Mat img = cv::imread(*fileNames.begin());
        auto bbox = *outputs.begin();
        for(const auto& item : bbox)
        {
            cv::rectangle(img,cv::Point(item.left,item.top),cv::Point(item.right,item.bot),cv::Scalar(0,0,255),3,8,0);
            cout << "class=" << item.classId << " prob=" << item.score*100 << endl;
            cout << "left=" << item.left << " right=" << item.right << " top=" << item.top << " bot=" << item.bot << endl;
        }
        cv::imwrite("result.jpg",img);
        cv::imshow("result",img);
        cv::waitKey(0);
    }
    gTensorFault->SetFaultState(0, 0, "ok");

    return 0;
}