philipp-mr / coilgen

Open source MR coil generator

The CoilGen Project is a community-based tool for the generation of coil Layouts within the MRI/NMR environment. It is based on a boundary element method and generates an interconnected non-overlapping wire-tracks on 3D support structures. The focus of this work is post processing.

The user must specify a target field (e.g., bz(x,y,z)=y for a constant gradient in the y-direction) and a surface mesh geometry (e.g., a cylinder defined in an .stl file). The code then generates a coil layout in the form of a non-overlapping, interconnected wire trace to achieve the desired field.

pyCoilGen: A Python Implementation

The author has supported the creation of pyCoilGen, a translation of this CoilGen Project to Python. The code is available on GitHub and documentation is available on ReadTheDocs.

A full description is given in the following publication: https://onlinelibrary.wiley.com/doi/10.1002/mrm.29294

Getting Started

Check the documentation to get started

Installation

1. Download and extract the file of the CoilGen repository
2. Run one of the examples in the folder "Examples"

The project requires MATLAB and optionally FastHenry2 for calculation of the coil inductance. The MATLAB version should not be older than 2020A.

Algorithm overview

License

See LICENSE.txt for more information.

(back to top)

Contact

Philipp Amrein, [email protected]

Project Link: [https://github.com/Philipp-MR/CoilGen]

(back to top)

Citation

For citation of this work, please refer to the following publication: https://onlinelibrary.wiley.com/doi/10.1002/mrm.29294 https://doi.org/10.1002/mrm.29294

(back to top)

coilgen's Issues

calculate_inductance_matrix seems run into infinite loop.

I try to add inductance into SF Optimization to get lower inductance. However, it seems run into infinite loop while calculate_inductance_matrix.

STL File Parameterization Error

Dear Dr. Amrein,

I am experiencing some issues with a PTC CREO STL file I created. The geometry is a hollow half-sphere similar to the "half_sphere_shape" STL included in the Geometry_Data directory, but without the side rims. The MATLAB error I'm getting is the following:

Index in position 2 exceeds
array bounds.

Error in parameterize_mesh (line 33)
is_new_node=[boundary_edges(:,1)' 0]==[0 boundary_edges(:,2)'];

Error in CoilGen (line 59)
coil_parts=parameterize_mesh(coil_parts,input);

Error in test (line 20)
coil_layouts(coil_ind).out=CoilGen(...

Could you please advise me on how to resolve this issue?
Thank you so much.

Startup problems

Dear Dr. Amrein,

I would like to use your software to design some imaging gradients. I downloaded the zip, extracted it and set the top level directory and subdirectories on the MATLAB path. I set the current directory to Examples. The ygradient_coil example complained that the "iteration_num_stream_function_refinement" parameter was invalid. I couldn't see it in the parse_input source so I removed it in the example code, which ran and produced a reasonable looking coil. So that's a small issue.

The coil is larger than I need, so I made a smaller cylinder stl in blender. I'm still working on getting the script to run. I want eventually to use the output of your software to produce a Gerber file, manufacture a flexible circuit 'board', and wrap that on a cylindrical form. Do you think the software is suitable for this application? Do you have any words of wisdom about how to do this?

Thanks very much.

The output of calculate_force_and_torque_matrix is resistance_matrix, not the force and torque matrix

The output of calculate_force_and_torque_matrix is resistance_matrix, not the force and torque matrix. And the code of calculate_force_and_torque_matrix nearly same as calculate_resistance_matrix.

Wire width?

Thanks for the great code. What is the width of the wire? Is there a way to change the width? Thank you.