some bug need to be fixed

Author: alexlib

pyptv/calibration_with_particles.ipynb

@@ -368,8 +368,7 @@
     "\n",
     "def _residuals(cal):\n",
     "    err = 0\n",
-    "    for i_frame in range(n_frames):\n",
-    "        # fig, ax = plt.subplots()\n",
+    "    for i_frame in range(n_frames):        \n",
     "        for i_point in range(len(points[i_frame])):\n",
     "            if ids[i_frame][i_point][i_cam] != -1:\n",
     "                tmp = convert_arr_metric_to_pixel(\n",
@@ -406,7 +405,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 54,
    "metadata": {},
    "outputs": [
     {
@@ -454,6 +453,85 @@
     "sol = minimize(error_function, x0, args=(cal), method='Nelder-Mead', tol=1e-4)\n"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 55,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "[ 1.29864389e+02  7.09995178e+01  6.06405459e+02 -1.18822423e-01\n",
+      "  2.34215634e-01  9.15198506e-03 -1.76412191e-01 -9.97652573e-02\n",
+      "  7.41781981e+01  3.41659939e-02  1.36738065e-01  3.02279045e-02\n",
+      "  2.17198498e-03 -1.57233580e-02  1.88667731e+01  2.95341468e-02]\n"
+     ]
+    }
+   ],
+   "source": [
+    "print(sol.x)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 58,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "80032.63451606502"
+      ]
+     },
+     "execution_count": 58,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "error_function(sol.x, cal)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [
+    {
+     "ename": "NameError",
+     "evalue": "name 'cals' is not defined",
+     "output_type": "error",
+     "traceback": [
+      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
+      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
+      "Cell \u001b[0;32mIn[1], line 3\u001b[0m\n\u001b[1;32m      1\u001b[0m i_frame \u001b[38;5;241m=\u001b[39m \u001b[38;5;241m0\u001b[39m  \n\u001b[1;32m      2\u001b[0m i_cam \u001b[38;5;241m=\u001b[39m \u001b[38;5;241m0\u001b[39m\n\u001b[0;32m----> 3\u001b[0m cal \u001b[38;5;241m=\u001b[39m \u001b[43mcals\u001b[49m[i_cam] \n\u001b[1;32m      4\u001b[0m array_to_calibration(sol\u001b[38;5;241m.\u001b[39mx, cal)\n\u001b[1;32m      6\u001b[0m fig, ax \u001b[38;5;241m=\u001b[39m plt\u001b[38;5;241m.\u001b[39msubplots()\n",
+      "\u001b[0;31mNameError\u001b[0m: name 'cals' is not defined"
+     ]
+    }
+   ],
+   "source": [
+    "i_frame = 0  \n",
+    "i_cam = 0\n",
+    "cal = cals[i_cam] \n",
+    "# array_to_calibration(sol.x, cal)\n",
+    "\n",
+    "fig, ax = plt.subplots()\n",
+    "for i_point in range(len(points[i_frame])):\n",
+    "    if ids[i_frame][i_point][i_cam] != -1:\n",
+    "        tmp = convert_arr_metric_to_pixel(\n",
+    "                    image_coordinates(np.atleast_2d(points[i_frame][i_point]), cal, cpar.get_multimedia_params()),\n",
+    "                    cpar\n",
+    "                )\n",
+    "        x, y  = targets[i_frame][i_cam][ids[i_frame][i_point][i_cam]].pos()\n",
+    "\n",
+    "\n",
+    "        ax.scatter(x, y, marker='+', color = 'blue')\n",
+    "        ax.scatter(tmp[0,0], tmp[0,1], marker='x', color = 'red')\n",
+    "\n",
+    "fig.show()           \n"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,