merge code for pooling figure

Author: mdeff

figure_pooling.ipynb

@@ -68,13 +68,6 @@
     "#     pl.savefig('test_fig1.pdf', bbox_inches='tight', transparent=True);"
    ]
   },
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": [
-    "# make figure"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -100,8 +93,7 @@
     "map_select_4[map_select_4!=hp.UNSEEN] = np.arange(n_filled, dtype=np.float64)/n_filled\n",
     "\n",
     "# map_select_2[map_select_2==hp.UNSEEN] = np.nan\n",
-    "# map_select_4[map_select_4==hp.UNSEEN] = np.nan\n",
-    "\n"
+    "# map_select_4[map_select_4==hp.UNSEEN] = np.nan"
    ]
   },
   {
@@ -129,7 +121,6 @@
     "filename_plot = os.path.join(pathfig,'figure_pooling1.svg')\n",
     "pl.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
     "\n",
-    "\n",
     "# pl.axes(ax[1])\n",
     "pl.figure()\n",
     "# newax = pl.gcf().add_axes([0.1, 0.4, 0.8, 0.8],frameon=False)\n",
@@ -149,7 +140,6 @@
     "# , sub=(ny, nx, 3)\n",
     "# pl.subplots_adjust(wspace=0.01, hspace=0.01)\n",
     "\n",
-    "\n",
     "# import matplotlib.pyplot as plt\n",
     "\n",
     "# newax = pl.gcf().add_axes([0.1, 0.1, 0.8, 0.8],frameon=False)\n",
@@ -163,7 +153,7 @@
     "# newax.axis('equal')\n",
     "# plt.show()\n",
     "filename_plot = os.path.join(pathfig,'figure_pooling.pdf')\n",
-    "pl.savefig(filename_plot, bbox_inches='tight'); print('saved {}'.format(filename_plot))\n"
+    "pl.savefig(filename_plot, bbox_inches='tight'); print('saved {}'.format(filename_plot))"
    ]
   },
   {
@@ -210,8 +200,24 @@
     "                   )\n",
     "\n",
     "svg_fig.save(os.path.join(pathfig, \"figure_pooling_svg.svg\"))\n",
-    "SVG(os.path.join(pathfig, 'figure_pooling_svg.svg'))\n",
-    "# svg_fig.save(\"figure_pooling_svg.pdf\")"
+    "SVG(os.path.join(pathfig, 'figure_pooling_svg.svg'))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "map_select_8 = hp.ud_grade(map_select_1, nside_out=8, order_in='NESTED', order_out='NESTED')\n",
+    "n_filled = float(np.count_nonzero(map_select_8!=hp.UNSEEN))\n",
+    "map_select_8[map_select_8!=hp.UNSEEN] = np.arange(n_filled, dtype=np.float64)/n_filled\n",
+    "\n",
+    "make_ball(map_select_8, cmap=plt.cm.RdYlBu)\n",
+    "filename_plot = os.path.join(pathfig,'figure_pooling4.svg')\n",
+    "plt.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
+    "filename_plot = os.path.join(pathfig,'figure_pooling4.pdf')\n",
+    "plt.savefig(filename_plot, bbox_inches='tight'); print('saved {}'.format(filename_plot))"
    ]
   }
  ],

figures.ipynb

@@ -1475,135 +1475,6 @@
     "plt.plot(x, np.sum(np.array(p[0:o])**2,axis=0))\n",
     "plt.plot(x, (o+0.5)/2*(x**0))"
    ]
