From 18a561783653b41b6eab530b700ec28bd0699cdb Mon Sep 17 00:00:00 2001 From: "Documenter.jl" Date: Thu, 10 Oct 2024 12:00:35 +0000 Subject: [PATCH] build based on 5c8123e --- dev/.documenter-siteinfo.json | 2 +- dev/basicReconstruction.html | 2 +- dev/datasetStore.html | 2 +- dev/index.html | 2 +- dev/matrixCompression.html | 2 +- dev/multiContrast.html | 2 +- dev/multiPatch.html | 2 +- dev/overview.html | 2 +- dev/parameters.html | 2 +- dev/recoResults.html | 2 +- 10 files changed, 10 insertions(+), 10 deletions(-) diff --git a/dev/.documenter-siteinfo.json b/dev/.documenter-siteinfo.json index 82fd406..295778d 100644 --- a/dev/.documenter-siteinfo.json +++ b/dev/.documenter-siteinfo.json @@ -1 +1 @@ -{"documenter":{"julia_version":"1.10.5","generation_timestamp":"2024-09-10T10:02:40","documenter_version":"1.7.0"}} \ No newline at end of file +{"documenter":{"julia_version":"1.11.0","generation_timestamp":"2024-10-10T12:00:31","documenter_version":"1.7.0"}} \ No newline at end of file diff --git a/dev/basicReconstruction.html b/dev/basicReconstruction.html index 4734143..6af41a1 100644 --- a/dev/basicReconstruction.html +++ b/dev/basicReconstruction.html @@ -16,4 +16,4 @@ spectralCleaning=true, fgFrames=1:10, bgCorrectionInternal=false, noiseFreqThresh=0.0, kargs...)

This function is responsible for loading the measurement data and potential background data that is subtracted from the measurements. For any frame to be reconstructed, the low level reconstruction routine is called.

Low Level Reconstruction

Finally, we have arrived at the low level reconstruction routine that has the signature

function reconstruction(S, u::Array; sparseTrafo = nothing,
                         lambd=0, progress=nothing, solver = "Kaczmarz",
-                        weights=nothing, kargs...)

One can see that it requires the system matrix S and the measurements u to be already loaded.

We note that S is typeless for a reason here. For a regular reconstruction one will basically feed in an Array{ComplexF32,2} in here, although more precisely it will be a Transposed version of that type if the Kaczmarz algorithm is being used for efficiency reasons.

However, in case that matrix compression is applied S will be of type SparseMatrixCSC. And for Multi-Patch Reconstruction S will be of type MultiPatchOperator. Hence, the solvers are implemented in a very generic way and require only certain functions to be implemented. The low level reconstruction method calls one of the solvers from RegularizedLeastSquares.jl.

+ weights=nothing, kargs...)

One can see that it requires the system matrix S and the measurements u to be already loaded.

We note that S is typeless for a reason here. For a regular reconstruction one will basically feed in an Array{ComplexF32,2} in here, although more precisely it will be a Transposed version of that type if the Kaczmarz algorithm is being used for efficiency reasons.

However, in case that matrix compression is applied S will be of type SparseMatrixCSC. And for Multi-Patch Reconstruction S will be of type MultiPatchOperator. Hence, the solvers are implemented in a very generic way and require only certain functions to be implemented. The low level reconstruction method calls one of the solvers from RegularizedLeastSquares.jl.

diff --git a/dev/datasetStore.html b/dev/datasetStore.html index e2ccd48..b55cb54 100644 --- a/dev/datasetStore.html +++ b/dev/datasetStore.html @@ -1,2 +1,2 @@ -- · MPI Reconstruction
GitHub

Layers

The reconstruction function has several layers starting from a high level over several middle layer to low layer functions. The most high level method has the following signature `$julia reconstruction(d::MDFDatasetStore, study::Study, exp::Experiment, recoParams::Dict{String,Any})$

+- · MPI Reconstruction
GitHub

Layers

The reconstruction function has several layers starting from a high level over several middle layer to low layer functions. The most high level method has the following signature `$julia reconstruction(d::MDFDatasetStore, study::Study, exp::Experiment, recoParams::Dict{String,Any})$

diff --git a/dev/index.html b/dev/index.html index 8f9331e..ea0d515 100644 --- a/dev/index.html +++ b/dev/index.html @@ -1,2 +1,2 @@ -Home · MPI Reconstruction
GitHub

