deploy: 9c3533c

LarsonLab · Nov 11, 2023 · 3f3a57e · 3f3a57e
1 parent efbb23a
commit 3f3a57e
Show file tree

Hide file tree

Showing 10 changed files with 61 additions and 53 deletions.
diff --git a/FOV and Resolution.html b/FOV and Resolution.html
diff --git a/_images/2a9e6a30cdf4e44bb913ba369e05776ffc6f98379cadde43502c1c4d5f326d41.png b/_images/2a9e6a30cdf4e44bb913ba369e05776ffc6f98379cadde43502c1c4d5f326d41.png
diff --git a/_images/3a2fc437b102fd207199631270b4ffd044cced66ee12d57212f8f268223da2fe.png b/_images/3a2fc437b102fd207199631270b4ffd044cced66ee12d57212f8f268223da2fe.png
diff --git a/...0bae952106681fc0bb0f06f2188c5a960754e.png → ...2e2e753e94a22c3b96d6343334b023bf8829a.png b/...0bae952106681fc0bb0f06f2188c5a960754e.png → ...2e2e753e94a22c3b96d6343334b023bf8829a.png
diff --git a/_images/4a80e1b5fc462784ae4b1d85c570bfa188869b427e7ff49274ff63a36cf6dc08.png b/_images/4a80e1b5fc462784ae4b1d85c570bfa188869b427e7ff49274ff63a36cf6dc08.png
diff --git a/_images/92ebb64b5ae058fa6dee8f06c18ea08381b27c2bb5ac926187082f4d196c58e9.png b/_images/92ebb64b5ae058fa6dee8f06c18ea08381b27c2bb5ac926187082f4d196c58e9.png
diff --git a/_images/9c5a73a49908a52a7c0d85e6b12641f21a2d0b345536f6a533be3e399a0d7fb1.png b/_images/9c5a73a49908a52a7c0d85e6b12641f21a2d0b345536f6a533be3e399a0d7fb1.png
diff --git a/_images/b7030ca9f2508ba6809e1f106121490f593d9ee006fb4141fe38522078613edd.png b/_images/b7030ca9f2508ba6809e1f106121490f593d9ee006fb4141fe38522078613edd.png
diff --git a/_sources/FOV and Resolution.ipynb b/_sources/FOV and Resolution.ipynb
@@ -67,14 +67,14 @@
     "% rectangular object to demonstrate FOV\n",
     "\n",
     "N = 256;\n",
-    "kx = [-N/2:N/2-1]/N;\n",
+    "kx = [-N/2+1:N/2]/N;\n",
     "N_rect = N/4;\n",
     "kdata = sinc(kx *N_rect).' * sinc(kx *N_rect);\n",
     "\n",
     "rect_recon = ifft2c(kdata);\n",
     "\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata)), [0 (1)])\n",
+    "imagesc(log(abs(kdata)), [0 max(log(abs(kdata(:)))))\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
     "imagesc(abs(rect_recon), [0 max(abs(rect_recon(:)))])\n",
@@ -88,23 +88,23 @@
     "\n",
     "figure\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata2)), [0 (1)])\n",
+    "imagesc(log(abs(kdata2)), [0 max(log(abs(kdata2(:))))])\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
     "imagesc(abs(rect_recon2), [0 max(abs(rect_recon2(:)))])\n",
     "colormap(gray), axis equal tight off\n",
     "title('Doubled \\Delta k, Half FOV')\n",
     "\n",
     "kdata3 = kdata;\n",
-    "kdata2(1:3:end,:) = 0;\n",
-    "kdata2(2:3:end,:) = 0;\n",
-    "kdata2(:,1:3:end) = 0;\n",
-    "kdata2(:,2:3:end) = 0;\n",
+    "kdata3(1:3:end,:) = 0;\n",
+    "kdata3(2:3:end,:) = 0;\n",
+    "kdata3(:,1:3:end) = 0;\n",
+    "kdata3(:,2:3:end) = 0;\n",
     "rect_recon3 = ifft2c(kdata3);\n",
     "\n",
     "figure\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata3)), [0 (1)])\n",
+    "imagesc(log(abs(kdata3)), [0 max(log(abs(kdata(:))))])\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
     "imagesc(abs(rect_recon3), [0 max(abs(rect_recon3(:)))])\n",
@@ -138,47 +138,47 @@
     "% rectangular object to demonstrate resolution\n",
     "\n",
     "N = 256;\n",
-    "kx = [-N/2:N/2-1]/N;\n",
+    "kx = [-N/2+1:N/2]/N;\n",
     "N_rect = N/2;\n",
     "kdata = sinc(kx *N_rect).' * sinc(kx *N_rect);\n",
     "\n",
     "rect_recon = ifft2c(kdata);\n",
     "\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata)), [0 (1)])\n",
+    "imagesc(log(abs(kdata)), [0 max(log(abs(kdata(:))))])\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
     "imagesc(abs(rect_recon), [0 max(abs(rect_recon(:)))])\n",
     "colormap(gray), axis equal tight off\n",
     "title('Full Resolution')\n",
     "\n",
-    "kdata2 = kdata;\n",
-    "kdata2([[1:N/4],[3*N/4+1:N]],:) = 0;\n",
-    "kdata2(:,[[1:N/4],[3*N/4+1:N]]) = 0;\n",
-    "rect_recon2 = ifft2c(kdata2);\n",
+    "kdata4 = kdata;\n",
+    "kdata4([[1:3*N/8],[5*N/8+1:N]],:) = 0;\n",
+    "kdata4(:,[[1:3*N/8],[5*N/8+1:N]]) = 0;\n",
+    "rect_recon4 = ifft2c(kdata4);\n",
     "\n",
     "figure\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata2)), [0 (1)])\n",
+    "imagesc(log(abs(kdata4)), [0 max(log(abs(kdata4(:))))])\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
-    "imagesc(abs(rect_recon2), [0 max(abs(rect_recon2(:)))])\n",
+    "imagesc(abs(rect_recon4), [0 max(abs(rect_recon4(:)))])\n",
     "colormap(gray), axis equal tight off\n",
-    "title('Half k_{max}, Double voxel size')\n",
+    "title('1/4 k_{max}, 4x voxel size')\n",
     "\n",
-    "kdata4 = kdata;\n",
-    "kdata4([[1:3*N/8],[5*N/8+1:N]],:) = 0;\n",
-    "kdata4(:,[[1:3*N/8],[5*N/8+1:N]]) = 0;\n",
-    "rect_recon4 = ifft2c(kdata4);\n",
+    "kdata8 = kdata;\n",
+    "kdata8([[1:7*N/16],[9*N/16+1:N]],:) = 0;\n",
+    "kdata8(:,[[1:7*N/16],[9*N/16+1:N]]) = 0;\n",
+    "rect_recon8 = ifft2c(kdata8);\n",
     "\n",
     "figure\n",
     "subplot(221)\n",
-    "imagesc(log(abs(kdata4)), [0 (1)])\n",
+    "imagesc(log(abs(kdata8)), [0 max(log(abs(kdata8(:))))])\n",
     "colormap(gray), axis equal tight off\n",
     "subplot(222)\n",
-    "imagesc(abs(rect_recon4), [0 max(abs(rect_recon4(:)))])\n",
+    "imagesc(abs(rect_recon8), [0 max(abs(rect_recon8(:)))])\n",
     "colormap(gray), axis equal tight off\n",
-    "title('1/4 k_{max}, 4x voxel size')\n",
+    "title('1/8 k_{max}, 8x voxel size')\n",
     "\n"
    ]
   },

diff --git a/searchindex.js b/searchindex.js