diff --git a/src/Reconstruction.jl b/src/Reconstruction.jl index 9eee13a..4aa3879 100644 --- a/src/Reconstruction.jl +++ b/src/Reconstruction.jl @@ -215,7 +215,7 @@ function reconstruction(S, bSF::Union{T,Vector{T}}, bMeas::MPIFile, freq::Array, sparseTrafo = nothing, loadasreal = false, maxload = 100, maskDFFOV=false, weightType=WeightingType.None, weightingLimit = 0, solver = "kaczmarz", spectralCleaning=true, spectralLeakageCorrection=spectralCleaning, - fgFrames=1:10, bgCorrectionInternal=false, + fgFrames=1:10, bgCorrectionInternal=false, tfCorrectionInternal=rxHasTransferFunction(bMeas), noiseFreqThresh=0.0, channelWeights=ones(3), numPeriodAverages=1, numPeriodGrouping=1, kargs...) where {T<:MPIFile} @@ -231,10 +231,10 @@ function reconstruction(S, bSF::Union{T,Vector{T}}, bMeas::MPIFile, freq::Array, #else uEmpty = getMeasurementsFD(emptyMeas, frequencies=freq, frames=bgFrames, numAverages=length(bgFrames), loadasreal = loadasreal,spectralLeakageCorrection=spectralLeakageCorrection, bgCorrection=bgCorrectionInternal, - numPeriodAverages=numPeriodAverages, numPeriodGrouping=numPeriodGrouping) + tfCorrection=tfCorrectionInternal, numPeriodAverages=numPeriodAverages, numPeriodGrouping=numPeriodGrouping) if bgFramesPost != nothing uEmptyPost = getMeasurementsFD(emptyMeas, false, frequencies=freq, frames=bgFramesPost, - numAverages = length(bgFramesPost), bgCorrection=bgCorrectionInternal, + numAverages = length(bgFramesPost), bgCorrection=bgCorrectionInternal, tfCorrection=tfCorrectionInternal, loadasreal = loadasreal, spectralLeakageCorrection=spectralLeakageCorrection, numPeriodAverages=numPeriodAverages, numPeriodGrouping=numPeriodGrouping) end @@ -264,7 +264,8 @@ function reconstruction(S, bSF::Union{T,Vector{T}}, bMeas::MPIFile, freq::Array, @debug "Loading measurements ..." u = getMeasurementsFD(bMeas, frequencies=freq, frames=partframes, numAverages=numAverages, loadasreal=loadasreal, spectralLeakageCorrection=spectralLeakageCorrection, - bgCorrection=bgCorrectionInternal, numPeriodAverages=numPeriodAverages, numPeriodGrouping=numPeriodGrouping) + bgCorrection=bgCorrectionInternal, tfCorrection=tfCorrectionInternal, + numPeriodAverages=numPeriodAverages, numPeriodGrouping=numPeriodGrouping) if emptyMeas!=nothing if bgFramesPost == nothing u = u .- uEmpty