MPIReco.jl

Julia package for the reconstruction of magnetic particle imaging (MPI) data

Introduction

This project provides functions for the reconstruction of MPI data. The project is implemented in the programming language Julia and contains algorithms for

Key features are

  • Frequency filtering for memory efficient reconstruction. Only frequencies used during reconstructions are loaded into memory.
  • Different solvers provided by the package RegularizedLeastSquares.jl
  • High-level until low-level reconstruction providing maximum flexibility for the user
  • Spectral leakage correction (implemented in MPIFiles.jl)

Installation

Start julia and open the package mode by entering ]. Then enter

add MPIReco

This will install the packages MPIReco.jl and all its dependencies. In particular this will install the core dependencies MPIFiles and RegularizedLeastSquares.

License / Terms of Usage

The source code of this project is licensed under the MIT license. This implies that you are free to use, share, and adapt it. However, please give appropriate credit by citing the project.

Contact

If you have problems using the software, find mistakes, or have general questions please use the issue tracker to contact us.

Contributors

+Home · MPI Reconstruction
GitHub

MPIReco.jl

Julia package for the reconstruction of magnetic particle imaging (MPI) data

Introduction

This project provides functions for the reconstruction of MPI data. The project is implemented in the programming language Julia and contains algorithms for

Key features are

  • Frequency filtering for memory efficient reconstruction. Only frequencies used during reconstructions are loaded into memory.
  • Different solvers provided by the package RegularizedLeastSquares.jl
  • High-level until low-level reconstruction providing maximum flexibility for the user
  • Spectral leakage correction (implemented in MPIFiles.jl)

Installation

Start julia and open the package mode by entering ]. Then enter

add MPIReco

This will install the packages MPIReco.jl and all its dependencies. In particular this will install the core dependencies MPIFiles and RegularizedLeastSquares.

License / Terms of Usage

The source code of this project is licensed under the MIT license. This implies that you are free to use, share, and adapt it. However, please give appropriate credit by citing the project.

Contact

If you have problems using the software, find mistakes, or have general questions please use the issue tracker to contact us.

Contributors

diff --git a/dev/matrixCompression.html b/dev/matrixCompression.html index f908566..c852e10 100644 --- a/dev/matrixCompression.html +++ b/dev/matrixCompression.html @@ -1,2 +1,2 @@ -Compression · MPI Reconstruction
GitHub

Matrix-Compression Techniques

The reconstruction can be accelerated by applying matrix compression. To this end, the system matrix S is transformed into a different domain by applying a basis transformation on the rows of the system matrix. In MPIReco.jl, matrix compression can be enabled by specifying sparseTrafo which can be "DCT-IV" or "FFT".

The transformations can be restricted to the drive-field field-of-view by setting useDFFoV = true. The compression factor that controls how many coefficients are dropped after application of the transformation is controlled by the parameter redFactor. For instance a reduction factor of redFactor = 0.01 will drop 99 % of the data.

+Compression · MPI Reconstruction
GitHub

Matrix-Compression Techniques

The reconstruction can be accelerated by applying matrix compression. To this end, the system matrix S is transformed into a different domain by applying a basis transformation on the rows of the system matrix. In MPIReco.jl, matrix compression can be enabled by specifying sparseTrafo which can be "DCT-IV" or "FFT".

The transformations can be restricted to the drive-field field-of-view by setting useDFFoV = true. The compression factor that controls how many coefficients are dropped after application of the transformation is controlled by the parameter redFactor. For instance a reduction factor of redFactor = 0.01 will drop 99 % of the data.

diff --git a/dev/multiContrast.html b/dev/multiContrast.html index a785003..313058c 100644 --- a/dev/multiContrast.html +++ b/dev/multiContrast.html @@ -2,4 +2,4 @@ Multi-Contrast · MPI Reconstruction
GitHub

Multi-Contrast Reconstruction

Until now we have discussed single-contrast reconstruction in which case the reconstructed image c has a singleton first dimension. To perform multi-contrast reconstruction one has to specify multiple system matrices

bSFa = MPIFile(filenameA)
 bSFb = MPIFile(filenameB)

and can then invoke

c = reconstruction([bSFa, bSFb], b;
                     SNRThresh=5, frames=1, minFreq=80e3,
-                    recChannels=1:2, iterations=1)

