for(int i=0; i < vsize; i++){

    int boxsize=10;

    int vispixi=int(visu[i])+1024;
    int vispixj=int(visv[i])+1024;

    int boxmini=vispixi-boxsize;
    int boxmaxi=vispixi+boxsize;
    int boxminj=vispixj-boxsize;
    int boxmaxj=vispixj+boxsize;

    for(int pix1=boxmini; pix1<=boxmaxi; pix1++){
        for(int pix2=boxminj; pix2<=boxmaxj; pix2++){
            realmodelcube[pix1][pix2][int(vischan[i])-1] += apeturearray[pix1-vispixi+boxsize][pix2-vispixj+boxsize][int(vischan[i])-1]*NoiseDiffSignal[i][0];
            immodelcube[pix1][pix2][int(vischan[i])-1] += apeturearray[pix1-vispixi+boxsize][pix2-vispixj+boxsize][int(vischan[i])-1]*NoiseDiffSignal[i][1];
        }
    }
}