-  },
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": [
-    "# Figure pooling"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "def make_ball(map_test1, cmap=plt.cm.gray_r, sub=None):\n",
-    "    cmap.set_under('w')\n",
-    "    cmap.set_bad('lightgray')\n",
-    "    dot_size=10\n",
-    "    rot = (0,30,345)\n",
-    "    vmin, vmax = -0.5, 1.5\n",
-    "    hp.visufunc.orthview(map=map_test1, half_sky=True, title='', rot=rot, cmap=cmap, cbar=False, hold=True, nest=True, min=vmin, max=vmax, notext=True, sub=sub);\n",
-    "    theta, phi = hp.pix2ang(hp.npix2nside(len(map_test1)), range(len(map_test1)), nest=True);\n",
-    "    hp.projscatter(theta, phi, c='k', s=dot_size);\n",
-    "    hp.graticule();\n",
-    "    hp.graticule(dmer=360,dpar=360,alpha=1, rot=(0,0,15), local=True);\n",
-    "    hp.graticule(dmer=360,dpar=360,alpha=1, rot=(0,0,195), local=True);      "
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "npix = hp.nside2npix(16)\n",
-    "map_data = np.arange(npix)\n",
-    "\n",
-    "\n",
-    "npix = hp.nside2npix(1)\n",
-    "map_select_1 = (np.arange(npix)==0).astype(np.float)\n",
-    "map_select_1[map_select_1==0] = hp.UNSEEN\n",
-    "map_select_2 = hp.ud_grade(map_select_1, nside_out=2, order_in='NESTED', order_out='NESTED')\n",
-    "map_select_4 = hp.ud_grade(map_select_1, nside_out=4, order_in='NESTED', order_out='NESTED')\n",
-    "map_select_8 = hp.ud_grade(map_select_1, nside_out=8, order_in='NESTED', order_out='NESTED')\n",
-    "\n",
-    "n_filled = float(np.count_nonzero(map_select_2!=hp.UNSEEN))\n",
-    "map_select_2[map_select_2!=hp.UNSEEN] = np.arange(n_filled, dtype=np.float64)/n_filled\n",
-    "n_filled = float(np.count_nonzero(map_select_4!=hp.UNSEEN))\n",
-    "map_select_4[map_select_4!=hp.UNSEEN] = np.arange(n_filled, dtype=np.float64)/n_filled\n",
-    "n_filled = float(np.count_nonzero(map_select_8!=hp.UNSEEN))\n",
-    "map_select_8[map_select_8!=hp.UNSEEN] = np.arange(n_filled, dtype=np.float64)/n_filled\n"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "nx, ny = 3,1\n",
-    "\n",
-    "\n",
-    "plt.figure()\n",
-    "\n",
-    "make_ball(map_select_1, cmap=plt.cm.RdYlBu)\n",
-    "filename_plot = os.path.join(pathfig,'figure_pooling1.svg')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
-    "\n",
-    "\n",
-    "\n",
-    "plt.figure()\n",
-    "make_ball(map_select_2, cmap=plt.cm.RdYlBu)\n",
-    "filename_plot = os.path.join(pathfig,'figure_pooling2.svg')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
-    "\n",
-    "\n",
-    "plt.figure()\n",
-    "\n",
-    "make_ball(map_select_4, cmap=plt.cm.RdYlBu)\n",
-    "filename_plot = os.path.join(pathfig,'figure_pooling3.svg')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
-    "\n"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "lw=3\n",
-    "point2=(377,87)\n",
-    "point3=(603,72)\n",
-    "\n",
-    "svg_fig = sc.Figure(\"22.5cm\", \"6.75cm\",  \n",
-    "                    sc.SVG(os.path.join(pathfig, \"figure_pooling3.svg\")).scale(0.83).move(0,0),\n",
-    "                    sc.SVG(os.path.join(pathfig, \"figure_pooling2.svg\")).scale(0.83).move(255, 0),\n",
-    "                    sc.SVG(os.path.join(pathfig, \"figure_pooling1.svg\")).scale(0.83).move(510, 0),\n",
-    "                    sc.Line( ((122,80), point2), width=lw, color='darkred'),\n",
-    "                    sc.Line( ((148,93), point2), width=lw, color='darkred'),\n",
-    "                    sc.Line( ((153,71), point2), width=lw, color='darkred'),\n",
-    "                    sc.Line( ((124,107), point2), width=lw, color='darkred'),\n",
-    "                    \n",
-    "                    sc.Line( ((384,87), point3), width=lw, color='darkblue'),\n",
-    "                    sc.Line( ((396,44), point3), width=lw, color='darkblue'),\n",
-    "                    sc.Line( ((348,55), point3), width=lw, color='darkblue'),\n",
-    "                    sc.Line( ((337,105), point3), width=lw, color='darkblue')\n",
-    "                   )\n",
-    "\n",
-    "svg_fig.save(os.path.join(pathfig, \"figure_pooling_svg.svg\"))\n",
-    "ipd.SVG(os.path.join(pathfig, 'figure_pooling_svg.svg'))\n",
-    "filename_plot = os.path.join(pathfig,'figure_pooling.pdf')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight'); print('saved {}'.format(filename_plot))"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "plt.figure()\n",
-    "\n",
-    "make_ball(map_select_8, cmap=plt.cm.RdYlBu)\n",
-    "filename_plot = os.path.join(pathfig,'figure_pooling4.svg')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight', transparent=True); print('saved {}'.format(filename_plot))\n",
-    "filename_plot = os.path.join(pathfig,'HEALpix.pdf')\n",
-    "plt.savefig(filename_plot, bbox_inches='tight'); print('saved {}'.format(filename_plot))"
-   ]
   }
  ],
  "metadata": {