reconstructionCUDA.cuh

#pragma once
#include <iostream>
#include "cuda_runtime.h"
#include "device_launch_parameters.h"
#include <stdio.h>
#include <string>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <complex>
#include <vector>
#include <cufft.h>
#include <chrono>
#include <thrust/host_vector.h>
#include <thrust/device_vector.h>
#include "gdal_alg.h";
#include "gdal_priv.h"
#include <gdal.h>


#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>



__global__ void Zhuan_Complex_kernel(float *PSF_1_gpu, cufftComplex *PSF_1_gpu_Complex, int total);

__global__ void PSF_unshort(float *PSF_1_gpu, unsigned short *PSF, int total);

__global__ void initial_kernel_1(float *ImgEst, float *Ratio, int total);

__global__ void gpuObjRecon_fuzhi(float *gpuObjRecon, int total);

__global__ void initial_kernel_3(float *gpuObjRecROI, int total);

__global__ void ImgExp_ge(unsigned short *Img_gpu, int BkgMean, float *ImgExp, int total);

__global__ void Ratio_fuzhi(float *Ratio, int total);

__global__ void OTF_mul_gpuObjRecon_Complex(cufftComplex *OTF, cufftComplex *gpuObjRecon_Complex, int total);

__global__ void ifftshift_real_max(cufftComplex *OTF, float *float_temp, int PSF_size_1, int PSF_size_2, int PSF_size_3);

__global__ void ifftshift(cufftComplex *OTF, float *float_temp, int PSF_size_1, int PSF_size_2, int PSF_size_3, cufftComplex *OTF_ifftshift);

__global__ void float_temp_sum(float *float_temp, float *ImgEst, int PSF_size_1, int PSF_size_2, int PSF_size_3);

__global__ void Ratio_fuzhi_2(float *ImgExp, float *ImgEst, float Tmp, int SNR, float *Ratio, int total);

__global__ void Ratio_Complex_ge(float *ImgExp, float *ImgEst, float Tmp, int SNR, cufftComplex *Ratio_Complex, int total);

__global__ void fftRatio_ge(cufftComplex *Ratio_Complex, cufftComplex *fftRatio, int PSF_size_1, int PSF_size_2, int PSF_size_3);

__global__ void fftceshi_gpu_fuzhi(cufftComplex *PSF_1_gpu_Complex, cufftComplex *fftceshi_gpu, int total);

__global__ void ifft2_divide(cufftComplex *OTF, int total, int scale);

__global__ void real_multiply(float *gpuObjRecon, float *float_temp, int total);

__global__ void fftRatio_mul_conjOTF(cufftComplex *fftRatio, cufftComplex *OTF, int total);

__global__ void cropReconImage_kernel(float *gpuObjRecon, float *gpuObjRecon_crop);