{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# What are LightCurve objects?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "`LightCurve` objects are data objects which encapsulate the brightness of a star over time. They provide a series of common operations, for example folding, binning, plotting, etc. There are a range of subclasses of `LightCurve` objects specific to telescopes, including `KeplerLightCurve` for Kepler and K2 data and `TessLightCurve` for TESS data.\n", "\n", "Although `lightkurve` was designed with Kepler, K2 and TESS in mind, these objects can be used for a range of astronomy data.\n", "\n", "You can create a `LightCurve` object from a `TargetPixelFile` object using Simple Aperture Photometry (see our tutorial for more information on Target Pixel Files [here](http://lightkurve.keplerscience.org/tutorials/1.02-target-pixel-files.html). Aperture Photometry is the simple act of summing up the values of all the pixels in a pre-defined aperture, as a function of time. By carefully choosing the shape of the aperture mask, you can avoid nearby contaminants or improve the strength of the specific signal you are trying to measure relative to the background.\n", "\n", "To demonstrate, lets create a `KeplerLightCurve` from a `KeplerTargetPixelFile`." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "INFO: Found cached file ./mastDownload/Kepler/kplr006922244_lc_Q111111111111111111/kplr006922244-2010078095331_lpd-targ.fits.gz with expected size 2009425. [astroquery.query]\n" ] } ], "source": [ "from lightkurve import KeplerTargetPixelFile, KeplerLightCurve\n", "#First we open a Target Pixel File from MAST, this one is already cached from our previous tutorial!\n", "tpf = KeplerTargetPixelFile.from_archive(6922244, quarter=4)\n", "\n", "#Then we convert the target pixel file into a light curve using the pipeline-defined aperture mask.\n", "lc = tpf.to_lightcurve(aperture_mask=tpf.pipeline_mask)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We've built a new `KeplerLightCurve` object called `lc`. Note in this case we've passed an **aperture_mask** to the `to_lightcurve` method. The default is to use the *Kepler* pipeline aperture. (You can pass your own aperture, which is a boolean `numpy` array.) By summing all the pixels in the aperture we have created a Simple Aperture Photometry (SAP) lightcurve.\n", "\n", "`KeplerLightCurve` has many useful functions that you can use. As with Target Pixel Files you can access the meta data very simply:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'Kepler'" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lc.mission" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "4" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lc.quarter" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And you still have access to time and flux attributes. In a light curve, there is only one flux point for every time stamp:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(array([352.37632485, 352.39675805, 352.43762445, ..., 442.16263546,\n", " 442.18306983, 442.2035041 ]),\n", " array([43689.15 , 43698.08 , 43694.105, ..., 43155.8 , 43148.465,\n", " 43151.562], dtype=float32))" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lc.time, lc.flux" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can also check the \"CDPP\" noise metric of the lightcurve using the built in method:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "407.35519723966718" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "lc.cdpp()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we can use the built in `plot` function on the `KeplerLightCurve` object to plot the time series. You can pass `plot` any keywords you would normally pass to `matplotlib.pyplot.plot`." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X14ZHV99/H3d2AtxFTBBrlWgrK3JSpQjem6IVKvItNW\nRCvVKmJRKA9SoqC1tDXojdruZV2hVmn1CkXkwWql6zNy04LGVrx7xxVYw8OqpBREoAjEytIhKLvM\n9/7j/E5yZnYeM3Nmzsx8XteVKzNnzkx+88vMfOf39P2ZuyMiIrJWuW4XQEREepsCiYiItESBRERE\nWqJAIiIiLVEgERGRliiQiIhISxRIRESkJQokIiLSEgUSERFpyd7dLkAnjIyM+CGHHNLtYtS0e/du\n9t57IP4dDVF9lFJ9rFJdlEqzPm6++eYldz+g3nkD8d845JBDuOmmm7pdjJqWlpYYGRnpdjEyQ/VR\nSvWxSnVRKs36MLN7GjlPXVsiItISBRIREWmJAomIiLREgURERFqiQCIiIi1RIBERkZYokIiISEsU\nSEREpCUKJE2anp5menq628UQEckMBRIREWmJAomIiLRkIHJttdPs7Gy3iyAikilqkYiISEsUSERE\npCUKJCIi0hIFEhERaYkCiYiItESBREREWpJaIDGzy8zsITO7vcrtZmZ/a2Z3mtmtZjaRuO1YM7sj\n3DaTOP4BM7vfzBbCz3FplV9ERBqTZovkCuDYGre/Ejg0/JwJzAKY2V7AJ8LthwFvMrPDEvf7qLuP\nh59r0yi4iIg0LrVA4u43AP9d45TjgU975DvAfma2HtgE3Onud7n7E8BV4VwREcmgbq5sPwi4N3H9\nvnCs0vHJxPVzzOxk4CbgXHf/WaUHN7MziVo6jI6OsrS01Mait9/OnTu7XYRMUX2UUn2sUl2UykJ9\n9FqKlFlgM+Dh90eA0yqd6O6XAJcAbNy40UdGRjpVxjXrhTJ2kuqjlOpjleqiVLfro5uB5H7g4MT1\n0XBsXZXjuPuD8UEz+yRwTfrFFBGRWroZSK4Gzjazq4i6rna6+wNm9jBwqJltIAogJwJ/AGBm6939\ngXD/1wIVZ4R1W/l+JUr0KCL9LLVAYmafA44GRszsPuD9RK0N3P1i4FrgOOBOYBk4Ndy228zOBq4D\n9gIuc/cd4WEvMLNxoq6tHwF/lFb522FxcZHl5WWmp6dXgsn09DSLi4uMjY1VDDBxEKoWfOrdLiLS\naakFEnd/U53bHXh7lduuJQo05cff0p7SNS6fzwMwNjZWcjwZGJLXIQogsWKxyPbt25mcnCSXyzE0\nNMTy8jILCwvk83nm5uYAmJmZYd26dTXLkgxCIiJZ0WuD7R2Vz+cpFAoALCwsUCwWARgeHq643e7U\n1NTKOZUUi8WVxwMoFApMTkYT0p761KdiZiwvL1MsFpmammJ8fHyPAFUoFFhcXKzZMlGrRUQ6SYGk\nhuXl5ZXLyQBRKBTYvn07EAWVuIVQK4jU89hjj5VcT7ZkAHK53MrjJ/9+fDvAxESUHCDZalFQEZG0\nKZDU0EhgiD/Ukx/o3SpLHFziy3G3HFAyTiMi0k4KJDVMTEyUfDj3mmQ3WnLcJh5rgWjsJ74cj9eI\niDRDgaSG2dnZ1FsanRLPHoPVsRYobcXk8/mVLjG1XkSkUQokdWzbtq2lYJLLraYzGx8fB1anBQMM\nDQ2VtBzSEo+51JIce5mammJ+fn7lMrByvRs01iOSXQokdUxPTzM8PLzywT8+Pr4ygys5AB5LDr5D\n7VlVSQsLCyuPH1+vNC5S/jeHh4dTCUTFYnGPAJpssZRLPk996IsMFgWSBg0NDVVcRJicVQWr4wyV\ngkWs0hqUM844g3Xr1pXcVim4JP9Wcnrw9PR06uM5yRZLLpdbKZeIDDYFkgZUW4UOUdcXrC0tSvKc\nLVu2lCReK7//9PT0SuCo9I2/PADFrZZKraZ2KBaLLCwsdCyYqHUjkl0WLTDvbxs3bvSbbrqp28Wo\naWlpqekMno12ISVbK+0OLHEgXUu5allLffQz1ccq1UWpNOvDzG529431zlOLZADEH+jJdSWFQoHh\n4WEgmtE1NDRUclu7xX97bm5OYygifUaBZIBUmtqbz+cZGhraYw1JnO4lnmhQrRUzOTm5MsGgPJ1L\n+SLIZEqYOIiJSO9TIOlh7fhGX20RYqWpvvFCxvLAEg/CJ3OTATVX/BcKhZVZYFoQKdLbFEikYcnA\nVSlANNslVigUSiYGJKcXb968uYWSikgnabA9I3pxADGtVf8TExPs2rWrZDr0oOvF10daVBelNNje\nY9IeJO61Qei0FkMmZ5hVqpNeqyeRfperf4pIZWNjYwwPD5ekgYnlcrmKx5sRp3VZWFhYCR7xVOZ4\ncaaIdJ9aJBnSa9+wy2d/xVOKk4Pm5Zt9xXumNLMKPw4oycH8OIVL+d9rhlo2Iu2hQNIEfeBUF8++\nKs/FFecmiy+X713fTNdYpXM7kfBSRGrTYHtG9PMAYr1tgeMANDs7W3e74moqrbBvtWxZ0s+vj2ap\nLkplYbBdYySSutnZ2Zof1MnFjPPz8yvdX90wPT1dM+FmWvcV6WXq2pKuqhRgZmdnWVpa4o1vfONK\n+pa1dmGV5xkDVtLBxI+ZDFzJnSRFpDFqkUhVWfiGHafvhygQTExM1J0Nlix3MjAUi0WKxSKFQmGP\nFfhTU1NMT0+zvLxMoVCo+7wr1U29lpdIv1KLpEm90qfeD5KzsaqtJykf4E/uRz89Pd1wS6Z8B8lq\nLZNuB1aRLFIgyRAFqeqqdYFVEucCa2WjrzgXWHm24uTumMmMyeWTBkQGiQKJVNWLH4jt3CUyzlYc\nd6UlU8IUi0WWl5f3SM2/ffv2klZLL9ahSLMUSJqkD4b2aXcLLPmhXslaN/Wqdp94vKX8cTVgL4NG\ngSRDFKTSk1wBn9wTpdZeK40qv39yXGZqagqonJZfpF8okEjXtDtwjo2N7dG1FW+glRyQj4/FgSXe\n5z6eIdbO7rF2bmssklWa/is9o9505PLAlMvlGBsb2yNty9zcXMmMsNnZWcbHx1cGytNYEKnZXtLP\n1CKRvhKPVyS38m2k5ZM8J5kPrLx1kpyt1cwiyfhxkq2ftSSb1OwwySLl2soI5Q8qtdb6aPcAfjKY\nVNqbvjzYNDOgPzEx0XA54/pI5iJr5v79RO+VUn2da8vMLjOzh8zs9iq3m5n9rZndaWa3mtlE4rZj\nzeyOcNtM4vgzzOzrZvYf4ff+aZVfelMaq8uHh4eZmJhgbm6u5mNPTEys5AprpHus2dld09PTGnOR\nTEqza+sK4OPAp6vc/krg0PAzCcwCk2a2F/AJ4LeB+4Abzexqd/8+MAPMufuWEGBmgHen+BxkwNUL\nSrUWSlbqGktqNQW+NveSrEitReLuNwD/XeOU44FPe+Q7wH5mth7YBNzp7ne5+xPAVeHc+D5XhstX\nAr+XTulFWjc7O1syVlOu2R0kywOHWieSFd0cbD8IuDdx/b5wrNLxeEnxge7+QLj8E+DAag9uZmcC\nZwKMjo6ytLTUpmKnY+fOnU3fZ2Ym6vXbsmVLu4vTdWupjyzasGEDt9122x7Hc7kchx9+eMOvy507\nd7Lvvvvy2GOPlRw/5phj2Lp1a1vK2iv65bXRLlmoj56dteXubmZVZwq4+yXAJRANtvfC4FyzZVy3\nbt2a7tcr+uF5XXrppSWJJOO0+PGU5Gaeo5ntceyxxx7j/PPPH7hB9354bbRTt+ujm+tI7gcOTlwf\nDceqHQd4MHR/EX4/1IFyirRkdnZ2JXCMj4+vDNo3++E/NzdXMY2+UrJIt3UzkFwNnBxmbx0J7Azd\nVjcCh5rZBjN7CnBiODe+zynh8inAVztd6CzR/he9IxlMWjU+Pl5yPTloPz09TT6f1wJI6ajUurbM\n7HPA0cCImd0HvB9YB+DuFwPXAscBdwLLwKnhtt1mdjZwHbAXcJm77wgPuwXYamanA/cAJ6RVfpEs\nGx4e3iOALC4urhyLsxDri4Z0ghYkZoQWWZVSfZSqVB/JtPblgSWWy+UYHx/vq4Ci10apnliQaGbP\nrHDseWstmHRPFrbOTcpaeXpNcmpxtTUprW7wJdKIRsZIvm1mK11IZnYu8OX0iiTSeb0Y1JrJ1dVr\nz016SyNjJEcDl5jZG4jWbfyAaNGgSF9ZXFzsuXGFRnN7le/cGOul5yrZVTeQuPsDZvYvwHlAEZhx\n99ZyO0hXZO1DI2vl6UXj4+Mr+8jH4jxf5V1ayetxAkqRdmhkjOQbRCvLjwBeBXzMzP467YL1g17s\nLhlU7Zye20mzs7PMz8+zbds2crkcuVxuJUDXShxZKBSUq0vappGurY+7+1fC5UfM7KVErRORvtLr\nLaR4nxRYfS75fL7mQPzk5GTJNsTJNPztTskv/auRrq2vlF3fDWxOrUQisiaVBt/n5uZKpglXUigU\nmJycZGJigsXFRZaXl0sCUD6fX9MmXDI46gYSM/sfIF5s8hSiRYUFd396mgUTKadvyGtTbY1JueQY\nSvL8QqGglorUVHeMxN1/2d2f5u5PA/YFfp9o7xAR6QHVcnQ1Y2FhYWVW2+LiIgsLCxr/kxVNpUjx\naBn8V8zs/USbSkkN+sYmWRKPoSwvLwPN7WeivU+klka6tl6XuJoDNgI/T61EIlUoMK9d+T7zi4uL\nK7PUGl35XigUSs6NWyjljy+Dp5EWye8mLu8GfsTqjoUi0mPKP/RrzeyqJZ5CXJ6NWAZPI2MkpyZ+\n3uruH3R37QMi0idaWT9TLBY1XiLVWyRm9nesztbag7u/I5USiUhHxS2UOBjEG2WVry2p1QXWi+ll\npH1qdW1lO++6ZI76y/tD+ZqReosbi8XimrrGpH/UCiSfDYsPRWQA1PsCMDY2tjJIn9xEK7Z9+3Ym\nJyfJ5XLMz8/XfKx8Pg80l8FYsqtWIPkuMAFRN5e7n9OZIkmvUkukv5XP/KrW1RWnXonXrVTbWKtQ\nKDA1NdV3G28NolqD7Za4fFTaBZH0KYmktMvs7Czbtm2reU6xWFzZWKv8dRe3ZuLb4xaK9KZagaT/\n9+DtcQoM0ivibq/p6emKQSOZhkV6T62ureeb2a1ELZPnhsuE6+7uL0y9dCKSaXGrZGpqquHV79UG\n5uPZYuU0iSP7agWSF3SsFAI0/4Zp9o2lN6KkZX5+nunp6T022SpXawpxnIW4XpeZZE/VQOLu93Sy\nINL79M2xuizWTbvLlMwMDFSc2dWIfD7P2NhYpupKamsqaeOgS/vDQG8c6XXlixtr7Skfb/dbHnDi\ny1kMvlKZAom0jd7wvSXt/1fcqqi0kDGXy5W0OsqnEzeaSFKywaLM8P1t48aNftNN2V6ov7S0xMjI\nSLeLkRmqj1K9XB/5fJ7l5eWVlkkul6u4dqTWTo7J7YB7uS7SkGZ9mNnN7r6x3nlVp/+a2W1mdmu1\nn/YWV6T90p4erenXjRkbG2NoaIhcLlc1iABMTExUfQxND862Wl1brw6/3x5+/0P4fVJ6xck29dmK\nNK9d75e4u2tmZoZ169btsdK+nX9LmlN31paZ/ba7vzhx04yZbUc7JErG6UOlt8SzvmqNjyTHW+KF\njbXS4MfdZRMTE3o9pKiRwXYzs6Pc/d/DlZfSwD4m0jx9q+ot+j+1X71xk0qzu+LAk9yTvnyAv9pi\nR2mPRgLJ6cBlZvb0cP0R4LT0ipRdg/bBocAm3TYxMdHwepRiscjU1NTK5aRCoaCMwymqG0jc/Wbg\nRXEgcfedqZdqQOkDWzqlV74kNNLdlVRrVf3y8jJDQ0PtKpok1A0kZnYg8FfAs9z9lWZ2GDDl7p9K\nvXQZ0+ybL+03qxZIyiCothalWY3mApPmNdK1dQVwOfDecH0R+CegbiAxs2OBi4C9gEvdfUvZ7fsD\nlwHPBX4OnObut4fb3gm8lShJ5Cfd/WPh+AfC8YfDw7zH3a9t4Hl0nPplJat67UvC2NhYWxYpJvdA\nSeq1+siaRgbNR9x9K1AECLsmPlnvTma2F/AJ4JXAYcCbQmsm6T3AQsgkfDJR0MHMjiAKFpuAFwGv\nNrNfTdzvo+4+Hn4yGUQgevHXmlFSLmvrErJWHhlcs7OzTExMMDw8vDKonhxcr6bS2pR4D5SFhQW2\nb9+uL3xt0EiL5DEz+xXC/iRmdiTQyDjJJuBOd78r3O8q4Hjg+4lzDgO2ALj7D83skNCV9gJgm7sv\nh/t+C3gdcEFDzyolyrbbXr3e9Sft08j/Kr7tjDPO4O6772ZsbGwl2/Dw8HDJ6nlYnfJba695aY9G\nAsm5wNVEe5L8O3AA8IYG7ncQcG/i+n1AeQ6EW4gCxLfNbBPwHGAUuB34YAhgjwPHAckcJ+eY2cnh\n2Lnu/rPyP25mZwJnAoyOjrK0tNRAkdvrhBNOAGDr1q11z925cyd33HEHQMNl3bVrV1Pnz8xES3+2\nbNlS58y1PX6zaj3+zp2tz+lIu/yd1I76yLJm/lfnnXceH/rQh9i1axdf+9rXVl7XseTre2lpiQ0b\nNrBjx46qgaNQKHDGGWc0/L7Imiy8NhrKtWVmewPPIxqvuMPddzVwn9cDx7r7GeH6W4BJdz87cc7T\niLqzXgzcBjwfeKu7L5jZ6cDbgMeAHcAv3P2PQ4tliaiFtBlY7+41pyO3K9dWs99wm5luuLS0xPnn\nn9/U4zerl76hK59SKdXHqrXURb29UnK5HPPz802XJQvvqSzk2mpk1tZ/Ahe6+8WJY9e4+6tr3A3g\nfuDgxPXRcGyFuz8KnBoe04C7gbvCbZ8iDOib2V8RtWhw9wcT5fgkcE2959AuzfalNjM+As2/GNPe\nCCtrsvCmjWWpLFJfvWnExWJR/9MWNDLYvgt4uZldbmZPCccOauB+NwKHmtmGcL8TibrIVpjZfonH\nPAO4IQQXzOyZ4fezibq//jFcX594iNcSdYNl0uzsbFMvSg1ui6QnHrDP5XIMDw/vcXs8AC/Na2SM\nZNnd32hmf040lvEGwsB7Le6+28zOBq4jmv57mbvvMLOzwu0XEw2qX2lmTtR9dXriIb4Yxkh2AW93\n90fC8QvMbDyU4UfAHzXyRNuh2RZG2pptIWXtG1fWytOMXiyz1E/BUiwW99ihsfzLnf73e2oo1xaA\nu18QkjVeDzyjkQcPU3OvLTt2ceLyPFDx09ndX1bl+Fsa+dtpSPsF1OzjZy2wpU1v4P7VrS8Vw8PD\ne8zoKhQKK1/SpqenS76wlb/nNHU40kggeV98wd2/YWavAE5Jr0iSlrTHYHq5hSGDaW5ujqmpqT0G\n4ZeXl2tutFV+n+Qix0F8/VcNJGb2fHf/IXC/mZWv6unYALdU1+wLNmtJ6wbxDSeVdfO1MD8/X7GL\nq5Lt27dXDDDFYnFlVlg+n8/Me6xTarVIziVaXf6RCrc5cEwqJRIRabN6reV4BXwraVji4NNqTrBe\nVGtjq7eG3y/vXHEkTc2OqWSt/1ddZ9Wpbtqj0pjJWsVrVwCGhoZKBvD7Ta2urdfVuqO7f6n9xZE0\nZW2ygEin1HttJmdotSM5ZKUNubL2xaydanVt/W6N2xxQIOlzaffzDtqCyl42KC2e2dnZlSCQy+VW\n9i+p1EpJbt9ba2A+1s+zLGt1bZ3ayYJI72v2w6afv6F1Wrxye3p6uu8/7NMWJ4AcHx/foy7j6cDL\ny8slr99KSSPL9fPrvZHpv5jZq4DDgX3iY+7+l2kVSnpTs2+U5eXllEqyNspGXF0vlnmt4pZDpecc\nH4tnQMbi1nulqcSxQqHA5OTkyqr6eMfG8pZKL9Z1I7m2LgaGgJcDlwKvB76bcrlkAGRt29Ne/8bY\nix9AWdRIPVbr9o3XktQaZykUCgwPD1MsFikUCmzfvn2lG61S91cvfAFpJNfWS939ZOBn7v4XwBRV\nVqPLYGt2I69mz+91zeZek94T/3/jDbgq5fSCPcdcisXiSndZL+bba6Rr6/Hwe9nMngX8FFhf43wZ\nUL2+8degLSJLUy98i05Tea6uRmaCxS2UOJjECxyTg/5Z1UggucbM9gMuBLYTzdi6NNVSiVSQtZX5\nIpWUB8/Z2dmaYyfl4u6uWNxaybK6gcTdN4eLXzSza4B93L37W3LJwGn2zTTo34q7qdfrvN2vnfHx\n8ZZXzcdfpLK4sLGRwfa9gFcBh8Tnmxnu/jfpFk2kVLPN+14fPJf+UWljrVwu19S+8fG4SnlAmpmZ\nYd26dV0NLo10bX0N+DnRVriNP2uRNlOXlnRKWh/KExMTK2lTktmC8/l8U6lZpqamVmZ57dixo+tj\nKHX3bDezW939hR0qTyratWd7mrQnd6leqI9Odp31Qn10Sr/WRaVB+UZbLblcruICylY1umd7I9N/\n/9nMfqcNZRLpK4uLi+o+6xG9sI11vBVwLJfLMT8/z7Zt20qOV1IsFtuSI2ytGgkk3wG+bGaPm9mj\nZvY/ZvZo2gUT6bRe+LCR/jc8PMzExATz8/NN33dycpJ8Pt/x13EjYyR/Q7QI8Tav1w8mIhVpBpk0\notbrI04SWS9BZDx9uJOvuUYCyb3A7Qoi0muURFJigxjA49QrndBIILkL+Dcz+2fgF/FBTf+VrFNg\nkH6RDIQTExMrr+16M72amV7cikYCyd3h5ynhR6QnNJvHK83zB/EbsaSj/LV0zDHH8Pjjj1cNGp3Y\nWqBmIAmLEX/Z3f801VKIpCDt3F/Nnq9xEknD1q1bGRkZKcnP1Wk1O9Dc/UngqA6VRaSvabqwpGl2\ndpbx8XGGh4dLfjqhka6tBTO7Gvg88Fh8UHu2izSnl1PmqzXVO7qRi6uRQLIPUer4YxLHtGe7SJPU\nGpG0dSvQN5L9V3u3i7RB1lokamV0T79t69xI9t9R4O9YHSv5NvBOd78vzYKJ9Ju039Rpfngo2Egt\njXRtXQ78I/CGcP3N4dhvp1UoEUmfgkP39FvdNxJIDnD3yxPXrzCzP06rQCIi9ahbrrYsDrb/1Mze\nDHwuXH8T0eC7iGRIlj5U9UFfW7/VTyOJWE4DTgB+AjwAvB7QALyI9C1lgm5OI7O27gFes5YHN7Nj\ngYuAvYBL3X1L2e37A5cBzyXahfE0d7893PZO4K2AAZ9094+F488A/olo698fASe4+8/WUj6RfpG1\nb7hplyMrz1MiVQOJmb2vxv3c3TfXeuCQXuUTRIPy9wE3mtnV7v79xGnvARbc/bVm9vxwft7MjiAK\nIpuAJ4B/MbNr3P1OYAaYc/ctZjYTrr+77jMVEWmQAmFzanVtPVbhB+B0Gvvg3gTc6e53ufsTwFXA\n8WXnHAZ8E8DdfwgcYmYHAi8Atrn7srvvBr4FvC7c53jgynD5SuD3GiiLiCSo60baqWqLxN0/El82\ns18G3kk0NnIV8JFq90s4iGgvk9h9QPmOLLcQBYhvm9km4DnAKHA78EEz+xXgceA4IN50/UB3fyBc\n/glwYANlEelrWjUv3VQv++8zgD8BTiL69j/R5vGILcBFZrYA3AZ8D3jS3X9gZh8GridqCS0AT5bf\n2d3dzCpuuGVmZwJnAoyOjrK0tNTGYrffzp07u12ETFF9lKpXHxs2bABo+HW+efPmps7PEr02SmWh\nPmqNkVxI1Fq4BPg1d6+9g8qe7gcOTlwfDcdWuPujhBlgZmZE+57cFW77FPCpcNtfEbVoAB40s/Xu\n/oCZrQceqvTH3f2SUHY2btzoIyMjTRa/83qhjJ2k+ihVqz7WrVtX95x+cv755wP9N9awVt3+v9dq\nkZxLtCPi/wbeG33OA9EsKnf3p9V57BuBQ81sA1EAORH4g+QJZrYfsBzGUM4AbgjBBTN7prs/ZGbP\nJgpoR4a7XQ2cQtSaOQX4aiNPVKSfNfuBms/nAZibm0ujODJgao2RtLTZr7vvNrOzgeuIpv9e5u47\nzOyscPvFRIPqV4buqR1EA/mxL4Yxkl3A2939kXB8C7DVzE4H7iFa4yIiA0QtkWxpZGX7mrn7tcC1\nZccuTlyeByqmRHX3l1U5/lMg38ZiigycZjMR91u22nbr9fK3KtVAIiLZNKgfeJKOlrqvRKQ3aR2J\ntJNaJCLSdb3eQur18rdKgURkAA36B5+0l7q2RKTtmu06U1dbb1OLRETaTilbBotaJCIi0hK1SESk\nrmbHVJpdMa8xm96mFomIiLREgUREJOOyPnlBgURE2i6fz68khmxE2h98mhWWLo2RiIhkXNbHkBRI\nRKTr0p4unPYHcdaSNna6HAokItJ22udksCiQiEjXNZvWPmuy0hLpFg22i4hIS9QiEZGuG/Rv9L1O\nLRIR6TmazpstapGISM9RUshsUYtERPqeWjDpUotERHpO2tOLs7YuJOsUSESk7ykgpEtdWyLS99S1\nlS61SESk58QJIRvt4mp2cF4tmOaoRSIiIi1Ri0RE+l6v5/7K+uC/AomI9BzN2soWBRIRkYzLekBT\nIBERKdPsB/egt2A02C4iIi1RIBER6bB+W9eiQCIiIi1JNZCY2bFmdoeZ3WlmMxVu39/Mvmxmt5rZ\nd83siMRt7zKzHWZ2u5l9zsz2Ccc/YGb3m9lC+DkuzecgItJus7OzfTWeklogMbO9gE8ArwQOA95k\nZoeVnfYeYMHdXwicDFwU7nsQ8A5go7sfAewFnJi430fdfTz8XJvWcxCRwdRvXU9pS3PW1ibgTne/\nC8DMrgKOB76fOOcwYAuAu//QzA4xswMTZdvXzHYBQ8B/pVhWEZEVaadU6bdZXml2bR0E3Ju4fl84\nlnQL8DoAM9sEPAcYdff7gb8Gfgw8AOx09+sT9zsndIddZmb7p/UEREQakc/nV/J/ZUGnW1TdXkey\nBbjIzBaA24DvAU+G4HA8sAF4BPi8mb3Z3T8DzAKbAQ+/PwKcVv7AZnYmcCbA6OgoS0tLHXg6a7dz\n585uFyFTVB+lVB+rOlEX7g7Q8OdGs+fv2rWrqfNrqVQf7Xz8RqQZSO4HDk5cHw3HVrj7o8CpAGZm\nwN3AXcArgLvd/eFw25eAlwKfcfcH4/ub2SeBayr9cXe/BLgEYOPGjT4yMtKeZ5WiXihjJ6k+Sqk+\nVqVdF8973vOa+jvf/OY3m3r8Sy+9tOky1VJeznXr1lU8npY0A8mNwKFmtoEogJwI/EHyBDPbD1h2\n9yeAM4Avuo4IAAAKTklEQVQb3P1RM/sxcKSZDQGPA3ngpnCf9e7+QHiI1wK3p/gcRGQA9frYRafL\nn1ogcffdZnY2cB3RrKvL3H2HmZ0Vbr8YeAFwpZk5sAM4Pdy2zcy+AGwHdhN1eV0SHvoCMxsn6tr6\nEfBHaT0HERGpL9UxkjA199qyYxcnLs8DY1Xu+37g/RWOv6XNxRQRkRZoZbuISIdlbZZXq7o9a0tE\nZOAsLy93uwhtpUAiItJhQ0ND3S5CWymQiIh0WLM7PMbdYFndMliBREQk47LeFaZAIiKScVnvClMg\nERHJuLGxiqskMkPTf0VE+kynkzYqkIiISEvUtSUiknHN5s7qdK4ttUhERPrMoO1HIiIibdbsDo+t\nUotERERaohaJiEif6fR0YbVIRESkJWqRiIj0Gc3aEhGRnqJAIiIiLVEgERGRliiQiIhISxRIRESk\nJQokIiLSEgUSERFpiQKJiIi0RIFERERaYu7e7TKkzsweBu7pdjnqGAGWul2IDFF9lFJ9rFJdlEqz\nPp7j7gfUO2kgAkkvMLOb3H1jt8uRFaqPUqqPVaqLUlmoD3VtiYhISxRIRESkJQok2XFJtwuQMaqP\nUqqPVaqLUl2vD42RiIhIS9QiERGRliiQdICZ7WNm3zWzW8xsh5n9ReK2c8zsh+H4BYnj55nZnWZ2\nh5m9ojslT0e1+jCzcTP7jpktmNlNZrYpcZ++rQ8AM9vLzL5nZteE688ws6+b2X+E3/snzu3ruoCK\n9XFheJ/camZfNrP9EucOXH0kjp9rZm5mI4ljna8Pd9dPyj+AAcPh8jpgG3Ak8HLgG8AvhdueGX4f\nBtwC/BKwAfhPYK9uP48O1Mf1wCvD8eOAfxuE+gjP8U+AfwSuCdcvAGbC5Rngw4NSF1Xq43eAvcPl\nDw96fYRjBwPXEa2RG+lmfahF0gEeKYSr68KPA9PAFnf/RTjvoXDO8cBV7v4Ld78buBPYRJ+oUR8O\nPC0cfzrwX+FyX9eHmY0CrwIuTRw+HrgyXL4S+L3E8b6tC6hcH+5+vbvvDle/A4yGywNZH8FHgT8n\net/EulIfCiQdEpqmC8BDwNfdfRswBrzMzLaZ2bfM7CXh9IOAexN3vy8c6xtV6uOPgQvN7F7gr4Hz\nwun9Xh8fI/pAKCaOHejuD4TLPwEODJf7vS6gcn0knQb8c7g8kPVhZscD97v7LWXndqU+FEg6xN2f\ndPdxom9Sm8zsCGBv4BlE3Tp/Bmw1M+tiMTumSn1MA+9y94OBdwGf6mYZO8HMXg085O43VzvHoz6L\ngZheWa8+zOy9wG7gsx0tWJdUqg8zGwLeA7yvawUrs3e3CzBo3P0RM/tX4FiibwtfCh8U3zWzIlHe\nnPuJ+j9jo+FY3ymrj1OAd4abPs9qU76f6+Mo4DVmdhywD/A0M/sM8KCZrXf3B8xsPVHLDfq7LqBK\nfbj7m83sD4FXA/nwnoEBrA/gH4jGP24J3ztHge1hckp36qPbg0iD8AMcAOwXLu8LfJvoDXEW8Jfh\n+BhRk9SAwykdMLuLPhpArFEfPwCODsfzwM3hcl/XR6JejmZ1cPlCSgfbLxikuqhQH8cC3wcOKDtn\nIOuj7PiPWB1s70p9qEXSGeuBK81sL6LuxK3ufo2ZPQW4zMxuB54ATvHo1bDDzLYSvXF2A2939ye7\nVfgUVKuPR4CLzGxv4OfAmQDu3u/1UckWoq7O04lm5ZwAA1sXAB8n+nD8evgW/h13P2uA66OibtWH\nVraLiEhLNNguIiItUSAREZGWKJCIiEhLFEhERKQlCiQiItISBRLJLDP7lZAJeMHMfmJm9yeu/78O\nl+VPzOz7IfvsnJk9p+z2p5nZfWb28cSxK8zs7kSZx8vu8xIz221mr08ce1fIiHy7mX3OzPYJxzeH\nv71gZteb2bPC8U2Jx7/FzF5b4zl8wcz+V4Xjf5gsd6vM7NfM7Ip2PZ5knwKJZJa7/9Tdxz1KpXIx\n8NH4uru/tMPF+R6w0d1fCHyBKDtv0mbghgr3+7NEmRfig2ENzYeJMh7Hxw4C3hH+zhHAXsCJ4eYL\n3f2FoS6uYTU9xu3h/HGiRXt/H9bhlDCzw4kWpt3V7BNvlrvfBoya2bPT/luSDQok0pPMrBB+Hx0S\nXn7VzO4ysy1mdpJF+53cZmbPDecdYGZfNLMbw89Rzfw9d/9Xd18OV5PZZzGzXydKqnh9pftWcQ7w\nRVZTn8T2BvYNwWCIkAHZ3R9NnPNUQu4td1/21ay4+1A9J9dJwFcTZT7VzBbN7LtEaTji478bkoh+\nz8y+YWYHmlnOon1RDgjn5Cza7+IAM3tDaD3dYmbJQPo1VoOg9DkFEukHLyJKN/MC4C3AmLtvIsrV\ndU445yKiFs1LgN9nz5TczTidkH3WzHLAR4A/rXLuh0KX1EfN7JfCfQ4CXgvMJk909/uJsh7/GHgA\n2OnuyRbLBy3KjHwSiYR9ZjZpZjuA24CzEoEl6Sjg5nD+euAvwrHfINrDIvZ/gSPd/cXAVcCfu3sR\n+Ez4uwC/Bdzi7g+HcrzC3V8EvCbxODcBL6tSJ9JnFEikH9zo7g94tK/Lf7LaMrgNOCRc/i3g4xal\nrr+aKBngcLN/yMzeDGwkyoUF8DbgWne/r8Lp5xHlUHsJUZbnd4fjHwPeHT6gk4+9P9F+EhuAZwFP\nDX8PAHd/r0eZkT8LnJ04vs3dDw9/57x4XKXMeuDhcHmSaNOwh939CeCfEueNAteZ2W1EGakPD8cv\nA04Ol08DLg+X/x24wszeStQVF3soPAcZAMq1Jf3gF4nLxcT1Iquv8RzRN+2fV3sQM7sceDHwX+5+\nXIXbfwt4L/CbIWgBTBHtKfM2YBh4ipkV3H3GV/cT+UV47LjVshG4KuSMGgGOM7PdRBt83R2+6WNm\nXwJeStQaSPoscC3w/uRBd/9B6PI7gqhFkPQ4UddXPX8H/I27X21mRwMfCI99r5k9aGbHEG2UdFI4\nfpaZTRJtvHSzmf26u/80/K3HG/h70gfUIpFBcT2r3VyUz6ACcPdTw6B4pSDyYuDvgdf46k6WuPtJ\n7v5sdz+EKFB82t1nwn3Wh99GtMPh7eE+G9z9kHCfLwBvc/evEHVpHWlmQ+E+eaKMyJjZoYniHA/8\nMBzfEA+uh5lkzyfKBlvuB8CvhsvbgN+0aFbcOuANifOezmra8VPKHuNSoqD2+TgRoJk9N7SI3kfU\n4olTmI/Fz1f6n1okMijeAXzCzG4let3fQDSu0qgLiVocnw8tiR+7+2tq34XPhgFqAxbq/T1332Zm\nXwC2E2Vu/R5wSbh5i5k9j6iVdU/isX4DmDGzXeG2t7n7UoWH/z9Eaci/4dEeJx8A5oFHQtliHwjP\n8WfAN4m62WJXE3VpXZ44dmEIcgbMEaUwB3h5+JsyAJT9V2QAmNm+wL8CR601rbiZbSSasFBzED1M\nKvgW8BtVBv6lz6hrS2QAuPvjRGMqa9q/28xmiKYrn9fA6c8m2pRLQWRAqEUiIiItUYtERERaokAi\nIiItUSAREZGWKJCIiEhLFEhERKQlCiQiItKS/w8CaeTrSANVTgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "There are a set of useful functions in `LightCurve` objects which you can use to work with the data. These include:\n", "* `flatten()`: Remove long term trends using a [Savitzky–Golay filter](https://en.wikipedia.org/wiki/Savitzky%E2%80%93Golay_filter)\n", "* `remove_outliers()`: Remove outliers using simple sigma clipping\n", "* `remove_nans()`: Remove infinite or NaN values (these can occur during thruster firings)\n", "* `fold()`: Fold the data at a particular period\n", "* `bin()`: Reduce the time resolution of the array, taking the average value in each bin.\n", "\n", "We can use these simply on a light curve object" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2UJXV95/H3p4dR0/YqaCOH0Oig0iqiTjrjwGhYH1oT\nQlRWNyYaXQmCHNqHaNasjrpmk3V1R9nEsNEz7ARQkpigxhgJiwFtE3F3hxEYGpgR6EyGGEATHLOM\nO/REZujv/lFVPdV37kPdvrfurb738zpnzvStp/ur362q7++h6leKCMzMzFZqpN8JMDOz1c2BxMzM\nOuJAYmZmHXEgMTOzjjiQmJlZRxxIzMysIw4kZmbWEQcSMzPriAOJmZl15Jh+J6AXxsfHY926df1O\nRlOHDx/mmGOG4ucoxPmxnPPjCOfFcmXmx6233rovIo5vtdxQ/Brr1q3jlltu6Xcymtq3bx/j4+P9\nTkZlOD+Wc34c4bxYrsz8kPTdIsu5acvMzDriQGJmZh1xIDEzs444kJiZWUccSMzMrCMOJGZm1hEH\nEjMz64gDiZmZdcSBpE0zMzPMzMz0OxlmZpXhQGJmZh1xIDEzs444kJiZWUccSMzMrCNDMfpvN23d\nurXfSTAzqxTXSMzMrCMOJGZm1hEHEjMz64gDiZmZdcSBxMzMOlJaIJF0paQHJe1qMF+S/rukPZLu\nkDSVm3e2pHvSeZtz0y+RdHe6/JclHVtW+s3MrJgyaySfBc5uMv/ngVPTfxcBWwEkrQE+nc4/DXij\npNPSdb4GnB4RzwfmgQ+UknIzMyustEASETcC/9xkkXOBP4zETcCxkk4ENgJ7ImJvRDwCXJ0uS0Tc\nEBGH0/VvAibKSr+ZmRXTzz6Sk4D7cp/vT6c1ml7rrcBXS0udmZkVsiqfbJf0IeAw8Lkmy1xE0mTG\nxMQE+/bt61HqVmb//v39TkKlOD+Wc34c4bxYrgr50c9A8gBwcu7zRDptbYPpAEj6VeBVwHRERKON\nR8Q2YBvAhg0bYnx8vGsJL8tqSGMvOT+Wc34c4bxYrt/50c+mrWuAt6R3b50J7I+I7wM3A6dKOkXS\nY4A3pMsi6WzgfcBrImKhXwm33vMLxcyqq7QaiaQ/BV4KjEu6H/hPJLUNIuIy4DrgHGAPsACcn847\nLOmdwPXAGuDKiNidbvZTwGOBr0kCuCkiLi5rH8zMrLXSAklEvLHF/ADe0WDedSSBpnb6M7uTOjMz\n6xY/2W5mZh1xIDEzs46sytt/bfj4hWJm1eUaiZmZdcSBxMzMOuJAYmZmHXEgMesSPzRpw8qBxMzM\nOuJAYkPLNQiz7nAgMbNKc8CvPgcS6zqf+Gb918vz0A8kmnWJH5ocPtmFuoq//fz8fM++y4HEVp0q\nn7y9MOz7b8UcOHCgZ9/lpi3rSL+asebn59m0aRPT09Mr3kZ2IV6NzXDz8/OrMt1V42bY7nCNpEO9\nLh1WpTTaz5Nvfn5+qbTVy1JXK734bWZmZti5c2dp2++V6elpFhYWWL9+fcv86vex3g1Fjo1uHT8z\nMzPMzc0t+1x2HjqQtCkrAc/OzvY5JY21e0B2+wJYbztlXmRXeqLMzMwwPz/P5ORkV9KQtUnng2y/\nL4Jl5Xs3tru4uNjTdvxOVaUQ18r8/DyLi4s9/U4HkoKyC8XCwgKjo6NLAaXeRajVAbfSA7LIha92\nmW4f/CutiZR5wciXvorI9mFubu6oE65oftVbbmGh/Lc/52sjZV+EG+VF/hirysW1WTo6TWOzY77I\nNludt/lCyEquJ7X6UUt3IGnDwsICi4uLSz/UyEjzLqb8AZA/GFsdVNk6nSpyAGfpWVhYKFSyzy6+\nIyMjrF+/vul2a7d14MABpqenj6rNtbPPMzMzR50oKyl95UtttRfkrP+htg+lVfqKpiP/u3TyO7cK\nXNPT0w0vKs2CRL3pednx0oksXUWPu3YVPd+KBuNmF/pW6cjOmbIC/+bNm1m7di3Q2zu18hxIOrC4\nuMjc3FzDC2pe7YFYe7FqpdUBmT9xFhYWmJubY3R0tOH2Gl1kWpVCs4tllo52m/gOHDjApk2b2L59\ne9352fd/5CMfqTu/3dpHfpuwvH8lk784ZPPbPSHbWb5eTaiI2pLx4uJi3cCcye/n9PR02xfBleR1\nu7JzaNOmTQ37S1ZauKoNeLUFu/yx0G5zZG0totE6RZqZsuNhbGxs2XVh06ZNAIWvL1lBtx8cSDq0\nuLjY8O6ffEk/+5GzpomRkZGlE2h0dJTPf/7zy9att63sIGlUGsxO/Nrl6gWtetvI0jY2NlZ327UH\nabNS6fz8/FKHaj5N2d/59GQn9BlnnAHA1NTU0rK1pfd2T5R853SjGmS9IFBbUs72B+qXSvMX7Xw+\n1utTa1QT6oaswFFbiChag8iO2U2bNi2ls/bYyfZ1586dS8dKFqiKXoTzGuVH/mJdLx356bX7XO94\nrd1+kSag2jRk52x+27V9Y/nrQf47GjVb5c/X0dHRtu/I27t3Lw8//HDD+b1odnQgKahZKTI7ifIX\nkOygyUrutevmPy8sLCyrntb7zvxFMLsQ5zX7jtqLSP7gLaL2xKm3D7WKnqT1lp2bm+PVr371spJY\n/kK+Uo3Se+DAgaW05APfzp07l0qF+W3UNpU0Oumb1fgazW9Xo23Ua/5rVHvO9//Vy6P8ha1ec+XI\nyMiyC2CzPpVse81k31V7d1pWUMrOlfz3ZE3OtTXn/D7kt1+77VZpavZbZU22k5OTy/KqyN11+Vrf\n4uLiUQXOVmmbnp5uGkRaNb93iwNJC1mputVFM//D1zaNtCoNLi4usnfvXp71rGctW6+2FJ+XLVfk\nYtRuKb4bfRCtZBe1eidJvoTaaXW9nVJ/s7TkZb/3yMhI28GtjOai/MW7WYEnf4HK1/qApvmcP55b\nBc3aC3a++Wjnzp1MTU01DX5ZE1ez32Lv3r0cPHhwKe/z51eji/eBAweWak/1tp2lqV5tPF9IbCRL\nw8LCQsM+pJ07dy41ReabquvtY23asoDQbjNfoybkbistXEm6UtKDknY1mC9J/13SHkl3SJrKzTtb\n0j3pvM256U+S9DVJf5v+f1xZ6e/UgQMHlg7OIhfCgwcPLh2A+XWLbL+I/IWgyMU1O0lbVbHrzS9y\nkGcl5GZB9sCBAx0HsaJ5tHPnzq58XzO1Nbt2S4uNfrds+vT0dOH079y5c6n0vHPnzpbrNfud8utm\nwSBfy2lHtq1Gv1tWAl/J75SlpVXtolaRmkV2PGc1o0bpW1hYYHp6uivNms1upoCjCwtlUkSUs2Hp\nXwMHgD+MiNPrzD8HeBdwDnAGcGlEnCFpDTAPvBK4H7gZeGNEfEfSJ4B/jogtaYA5LiLe3yotGzZs\niFtuuWVF+9GoWacMIyMjpX5XvrRV5AI7NTVV6CSamppaFjxaHeCZbu3vyMjIspJXvtSW9buUqch+\njI2N1c2THTt2FPqOffv28cu//Mt1tzE2Nsbs7Gzp+zo2NtZxDbHsY7wb359vmi7j4c/sPFxp02YW\nIFqlreix1YykWyNiQ6vlSmvaiogbJa1rssi5JEEmgJskHSvpRGAdsCci9gJIujpd9jvp/y9N178K\n+BugZSDpRC8P+l59Vzul9CLm5uaWquxFgwh0vr9Z23xtf0W+hN4LRfajXp50q/26F8+vQHf6dPoZ\nRIp+f7PmqW7odLtVHNmgn30kJwH35T7fn06rNz0rap0QEd9P//5H4IQyE9irC1GvlPWgUtZJWOZ3\n1JN1HPf7YbiVKnJbZxH1br6wlct3epe1/bLV6+sp06rtbI+IkNSwXU7SRcBFABMTE+zbt28l33HU\ntH5XzatqcXGxJ81ImZGRkaVnTfK/7Uqf0eiHQ4cOFT4u9+/fX/d4zG/Lume1HEONnHLKKSu65q1U\nPwPJA8DJuc8T6bS1DaYD/JOkEyPi+2kz2IONNh4R24BtkPSRjI+Pt53AgwcPHjVtdHR02ZPtq/2A\nW4kq7Pfo6Cj1ftN+p6sda9eurbsPjXzjG99YNhpBvva3e/fuuusU7efqh270udjRxsbG2j62OtXP\nYeSvAd6S3r11JrA/bba6GThV0imSHgO8IV02W+e89O/zgK/0OtGTk5Ps2LGDqamprjVNrCbZ0ChF\n2vfLuIe9m1X2Mqr/Rfd5x44dHTXJ1T7cVu9iPDIywtatW3vezFHU5OQk27dvZ2pqqmfPO+SNjIx0\n9L1VzNeRkZGjxtnryfeWtWFJfwpsB54l6X5JF0i6WNLF6SLXAXuBPcAfAG8HiIjDwDuB64G7gC9E\nRFbc2gK8UtLfAq9IP5emWUlp69atK7oQjI2NMTY2tnQQt3sgl7l8q2XHxsaWnh4ucn96vUDb7ORr\ndUEZGRlhdnaWqamppoPbNVs///3dGPW31ujo6NJv3CwdZRsbG6t0QSefB1u3bi0lrWNjYw1vgW02\nVlwztb9ddi6vJA21y3WqH8F46bvL2nBEvDEiToyItRExERFXRMRlEXFZOj8i4h0R8YyIeF5E3JJb\n97qImEznfTQ3/YcRMR0Rp0bEKyLin8tKPxz9446NjR0VPNr58fK3FY6Oji79K6JR0MmnsfagHhsb\na3rBzwe1LEg0O6BnZ2eZnZ0tHEC3bt16VHpmZ2dXdNJkd2hl262Xhmb35me3CGffn/2WKwnMjdbJ\nAt3k5GTDIJXfj060qmnU/k4rKbRMTU01/I6RkZGOnlOoHVer3v60uki30qygMDo6uhTAin5HFnzy\n6Zqdna0bkLL52bAxrb5jcnLyqP1vN4/buZ50m9+Q2ETtgVjvwFy/fv2yH7vRwZ9dZFp9Rz1Fm9Fm\nZ2eXmgqmpqYKD6g4Ojpa6AJYq0hAyE7UfJW73ndkF/bt27cftd2pqamOT5D8+s32M5+e2hpSFmwb\n1Tryd5H14k6yRr9vbboaXezqyQoV2e8xOTl5VE0OjlyIs9+23nHfboGh9jfJjucdO3YU2lZ+mfyz\nTa0u4tk53OyindXwtm7dyuzs7LLjqVFQLxqgsrzOattZ/me/WasabrZMdlz34y7GVXvXVq9l7c21\nsmn5AQrryZemM/WaYmpLrNmBMTMzs/TMRO0ghLUHbLMDKVs2Cx75W5zz+1I7ZES9C282tlCjh+Ty\n+54/wLP9yd9hld9+7bhlnZwYU1NTR42Nld/e9u3bl/Y3K3Fmw5jk01lv32sfeCySztrnXtpV5Dsa\nbX90dLRu53bW6Q0cNTxJ/hUI2W+dLxTkg2f2lHx2DC8sLNS9MaP2AdaiareVf8gzCwK1w/PnR/rN\n73uW/qKvW25VKMvOhSwfa4/5enmf5VOjAmZtHjUbrqlfASTjGkkLtaWDRiYnJ1m/fv1SKSq7YDdr\nHshkJd9WzR61B0v++5o1YeVL1rXV36y5qtl3NmvOalS7yGSl2nqydNQL0s973vM6LtG2mp6fn+V9\n1tyRXyfLsywv8vOyC1i90me9Gk07TYMr0azAk9VMsmNtamqKHTt21K151wbe7HfMl8zz87K/s1rb\n5OTksuOs3T7B2uVqtwf1H+yrd0HN0l5biyjyvVCsZpXl7/bt25eOn/x31Et/O4WKmZmZpZsTaj3+\n8Y+vu8+9DCyukTSRL421+lFq5xd5WC4/7+UvfzkHDx5cNkJqbckqky9BFpUvabY6eLOaSqvvqFeD\nqdeE1ywI1evX2Lp169KQIO1Y6ZAotS8eq5feRr8FLO/7qt1uvrZTRuf+Sm7FrldYyReUmh2zrUq+\ntTXLfI2g1fr5ZWvzKl9zuOeee5AELH8AtsiFs1FhrdGx3Ohi3+57ePI1tlZa5XH2m2c1u14+L9KI\nA0lFPP3pTz9qaOxGaqvlreRPkqLrtfMd+cBTtC2+SJNCty68Rfah2cWtk21nI73Oz8+3/TKzIvLP\nNRVRL9+LNu+sdPuZos0vrYLNhRdeyL333tuwubLddNV+d3axb9Ts1ErR37dR4aOVRqMj95MDyRBp\n58JchWFHqpCG1SDfJ1H0poSV5G2nz710EtDztmzZwoc//OEVpaNIwMlqtWXUIPN9jr3Iz15xIClJ\nuz/yli1bCj+JWlu6qsIB1W4tKVNm2qvywFht8043rTTf6+nmb1F7a29Zyth2WSM3dCutVTjfazmQ\n2EDJ+o9q75rpt9p+mG5udxiVud9Z7a5Kqv47O5C0UPUfsCqqkk9ZB36VgoitLivpFxl2LQOJpKdE\nxIM1054VEfeUlyyzlaly8Khy2sw6UaRG8i1JH46ILwBIei9wAXBaqSmzrvOFzMzKUCSQvBTYJun1\nJC+SugvYWGairJhhCwxVurnAzI5o+ZhpOrT7XwGbSF6De1VE9O41eGZmVmlF+ki+DnwPOJ3khVNX\nSLoxIn6j7MSZmVn1FRn45lMR8ZaIeCgi7gReBOwvOV2VNDMz43djm5nVaFkjiYi/qPl8GPhIaSky\nM+sy96+Vq0jT1v8DIv34GJJ3qh+IiCeWmTAbfO2e3L4IWK848LSnSI3kX2V/Kxly81zgzDITZWbD\nxRfu1a2t95Gkr8f9C+DnSkrPUHMfjJmtRkWatl6X+zgCbAD+pbQUVZhLS2ZmRyvyQOKrc38fBv6e\npHnLrCMOzNYrPtbKVaSP5PxeJMTMrCoceNrTMJBI+n2O3K11lIj4tVJSZGY9UaUO7iqkwVauWY3k\nlk43Luls4FJgDXB5RGypmX8ccCXwDJJ+l7dGxK503ruBtwEC/iAifi+dvh64DHgcSVPb2yPi252m\ntQqqdjJV6UJjZtXVLJB8Ln34cEUkrQE+DbwSuB+4WdI1EfGd3GIfBOYi4rWSnp0uPy3pdJIgshF4\nBPgrSddGxB7gE8BvR8RXJZ2Tfn7pStNpZmadaXb771IpP23matdGYE9E7I2IR4CrObqT/jTgGwAR\ncTewTtIJwHOAHRGxkAazbwLZ3WMBPCH9+4kk44CZmVmfNKuRKPf3i1ew7ZOA+3Kf7wfOqFnmdpIA\n8S1JG4GnARPALuCjkp4MHATO4UhT23uA6yX9N5JA+KK6iZcuAi4CmJiYYN++fSvYhd7Zv796w5cd\nOnQIoC95V8X86Kcy8qOfv28nBuHY2Lx5MwBbtmxpsWRrVciPZoGkYUd7F20BLpU0B9wJ3AY8GhF3\nSfo4cAPwMDAHPJquMwP8ekR8SdIvAVcArzgq8RHbgG0AGzZsiPHx8dJ3plNVS+PatWuB/qWravnR\nb93Oj8svv7yr2+ul1X5sdPvc6nd+NAskz5Z0B0nN5Bnp36SfIyKe32LbD5AMO5+ZSKctiYgfAefD\n0vAr9wJ703lXkAQJJH2MpEYDcB7w7vTvLwKr92youGHqZPeNBWYr1yyQPKfDbd8MnCrpFJIA8gbg\nV/ILSDoWWEj7UC4EbkyDy9K74iU9laT5Kxvf63vAS4C/AV4O/G2H6TSzPnMgX90aBpKI+G4nG46I\nw5LeCVxPcvvvlRGxW9LF6fzLSILVVZIC2E3yLvjMl9I+kkPAOyLioXT620iaw44huWX4ok7S2S4f\n8INrfn6emZkZ/7ZWukE7xooMkbJiEXEdcF3NtMtyf28HJhuse1aD6f8L+OkuJtPMzDrQ1ui/ZmZm\ntRxIzEiaGiYn61aOzayFZmNt3UnzsbZa3bU1kAatbdOO8G/bP8771a1ZH8mr0v/fkf7/R+n/byov\nOWZmg2/QbtppedeWpFdGxE/lZm2WtBPYXHbihs2gHVxmNhyK9JFI0otzH15UcD0zMxsCRW7/vQC4\nUtIT088PAW8tL0lmZraaFHlD4q3AC7JAEhH9HyHMzGwVG7Tm65ZNVJJOkHQFcHVE7Jd0mqQLWq1n\nZtU1MzOz1Cdn1qkiTVufBT4DfCj9PA98nnRAReueQSulmNlwKNJpPh4RXwAWIRlDiyNDupuZ2ZAr\nEkgeTgdPDABJZwLuJzEzM6BY09Z7gWtI3knyv4HjgdeXmiozsy7yM1rlKnTXlqSXAM8ieanVPRFx\nqPSUmVlpfEEdbL0OnEXu2vo74MKI2B0RuyLikKRre5A2G3C+c8hsMBTpIzkEvEzSZyQ9Jp12Uolp\nMqvLgcesmor0kSxExC9Leh/wLUmvp8mowGZmVeOmvHIVCSQCiIhPpIM13gA8qdRUmZnZqlEkkPxm\n9kdEfF3SzwHnlZckMzPrRK9rYM1ebPXsiLgbeEDSVM1sd7Zbx9zcYDYYmtVI3gu8DfidOvMCeHkp\nKTJrYDUHHj/HsLr492pPsxdbvS39/2W9S46Zma02zZq2XtdsxYj481Ybl3Q2cCmwBrg8IrbUzD8O\nuBJ4BvAvwFsjYlc6790kNSIBfxARv5db710krwB+FPifEfG+VmkxM7NyNGvaenWTeQE0DSSS1gCf\nBl4J3A/cLOmaiPhObrEPAnMR8VpJz06Xn5Z0OkkQ2Qg8AvyVpGsjYo+klwHnAi+IiB9LekqLfTQz\nsxI1a9o6v8NtbwT2RMReAElXkwSAfCA5DdiSft/dktZJOgF4DrAjIhbSdb8JvA74BDADbImIH6fr\nPdhhOlctt+OaWRUUuf0XSb8APBd4XDYtIv5zi9VOAu7Lfb4fOKNmmdtJAsS3JG0EngZMALuAj6aj\nDh8EzgFuSdeZBM6S9FGS5rDfiIibi+yHWb842K8u/r3a0zKQSLoMGAVeBlwO/CLw7S59/xbgUklz\nwJ3AbcCjEXGXpI+TPPz4MDDHkXegHEPyQOSZwAuBL0h6ekQse9pe0kXARQATExPs27evS0kux/79\n7Y/Mf+hQMnZm1fdtJVaSH4NsteXH5s2bAdiyZUuLJdu32vKibFXIjyI1khdFxPMl3RERvy3pd4Cv\nFljvAeDk3OeJdNqSiPgRcD6AJAH3AnvTeVeQvoVR0sdIajSk//95Gji+LWkRGAd+ULPtbcA2gA0b\nNsT4+HiBJPdXu2lcu3btitZbLQZ1v1ZqNeVH2cfmasqLXuh3fhQJJAfT/xck/STwQ+DEAuvdDJwq\n6RSSAPIG4FfyC0g6lmQsr0eAC4Eb0+CCpKdExIOSnkrS/HVmutpfkNSO/lrSJPAYYPCK5AW4+m1m\nVVAkkFybXvAvAXaS3LF1eauVIuKwpHcC15Pc/ntlROyWdHE6/zKSTvWrJAWwG7ggt4kvpX0kh4B3\nRMRD6fQrgSsl7SK5o+u82mYtMzPrHbVzDZb0WOBxEdH/Rrk2bNiwIW655ZbWC/bRvn37Sq+erqa7\nvHqRH6tJv/OjSsdOv/OiasrMD0m3RsSGVssV6WxfA/wCsC5bXhIR8budJtLMzFa/Ik1bf0lym+2d\nwGK5yTEzs1pVqhHWUySQTETE80tPiZmZrUpFAslXJf1sRNxQemrMrJKqWhK2aigSSG4CvixphOQO\nKgEREU8oNWXWdb4YmJWj6k1PZSsSSH4X2ATc6dtszY4Y9ouHWaZIILkP2OUgYmbWH1UvrBQJJHuB\nv5H0VeDH2UTf/mtmZlAskNyb/ntM+s/MzHKqXmMoW9NAkj6M+K8i4jd6lB4zM1tlmgaSiHhU0ot7\nlRizbiq7M3zYS6H95BsdqqVI09acpGuAL5K8GwQo9s52Gy4+uc2GU5FA8jiSoeNfnpvW8p3tZmY2\nHFoGki68u93MzAZYkdF/J4DfB7K+km8B746I+xuvZeCmHrOy+JyqliJNW58B/gR4ffr5zem0V5aV\nKLNu8MXGrDeKBJLjI+Izuc+flfSeshJkq5cv3GbDaaTAMj+U9GZJa9J/bybpfDczMytUI3krSR/J\nJ0nu1vo/gDvgC3AJ3XrF/XHWT0Xu2vou8JoepGXg+OQ2s2HQMJBI+s0m60VEfKSE9JiZ2SrTrEby\ncJ1pjwcuAJ4MOJCYmVnjzvaI+J3sH7AN+AmSvpGrgacX2biksyXdI2mPpM115h8n6cuS7pD0bUmn\n5+a9W9IuSbvr3SUm6b2SQtJ4kbSYDbKs+TRrTjXrpaZ3bUl6kqT/AtxBUnuZioj3R8SDrTacjhz8\naeDngdOAN0o6rWaxDwJzEfF84C3Apem6pwNvAzYCLwBeJemZuW2fDPws8A+F9tLMzErTrI/kEuB1\nJLWR50XEgTa3vRHYExF70+1dDZwLfCe3zGnAFoCIuFvSOkknAM8BdkTEQrruN9O0fCJd75PA+4Cv\ntJmmnnInu1k5fCNLtTSrkbwX+EngPwLfk/Sj9N//k/SjAts+ieQ1vZn702l5t5MECCRtBJ4GTAC7\ngLMkPVnSKHAOcHK63LnAAxFxe4E0WA/NzMy4acVsCDWskUREkYcVO7UFuFTSHHAncBvwaETcJenj\nwA0knf5zwKNpUPkgSbNWU5IuAi4CmJiYYN++fSXtQnfs37+/30no2KFDhwC6kteDkB/dVCQ/upn/\nVbZ///6h2dciqnCuFHkgcaUeIK1FpCbSaUsi4kekDzdKEskrffem864ArkjnfYykRvMM4BTg9mRx\nJoCdkjZGxD/WbHsbSbMcGzZsiPHx6vfJr4Y0NrN27Vqge/vR6XbKbP7oR9NKq/y4/PLLe5SS/uv2\nsdapfje19TsfygwkNwOnSjqFJIC8AfiV/AKSjgUWIuIR4ELgxjS4IOkpEfGgpKeSNH+dGREPAU/J\nrf/3wIaIcLHEzKxPSgskEXFY0juB64E1wJURsVvSxen8y0g61a+SFMBukmdUMl+S9GTgEPCONIiY\nmbmTvWLKrJEQEdcB19VMuyz393ZgssG6ZxXY/roOk2hd5JPbbDj1okPdKsJ3VZkNh16f66XWSMz6\nqcwakmtfzfW787nXyt7PquenayRmZtYRBxIzM+uIA4mZmXXEfSRDpKrtq1Z9VW+jt+V6/Ts5kKxi\nPrmtqnxMdlfV89NNW2Zm1hEHErMe8XM8NqgcSMzMrCPuI2lD1fokqpIOG3w+1qwZ10hs1XDTkFk1\nuUZSIVWr8Vh3+Xe1QeUaiZmZdcSBxMz6ruxmSzeLlstNW21w00R3uSnPbDA4kNiq4YDTXQ7k1i0O\nJBXiE9qsGhxk2+M+EjOzHhu0PhvXSNowbKWUYdtf65+yjzEfw+VyILG+8cltNhgcSMyGlAN5Y86b\n9riPxMzMOlJqIJF0tqR7JO2RtLnO/OMkfVnSHZK+Len03Lx3S9olabek9+SmXyLp7nSdL0s6tsx9\nsOrwQ2t0Xu2EAAAPZUlEQVSDa9jyfuvWrQNV6ymtaUvSGuDTwCuB+4GbJV0TEd/JLfZBYC4iXivp\n2eny02lAeRuwEXgE+CtJ10bEHuBrwAci4rCkjwMfAN5f1n7kDdIPX8Sw7a+ZrUyZNZKNwJ6I2BsR\njwBXA+fWLHMa8A2AiLgbWCfpBOA5wI6IWIiIw8A3gdely92QTgO4CZgocR9siMzPzzM/P9/vZJh1\nrNc1vDI7208C7st9vh84o2aZ20kCxLckbQSeRhIYdgEflfRk4CBwDnBLne94K/D5el8u6SLgIoCJ\niQn27du38j3pgf379/c7CT23eXPS2rlly5aj5tXLj0OHDgGU9ltGRKnb78SgHx/t/LaDnhft6se5\nUqvfd21tAS6VNAfcCdwGPBoRd6XNVjcADwNzwKP5FSV9CDgMfK7ehiNiG7ANYMOGDTE+Pl7aTnTL\nhz/8YWB4mpTWrl0LQKPfpnb65ZdfXmp6JDVNT79VNV3d0OpYqFV2XpT9DFW3t1+bH+3mZ6fKDCQP\nACfnPk+k05ZExI+A8wGUnMX3AnvTeVcAV6TzPkZSoyH9/KvAq4DpyIqR1nXD9kDi5ORkv5MwtFb7\nMTZs50qtMgPJzcCpkk4hCSBvAH4lv0B6x9VC2odyIXBjGlyQ9JSIeFDSU0mav85Mp58NvA94SUQs\nlJj+owz7wTLo/LuarUxpgSS9q+qdwPXAGuDKiNgt6eJ0/mUknepXSQpgN3BBbhNfSvtIDgHviIiH\n0umfAh4LfC1tirgpIi4uaz/MusGFEBtkpfaRRMR1wHU10y7L/b0dqNueEBFnNZj+zG6m0frHF9XB\n5cDZ3KDlS7872y2n3YOraidr1dIzbNrJ/2H7rdrd39WeL72+jd2BpA2r/eAatpPJ+mfYnscZ9nPF\nY22VaNiGfagaD6kyuKqW91VLT6+5RmIDq0ql4tVeYl1tz1UMu17fyu5AUqKyL2RVulBC9dJTtnb3\nt0oXy7LTUPaFbNiOtXa5j8QKq9oDdFVLz7CpQoDKVCktUL30lM01kgpb7dX71X4ytZs/ZZ9MDpzW\nK9PT0wDMzs72OSX1OZCUyBea5toNDO1W18sOnO2mp0rNMVVqZoP286bdc6tqhbSq5X+nHEjaUPaF\nrEoXmipaWOjpiDgDrWrH2mrvUym7UNTusd/rAOVAUqJhKwU5cDZX5sVytZdwq1bbrFpgqDoHkhKV\nXV0v22q/OLWr7P0dtsBpw8OBpA2+0HdXu6Wy0dHRklJi3dZu53DVzq2qWb9+fb+T0JQDifVNu4Gh\nahebdtNTpfS3e/fPai2sZNrd37Jrj4NWKHIgaUPZJ9NqP1mr1gczTKrWmbzatVtbbjdQValQ0Q0O\nJCUq+xbFqnVQVq1PaJgCW9Vuh21Xu+ko+7mKsmsMVcn3bnEgKVHVDpaqXTyqpuxSfZn57o586ycH\nklWs3VJo1R6gq1oNybqnqk9gWzkcSMzMagxaH0bZHEisa3zy9U/ZeV+15tBh29+qpaeWA0mFtHuw\nVO3gqlp62r3YDFNzTNV+q3at9vQPGgeSIeIL6+DyhdX6yYFkiAzbxWbY9tesX0p9Z7uksyXdI2mP\npM115h8n6cuS7pD0bUmn5+a9W9IuSbslvSc3/UmSvibpb9P/jytzH8zMrLnSAomkNcCngZ8HTgPe\nKOm0msU+CMxFxPOBtwCXpuueDrwN2Ai8AHiVpGem62wGZiPiVGA2/WxmZn1SZo1kI7AnIvZGxCPA\n1cC5NcucBnwDICLuBtZJOgF4DrAjIhYi4jDwTeB16TrnAlelf18F/JsS98HMzFooM5CcBNyX+3x/\nOi3vdtIAIWkj8DRgAtgFnCXpyZJGgXOAk9N1ToiI76d//yNwQjnJNzOzIvrd2b4FuFTSHHAncBvw\naETcJenjwA3Aw8Ac8GjtyhERkqLehiVdBFwEMDExwb59+0rahe7Yv39/v5NQKc6P5ZwfRzgvlqtC\nfpQZSB7gSC0CkprGA/kFIuJHwPkAkgTcC+xN510BXJHO+xhJjQbgnySdGBHfl3Qi8GC9L4+IbcA2\ngA0bNsT4+HiXdqs8qyGNveT8WM75cYTzYrl+50eZTVs3A6dKOkXSY4A3ANfkF5B0bDoP4ELgxjS4\nIOkp6f9PJWn++pN0uWuA89K/zwO+UuI+mJlZC6XVSCLisKR3AtcDa4ArI2K3pIvT+ZeRdKpflTZP\n7QYuyG3iS5KeDBwC3hERD6XTtwBfkHQB8F3gl8raBzMza63UPpKIuA64rmbaZbm/twN1H7eOiLMa\nTP8hMN3FZJqZWQcUUbeveqBI+gFJ7aXKxoFq3xHQW86P5ZwfRzgvliszP54WEce3WmgoAslqIOmW\niNjQ73RUhfNjOefHEc6L5aqQH6UOkWJmZoPPgcTMzDriQFId2/qdgIpxfizn/DjCebFc3/PDfSRm\nZtYR10jMzKwjDiQ9IOlx6ftWbk/fr/LbuXnvknR3Ov0TuekfSN/jco+kn+tPysvRKD8krZd0k6Q5\nSbekA3lm6wxsfkDy2gVJt0m6Nv3c8L07g54XUDc/LknPkzvSdxgdm1t26PIjN/29kkLSeG5a7/Mj\nIvyv5H+AgLH077XADuBM4GXA14HHpvOekv5/GsnIyI8FTgH+DljT7/3oQX7cAPx8Ov0c4G+GIT/S\nffz3JMMAXZt+/gSwOf17M/DxYcmLBvnxs8Ax6d8fH/b8SKedTDJyyHeB8X7mh2skPRCJA+nHtem/\nAGaALRHx43S5bADKc4GrI+LHEXEvsIfk/S4DoUl+BPCEdPoTge+lfw90fkiaAH4BuDw3udF7dwY6\nL6B+fkTEDZG8mwjgJpJBYGFI8yP1SeB9JOdNpi/54UDSI2nVdI5ktOKvRcQOkuFhzpK0Q9I3Jb0w\nXbzIu1xWtQb58R7gEkn3Af8N+EC6+KDnx++RXBAWc9MavXdn0PMC6udH3luBr6Z/D2V+SDoXeCAi\nbq9Zti/54UDSIxHxaESsJylJbUxfJ3wM8CSSZp3/QDIYpfqYzJ5pkB8zwK9HxMnAr5O+RmCQSXoV\n8GBE3NpomUjaLIbi9spW+SHpQ8Bh4HM9TVif1MuP9GV/HwR+s28Jq9HvF1sNnYh4SNJfA2eTlBb+\nPL1QfFvSIsm4OS3f5TIoavLjPODd6awvcqQqP8j58WLgNZLOAR4HPEHSH9P4vTuDnBfQID8i4s2S\nfhV4FTCdnjMwhPkB/BFJ/8ftablzAtiZ3pzSn/zodyfSMPwDjgeOTf/+CeBbJCfExcB/TqdPklRJ\nBTyX5R1mexmgDsQm+XEX8NJ0+jRwa/r3QOdHLl9eypHO5UtY3tn+iWHKizr5cTbwHeD4mmWGMj9q\npv89Rzrb+5IfrpH0xokk711ZQ9Kc+IWIuDZ9qdeVknYBjwDnRXI07Jb0BZIT5zDJ+1iOetXwKtYo\nPx4iefXyMcC/kL4qOZL32AxyftRT9707Q5oXAJ8iuTh+LS2F3xQRFw9xftTVr/zwk+1mZtYRd7ab\nmVlHHEjMzKwjDiRmZtYRBxIzM+uIA4mZmXXEgcQqS9KT05GA5yT9o6QHcp//T4/T8u8lfScdfXZW\n0tNq5j9B0v2SPpWb9llJ9+bSvL5mnRdKOizpF3PTfj0dEXmXpD+V9Lh0+kfS756TdIOkn0ynb8xt\n/3ZJr22yD38m6el1pv9qPt2dkvQ8SZ/t1vas+hxIrLIi4ocRsT6SoVQuAz6ZfY6IF/U4ObcBGyLi\n+cCfkYzOm/cR4MY66/2HXJrnsonpMzQfJxnxOJt2EvBr6fecDqwB3pDOviQinp/mxbUcGR5jV7r8\nepKH9v5H+hzOMpKeS/Jg2t52d7xdEXEnMCHpqWV/l1WDA4mtSpIOpP+/NB3w8iuS9kraIulNSt53\ncqekZ6TLHS/pS5JuTv+9uJ3vi4i/joiF9GN+9Fkk/TTJoIo31Fu3gXcBX+LI0CeZY4CfSIPBKOkI\nyBHxo9wyjycdeysiFuLIqLiPo/GYXG8CvpJL8/mS5iV9m2QYjmz6q9NBRG+T9HVJJ0gaUfJelOPT\nZUaUvO/ieEmvT2tPt0vKB9K/5EgQtAHnQGKD4AUkw808B/h3wGREbCQZq+td6TKXktRoXgj8W44e\nkrsdF5COPitpBPgd4DcaLPtf0yapT0p6bLrOScBrga35BSPiAZJRj/8B+D6wPyLyNZaPKhkZ+U3k\nBuyTdIak3cCdwMW5wJL3YuDWdPkTgd9Op/0MyTssMv8LODMifgq4GnhfRCwCf5x+L8ArgNsj4gdp\nOn4uIl4AvCa3nVuAsxrkiQ0YBxIbBDdHxPcjea/L33GkZnAnsC79+xXAp5QMXX8NyWCAY+1+kaQ3\nAxtIxsICeDtwXUTcX2fxD5CMofZCklGe359O/z3g/ekFOr/t40jeJ3EK8JPA49PvAyAiPhTJyMif\nA96Zm74jIp6bfs8Hsn6VGicCP0j/PoPkpWE/iIhHgM/nlpsArpd0J8mI1M9Np18JvCX9+63AZ9K/\n/zfwWUlvI2mKyzyY7oMNAY+1ZYPgx7m/F3OfFzlyjI+QlLT/pdFGJH0G+CngexFxTp35rwA+BLwk\nDVoAm0jeKfN2YAx4jKQDEbE5jrxP5MfptrNaywbg6nTMqHHgHEmHSV7wdW9a0kfSnwMvIqkN5H0O\nuA74T/mJEXFX2uR3OkmNIO8gSdNXK78P/G5EXCPppcBvpdu+T9I/SXo5yYuS3pROv1jSGSQvXrpV\n0k9HxA/T7zpY4PtsALhGYsPiBo40c1F7BxVARJyfdorXCyI/BfwP4DVx5E2WRMSbIuKpEbGOJFD8\nYURsTtc5Mf1fJG843JWuc0pErEvX+TPg7RHxFyRNWmdKGk3XmSYZERlJp+aScy5wdzr9lKxzPb2T\n7Nkko8HWugt4Zvr3DuAlSu6KWwu8PrfcEzky7Ph5Ndu4nCSofTEbCFDSM9Ia0W+S1HiyIcwns/21\nwecaiQ2LXwM+LekOkuP+RpJ+laIuIalxfDGtSfxDRLym+Sp8Lu2gFjDX6vsiYoekPwN2kozcehuw\nLZ29RdKzSGpZ381t62eAzZIOpfPeHhH76mz+f5IMQ/71SN5x8lvAduChNG2Z30r38f8C3yBpZstc\nQ9Kk9ZnctEvSICdglmQIc4CXpd9pQ8Cj/5oNAUk/Afw18OKVDisuaQPJDQtNO9HTmwq+CfxMg45/\nGzBu2jIbAhFxkKRPZUXv75a0meR25Q8UWPypJC/lchAZEq6RmJlZR1wjMTOzjjiQmJlZRxxIzMys\nIw4kZmbWEQcSMzPriAOJmZl15P8DpQvenX75aXgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "flat_lc = lc.flatten(window_length=401)\n", "flat_lc.plot();" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucXWV97/HPd2KoTKcKEuUggwaQqQaEOI0BpLy8RFug\nVqqn3qqVIpDDiC32aG3Uc6w9HD1R6oVWzqQRUHpqS7WoRRsLGq3YGsIlDJAUGGNSBdRibIkOk0rC\n/M4fa63Jysq+rJk9a/bsvb/v12tes/e6/p59Wb/1PGvt51FEYGZmNlt97Q7AzMw6mxOJmZm1xInE\nzMxa4kRiZmYtcSIxM7OWOJGYmVlLnEjMzKwlTiRmZtYSJxIzM2vJE9odwHxYsmRJLF26tN1hzNi+\nfft4whN64i0Ceq+84DL3ik4t8x133LErIp7abLnOK9ksLF26lNtvv73dYczYrl27WLJkSbvDmDe9\nVl5wmXtFp5ZZ0nfLLOemLTMza4kTiZmZtcSJxMzMWuJEYmZmLXEiMTOzljiRmJlZS5xIzMysJU4k\nZmbWEicSM7MFYmRkhJGRkXaHMWNOJGZm1hInEjMza4kTiZmZtcSJxMzMWtITvf+amXWC0dHRdocw\nK66RmJlZS5xIzMysJU4kZmbWEicSMzNriROJmZm1pLJEIukaSQ9L2lpnviT9qaTtku6WNJybd5ak\n+9N5a3LTL5d0X7r85yUdVlX8ZmZWTpU1kk8BZzWYfzZwQvq3GhgFkLQIuDKdvwx4vaRl6TpfAU6K\niJOBceBdlURuZmalVZZIIuJm4N8bLHIu8BeRuAU4TNJRwEpge0TsiIjHgOvSZYmImyJiX7r+LcBg\nVfGbmVk57bxGcjTwQO75g+m0etOL3gx8ubLozMyslI78Zbuk9wD7gE83WGY1SZMZg4OD7Nq1a56i\nmzu7d+9udwjzqtfKCy5zr+j2MrczkTwEHJN7PphOW1xnOgCSfgd4ObAqIqLexiNiPbAeYMWKFbFk\nyZI5C3w+dWrcs9Vr5QWXuVd0c5nb2bR1A/Cm9O6t04DdEfED4DbgBEnHSjoEeF26LJLOAt4JvCIi\nJtsVuJl1hk4dKKrTVFYjkfTXwIuAJZIeBP6IpLZBRKwDNgDnANuBSeD8dN4+SW8FbgQWAddExLZ0\nsx8Hfg74iiSAWyLi4qrKYGZmzVWWSCLi9U3mB3BJnXkbSBJNcfqz5iY6MzObK/5lu5mZtcSJxMzM\nWtKRt/+amZXRqQNFdRrXSMzMrCVOJGZm1hInEjMza4kTibWNfyxm1h2cSMzMrCVOJGZdwLU7aycn\nErMK+QBvvcCJZAHywcfMWjHfxxD/ILGk7E3xD5zmjl/L3uDvTvdzjWQGxsfHK8nyroF0Jr9vthCN\njIwwPj4+r/t0jWQWigePhXqm1cqZYKtnke08C+2UM+CRkRH27t3LVVdd1fK25qOstV7XTnmtrVpO\nJHOkii/UmjVrWLx4celtlolhruLMznqGhoYO2lZ+XmbVqlUAbNy4se42i+VttI+qzOT1KZ711Vq3\nlbgbxTLT93GuT346tSa2kBJfvVgWUoxluWlrARodHZ3RAaLRl3o2zXEzrRo3imG2zT/5GIplyLbZ\nyrZrrTcXTQIzfb137NjRloNyrdegV5vqFnK5671Pq1atmp5ea5n5PPkC10hKGx0dbXiwHBsbo7+/\nv6V9FM9Etm3b1mjxg9Yt1gJqbbvevvLTxsfHmZycrLl8LY0OvmXO2jM7duxAUtMDer159ZJDXr0v\nV3Gfzc4Ks/nNXu+5bgbKx1lv37XiKLPcXHyGa+276gNas5pbo5pzvW0B814bbqQYV7Nl5psTyQzk\nm1waqXXwqvUBLbPdqampAw7w2Qc7+4JkB/z+/n4mJydLnVEXD5jZ/suejTdaLotrYmKi7oGp7Fl7\nvmy1TExM1IxrcnKSqakpVq1aNf06F2POYsi/lo3KlJc/sOS3UzQ+Pj7dpFc86GfbPO644w6aVnw/\nGjULNoqt+BlZvnx50/VrxTFXsu1edtllpfZRpjl0rhTfq1onUzB31w6zfQ4NDZVOAGW/2+1Ifk4k\nJeU/QPnE0Kimkl83nzyyA37+4JKdEU5NTXHqqafS19fH1NQUwPQBeXJykrGxsQP2ly2T/9DnY1u1\natUB+8p/QaamptiyZQurVq1q+mUdGxsD9h+M6n3Jsn1l+5mcnJxOBFm580kmO+gD9PX1TceY6e/v\nn449K1dx39kBIL+vTD65NJJPSlls+fcnXwtYtWoVExMT0/GOjY2xatWq6X3USvBbtmyhr6+v5sF8\nJrJtZ+WsVTvJn3gUywXJ+97X1zcdW2ZqauqAMpx++unTy27atGn685S93pAc5LOmllonR9kBukzt\nqRh/LbNNLtl7lP8s1XqfZxPTXKj3Ps60qTX7bs33dRYlQ6d3txUrVsTtt9/e0jbyX8687KCZPS6e\n4eYPqtlBJPtC5z/EWRKpJTtg1ZufNzAwUDe2TPEAU2t/U1NTDAwMHHCgL24/kz+w5Q+w+ZizA1e2\nTLat4vZrxVFmWqa/v396+/mDZTGOYmKpFUPxdS/Gm20r/1oW18k/r7W9Y489lm3bth0QZ7bd7D3b\nsmULAJs3b55OYvVek3qPs/3Vet/zy9V7T4aHhxkdHZ1OLvnp+c9Ttn7+5CEvXz7Yn4jgwJp5MRlD\n7c9z/vvW6MaPLKYsARUTSfF7na/F1fpOz/QAvWvXLl772tce8NrUO15k285OToqf1/xno3gimz9h\nmYvaiaQ7ImJFs+VcI5mh4hcj/wHMzu7hwANq9iXLPtD5WkSZdukyCSQfX7Z89qGamppibGyM5cuX\nTx+Uyuyv1gGl3kEzP71WvMVlittqFEeZafXOsIv7qDWt7P6Lr0e97dd7Xuu13LZtW83XLzuzzNdg\nigfxZvuoFX8t+eXqJfYtW7Zw6qmnHjS92ASUrT8xMcHAwMBBiTZfvlo10Fqx1appFmOYmJiY/v7l\nT5qyz3tfXx+Tk5MHNEPmX+fsefGkJJOd6PX19R3QNDvTg3St8uQTXX7b2bShoaHpGDNZAsnKW6Zl\npEqukZSUTxB5+TOw/JkwHPxFbnRGWGv5uVbrbN4WtnqfmW6Rr/1C49pq/vM7PDw8PT1fm8+vly2f\nrznWqqWWjSmbV6wRlJGvkeRjbbT94klfo9gHBgbYuHEjp59+OsAB5c1qVbNRtkZS2e2/kq6R9LCk\nrXXmS9KfStou6W5Jw7l5Z0m6P523Jjf9KZK+Iunb6f/Dq4o/k2X9emdzExMTpc9y6y0zNTU1Lwd4\nJ5HO081JBPZ/JicmJg74jmXfq3rGxsbYsmXLQU3C+YNsfnp24M++q422XdxeLbO9TTy/7YmJiYOa\ngvPX1Gqt2+g4lNVY88ej+frOV/k7kk8BZzWYfzZwQvq3GhgFkLQIuDKdvwx4vaRl6TprgI0RcQKw\nMX1emZGREbZs2TL9wWt2wTYzX4nBrJs0+87UarZr1JSYn1ZsGiqjVsKZyXGgqJWmxfyy9ea385hT\nWSKJiJuBf2+wyLnAX0TiFuAwSUcBK4HtEbEjIh4DrkuXzda5Nn18LfAb1URfm5ODWWeay+/u1NTU\nrC5il4mhU48x7fxl+9HAA7nnD6bT6k0HODIifpA+/iFwZJUBznfHZ2bWGXxsOFDH3rUVESGp7p0C\nklaTNJkxODjIrl27ZrOP2QdoZl1rYmKCCy+8kLVr15Zafvfu3RVHVN9ll102q+PfTLQzkTwEHJN7\nPphOW1xnOsC/SToqIn6QNoM9XG/jEbEeWA/JXVtLliyZcYB79uyZ8Tpm1ht27tzJbI4r820+Ymxn\n09YNwJvSu7dOA3anzVa3ASdIOlbSIcDr0mWzdc5LH58H/F2VAc5lv0O9JP+jxfwdKZ2gr6+v5o8u\nzVrxmte8pt0hVKrK23//GtgE/KKkByVdIOliSReni2wAdgDbgU8AbwGIiH3AW4EbgXuBz0RE1nvh\nWuBlkr4NvDR9XpmFcOtlmQPxXB/46m2v0X76+voOiHUuEkg7DurFLjNmWw4no+4327u3ulGVd229\nPiKOiojFETEYEVdHxLqIWJfOj4i4JCKOj4jnRsTtuXU3RMRQOu/9uek/johVEXFCRLw0IhrdFday\nhXA2Xa87kIGBAQYGBg56XlZ24O/r62N4eHh63WKZs+d9fX1s3LjxoH0MDAwwPDzMpk2bDqjB9ff3\nMzAwcMAvs/PJJttvI/lfqjdKKtm8ZuXPlismvfz87G6cbFv9/f3Tyzfaf1HxINMsyeZf/+Hh4Zqv\nVT72mSbubH5xO8UY6r3/+TjK7rPbucViv97+JLRgtmectb6sefUOcNlfdoAeGhpi48aNB3UCmB1Q\niwed/Layriv6+/unf/U6NDR00C+H82XMf2nyB9b82Xu2jVplHh4ent5GlmBGR0cPOFjm181+qZtP\nSENDQ2zevPmgxJfvt6kYd/Y65JPa8uXL2bRp00HvRf4XwENDQ9Ply16rrHzZAXXz5s1s3rx5OsZs\n39nrmz0eHh5m+fLl0+suX7685gE9XxbYfxKR75tp+fLlLF++/KD3I7/9fNmz+dn/0dFRNm7cOJ34\n80k4K2sxMWflzV6/bF69Dig7McHMJuaZdEQ5F/utdwI02+XmUsfetTUfavUTBEx/6er1W1WrM79s\nerZdOLCb9Hznjdk+IoJHH330gGVqdRMO+z/U+dsSix3O5fvuKSp2fpffZjY/P61Zl/r1OpXLd5yX\nxZgtl3UjUbZJMXud89urVY568dbq3bdsWfLvQ63XpNhhXr2hBe6//34kHdRxYD5R1CtHsfPNYqeY\nWR9NtTqprFW+4uNa8cKBPe8WOxas12FkNg3K/VaiHd35VNkdTb0bd7LPf74fvuzEIDu+5Lthyjq5\nzPchVnyd8l2szFdCcSJpoHgAyd6Y4pct3wNq9iGAAzuMyx/AGx3ksl47a8VS5gdQ9cbfKKPegXIm\n6+e7GM9PLys7oBaTWXH42rJxNprfKBnWGpCqrHod6BX7PBodHeXCCy9k586ddeOq9z7W6/U236tt\nfl+13pdavd7W2k+tuPLP8z3o5jsQzL4vxSagYl9YmUbJBzho2fy0Momn1jJZrStLtPlk0mybM/mO\nHXrooTz66KP09SVd8tfqdql4sph//4ufnSzuWt/3dnTi6ETSQPEglHWlnb1B2ZkCUPOLVOsaQKM3\nuPhFvf/++5t2ulamej2XY0I3WyffTfdsFLdfb3/N4pjpL48bLd9sDPYy0+pte+3atQ1vz2xWjuL8\nesvWGr+j2eiN2dgjjcb+qFfzhIOb6rLvyqGHHjpdC8uSXraPrIfhWs1mxRaAYkeMtVoA8tOLJ3j1\nylMceydf86uiptToc9foZKJRjXe+OZGUVOvNLo6l0OzMttG0WvN37dpV6h7wWomi7IGx3aqIpUwt\npOoY5sNCizv7DuRHnqzl/vvvn35c6y65egfqrLaQNS2Pj48f8Di/jfx1pXpj/eSv8eVrwbUuord6\nR2LWtJWPC2ZXoy5b25jPz4cTSYs64WA9E3MVd6eWfyGZz9ewyn0Vz5qPO+44Fi9eDBxco6p1PadW\nk1rxZKnWmXs2P2vWy9ar1RV8cdvZturte6adQFZdi2n3982JZBYWSvJYCHG0+wNsC9dsmvPytYuZ\nHCzrNd9m22y2Tr1rcPW2n9Us8gNlNZLVtGY6jkmncCIxsxkpc02m0XrN+n0qc31rtgfiuboIPdNb\nf7PrQt2WQDJOJGZWU6OL7Atdo2RX5mA+27v9elXTRCLpaRHxcGHaL0bE/fXWMbPe1OoZd7eeseev\nC3WjMjWSb0r6nxHxGQBJbwcuIBm90MzMmmh2m3enK5NIXgSsl/RqkoGk7iUZxdDMzKx5X1tp1+7/\nAJwOLAWujYj2d4trZmYLQplrJF8Fvg+cRDLg1NWSbo6Id1QdnJmZLXxlmrY+HhFfSB8/IukFwLsq\njMnMbMHzuO37lWna+kLh+b6IuKy6kMzMusuaNWu6+nbhMk1bPwUifXoIyZjqExHx5CoDMzPrFjt2\n7EBSu8OoTNNEEhG/kD1W8kqcC5xWZVBmZgudh9rdb0bdWKbD434B+NWK4jEz6zp79uzp6sRTpmnr\nVbmnfcAK4D8ri8jMrAN4zPb9yty19eu5x/uAfyVp3jIzsxJOPPHE3u4iJSLOn49AzMw6yUx7AO5m\ndROJpD9j/91aB4mI36skIjOzDjCTDiZ7+a6t21vduKSzgCuARcBVEbG2MP9w4BrgeJLrLm+OiK3p\nvEuBiwABn4iIj6XTlwPrgCeSNLW9JSJubTVWM7OZqDXOfT293PvvpyNi32w3LGkRcCXwMuBB4DZJ\nN0TEv+QWezcwFhGvlPTsdPlVkk4iSSIrgceAf5D0pYjYDnwI+OOI+LKkc9LnL5ptnGZm1ppGt/9O\nn+WnzVwztRLYHhE7IuIx4DoOvki/DPgaQETcByyVdCTwHGBzREymyewbQHb3WABPSh8/maQfMDMz\na5NGNZJ8g94Zs9j20cADuecPAqcWlrmLJEF8U9JK4JnAILAVeL+kI4A9wDnsb2p7G3CjpD8hSYQv\nqBm8tBpYDTA4ONh0eM+FaPfu3e0OYV71WnnBZe5ke/fuBZoPHQzwne98B0kdeRwqo1EiqXuhfQ6t\nBa6QNAbcA9wJPB4R90r6IHAT8CgwBjyerjMC/H5EXC/pNcDVwEsPCj5iPbAeYMWKFdGpg8p0atyz\n1WvlBZe5U2XXPMqU5fjjj2fx4sVdUe5aGiWSZ0u6m6Rmcnz6mPR5RMTJTbb9EEm385nBdNq0iPgJ\ncD5Md7+yE9iRzruaJEkg6QMkNRqA84BL08efBa5qEoeZmVWoUSJ5Tovbvg04QdKxJAnkdcBv5ReQ\ndBgwmV5DuRC4OU0u02PFS3oGSfNX1r/X94EXAv8IvAT4dotxmplVqmdv/42I77ay4YjYJ+mtwI0k\nt/9eExHbJF2czl9HkqyulRTANpKx4DPXp9dI9gKXRMQj6fSLSJrDnkByy/DqVuI0M6vanj172h1C\npcp0kTJrEbEB2FCYti73eBNQ8+ehEXFmnen/BPzSHIZpZmYtmFHvv2ZmZkWV1kjMzLrVTLpI6XaN\n+tq6h8Z9bTW7a8vMzIAvfvGLXXvrLzSukbw8/X9J+v//pf/fUF04ZmadYSZ9bXW7pndtSXpZRDwv\nN2uNpC3AmqqDMzNbqMbHx9sdwoJR5mK7JJ2Re/KCkuuZmVkPKHOx/QLgGklPTp8/Ary5upDMzKyT\nlBkh8Q7glCyRRER39LhmZtYCj5C4X9NEknbr/gHg6RFxtqRlwOlpX1hmZj3J10j2K3Ot41Mk3Zw8\nPX0+TtKVu5lZzxoaGnKtJFUmkSyJiM8AU5D0ocX+Lt3NzKzHlUkkj6adJwaApNMAXycxMzOg3F1b\nbwduIBmT5J+BpwKvrjQqM7MFztdI9mtaI0nv2nohyZC2/w04MSLuqjowM7OFLLs+kv3CvZc1TSSS\nvgNcGBHbImJrROyV9KV5iM3MbMEaHx9ncnKy3WEsCGWukewFXizpk5IOSacdXWFMZmZdY2RkhDVr\nurtHqTKJZDIiXgvcC3wzHfq2bq/AZma9YGhoiP7+/naHsSCUudgugIj4UNpZ403AUyqNyszMOkaZ\nRPLe7EFEfFXSrwLnVReSmVlnGBoacjfyNB7Y6tkRcR/wkKThwmxfbDeznuYEsl+jGsnbgYuAD9eY\nF8BLKonIzKxDeHCrRKOBrS5K/794/sIxM7NO06hp61WNVoyIzzXbuKSzgCuARcBVEbG2MP9w4Brg\neOA/gTdHxNZ03qUkNSIBn4iIj+XW+12SIYAfB/4+It7ZLBYzM6tGo6atX28wL4CGiUTSIuBK4GXA\ng8Btkm6IiH/JLfZuYCwiXinp2enyqySdRJJEVgKPAf8g6UsRsV3Si4FzgVMi4meSntakjGZmVqFG\nTVvnt7jtlcD2iNgBIOk6kgSQTyTLgLXp/u6TtDQd/+Q5wOaImEzX/QbwKuBDwAiwNiJ+lq73cItx\nmplVZnx8nIju/uldmdt/kfRrwInAE7NpEfG/mqx2NPBA7vmDwKmFZe4iSRDflLQSeCYwCGwF3p/2\nOrwHOAe4PV1nCDhT0vtJmsPeERG3lSmHmdlcKnORfWhoiL17985DNO1TZoTEdUA/8GLgKuA3gVvn\naP9rgSskjQH3AHcCj0fEvZI+SPLjx0eBMfaPgfIEkh9EngY8H/iMpOOikPIlrQZWAwwODrJr1645\nCnn+7N7dW73191p5wWXudFnXJ2vXrq27zN69e9m3b19HHoPKKlMjeUFEnCzp7oj4Y0kfBr5cYr2H\ngGNyzwfTadMi4ifA+QCSBOwEdqTzrgauTud9gKRGQ/r/c2niuFXSFLAE+FFh2+uB9QArVqyIJUuW\nlAh54enUuGer18oLLnMnW7x4MdC4PGWW6XRlEsme9P+kpKcDPwaOKrHebcAJko4lSSCvA34rv4Ck\nw0j68noMuBC4OU0uSHpaRDyc9u31KpIaCMAXSGpHX5c0BBwCdG+qN7OONjo62tW1ESiXSL6UHvAv\nB7aQ3LF1VbOVImKfpLeSjPe+CLgmIrZJujidv47kovq1kgLYBlyQ28T16TWSvcAlEfFIOv0a4BpJ\nW0nu6Dqv2KxlZmbzp2kiiYjL0ofXp+OQPDEiSjVyRsQGYENh2rrc400kF89rrXtmnemPAW8ss38z\ns3YbGRlh7969XHVV0/PvjlXmYvsi4NeApdnykoiIj1QbmpmZdYIyTVtfJLnN9h5gqtpwzMy6i39H\nkhiMiJMrj8TMzDpSmUTyZUm/EhE3VR6NmVmX8Q8SE7cAn5fUR3IHlYCIiCdVGpmZmXWEMonkI8Dp\nwD2+zdbMbGZ8jSTxALDVScTMbObctJXYAfyjpC8DP8sm+vZfMzODcolkZ/p3SPpnZmZ4iN1Mw0SS\n/hjxFyLiHfMUj5mZdZi+RjMj4nHgjHmKxcyso4yMjDAyMtLuMNquTNPWmKQbgM+SjA0ClBuz3cys\n1/murcQTSbqOf0luWtMx283MrDeU6f231bHbzcysizW8RgIgaVDS5yU9nP5dL2lwPoIzM7OFr0zT\n1ieBvwJenT5/YzrtZVUFZWbWCXz7b6JMInlqRHwy9/xTkt5WVUBmZt2kF37Z3rRpC/ixpDdKWpT+\nvZHk4ruZmVmpGsmbgT8DPkpyt9a3AF+AN7Oel/2GpFET1+joKLt27ZqvkNqizF1b3wVeMQ+xmJl1\nlPHx8XaHsCDUTSSS3ttgvYiIyyqIx8zMOkyjGsmjNab9PHABcATgRGJmZvUvtkfEh7M/YD1wKMm1\nkeuA48psXNJZku6XtF3SmhrzD09/o3K3pFslnZSbd6mkrZK21bpLTNLbJYWkJWViMTOba0NDQwwN\nDTVcZmRkhDVrDjr8dZWGd21Jeoqk/w3cTVJ7GY6IP4yIh5ttOO05+ErgbGAZ8HpJywqLvRsYi4iT\ngTcBV6TrngRcBKwETgFeLulZuW0fA/wK8L1SpTQzs8rUTSSSLgduA34KPDci3hcR/zGDba8EtkfE\njoh4jKQmc25hmWXA1wAi4j5gqaQjgecAmyNiMiL2Ad8AXpVb76PAO0nuIjMzszZqVCN5O/B04H8A\n35f0k/Tvp5J+UmLbR5MM05t5MJ2WdxdpgpC0EngmMAhsBc6UdISkfuAc4Jh0uXOBhyLirhIxmJlV\nZnx83Hdu0eBie0SU+bFiq9YCV0gaA+4B7gQej4h7JX0QuInkov8Y8HiaVN5N0qzVkKTVwGqAwcHB\njryPe/fu3e0OYV71WnnBZe50WffwjY4ve/fuZd++fR15DCqrzA8SZ+sh0lpEajCdNi0ifkL640ZJ\nIhnSd0c672rg6nTeB0hqNMcDxwJ3JYszCGyRtDIifljY9nqSmwRYsWJFLFnSmdfkOzXu2eq18oLL\n3Mn27NkDNC7Pzp07iYiuKXMtVSaS24ATJB1LkkBeB/xWfgFJhwGT6TWUC4Gb0+SCpKdFxMOSnkHS\n/HVaRDwCPC23/r8CKyKie1O9mdkCV1kiiYh9kt4K3AgsAq6JiG2SLk7nryO5qH6tpAC2kfxGJXO9\npCOAvcAlaRIxM1sw+vv72x3CglBljYSI2ABsKExbl3u8Cah5E3ZEnFli+0tbDNHMbNaa/YakV1Sa\nSMzMetnIyAgTExP09c3HvUvt092lMzNro/Hxcfr6+jjxxBPbHUqlnEjMzKwlTiRmZtYSJxIzM2uJ\nL7abmc1Ss+5Rsru6Lrusu0fdcCIxM6tINgRvN3ePAm7aMjOzFjmRmJlVZGRkhJGRkXaHUTknEjMz\na4kTiZnZLG3cuJGhoaGeqHU04kRiZjZLIyMjHtgKJxIzM2uRE4mZmbXEvyMxM2vB5ORk3eatXmn2\nciIxM2tBf39/3XFJJicn5zma9nDTlpnZLI2OjjYc3Gpqaoqpqal5jKg9nEjMzGap2V1bfX19XT+o\nFbhpy8ysMr1QGwHXSMzMrEWukZiZVaQXmrXANRIzs1lrdrG9v7+f/v7+eYyoPZxIzMysJZUmEkln\nSbpf0nZJa2rMP1zS5yXdLelWSSfl5l0qaaukbZLelpt+uaT70nU+L+mwKstgZtbI+Ph43Tu3Jicn\ne+K3JJUlEkmLgCuBs4FlwOslLSss9m5gLCJOBt4EXJGuexJwEbASOAV4uaRnpet8BTgpXWcceFdV\nZTAza2ZoaKhu85abtlq3EtgeETsi4jHgOuDcwjLLgK8BRMR9wFJJRwLPATZHxGRE7AO+AbwqXe6m\ndBrALcBghWUwM2vINZJq79o6Gngg9/xB4NTCMneRJIhvSloJPJMkMWwF3i/pCGAPcA5we419vBn4\nm1o7l7QaWA0wODjYkWMm7969u90hzKteKy+4zN0gSxSNjjHdVuaidt/+uxa4QtIYcA9wJ/B4RNwr\n6YPATcCjwBjweH5FSe8B9gGfrrXhiFgPrAdYsWJFLFmypLJCVKlT456tXisvuMydLvvRYaMyPfnJ\nT+6qMhdVmUgeAo7JPR9Mp02LiJ8A5wNIErAT2JHOuxq4Op33AZIaDenz3wFeDqyKiKisBGZmszQy\nMsLU1FS1nUv6AAALxElEQVRP/JakyhLeBpwg6VhJhwCvA27ILyDpsHQewIXAzWlyQdLT0v/PIGn+\n+qv0+VnAO4FXRET3Nz6aWUeoNdxuX19fT1xsr6xGEhH7JL0VuBFYBFwTEdskXZzOX0dyUf1aSQFs\nAy7IbeL69BrJXuCSiHgknf5x4OeArySVGG6JiIurKoeZ2WyMj4/3TF9blV4jiYgNwIbCtHW5x5uA\nmvfNRcSZdaY/q9Z0M7N2GBgYqHlnVi/crZXp/sY7M7MKTU5OMjU1xfj4+AHNW1ltpFEXKt3CicTM\nrAXNEsbo6Oh8htMWTiRmZnNgbGys3SG0jROJmdkcKTZv9QonEjOzFmS/E5mammJycpKxsbGeSyZO\nJGZmLVi+fPn04174zUgtTiRmZi3IX0zP7uDqteslTiRmZnMku4Mr+z8wMNDOcOaNE4mZmbXEicTM\nzFriRGJmZi1xIjEzq0gvdI8CTiRmZi0bHh5udwht5URiZtaiemO29wonEjOzFtVqwhoeHu6JDhvB\nicTMrBK9VEtxIjEzq0CvXGgHJxIzM2uRE4mZWYtGR0enewHuRb1bcjOzOdTf309fX9/0Xy95QrsD\nMDPrBsVrIr1yxxY4kZiZzYleShxFlda/JJ0l6X5J2yWtqTH/cEmfl3S3pFslnZSbd6mkrZK2SXpb\nbvpTJH1F0rfT/4dXWQYzM2usskQiaRFwJXA2sAx4vaRlhcXeDYxFxMnAm4Ar0nVPAi4CVgKnAC+X\n9Kx0nTXAxog4AdiYPjczszapskayEtgeETsi4jHgOuDcwjLLgK8BRMR9wFJJRwLPATZHxGRE7AO+\nAbwqXedc4Nr08bXAb1RYBjMza6LKRHI08EDu+YPptLy7SBOEpJXAM4FBYCtwpqQjJPUD5wDHpOsc\nGRE/SB//EDiymvDNzKyMdl9sXwtcIWkMuAe4E3g8Iu6V9EHgJuBRYAx4vLhyRISkqLVhSauB1QCD\ng4Ps2rWroiJUZ/fu3e0OYV71WnnBZe4V3V7mKhPJQ+yvRUBS03gov0BE/AQ4H0CSgJ3AjnTe1cDV\n6bwPkNRoAP5N0lER8QNJRwEP19p5RKwH1gOsWLEilixZMkfFml+dGvds9Vp5wWXuFd1c5iqbtm4D\nTpB0rKRDgNcBN+QXkHRYOg/gQuDmNLkg6Wnp/2eQNH/9VbrcDcB56ePzgL+rsAxmZtZEZTWSiNgn\n6a3AjcAi4JqI2Cbp4nT+OpKL6temzVPbgAtym7he0hHAXuCSiHgknb4W+IykC4DvAq+pqgxmZtZc\npddIImIDsKEwbV3u8SagZheZEXFmnek/BlbNYZhmZtYCRdS8Vt1VJP2IpPbSaZYAnXeXwOz1WnnB\nZe4VnVrmZ0bEU5st1BOJpFNJuj0iVrQ7jvnSa+UFl7lXdHuZe6uLSjMzm3NOJGZm1hInkoVtfbsD\nmGe9Vl5wmXtFV5fZ10jMzKwlrpGYmVlLnEgWiJmMsyJpkaQ7JX1pPmOca2XKLOkYSV+X9C/p2DSX\ntiPWVpUYm0eS/jSdf7ek4XbEOZdKlPkNaVnvkfQtSae0I8651KzMueWeL2mfpN+cz/iq4kSycMxk\nnJVLgXvnJapqlSnzPuDtEbEMOA24pMa4NgtaybF5zgZOSP9WAx093F7JMu8EXhgRzwUuo8OvI5Qs\nc7Zc1iltV3AiWThKjbMiaRD4NeCqeYqrSk3LHBE/iIgt6eOfkiTQ4nAEC12ZsXnOBf4iErcAh6Wd\nknaqpmWOiG9FxH+kT28h6di1k5V5nwF+F7ieOh3OdiInkoWj7DgrHwPeCUzNS1TVmtHYMpKWAs8D\nNlcb1pwrMzZPmWU6yUzLcwHw5Uojql7TMks6GnglHV7jLGr3eCQ9RdJXgf9SY9Z78k/qjbMi6eXA\nwxFxh6QXVRPl3Gq1zLntDJCcxb0t6yHauoOkF5Mkkl9udyzz4GPAH0bEVDJyRndwIplHEfHSevMk\nlRln5QzgFZLOAZ4IPEnSX0bEGysKuWVzUGYkLSZJIp+OiM9VFGqVmo7NU3KZTlKqPJJOJmmmPTvt\nkLWTlSnzCuC6NIksAc6RtC8ivjA/IVbDTVsLR9NxViLiXRExGBFLScZ3+dpCTiIlNC1zOuDZ1cC9\nEfGReYxtLjUdmyd9/qb07q3TgN25Zr9OVGY8omcAnwN+OyLG2xDjXGta5og4NiKWpt/hvwXe0ulJ\nBJxIFpK1wMskfRt4afocSU+XtKHhmp2rTJnPAH4beImksfTvnPaEOzsRsQ/Ixua5F/hMNjZPNj4P\nyXALO4DtwCeAt7Ql2DlSsszvBY4A/m/6vt7epnDnRMkydyX/st3MzFriGomZmbXEicTMzFriRGJm\nZi1xIjEzs5Y4kZiZWUucSGzBknRE7pbfH0p6KPf8W/Mcy39PeyC+W9JGSc8szH+SpAclfTw37VOS\nduZiXl5Y56AeYCX9ftrL8VZJfy3pien0y9J9j0m6SdLT0+krc9u/S9IrG5ThbyUdV2P67+TjbpWk\n50r61FxtzxY+JxJbsCLixxGxPCKWA+uAj2bPI+IF8xzOncCKiDiZ5IdkHyrMvwy4ucZ6f5CLeSyb\nWKsH2LQfpt9L93MSsIjkR20Al0fEyelr8SWS32AAbE2XXw6cBfy5pIN6rJB0IrAoInbMtOAzFRH3\nAIPpDw6tBziRWEeSNJH+f5Gkb0j6O0k7JK1VMs7FrUrGuTg+Xe6pkq6XdFv6d8ZM9hcRX4+IyfTp\nAT3VSvolkg4nZ9IteL0eYJ8AHJomg37g++n+8/2L/TwQ6fTJ9IdwkHSbU++HYW8g13OApPMljUu6\nleRHn9n0X5e0Wcl4N1+VdKSkPiVjxjw1XaZPyXgbT5X06rT2dJekfCL9IvuToHU5JxLrBqcAFwPP\nIfkV/FBErCTpw+l302WuIKnRPB/4r7TWDf90T7WS+oAPA++os+z/SZukPirp59J1avYAGxEPAX8C\nfA/4AUk3Kfkay/slPUCSFN6bm36qpG3APcDFucSSdwZwR7r8UcAfp9N+mWTsjMw/AadFxPNIukF/\nZ0RMAX+Z7heSXgjuiogfpXH8akScArwit53bgTPrvCbWZZxIrBvclo5b8jPgO+yvGdwDLE0fvxT4\nuKQxkv6PnqSkR+EZkfRGko73Lk8nvQXYEBEP1lj8XcAQ8HzgKcAfptOne4AtbPtwkvErjgWeDvx8\nuj8AIuI9EXEM8GmSrjiy6Zsj4sR0P+/KrqsUHAX8KH18KvCPEfGjdNyMv8ktNwjcKOke4A+AE9Pp\n1wBvSh+/Gfhk+vifgU9JuoikKS7zcFoG6wHu/de6wc9yj6dyz6fY/xnvIznT/s96G5H0SZLxTr4f\nEQf15yXppSTd378wTVoApwNnSnoLMAAcImkiItbkOl38WbrtrNZSswdYYDGwMz3TR9LngBeQ1Aby\nPk3SN9cf5SdGxL1pk99JJDWCvD0kTV/N/BnwkYi4QclQBe9Lt/2Akt6aX0IygNMb0ukXSzqVZLC1\nOyT9UtqL7xPTfVoPcI3EesVN7G/mongHFUBEnJ9eFK+VRJ4H/Dnwioh4OLfOGyLiGWlvru8gGeVw\nTbrOUel/kYz+uDVdp14PsN8DTpPUn66zinRIZUkn5MI5F7gvnX5sdnE9vZPs2cC/1ij/vcCz0seb\ngRcquStuMfDq3HJPZn/X5+dxoKtIktpnI+LxdJ/HpzWi95LUeLJu1Iey8lr3c43EesXvAVdKupvk\nc38zyXWVsi4nqXF8Nq1JfC8iXtF4FT6dXqAWMNZsfxGxWdLfAltIxqq/k/3jmK+V9Isktazv5rb1\ny8AaSXvTeW+JiF01Nv/3wIuAr6bjv7wP2AQ8ksaWeV9axv8AvkbSzJa5gaRJ65O5aZenSU7ARuCu\ndPqL031aD3Dvv2Y9QNKhwNeBM7LaxCy2sYLkhoWGF9HTmwq+AfxynQv/1mXctGXWAyJiD8k1lVmN\nAy9pDcntyu8qsfgzgDVOIr3DNRIzM2uJayRmZtYSJxIzM2uJE4mZmbXEicTMzFriRGJmZi1xIjEz\ns5b8f4M3fCAqTad3AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "folded_lc = flat_lc.fold(period=3.5225)\n", "folded_lc.plot();" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYHVWZ7/HvjwAD7Si3AAdoIInSKiBEjQFEjkJEARWO\nnCODA8KgwEMLCjMwGvQ8jB4etYXxwiinkRNA5gyKKKIZThQ0KnjhLg1JBNpMMtzEgTgShEZI4D1/\n1NpNZWfv3ZVU1969u3+f5+mn67Kq6l17766316raqxQRmJmZbaxNOh2AmZl1NycSMzMrxYnEzMxK\ncSIxM7NSnEjMzKwUJxIzMyvFicTMzEpxIjEzs1KcSMzMrJRNOx1AO0yfPj1mzJjR6TA2yNq1a9l0\n0ynx9oxynacG17l73HXXXasiYvuxynVfzTbCjBkzuPPOOzsdxgZZtWoV06dP73QYbeU6Tw2uc/eQ\n9GCRcu7aMjOzUpxIzMysFCcSMzMrxYnEzMxKcSIxM7NSnEjMzKwUJxIzMyvFicTMzEpxIjHL6e/v\np7+/v9NhmHUVJxIrrdnJt9VJuegJe2P23Q79/f3MmzdvwiSd8Xitx+t4tvG69XWtLJFIulzS45KW\nNlkvSf8kabmkeyW9IbfuMEkPpHXzc8u3lfQjSb9Nv7epKv5O2NAPUb78xp6YG803OkEWOXG2KlMf\n6/Dw8Jhxzps3r2W54eHhUsmoUbn66UbbjhV/keM122/R939DX8/67cb6HEwWk7VeE40iopodS/8V\neBr454jYu8H6I4CPAEcA+wEXRcR+kqYBw8ChwCPAHcD7I+I3ki4A/jMiBlKC2SYiPj5WLHPmzImJ\nNNZW7YM9ODjYdPkhhxyCJBYvXsy8efMYGRlh9uzZAOudNPr6+kan6/fZ7Nj1BgcHm56QavsfGhqi\np6eHxYsXj5bt6+sbPWb99vm48nHXx1uLac2aNaxcuXJ0fW3/9dsvXrx4nbrkjz/Wseq3q9+mr6+P\n4eHh9V7v+nrWb9/s9WtW19p+Z86cyYIFC9aLv36f9a9l/XvRKp5a7M32Va/ZNo3232h6LO0cd2pD\n4qpSF4+1dVdEzBmrXGWDNkbEzZJmtChyFFmSCeBWSVtL2gmYASyPiBUAkq5OZX+Tfr8tbX8l8DNg\nzERShUZ/9EXKj1VueHiYefPmNTxOs5NUfnn+ZFafhGrHrq0fGRmhp6eHvr6+9barPynmTyq1+EZG\nRkbL1J+khoaG1jkxN6trbbuRkRG23HJLnn32WYaHh9dJFvXx1p/canXs6ekZLQfQ09OzTkz5uPPz\n+dczX8/8tvlYG732/f39DA0NAYwmoKGhodETfu14tWPPnj2bvr4+1qxZM7p9/nWpX1ZbXqtn7X2r\nvTeNtq1/3fPzY32W6lt7Q0NDHHDAAescdyz1r1fttV2zZs1o8qxapxPIeCp6DtnQc9N46OQ1kl2A\nh3Pzj6RlzZYD7BgRj6Xp3wM7Vh0kNO7+KdI106yLB9b977hRd9CsWbNGT8J9fX3ccssto+sWL148\neqIF1vuvfWRkhJGRkXVOXPUnhr6+vtETXH5Zfl99fX0sXrx4NNbBwcHR9X19faMnzEYf3NqJsn6b\n2v5qy2uvY09PD7NmzWq4z1pcs2fPXqfetXK1OtbK9fT0jJatHSP/ftVOxPWvXb6etW3zZVq11oaH\nh0ePW6vf7NmzR0+8tXL5YwOsWLFivS68Wvn6ZAaM1qv2OsybN2+dz0/+WPnXO7++WfyNjg/Z563+\ns1KvVmfonu6kVl3D7bqu1OjvP18u//7WPg/Nup5ry4aHh0f/kWuXyrq2AFKL5PomXVvXAwMR8Ys0\nv5isdTEDOCwiTk7LPwDsFxFnSHoyIrbO7eOPEdHwOomkU4FTAXp7e9949913b3Q95s+f33D5wMDA\n6PoVK1Ywa9YsBgYGRuefffZZttxyS6655pp1ykGWKACWLVsGwJZbbjm6PcDq1avZaqutRo9d22/9\ncfPz9bE2iil/jEZ1bLTvRsdptt2Gym/fqM4bEmeRY2zM+rGU2f6cc87hoYceGv085PeTf8/q1+XX\nA+uUqbch72Gzz0/9MWvHqz9+bb62LP/Zq9ltt91Gn81R/z7mj1v7G9prr70Kvb/5+jZyzDHHrLO/\n+mMCo3+ztfps7Gei3urVq/nc5z63XpytXs+BgQGOOeaYdbaZNWvW6Dljr732anis+n2Usf322xfq\n2upkIvka8LOI+Gaaf4Cs22oG8KmIeGdafi5ARHyuViYiHkvdYD+LiFePFUeZayT1XSvwUtcF0LBr\nBVinr71+Oaz/H2B+f1BNn+rG9he3q5+5SJ0nSp/3eMnXuWzdWl3/2th9Ndq26HWh+q7BmpkzZ7Jy\n5cp1ulbry9T//dR3O+bluzKbdbvV9gWs080LL3Vx9vX1jXZFNuuSzV8r2xC19znflVuLpVHda8vq\np+tfg/zrU9+NCOU/Ux2/RlLAQuCMdA1kP2B1ShBPAHtImgk8ChwL/HVumxOBgfT7++0MuFnzPv/h\nrX0QG12AHs83eGNs7LEm0kl7IsUy3srWbTxfm1b7yq9rdKNDbXmzz/jJJ58MNP57anSinjdvHkND\nQ+uchHt6etZJILVt6/95q+8WbBZvq5tG8mpdxGVe63ycrWJotq5RnI3OMe1U5V1b3yRrYUwH/gP4\nB2AzgIi4RJKArwKHASPASRFxZ9r2CODLwDTg8oj4TFq+HXANsBvwIHBMRPznWLFsbItkY07245Ug\nuvUujzJc5+5X5PN/8skns9lmm613QmzUmmjUCml0B1lN/Q0U9dfUmh0vv49mdwM2uoutqCrf5yr/\nKe14iyQi3j/G+gBOb7JuEbCowfI/APPW36IzGr2Bk/k/ZrOxFPn8DwwMrHdSrU8OjfZZ//fW6Fj1\nt4ZvTLxjtcYm2o0EE+GcU+k1komiqu+RVPmfwGT7T7UI13lqaHedG/2djtf1qKLbd+v73PEWyVQw\nEf4TMLPWNrTVYRvOicTMbAM5Ea3LgzaamVkpTiRmZlaKE4mZWRt0y9AxG8OJxMzMSvHFdjOzNpjM\nF+jdIjEzs1KcSMzMrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxYnEzMxKcSIxM7NSnEjM\nzKwUJxIzMyvFicTMzEpxIjEzs1KcSMzMrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxYnE\nzMxKcSIxM7NSnEjMzDZCf38//f39nQ5jQnAiMTOzUjbtdABmZt1ocHCw0yFMGJW2SCQdJukBScsl\nzW+wfhtJ10m6V9LtkvbOrTtT0lJJyySdlVs+W9KtkoYk3SlpbpV1MDOz1ipLJJKmARcDhwN7Au+X\ntGddsU8AQxGxD3ACcFHadm/gFGAusC/wbkmvSttcAHw6ImYD56V5MzPrkCpbJHOB5RGxIiKeB64G\njqorsyfwE4CIuB+YIWlH4LXAbRExEhFrgZuAo9M2AbwiTW8F/K7COpiZ2RiqTCS7AA/n5h9Jy/Lu\nISWI1EW1O9ALLAUOkrSdpB7gCGDXtM1ZwIWSHgb+ETi3shqYmdmYOn2xfQC4SNIQsAS4G3ghIu6T\n9HngRuAZYAh4IW3TD/xtRFwr6RjgMuDt9TuWdCpwKkBvby+rVq2qvDLjafXq1Z0Ooe1c56nBdZ58\nqkwkj/JSKwKylsaj+QIR8RRwEoAkASuBFWndZWRJAkmfJWvRAJwInJmmvw0saHTwiLgUuBRgzpw5\nMX369NIVardujLks13lqcJ0nlyq7tu4A9pA0U9LmwLHAwnwBSVundQAnAzen5IKkHdLv3ci6v76R\nyv0OeGuaPgT4bYV1MDOzMVTWIomItZLOAG4ApgGXR8QySael9ZeQXVS/UlIAy4AP5XZxraTtgDXA\n6RHxZFp+Cll32KbAn0ndV2Zm1hmVXiOJiEXAorpll+SmbwH6mmx7UJPlvwDeOI5hmplZCR4ixcys\nYvPnz5/U43I5kZiZWSmdvv3XzGzSGxgY8F1bZmZmzTiRmJlZKU4kZmZWihOJmZmV4kRiZmaljJlI\nakOV1C17dTXhmJlZtynSIvl5GmUXAElnA9dVF5KZmXWTIt8jeRtwqaT3ATsC95E9tMrMzGzsFklE\nPAb8EDgAmAFcGRFPVxyXmZl1iTFbJJJ+TDZ0+95kzxe5TNLNEXFO1cGZmdnEV+QayVcj4oSIeDIi\nlgBvBib3477MzKywMVskEfG9uvm1wPmVRWRmZl2lSNfWn4BIs5sDmwFPR8RWVQZmZmbdoUiL5OW1\n6fRc9aOA/asMyszMuscGfbM9Mt8D3llRPGZm1mWKdG0dnZvdBJhD9qx0MzOzQl9IfE9uei3w72Td\nW2ZmZoWukZzUjkDMzKw7NU0kkr7CS3drrSciPlpJRGZm1lVatUjubFsUZmbWtVolkqvSlw/NzMya\nanX77+21idTNZWZmtp5WiUS56QOrDsTMzLpTq0TS9EK7mZlZTatrJK+RdC9Zy+SVaZo0HxGxT+XR\nmZnZhNcqkby2bVGYmVnXappIIuLBdgZiZmbdaYMGbTQzM6vnRGJmZqVUmkgkHSbpAUnLJc1vsH4b\nSddJulfS7ZL2zq07U9JSScsknVW33Uck3Z/WXVBlHczMrLVWY20tofVYWy3v2pI0DbgYOBR4BLhD\n0sKI+E2u2CeAoYh4r6TXpPLzUkI5BZgLPA/8UNL1EbFc0sFkow/vGxHPSdqhUE3NzKwSre7aenf6\nfXr6/X/T7+MK7nsusDwiVgBIuposAeQTyZ7AAEBE3C9phqQdye4Yuy0iRtK2NwFHAxcA/cBARDyX\ntnu8YDxmZlaBpl1bEfFgunPr0Ij4WEQsST/zgXcU2PcuwMO5+UfSsrx7yBIEkuYCuwO9wFLgIEnb\nSeoBjgB2Tdv0pXW3SbpJ0psKxGJmZhUp8mArSTowIn6ZZt7M+F1bGQAukjQELAHuBl6IiPskfR64\nEXgGGAJeyMW8Ldlz498EXCNpVkSs0w0n6VTgVIDe3l5WrVo1TiG3x+rVqzsdQtu5zlOD6zz5FEkk\nHwIul7RVmn8S+GCB7R7lpVYEZC2NR/MFIuIp4CTIshWwEliR1l0GXJbWfZasRUP6/d2UOG6X9CIw\nHXiibt+XApcCzJkzJ6ZPn14g5ImlG2Muy3WeGiZLnfv7+wEYHBwcs+xkqXMjRZ6QeBewby2RRETR\n1HoHsIekmWQJ5Fjgr/MFJG0NjETE88DJwM0puSBph4h4XNJuZN1f+6fNvgccDPxUUh+wOdBdzQ0z\ns0lkzESSLn5/Ftg5Ig6XtCdwQGoxNBURayWdAdwATAMuj4hlkk5L6y8hu6h+paQAlpG1fmqulbQd\nsAY4PSKeTMsvJ2shLSW7o+vE+m4tM7N2KNISmQqKdG19HbgC+GSaHwa+Rep2aiUiFgGL6pZdkpu+\nhezieaNtD2qy/Hng+AJxm5lZGxS5aD49Iq4BXoSspcFLF77NzGyKK5JInkldTAEgaX9gct+CYGZm\nhRXp2jobWEj2TJJfAtsD76s0KjMz6xqF7tqS9Fbg1WQPtXogItZUHpmZmXWFMbu2JP0bcHJELIuI\npRGxRtL1bYjNzMy6QJFrJGuAgyVdIWnztKx+qBMzM5uiiiSSkYj4K+A+4OfpC4L+3oaZmQEFx9oC\niIgLJP2abPyrbSuNyszMukaRRHJebSIifizpncCJ1YVkZmbdpNWDrV4TEfcDj0p6Q91qX2w3MzOg\ndYvkbLKnFH6hwboADqkkIjMz6ypNE0lEnJJ+H9y+cMzMrNu06to6utWGEfHd8Q/HzMy6Tauurfe0\nWBeAE4mZmbXs2jqpnYGYmVl3KnL7L5LeBewFbFFbFhH/q6qgzMysexQZa+sS4K+Aj5B9OfF9wO4V\nx2VmZl2iyBApb46IE4A/RsSngQNo8lRDMzObeookkmfT7xFJO5MN4rhTdSGZmVk3KXKN5HpJWwMX\nAr8mu2NrQaVRmZlNEv39/axZs4YFCybvabPIg63OT5PXpueQbBERftSumZkBBRKJpGnAu4AZtfKS\niIgvVhuamVn3GxwcZNWqVZ0Oo1JFurb+FfgzsAR4sdpwzMys2xRJJL0RsU/lkZiZWVcqctfWDyS9\no/JIzMysKxVpkdwKXCdpE7JbfwVERLyi0sjMzKwrFEkkXyT7EuKSiPCz2s3MbB1FurYeBpY6iZiZ\nWSNFWiQrgJ9J+gHwXG2hb/81MzMolkhWpp/N04+ZmdmolokkfRnx5RFxTpviMTOzLtPyGklEvAAc\n2KZYzMysCxW52D4kaaGkD0g6uvZTZOeSDpP0gKTlkuY3WL+NpOsk3Svpdkl759adKWmppGWSzmqw\n7dmSQtL0IrGYmVk1ilwj2QL4A3BIbtmYz2xP3WIXA4cCjwB3SFoYEb/JFfsEMBQR75X0mlR+Xkoo\npwBzgeeBH0q6PiKWp33vCrwDeKhA/GZmVqEio/9u7LPb5wLLI2IFgKSrgaOAfCLZExhIx7lf0gxJ\nOwKvBW6LiJG07U3A0cAFabsvAR8Dvr+RsZmZ2Tgp8qjd3tT99Hj6uVZSb4F970L2HZSaR9KyvHvI\nEgSS5pI9wrcXWAocJGk7ST3AEcCuqdxRwKMRcU+BGMzMrGJFurauAL5B9qx2gOPTskPH4fgDwEWS\nhshGF74beCEi7pP0eeBG4BlgCHghJZVPkHVrtSTpVOBUgN7e3q4bxnn16qn3yBfXeWpwnSefIolk\n+4i4Ijf/9UYXvxt4lNSKSHrTslER8RRwEoAkkX1fZUVadxlwWVr3WbIWzSuBmcA9WXF6gV9LmhsR\nv6/b96XApQBz5syJ6dO775p8N8Zclus8NbjOk0uRu7b+IOl4SdPSz/FkF9/Hcgewh6SZkjYHjgUW\n5gtI2jqtAzgZuDklFyTtkH7vRtb99Y2IWBIRO0TEjIiYQZZc3lCfRMzMrH2KtEg+CHyF7AJ3AL8i\ntSJaiYi1ks4AbgCmAZdHxDJJp6X1l5BdVL9SUgDLgA/ldnGtpO3IRhw+PSKeLF4tMzNrlyJ3bT0I\nHLkxO4+IRcCiumWX5KZvAfqabHtQgf3P2Ji4zMxs/DRNJJLOa7FdRMT5FcRjZmZdplWL5JkGy15G\n1v20HeBEYmZmzRNJRHyhNi3p5cCZZNdGrga+0Gw7MzObWsYa/Xdb4O+A44Arye6Q+mM7AjMzs+7Q\n6hrJhWS33V4KvC4inm5bVGZm1jVafY/kbGBn4H8Cv5P0VPr5k6Sn2hOemZlNdK2ukRT5sqKZmU1x\nThZmZiX09/fT39/f6TA6yonEzMxKKTJEipmZNTE4ONjpEDrOLRIzMyvFicTMzEpxIjEzs1KcSMzM\nrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxd9sNzOrUH9/P2vWrGHBggWdDqUybpGYmVkp\nbpGYmVVocHCQVatWdTqMSrlFYmZmpTiRmJlZKU4kZmZWihOJmZmV4kRiZmalOJGYmVkpTiRmZlaK\nE4mZmZXiRGJmZqU4kZiZWSlOJGZmVkqliUTSYZIekLRc0vwG67eRdJ2keyXdLmnv3LozJS2VtEzS\nWbnlF0q6P21znaStq6yDmZm1VlkikTQNuBg4HNgTeL+kPeuKfQIYioh9gBOAi9K2ewOnAHOBfYF3\nS3pV2uZHwN5pm2Hg3KrqYGZmY6uyRTIXWB4RKyLieeBq4Ki6MnsCPwGIiPuBGZJ2BF4L3BYRIxGx\nFrgJODqVuzEtA7gV6K2wDmZmNoYqh5HfBXg4N/8IsF9dmXvIEsTPJc0FdidLDEuBz0jaDngWOAK4\ns8ExPgh8q9HBJZ0KnArQ29vbdcM4r169utMhtJ3rPDW4zpNPp59HMgBcJGkIWALcDbwQEfdJ+jxw\nI/AMMAS8kN9Q0ieBtcBVjXYcEZcClwLMmTMnpk+fXlklqtKNMZflOk8NrvPkUmUieRTYNTffm5aN\nioingJMAJAlYCaxI6y4DLkvrPkvWoiHN/w3wbmBeRERlNTAzG0N/fz+QPcBqqqryGskdwB6SZkra\nHDgWWJgvIGnrtA7gZODmlFyQtEP6vRtZ99c30vxhwMeAIyNipML4zcysgMpaJBGxVtIZwA3ANODy\niFgm6bS0/hKyi+pXSgpgGfCh3C6uTddI1gCnR8STaflXgb8AfpQ1Yrg1Ik6rqh5mZq1M5ZZITaXX\nSCJiEbCobtkluelbgL4m2x7UZPmrGi03M7PO8DfbzcysFCcSM7OK9Pf3j16Mn8ycSMzMrJROf4/E\nzGzSql2I77YvRG8ot0jMzKwUJxIzMyvFicTMzEpxIjEzq8BUuWMLnEjMzKwk37VlZlaBqTR0ilsk\nZmYlTaVurEacSMzMrBR3bZmZlTSVurEacYvEzMxKcSIxMxsHU/k6iROJmZmV4mskZmbjIH+dZKo9\nx90tEjOzceCuLTMzGxdTrTUCTiRmZuOiljiGh4c7HEn7OZGYmY2jvr4+gCnVzeVEYmY2TgYHB6dU\nl1aN79oyMxtnUy2ZuEViZmalOJGYmVkpTiRmZlaKE4mZmZXiRGJmZqU4kZiZWSlOJGZmVooTiZmZ\nleJEYmZmpSgiOh1D5SQ9ATzY6Tg20HRgVaeDaDPXeWpwnbvH7hGx/ViFpkQi6UaS7oyIOZ2Oo51c\n56nBdZ583LVlZmalOJGYmVkpTiQT16WdDqADXOepwXWeZHyNxMzMSnGLxMzMSnEimSAkbSvpR5J+\nm35v06LsNEl3S7q+nTGOtyJ1lrSrpJ9K+o2kZZLO7ESsZUk6TNIDkpZLmt9gvST9U1p/r6Q3dCLO\n8VSgzselui6R9CtJ+3YizvE0Vp1z5d4kaa2k/9HO+KriRDJxzAcWR8QewOI038yZwH1tiapaReq8\nFjg7IvYE9gdOl7RnG2MsTdI04GLgcGBP4P0N6nA4sEf6ORXo6kfsFazzSuCtEfE64Hy6/DpCwTrX\nyn0euLG9EVbHiWTiOAq4Mk1fCfy3RoUk9QLvAha0Ka4qjVnniHgsIn6dpv9ElkB3aVuE42MusDwi\nVkTE88DVZHXPOwr458jcCmwtaad2BzqOxqxzRPwqIv6YZm8Fetsc43gr8j4DfAS4Fni8ncFVyYlk\n4tgxIh5L078HdmxS7svAx4AX2xJVtYrWGQBJM4DXA7dVG9a42wV4ODf/COsnwyJlusmG1udDwA8q\njah6Y9ZZ0i7Ae+nyFme9TTsdwFQi6cfAf2mw6pP5mYgISevdTifp3cDjEXGXpLdVE+X4Klvn3H7+\nkuy/uLMi4qnxjdI6SdLBZInkLZ2OpQ2+DHw8Il6U1OlYxo0TSRtFxNubrZP0H5J2iojHUpdGo2bv\ngcCRko4AtgBeIelfIuL4ikIubRzqjKTNyJLIVRHx3YpCrdKjwK65+d60bEPLdJNC9ZG0D1k37eER\n8Yc2xVaVInWeA1ydksh04AhJayPie+0JsRru2po4FgInpukTge/XF4iIcyOiNyJmAMcCP5nISaSA\nMeus7C/uMuC+iPhiG2MbT3cAe0iaKWlzsvduYV2ZhcAJ6e6t/YHVuW6/bjRmnSXtBnwX+EBEDHcg\nxvE2Zp0jYmZEzEh/w98BPtztSQScSCaSAeBQSb8F3p7mkbSzpEUdjaw6Rep8IPAB4BBJQ+nniM6E\nu3EiYi1wBnAD2c0C10TEMkmnSTotFVsErACWA/8H+HBHgh0nBet8HrAd8L/T+3pnh8IdFwXrPCn5\nm+1mZlaKWyRmZlaKE4mZmZXiRGJmZqU4kZiZWSlOJGZmVooTiU1YkrbL3fL7e0mP5uZ/1eZY/i6N\nQHyvpMWSdq9b/wpJj0j6am7Z1yWtzMU8u26b9UaAlfS3aZTjpZK+KWmLtPz8dOwhSTdK2jktn5vb\n/z2S3tuiDt+RNKvB8r/Jx12WpNdJ+vp47c8mPicSm7Ai4g8RMTsiZgOXAF+qzUfEm9sczt3AnIjY\nh+yLZBfUrT8fuLnBdn+fi3motrDRCLBpHKaPpuPsDUwj+1IbwIURsU96La4n+w4GwNJUfjZwGPA1\nSeuNWCFpL2BaRKzY0IpvqIhYAvSmLxzaFOBEYl1J0tPp99sk3STp+5JWSBpQ9pyL25U95+KVqdz2\nkq6VdEf6OXBDjhcRP42IkTS7zki1kt5INuDkhgwL3mwE2E2BLVMy6AF+l46fH1/sZUCk5SPpi3CQ\nDZvT7Ithx5EbOUDSSZKGJd1O9qXP2vL3SLpN2fNufixpR0mbKHtmzPapzCbKnrexvaT3pdbTPZLy\nifRfeSkJ2iTnRGKTwb7AacBryb4F3xcRc8nGcPpIKnMRWYvmTcB/p9ww/KMj1UraBPgCcE6Tsp9L\nXVJfkvQXaZuGI8BGxKPAPwIPAY+RDZOSb7F8RtLDZEnhvNzy/SQtA5YAp+USS96BwF2p/E7Ap9Oy\nt5A9O6PmF8D+EfF6smHQPxYRLwL/ko4L2SgE90TEEymOd0bEvsCRuf3cCRzU5DWxScaJxCaDO9Jz\nS54D/o2XWgZLgBlp+u3AVyUNkY1/9AplIwpvEEnHkw28d2Fa9GFgUUQ80qD4uUAf8CZgW+Djafno\nCLB1+96G7PkVM4GdgZel4wEQEZ+MiF2Bq8iG4qgtvy0i9krHObd2XaXOTsATaXo/4GcR8UR6bsa3\ncuV6gRskLQH+HtgrLb8cOCFNfxC4Ik3/Evi6pFPIuuJqHk91sCnAo//aZPBcbvrF3PyLvPQZ34Ts\nP+0/N9uJpCvInnfyu4hYbzwvSW8nG/7+rSlpARwAHCTpw8BfAptLejoi5ucGXXwu7bvWamk4Aiyw\nGbAy/aePpO8CbyZrDeRdRTY21z/kF0bEfanLb2+yFkHes2RdX2P5CvDFiFio7FEFn0r7fljZaM2H\nkD3A6bi0/DRJ+5E9bO0uSW9Mo/hukY5pU4BbJDZV3MhL3VzU30EFEBEnpYvijZLI64GvAUdGxOO5\nbY6LiN3SaK7nkD3lcH7aZqf0W2RPf1yatmk2AuxDwP6SetI280iPVJa0Ry6co4D70/KZtYvr6U6y\n1wD/3qD+9wGvStO3AW9VdlfcZsD7cuW24qWhz09kXQvIktq3I+KFdMxXphbReWQtntow6n21+trk\n5xaJTRXcNqu0AAAA8klEQVQfBS6WdC/Z5/5msusqRV1I1uL4dmpJPBQRR7behKvSBWoBQ2MdLyJu\nk/Qd4Ndkz6q/m5eeYz4g6dVkrawHc/t6CzBf0pq07sMRsarB7v8f8Dbgx+n5L58CbgGeTLHVfCrV\n8Y/AT8i62WoWknVpXZFbdmFKcgIWA/ek5QenY9oU4NF/zaYASVsCPwUOrLUmNmIfc8huWGh5ET3d\nVHAT8JYmF/5tknHXltkUEBHPkl1T2ajnwEuaT3a78rkFiu8GzHcSmTrcIjEzs1LcIjEzs1KcSMzM\nrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMyslP8PnBK2S4SKLBgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "binned_lc = folded_lc.bin(binsize=10)\n", "binned_lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Or we can do these all in a single (long) line!" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XuYHVWZ7/HvjwAD7Si3AAdoIInSKiBEjQFEjkJEARWO\nnCODA8KgwEMLCjMwGvQ8jB4etYXxwiinkRNA5gyKKKIZThQ0KnjhLg1JBNpMMtzEgTgShEZI4D1/\n1NpNZWfv3ZVU1969u3+f5+mn67Kq6l17766316raqxQRmJmZbaxNOh2AmZl1NycSMzMrxYnEzMxK\ncSIxM7NSnEjMzKwUJxIzMyvFicTMzEpxIjEzs1KcSMzMrJRNOx1AO0yfPj1mzJjR6TA2yNq1a9l0\n0ynx9oxynacG17l73HXXXasiYvuxynVfzTbCjBkzuPPOOzsdxgZZtWoV06dP73QYbeU6Tw2uc/eQ\n9GCRcu7aMjOzUpxIzMysFCcSMzMrxYnEzMxKcSIxM7NSnEjMzKwUJxIzMyvFicTMzEpxIjHL6e/v\np7+/v9NhmHUVJxIrrdnJt9VJuegJe2P23Q79/f3MmzdvwiSd8Xitx+t4tvG69XWtLJFIulzS45KW\nNlkvSf8kabmkeyW9IbfuMEkPpHXzc8u3lfQjSb9Nv7epKv5O2NAPUb78xp6YG803OkEWOXG2KlMf\n6/Dw8Jhxzps3r2W54eHhUsmoUbn66UbbjhV/keM122/R939DX8/67cb6HEwWk7VeE40iopodS/8V\neBr454jYu8H6I4CPAEcA+wEXRcR+kqYBw8ChwCPAHcD7I+I3ki4A/jMiBlKC2SYiPj5WLHPmzImJ\nNNZW7YM9ODjYdPkhhxyCJBYvXsy8efMYGRlh9uzZAOudNPr6+kan6/fZ7Nj1BgcHm56QavsfGhqi\np6eHxYsXj5bt6+sbPWb99vm48nHXx1uLac2aNaxcuXJ0fW3/9dsvXrx4nbrkjz/Wseq3q9+mr6+P\n4eHh9V7v+nrWb9/s9WtW19p+Z86cyYIFC9aLv36f9a9l/XvRKp5a7M32Va/ZNo3232h6LO0cd2pD\n4qpSF4+1dVdEzBmrXGWDNkbEzZJmtChyFFmSCeBWSVtL2gmYASyPiBUAkq5OZX+Tfr8tbX8l8DNg\nzERShUZ/9EXKj1VueHiYefPmNTxOs5NUfnn+ZFafhGrHrq0fGRmhp6eHvr6+9barPynmTyq1+EZG\nRkbL1J+khoaG1jkxN6trbbuRkRG23HJLnn32WYaHh9dJFvXx1p/canXs6ekZLQfQ09OzTkz5uPPz\n+dczX8/8tvlYG732/f39DA0NAYwmoKGhodETfu14tWPPnj2bvr4+1qxZM7p9/nWpX1ZbXqtn7X2r\nvTeNtq1/3fPzY32W6lt7Q0NDHHDAAescdyz1r1fttV2zZs1o8qxapxPIeCp6DtnQc9N46OQ1kl2A\nh3Pzj6RlzZYD7BgRj6Xp3wM7Vh0kNO7+KdI106yLB9b977hRd9CsWbNGT8J9fX3ccssto+sWL148\neqIF1vuvfWRkhJGRkXVOXPUnhr6+vtETXH5Zfl99fX0sXrx4NNbBwcHR9X19faMnzEYf3NqJsn6b\n2v5qy2uvY09PD7NmzWq4z1pcs2fPXqfetXK1OtbK9fT0jJatHSP/ftVOxPWvXb6etW3zZVq11oaH\nh0ePW6vf7NmzR0+8tXL5YwOsWLFivS68Wvn6ZAaM1qv2OsybN2+dz0/+WPnXO7++WfyNjg/Z563+\ns1KvVmfonu6kVl3D7bqu1OjvP18u//7WPg/Nup5ry4aHh0f/kWuXyrq2AFKL5PomXVvXAwMR8Ys0\nv5isdTEDOCwiTk7LPwDsFxFnSHoyIrbO7eOPEdHwOomkU4FTAXp7e9949913b3Q95s+f33D5wMDA\n6PoVK1Ywa9YsBgYGRuefffZZttxyS6655pp1ykGWKACWLVsGwJZbbjm6PcDq1avZaqutRo9d22/9\ncfPz9bE2iil/jEZ1bLTvRsdptt2Gym/fqM4bEmeRY2zM+rGU2f6cc87hoYceGv085PeTf8/q1+XX\nA+uUqbch72Gzz0/9MWvHqz9+bb62LP/Zq9ltt91Gn81R/z7mj1v7G9prr70Kvb/5+jZyzDHHrLO/\n+mMCo3+ztfps7Gei3urVq/nc5z63XpytXs+BgQGOOeaYdbaZNWvW6Dljr732anis+n2Usf322xfq\n2upkIvka8LOI+Gaaf4Cs22oG8KmIeGdafi5ARHyuViYiHkvdYD+LiFePFUeZayT1XSvwUtcF0LBr\nBVinr71+Oaz/H2B+f1BNn+rG9he3q5+5SJ0nSp/3eMnXuWzdWl3/2th9Ndq26HWh+q7BmpkzZ7Jy\n5cp1ulbry9T//dR3O+bluzKbdbvV9gWs080LL3Vx9vX1jXZFNuuSzV8r2xC19znflVuLpVHda8vq\np+tfg/zrU9+NCOU/Ux2/RlLAQuCMdA1kP2B1ShBPAHtImgk8ChwL/HVumxOBgfT7++0MuFnzPv/h\nrX0QG12AHs83eGNs7LEm0kl7IsUy3srWbTxfm1b7yq9rdKNDbXmzz/jJJ58MNP57anSinjdvHkND\nQ+uchHt6etZJILVt6/95q+8WbBZvq5tG8mpdxGVe63ycrWJotq5RnI3OMe1U5V1b3yRrYUwH/gP4\nB2AzgIi4RJKArwKHASPASRFxZ9r2CODLwDTg8oj4TFq+HXANsBvwIHBMRPznWLFsbItkY07245Ug\nuvUujzJc5+5X5PN/8skns9lmm613QmzUmmjUCml0B1lN/Q0U9dfUmh0vv49mdwM2uoutqCrf5yr/\nKe14iyQi3j/G+gBOb7JuEbCowfI/APPW36IzGr2Bk/k/ZrOxFPn8DwwMrHdSrU8OjfZZ//fW6Fj1\nt4ZvTLxjtcYm2o0EE+GcU+k1komiqu+RVPmfwGT7T7UI13lqaHedG/2djtf1qKLbd+v73PEWyVQw\nEf4TMLPWNrTVYRvOicTMbAM5Ea3LgzaamVkpTiRmZlaKE4mZWRt0y9AxG8OJxMzMSvHFdjOzNpjM\nF+jdIjEzs1KcSMzMrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxYnEzMxKcSIxM7NSnEjM\nzKwUJxIzMyvFicTMzEpxIjEzs1KcSMzMrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxYnE\nzMxKcSIxM7NSnEjMzDZCf38//f39nQ5jQnAiMTOzUjbtdABmZt1ocHCw0yFMGJW2SCQdJukBScsl\nzW+wfhtJ10m6V9LtkvbOrTtT0lJJyySdlVs+W9KtkoYk3SlpbpV1MDOz1ipLJJKmARcDhwN7Au+X\ntGddsU8AQxGxD3ACcFHadm/gFGAusC/wbkmvSttcAHw6ImYD56V5MzPrkCpbJHOB5RGxIiKeB64G\njqorsyfwE4CIuB+YIWlH4LXAbRExEhFrgZuAo9M2AbwiTW8F/K7COpiZ2RiqTCS7AA/n5h9Jy/Lu\nISWI1EW1O9ALLAUOkrSdpB7gCGDXtM1ZwIWSHgb+ETi3shqYmdmYOn2xfQC4SNIQsAS4G3ghIu6T\n9HngRuAZYAh4IW3TD/xtRFwr6RjgMuDt9TuWdCpwKkBvby+rVq2qvDLjafXq1Z0Ooe1c56nBdZ58\nqkwkj/JSKwKylsaj+QIR8RRwEoAkASuBFWndZWRJAkmfJWvRAJwInJmmvw0saHTwiLgUuBRgzpw5\nMX369NIVardujLks13lqcJ0nlyq7tu4A9pA0U9LmwLHAwnwBSVundQAnAzen5IKkHdLv3ci6v76R\nyv0OeGuaPgT4bYV1MDOzMVTWIomItZLOAG4ApgGXR8QySael9ZeQXVS/UlIAy4AP5XZxraTtgDXA\n6RHxZFp+Cll32KbAn0ndV2Zm1hmVXiOJiEXAorpll+SmbwH6mmx7UJPlvwDeOI5hmplZCR4ixcys\nYvPnz5/U43I5kZiZWSmdvv3XzGzSGxgY8F1bZmZmzTiRmJlZKU4kZmZWihOJmZmV4kRiZmaljJlI\nakOV1C17dTXhmJlZtynSIvl5GmUXAElnA9dVF5KZmXWTIt8jeRtwqaT3ATsC95E9tMrMzGzsFklE\nPAb8EDgAmAFcGRFPVxyXmZl1iTFbJJJ+TDZ0+95kzxe5TNLNEXFO1cGZmdnEV+QayVcj4oSIeDIi\nlgBvBib3477MzKywMVskEfG9uvm1wPmVRWRmZl2lSNfWn4BIs5sDmwFPR8RWVQZmZmbdoUiL5OW1\n6fRc9aOA/asMyszMuscGfbM9Mt8D3llRPGZm1mWKdG0dnZvdBJhD9qx0MzOzQl9IfE9uei3w72Td\nW2ZmZoWukZzUjkDMzKw7NU0kkr7CS3drrSciPlpJRGZm1lVatUjubFsUZmbWtVolkqvSlw/NzMya\nanX77+21idTNZWZmtp5WiUS56QOrDsTMzLpTq0TS9EK7mZlZTatrJK+RdC9Zy+SVaZo0HxGxT+XR\nmZnZhNcqkby2bVGYmVnXappIIuLBdgZiZmbdaYMGbTQzM6vnRGJmZqVUmkgkHSbpAUnLJc1vsH4b\nSddJulfS7ZL2zq07U9JSScsknVW33Uck3Z/WXVBlHczMrLVWY20tofVYWy3v2pI0DbgYOBR4BLhD\n0sKI+E2u2CeAoYh4r6TXpPLzUkI5BZgLPA/8UNL1EbFc0sFkow/vGxHPSdqhUE3NzKwSre7aenf6\nfXr6/X/T7+MK7nsusDwiVgBIuposAeQTyZ7AAEBE3C9phqQdye4Yuy0iRtK2NwFHAxcA/cBARDyX\ntnu8YDxmZlaBpl1bEfFgunPr0Ij4WEQsST/zgXcU2PcuwMO5+UfSsrx7yBIEkuYCuwO9wFLgIEnb\nSeoBjgB2Tdv0pXW3SbpJ0psKxGJmZhUp8mArSTowIn6ZZt7M+F1bGQAukjQELAHuBl6IiPskfR64\nEXgGGAJeyMW8Ldlz498EXCNpVkSs0w0n6VTgVIDe3l5WrVo1TiG3x+rVqzsdQtu5zlOD6zz5FEkk\nHwIul7RVmn8S+GCB7R7lpVYEZC2NR/MFIuIp4CTIshWwEliR1l0GXJbWfZasRUP6/d2UOG6X9CIw\nHXiibt+XApcCzJkzJ6ZPn14g5ImlG2Muy3WeGiZLnfv7+wEYHBwcs+xkqXMjRZ6QeBewby2RRETR\n1HoHsIekmWQJ5Fjgr/MFJG0NjETE88DJwM0puSBph4h4XNJuZN1f+6fNvgccDPxUUh+wOdBdzQ0z\ns0lkzESSLn5/Ftg5Ig6XtCdwQGoxNBURayWdAdwATAMuj4hlkk5L6y8hu6h+paQAlpG1fmqulbQd\nsAY4PSKeTMsvJ2shLSW7o+vE+m4tM7N2KNISmQqKdG19HbgC+GSaHwa+Rep2aiUiFgGL6pZdkpu+\nhezieaNtD2qy/Hng+AJxm5lZGxS5aD49Iq4BXoSspcFLF77NzGyKK5JInkldTAEgaX9gct+CYGZm\nhRXp2jobWEj2TJJfAtsD76s0KjMz6xqF7tqS9Fbg1WQPtXogItZUHpmZmXWFMbu2JP0bcHJELIuI\npRGxRtL1bYjNzMy6QJFrJGuAgyVdIWnztKx+qBMzM5uiiiSSkYj4K+A+4OfpC4L+3oaZmQEFx9oC\niIgLJP2abPyrbSuNyszMukaRRHJebSIifizpncCJ1YVkZmbdpNWDrV4TEfcDj0p6Q91qX2w3MzOg\ndYvkbLKnFH6hwboADqkkIjMz6ypNE0lEnJJ+H9y+cMzMrNu06to6utWGEfHd8Q/HzMy6Tauurfe0\nWBeAE4mZmbXs2jqpnYGYmVl3KnL7L5LeBewFbFFbFhH/q6qgzMysexQZa+sS4K+Aj5B9OfF9wO4V\nx2VmZl2iyBApb46IE4A/RsSngQNo8lRDMzObeookkmfT7xFJO5MN4rhTdSGZmVk3KXKN5HpJWwMX\nAr8mu2NrQaVRmZlNEv39/axZs4YFCybvabPIg63OT5PXpueQbBERftSumZkBBRKJpGnAu4AZtfKS\niIgvVhuamVn3GxwcZNWqVZ0Oo1JFurb+FfgzsAR4sdpwzMys2xRJJL0RsU/lkZiZWVcqctfWDyS9\no/JIzMysKxVpkdwKXCdpE7JbfwVERLyi0sjMzKwrFEkkXyT7EuKSiPCz2s3MbB1FurYeBpY6iZiZ\nWSNFWiQrgJ9J+gHwXG2hb/81MzMolkhWpp/N04+ZmdmolokkfRnx5RFxTpviMTOzLtPyGklEvAAc\n2KZYzMysCxW52D4kaaGkD0g6uvZTZOeSDpP0gKTlkuY3WL+NpOsk3Svpdkl759adKWmppGWSzmqw\n7dmSQtL0IrGYmVk1ilwj2QL4A3BIbtmYz2xP3WIXA4cCjwB3SFoYEb/JFfsEMBQR75X0mlR+Xkoo\npwBzgeeBH0q6PiKWp33vCrwDeKhA/GZmVqEio/9u7LPb5wLLI2IFgKSrgaOAfCLZExhIx7lf0gxJ\nOwKvBW6LiJG07U3A0cAFabsvAR8Dvr+RsZmZ2Tgp8qjd3tT99Hj6uVZSb4F970L2HZSaR9KyvHvI\nEgSS5pI9wrcXWAocJGk7ST3AEcCuqdxRwKMRcU+BGMzMrGJFurauAL5B9qx2gOPTskPH4fgDwEWS\nhshGF74beCEi7pP0eeBG4BlgCHghJZVPkHVrtSTpVOBUgN7e3q4bxnn16qn3yBfXeWpwnSefIolk\n+4i4Ijf/9UYXvxt4lNSKSHrTslER8RRwEoAkkX1fZUVadxlwWVr3WbIWzSuBmcA9WXF6gV9LmhsR\nv6/b96XApQBz5syJ6dO775p8N8Zclus8NbjOk0uRu7b+IOl4SdPSz/FkF9/Hcgewh6SZkjYHjgUW\n5gtI2jqtAzgZuDklFyTtkH7vRtb99Y2IWBIRO0TEjIiYQZZc3lCfRMzMrH2KtEg+CHyF7AJ3AL8i\ntSJaiYi1ks4AbgCmAZdHxDJJp6X1l5BdVL9SUgDLgA/ldnGtpO3IRhw+PSKeLF4tMzNrlyJ3bT0I\nHLkxO4+IRcCiumWX5KZvAfqabHtQgf3P2Ji4zMxs/DRNJJLOa7FdRMT5FcRjZmZdplWL5JkGy15G\n1v20HeBEYmZmzRNJRHyhNi3p5cCZZNdGrga+0Gw7MzObWsYa/Xdb4O+A44Arye6Q+mM7AjMzs+7Q\n6hrJhWS33V4KvC4inm5bVGZm1jVafY/kbGBn4H8Cv5P0VPr5k6Sn2hOemZlNdK2ukRT5sqKZmU1x\nThZmZiX09/fT39/f6TA6yonEzMxKKTJEipmZNTE4ONjpEDrOLRIzMyvFicTMzEpxIjEzs1KcSMzM\nrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMysFCcSMzMrxd9sNzOrUH9/P2vWrGHBggWdDqUybpGYmVkp\nbpGYmVVocHCQVatWdTqMSrlFYmZmpTiRmJlZKU4kZmZWihOJmZmV4kRiZmalOJGYmVkpTiRmZlaK\nE4mZmZXiRGJmZqU4kZiZWSlOJGZmVkqliUTSYZIekLRc0vwG67eRdJ2keyXdLmnv3LozJS2VtEzS\nWbnlF0q6P21znaStq6yDmZm1VlkikTQNuBg4HNgTeL+kPeuKfQIYioh9gBOAi9K2ewOnAHOBfYF3\nS3pV2uZHwN5pm2Hg3KrqYGZmY6uyRTIXWB4RKyLieeBq4Ki6MnsCPwGIiPuBGZJ2BF4L3BYRIxGx\nFrgJODqVuzEtA7gV6K2wDmZmNoYqh5HfBXg4N/8IsF9dmXvIEsTPJc0FdidLDEuBz0jaDngWOAK4\ns8ExPgh8q9HBJZ0KnArQ29vbdcM4r169utMhtJ3rPDW4zpNPp59HMgBcJGkIWALcDbwQEfdJ+jxw\nI/AMMAS8kN9Q0ieBtcBVjXYcEZcClwLMmTMnpk+fXlklqtKNMZflOk8NrvPkUmUieRTYNTffm5aN\nioingJMAJAlYCaxI6y4DLkvrPkvWoiHN/w3wbmBeRERlNTAzG0N/fz+QPcBqqqryGskdwB6SZkra\nHDgWWJgvIGnrtA7gZODmlFyQtEP6vRtZ99c30vxhwMeAIyNipML4zcysgMpaJBGxVtIZwA3ANODy\niFgm6bS0/hKyi+pXSgpgGfCh3C6uTddI1gCnR8STaflXgb8AfpQ1Yrg1Ik6rqh5mZq1M5ZZITaXX\nSCJiEbCobtkluelbgL4m2x7UZPmrGi03M7PO8DfbzcysFCcSM7OK9Pf3j16Mn8ycSMzMrJROf4/E\nzGzSql2I77YvRG8ot0jMzKwUJxIzMyvFicTMzEpxIjEzq8BUuWMLnEjMzKwk37VlZlaBqTR0ilsk\nZmYlTaVurEacSMzMrBR3bZmZlTSVurEacYvEzMxKcSIxMxsHU/k6iROJmZmV4mskZmbjIH+dZKo9\nx90tEjOzceCuLTMzGxdTrTUCTiRmZuOiljiGh4c7HEn7OZGYmY2jvr4+gCnVzeVEYmY2TgYHB6dU\nl1aN79oyMxtnUy2ZuEViZmalOJGYmVkpTiRmZlaKE4mZmZXiRGJmZqU4kZiZWSlOJGZmVooTiZmZ\nleJEYmZmpSgiOh1D5SQ9ATzY6Tg20HRgVaeDaDPXeWpwnbvH7hGx/ViFpkQi6UaS7oyIOZ2Oo51c\n56nBdZ583LVlZmalOJGYmVkpTiQT16WdDqADXOepwXWeZHyNxMzMSnGLxMzMSnEimSAkbSvpR5J+\nm35v06LsNEl3S7q+nTGOtyJ1lrSrpJ9K+o2kZZLO7ESsZUk6TNIDkpZLmt9gvST9U1p/r6Q3dCLO\n8VSgzselui6R9CtJ+3YizvE0Vp1z5d4kaa2k/9HO+KriRDJxzAcWR8QewOI038yZwH1tiapaReq8\nFjg7IvYE9gdOl7RnG2MsTdI04GLgcGBP4P0N6nA4sEf6ORXo6kfsFazzSuCtEfE64Hy6/DpCwTrX\nyn0euLG9EVbHiWTiOAq4Mk1fCfy3RoUk9QLvAha0Ka4qjVnniHgsIn6dpv9ElkB3aVuE42MusDwi\nVkTE88DVZHXPOwr458jcCmwtaad2BzqOxqxzRPwqIv6YZm8Fetsc43gr8j4DfAS4Fni8ncFVyYlk\n4tgxIh5L078HdmxS7svAx4AX2xJVtYrWGQBJM4DXA7dVG9a42wV4ODf/COsnwyJlusmG1udDwA8q\njah6Y9ZZ0i7Ae+nyFme9TTsdwFQi6cfAf2mw6pP5mYgISevdTifp3cDjEXGXpLdVE+X4Klvn3H7+\nkuy/uLMi4qnxjdI6SdLBZInkLZ2OpQ2+DHw8Il6U1OlYxo0TSRtFxNubrZP0H5J2iojHUpdGo2bv\ngcCRko4AtgBeIelfIuL4ikIubRzqjKTNyJLIVRHx3YpCrdKjwK65+d60bEPLdJNC9ZG0D1k37eER\n8Yc2xVaVInWeA1ydksh04AhJayPie+0JsRru2po4FgInpukTge/XF4iIcyOiNyJmAMcCP5nISaSA\nMeus7C/uMuC+iPhiG2MbT3cAe0iaKWlzsvduYV2ZhcAJ6e6t/YHVuW6/bjRmnSXtBnwX+EBEDHcg\nxvE2Zp0jYmZEzEh/w98BPtztSQScSCaSAeBQSb8F3p7mkbSzpEUdjaw6Rep8IPAB4BBJQ+nniM6E\nu3EiYi1wBnAD2c0C10TEMkmnSTotFVsErACWA/8H+HBHgh0nBet8HrAd8L/T+3pnh8IdFwXrPCn5\nm+1mZlaKWyRmZlaKE4mZmZXiRGJmZqU4kZiZWSlOJGZmVooTiU1YkrbL3fL7e0mP5uZ/1eZY/i6N\nQHyvpMWSdq9b/wpJj0j6am7Z1yWtzMU8u26b9UaAlfS3aZTjpZK+KWmLtPz8dOwhSTdK2jktn5vb\n/z2S3tuiDt+RNKvB8r/Jx12WpNdJ+vp47c8mPicSm7Ai4g8RMTsiZgOXAF+qzUfEm9sczt3AnIjY\nh+yLZBfUrT8fuLnBdn+fi3motrDRCLBpHKaPpuPsDUwj+1IbwIURsU96La4n+w4GwNJUfjZwGPA1\nSeuNWCFpL2BaRKzY0IpvqIhYAvSmLxzaFOBEYl1J0tPp99sk3STp+5JWSBpQ9pyL25U95+KVqdz2\nkq6VdEf6OXBDjhcRP42IkTS7zki1kt5INuDkhgwL3mwE2E2BLVMy6AF+l46fH1/sZUCk5SPpi3CQ\nDZvT7Ithx5EbOUDSSZKGJd1O9qXP2vL3SLpN2fNufixpR0mbKHtmzPapzCbKnrexvaT3pdbTPZLy\nifRfeSkJ2iTnRGKTwb7AacBryb4F3xcRc8nGcPpIKnMRWYvmTcB/p9ww/KMj1UraBPgCcE6Tsp9L\nXVJfkvQXaZuGI8BGxKPAPwIPAY+RDZOSb7F8RtLDZEnhvNzy/SQtA5YAp+USS96BwF2p/E7Ap9Oy\nt5A9O6PmF8D+EfF6smHQPxYRLwL/ko4L2SgE90TEEymOd0bEvsCRuf3cCRzU5DWxScaJxCaDO9Jz\nS54D/o2XWgZLgBlp+u3AVyUNkY1/9AplIwpvEEnHkw28d2Fa9GFgUUQ80qD4uUAf8CZgW+Djafno\nCLB1+96G7PkVM4GdgZel4wEQEZ+MiF2Bq8iG4qgtvy0i9krHObd2XaXOTsATaXo/4GcR8UR6bsa3\ncuV6gRskLQH+HtgrLb8cOCFNfxC4Ik3/Evi6pFPIuuJqHk91sCnAo//aZPBcbvrF3PyLvPQZ34Ts\nP+0/N9uJpCvInnfyu4hYbzwvSW8nG/7+rSlpARwAHCTpw8BfAptLejoi5ucGXXwu7bvWamk4Aiyw\nGbAy/aePpO8CbyZrDeRdRTY21z/kF0bEfanLb2+yFkHes2RdX2P5CvDFiFio7FEFn0r7fljZaM2H\nkD3A6bi0/DRJ+5E9bO0uSW9Mo/hukY5pU4BbJDZV3MhL3VzU30EFEBEnpYvijZLI64GvAUdGxOO5\nbY6LiN3SaK7nkD3lcH7aZqf0W2RPf1yatmk2AuxDwP6SetI280iPVJa0Ry6co4D70/KZtYvr6U6y\n1wD/3qD+9wGvStO3AW9VdlfcZsD7cuW24qWhz09kXQvIktq3I+KFdMxXphbReWQtntow6n21+trk\n5xaJTRXcNqu0AAAA8klEQVQfBS6WdC/Z5/5msusqRV1I1uL4dmpJPBQRR7behKvSBWoBQ2MdLyJu\nk/Qd4Ndkz6q/m5eeYz4g6dVkrawHc/t6CzBf0pq07sMRsarB7v8f8Dbgx+n5L58CbgGeTLHVfCrV\n8Y/AT8i62WoWknVpXZFbdmFKcgIWA/ek5QenY9oU4NF/zaYASVsCPwUOrLUmNmIfc8huWGh5ET3d\nVHAT8JYmF/5tknHXltkUEBHPkl1T2ajnwEuaT3a78rkFiu8GzHcSmTrcIjEzs1LcIjEzs1KcSMzM\nrBQnEjMzK8WJxMzMSnEiMTOzUpxIzMyslP8PnBK2S4SKLBgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.remove_nans().flatten(window_length=401).fold(period=3.5225).bin(binsize=10).plot();" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.6" } }, "nbformat": 4, "nbformat_minor": 2 }