Merge pull request #97 from geoscixyz/gpr_add_notebooks

Gpr add notebooks

Author: sgkang

Notebooks/GPR_Attenuation.ipynb

@@ -99,21 +99,21 @@
  "metadata": {
   "anaconda-cloud": {},
   "kernelspec": {
-   "display_name": "Python [conda root]",
+   "display_name": "Python 3",
    "language": "python",
-   "name": "conda-root-py"
+   "name": "python3"
   },
   "language_info": {
    "codemirror_mode": {
     "name": "ipython",
-    "version": 2
+    "version": 3
    },
    "file_extension": ".py",
    "mimetype": "text/x-python",
    "name": "python",
    "nbconvert_exporter": "python",
-   "pygments_lexer": "ipython2",
-   "version": "2.7.13"
+   "pygments_lexer": "ipython3",
+   "version": "3.6.6"
   },
   "widgets": {
    "state": {

Notebooks/GPR_Lab6_FitData.ipynb

@@ -28,10 +28,8 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
-   "metadata": {
-    "collapsed": true
-   },
+   "execution_count": 7,
+   "metadata": {},
    "outputs": [],
    "source": [
     "%matplotlib inline\n",
@@ -59,25 +57,30 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 8,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "94b533cad0a54e598df4e27ade6ebe50"
-      }
+       "model_id": "0fa2798273f4471990b12b305ecbb46e",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "interactive(children=(IntSlider(value=50, description='epsr'), FloatSlider(value=1.0, description='h', max=2.0…"
+      ]
      },
      "metadata": {},
      "output_type": "display_data"
     },
     {
      "data": {
       "text/plain": [
-       "<function gpgLabs.GPR.GPRlab1.PipeWidgetFcn>"
+       "<function gpgLabs.GPR.GPRlab1.PipeWidgetFcn(epsr, h, xc, r, imgcmp)>"
       ]
      },
-     "execution_count": 2,
+     "execution_count": 8,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -108,27 +111,32 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 9,
    "metadata": {
-    "scrolled": true
+    "scrolled": false
    },
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "9b6c195504c14ea38c09ac60dd1cc3c6"
-      }
+       "model_id": "ac5d5481f7994f93b54e634da9f0e920",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "interactive(children=(IntSlider(value=50, description='epsr'), FloatSlider(value=1.0, description='h', max=2.0…"
+      ]
      },
      "metadata": {},
      "output_type": "display_data"
     },
     {
      "data": {
       "text/plain": [
-       "<function gpgLabs.GPR.GPRlab1.WallWidgetFcn>"
+       "<function gpgLabs.GPR.GPRlab1.WallWidgetFcn(epsr, h, x1, x2, imgcmp)>"
       ]
      },
-     "execution_count": 4,
+     "execution_count": 9,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -142,31 +150,29 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": []
   }
  ],
  "metadata": {
   "anaconda-cloud": {},
   "kernelspec": {
-   "display_name": "Python [conda root]",
+   "display_name": "Python 3",
    "language": "python",
-   "name": "conda-root-py"
+   "name": "python3"
   },
   "language_info": {
    "codemirror_mode": {
     "name": "ipython",
-    "version": 2
+    "version": 3
    },
    "file_extension": ".py",
    "mimetype": "text/x-python",
    "name": "python",
    "nbconvert_exporter": "python",
-   "pygments_lexer": "ipython2",
-   "version": "2.7.13"
+   "pygments_lexer": "ipython3",
+   "version": "3.6.6"
   },
   "widgets": {
    "state": {

Notebooks/GPR_TBL4_DOI_Resolution.ipynb

@@ -22,7 +22,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": 4,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -54,13 +54,36 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 12,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "1.3416407864998738"
+      ]
+     },
+     "execution_count": 12,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "fc = 250*1e6\n",
+    "d = 6\n",
+    "v = 3*1e8 / np.sqrt(4)\n",
+    "np.sqrt(v*d / (2*fc))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 14,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "939fe164aae44aed869231a9df22d4e9",
+       "model_id": "21d13f084aaa48768dba032d4db5f9fe",
        "version_major": 2,
        "version_minor": 0
       },
@@ -77,7 +100,7 @@
        "<function gpgLabs.GPR.GPR_zero_offset.GPRWidgetWaveRegime(sig, epsr, fc, x1, d1, R1, x2, d2, R2)>"
       ]
      },
-     "execution_count": 2,
+     "execution_count": 14,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -107,13 +130,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 15,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "application/vnd.jupyter.widget-view+json": {
-       "model_id": "89228aede5c54d3d8423e8244b46fa01",
+       "model_id": "fbc99ffb253245358914d26a9d865717",
        "version_major": 2,
        "version_minor": 0
       },
@@ -130,14 +153,21 @@
        "<function gpgLabs.GPR.Attenuation.WaveVelandSkindWidgetTBL(epsr, sigma)>"
       ]
      },
-     "execution_count": 3,
+     "execution_count": 15,
      "metadata": {},
      "output_type": "execute_result"
     }
    ],
    "source": [
     "AttenuationWidgetTBL()"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {

Notebooks/index.ipynb

@@ -38,8 +38,8 @@
     "- <a href=\"Seis_VerticalResolution.ipynb\">Seis_VerticalResolution.ipynb</a> - Vertical resolution in reflection\n",
     "\n",
     "### Ground Penetrating Radar\n",
-    "- <a href=\"GPR_Attenuation.ipynb\">GPR_Attenuation.ipynb</a> - EM attenuation in GPR\n",
-    "- <a href=\"GPR_FitData.ipynb\">GPR_FitData.ipynb</a> - Fit field data\n",
+    "- <a href=\"GPR_TBL4_DOI_Resolution.ipynb\">GPR_TBL4_DOI_Resolution.ipynb</a> - Horizontal resolution + Probing distance of GPR\n",
+    "- <a href=\"GPR_Lab6_FitData.ipynb\">GPR_Lab6_FitData.ipynb</a> - Fit field GPR data\n",
     "\n",
     "### Electromagnetics\n",
     "- <a href=\"EM_ThreeLoopModel.ipynb\">EM_ThreeLoopModel.ipynb</a> - EM induction explained by a 3-loop circuit model\n",

gpgLabs/Mag/MagDipoleApp.py

@@ -345,6 +345,7 @@ def plot_map(self):
             ax2.remove()
         else:
             ax1.plot(self.xy_profile[:, 0], self.xy_profile[:, 1], 'w')
+
     def get_half_width(self, n_points=200):
         ind_max = np.argmax(abs(self.data_profile))
         A_half = self.data_profile[ind_max] / 2.