Now one can access the first and second channel by c[1,:,:,:] and c[2,:,:,:].

+ recChannels=1:2, iterations=1)

Now one can access the first and second channel by c[1,:,:,:] and c[2,:,:,:].

diff --git a/dev/multiPatch.html b/dev/multiPatch.html index 715b57b..90d7977 100644 --- a/dev/multiPatch.html +++ b/dev/multiPatch.html @@ -18,4 +18,4 @@ c4 = reconstruction(bSFs, b; SNRThresh=5, frames=1, minFreq=80e3, recChannels=1:2,iterations=1, spectralLeakageCorrection=false, mapping=mapping, systemMatrices = S, SFGridCenter=SFGridCenter, - FFPos=FFPos, FFPosSF=FFPos) + FFPos=FFPos, FFPosSF=FFPos) diff --git a/dev/overview.html b/dev/overview.html index 99cf4d8..2b4f167 100644 --- a/dev/overview.html +++ b/dev/overview.html @@ -12,4 +12,4 @@ recChannels=1:2, iterations=1, spectralLeakageCorrection=true) -

Lets go through that script step by step. First, we create handles for the system matrix and the measurement data. Both are of the type MPIFile which is an abstract type that can for instance be an MDFFile or a BrukerFile.

Using the handles to the MPI datasets we can call the reconstruction function that has various variants depending on the types that are passed to it. Here, we exploit the multiple dispatch mechanism of julia. In addition to the file handles we also apply several reconstruction parameters using keyword arguments. In this case, we set the SNR threshold to 5 implying that only matrix rows with an SNR above 5 are used during reconstruction. The parameter frame decides which frame of the measured data should be reconstructed.

The object c is of type ImageMeta and contains not only the reconstructed data but also several metadata such as the reconstruction parameters being used. More details on the return type are discussed in the Reconstruction Results

Data Storage

One can store the reconstruction result into an MDF file by calling

saveRecoData("filename.mdf", c)

In order to load the data one calls

c = loadRecoData("filename.mdf", c)

We will next take a closer look at different forms of the reconstruction routine.

+

Lets go through that script step by step. First, we create handles for the system matrix and the measurement data. Both are of the type MPIFile which is an abstract type that can for instance be an MDFFile or a BrukerFile.

Using the handles to the MPI datasets we can call the reconstruction function that has various variants depending on the types that are passed to it. Here, we exploit the multiple dispatch mechanism of julia. In addition to the file handles we also apply several reconstruction parameters using keyword arguments. In this case, we set the SNR threshold to 5 implying that only matrix rows with an SNR above 5 are used during reconstruction. The parameter frame decides which frame of the measured data should be reconstructed.

The object c is of type ImageMeta and contains not only the reconstructed data but also several metadata such as the reconstruction parameters being used. More details on the return type are discussed in the Reconstruction Results

Data Storage

One can store the reconstruction result into an MDF file by calling

saveRecoData("filename.mdf", c)

In order to load the data one calls

c = loadRecoData("filename.mdf", c)

We will next take a closer look at different forms of the reconstruction routine.

diff --git a/dev/parameters.html b/dev/parameters.html index c55c0b6..a03cf00 100644 --- a/dev/parameters.html +++ b/dev/parameters.html @@ -1,2 +1,2 @@ -Parameters · MPI Reconstruction
GitHub
+Parameters · MPI Reconstruction
GitHub
diff --git a/dev/recoResults.html b/dev/recoResults.html index 7078c8f..b30055c 100644 --- a/dev/recoResults.html +++ b/dev/recoResults.html @@ -1,2 +1,2 @@ -Results · MPI Reconstruction
GitHub

Reconstruction Results

The object c is of type ImageMeta and contains not only the reconstructed data but also several metadata such as the reconstruction parameters being used. c has in total 5 dimensions. The first dimension encodes multi-spectral channels. Dimensions 2-4 encode the three spatial dimensions. The last dimension contains the number of frames being stored in c.

+Results · MPI Reconstruction
GitHub

Reconstruction Results

The object c is of type ImageMeta and contains not only the reconstructed data but also several metadata such as the reconstruction parameters being used. c has in total 5 dimensions. The first dimension encodes multi-spectral channels. Dimensions 2-4 encode the three spatial dimensions. The last dimension contains the number of frames being stored in c.