ηT=0.9;
DC=-100;
ηM=0.53;
ηLe=0.98;

APDdata={Import["E:\\Imaging\\Fig3\\Fig3bAPD.csv",{"Data",All}]};
CAMdata={Import["E:\\Imaging\\Fig3\\Fig3bCCD.csv",{"Data",All}]};
APD=Take[APDdata,1,{3,Length[APDdata[[1]]]}];
CAM=Take[CAMdata,1,{3,Length[CAMdata[[1]]]}];

APD1=Table[Table[If[APD[[1,j,i]]<1,APD[[1,j,i-3]],APD[[1,j,i]]],{i,5,Length[APD[[1,1]]]}],{j,1,Length[APD[[1]]]}];
CAM1=Table[Table[CAM[[1,j,i]],{i,5,Length[CAM[[1,1]]]}],{j,1,Length[CAM[[1]]]}];
η=Table[Mean[Table[((CAM1[[j,i]]*ηT)/(APD1[[j,i]]+DC)),{i,1,Length[APD1[[1]]]}]],{j,1,Length[CAM1]}];
ηC=Table[{(APD[[1,j,4]]-17.3)*1000,(APD[[1,j,3]]-23.45)*1000,Mean[Table[(((CAM1[[j,i]]*ηT)/(APD1[[j,i]]+DC))/(ηM*ηLe)),{i,1,Length[APD1[[1]]]}]]},{j,1,Length[CAM1]}];
Direct=Table[{(APD[[1,j,4]]-17.3)*1000,(APD[[1,j,3]]-23.45)*1000,Mean[Table[CAM1[[j,i]],{i,1,Length[APD1[[1]]]}]]},{j,1,Length[CAM1]}];

MatrixForm[Partition[η,APD[[1,Length[APD[[1]]],1]]]];
Imagedata=Partition[η,APD[[1,Length[APD[[1]]],1]]];
Image1=ListDensityPlot[Imagedata,ColorFunction->"GreenPinkTones",PlotLegends->Automatic];
ListDensityPlot[ηC,ColorFunction->"Rainbow",PlotLegends->Automatic,AspectRatio->Automatic,FrameLabel->{"X Position (μm)","Y Position (μm)"},FrameStyle->{Thick,Directive[Thick,Thick]},LabelStyle->{FontSize->19,Black,FontFamily->"Helvetica"},FrameTicksStyle->Black,PlotLegends->Placed[BarLegend[Automatic,LegendLabel->"η",LabelStyle->{FontFamily->"Helvetica",FontSize->19,Black}],Right],InterpolationOrder->0,CoordinatesToolOptions->{"DisplayFunction"->Function[pt,Row[{"point=",NumberForm[pt,6]}]],"CopiedValueFunction"->Function[pt,pt]}]