|
41 | 41 | ], |
42 | 42 | "metadata": {} |
43 | 43 | }, |
| 44 | + { |
| 45 | + "outputs": [], |
| 46 | + "cell_type": "code", |
| 47 | + "source": [ |
| 48 | + "using SEQ_BRUKER_a_MP2RAGE_CS_360.MRIReco\n", |
| 49 | + "using SEQ_BRUKER_a_MP2RAGE_CS_360.MRIReco.RegularizedLeastSquares" |
| 50 | + ], |
| 51 | + "metadata": {}, |
| 52 | + "execution_count": null |
| 53 | + }, |
44 | 54 | { |
45 | 55 | "cell_type": "markdown", |
46 | 56 | "source": [ |
47 | | - "## Loading Package" |
| 57 | + "## Download the datasets\n", |
| 58 | + "if you run the literate example offline change the following line by : `MP2_artifacts = artifact\"MP2RAGE_data\"" |
48 | 59 | ], |
49 | 60 | "metadata": {} |
50 | 61 | }, |
51 | 62 | { |
52 | 63 | "outputs": [], |
53 | 64 | "cell_type": "code", |
54 | 65 | "source": [ |
55 | | - "using LazyArtifacts # loading data\n", |
56 | | - "using SEQ_BRUKER_a_MP2RAGE_CS_360\n", |
57 | | - "using CairoMakie # plotting\n", |
58 | | - "\n", |
59 | 66 | "datadir = Main.MP2_artifacts\n", |
60 | 67 | "@info \"The test data is located at $datadir.\"" |
61 | 68 | ], |
|
86 | 93 | ], |
87 | 94 | "metadata": {} |
88 | 95 | }, |
89 | | - { |
90 | | - "outputs": [], |
91 | | - "cell_type": "code", |
92 | | - "source": [ |
93 | | - "using SEQ_BRUKER_a_MP2RAGE_CS_360.MRIReco\n", |
94 | | - "using SEQ_BRUKER_a_MP2RAGE_CS_360.MRIReco.RegularizedLeastSquares" |
95 | | - ], |
96 | | - "metadata": {}, |
97 | | - "execution_count": null |
98 | | - }, |
99 | 96 | { |
100 | 97 | "cell_type": "markdown", |
101 | 98 | "source": [ |
|
121 | 118 | { |
122 | 119 | "cell_type": "markdown", |
123 | 120 | "source": [ |
124 | | - "for comparison purpose let's perform the undersampled reconstruction (without the paramCS keyword)" |
| 121 | + "for comparison purposes let's perform the undersampled reconstruction (without the paramCS keyword)" |
125 | 122 | ], |
126 | 123 | "metadata": {} |
127 | 124 | }, |
|
148 | 145 | "begin\n", |
149 | 146 | " f = Figure(size=(500,400))\n", |
150 | 147 | " ax=Axis(f[1,1],title=\"TI₁ undersampled\")\n", |
151 | | - " h=heatmap!(ax,abs.(d_under[\"im_reco\"][:,:,60,1,1,1]),colormap=:grays)\n", |
| 148 | + " h=heatmap!(ax,abs.(d_under[\"im_reco\"][:,:,60,1,1]),colormap=:grays)\n", |
152 | 149 | "\n", |
153 | 150 | " ax=Axis(f[1,2],title=\"TI₁ CS\")\n", |
154 | | - " h=heatmap!(ax,abs.(d[\"im_reco\"][:,:,60,1,1,1]),colormap=:grays)\n", |
| 151 | + " h=heatmap!(ax,abs.(d[\"im_reco\"][:,:,60,1,1]),colormap=:grays)\n", |
155 | 152 | "\n", |
156 | 153 | "\n", |
157 | | - " ax=Axis(f[2,1],title=\"UNIT1 undersampled\")\n", |
158 | | - " h=heatmap!(ax,d_under[\"T1map\"][:,:,60,1,1],colorrange = (500,2000))\n", |
| 154 | + " ax=Axis(f[2,1],title=\"T₁ map undersampled\")\n", |
| 155 | + " h=heatmap!(ax,d_under[\"T1map\"][:,:,60,1],colorrange = (500,2000))\n", |
159 | 156 | "\n", |
160 | | - " ax=Axis(f[2,2],title=\"UNIT1 CS\")\n", |
161 | | - " h=heatmap!(ax,d[\"T1map\"][:,:,60,1,1],colorrange = (500,2000))\n", |
| 157 | + " ax=Axis(f[2,2],title=\"T₁ map CS\")\n", |
| 158 | + " h=heatmap!(ax,d[\"T1map\"][:,:,60,1],colorrange = (500,2000))\n", |
162 | 159 | "\n", |
163 | 160 | " for ax in f.content # hide decoration befor adding colorbar\n", |
164 | 161 | " hidedecorations!(ax)\n", |
|
0 commit comments