d007-Circle_and_Cycloid.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import cv2\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x1133a40f0>"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQsAAAD8CAYAAABgtYFHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAElZJREFUeJzt3X2sHNV9xvHvUwMmbWiMiWu5AQmSuI1I1TiYglHyBwXR\nGLeqiRQho6pYyJJTCSSiRG1wK7WJ1EqJ1MQtUovqFIqp0gAlibAQLXUMUpU/eLvgEGOHYApRsAzO\nCziJkNza/PrHnrWHy32Z3ZnZOTP7fK5Wd3dm9u4Z353H58ycuz9FBGZmi/mlthtgZt3gsDCzUhwW\nZlaKw8LMSnFYmFkpDgszK6WxsJC0XtJzkg5KuqWp1zGzyVAT8ywkLQG+D1wFvAw8AVwXEftrfzEz\nm4imehaXAAcj4n8i4n+Bu4GNDb2WmU3AaQ393PcAPyw8fhm4dL6NJXkaaQvWspYZZljL2pGfO8NM\nAy2yhv04IlaM++SmwmJRkrYCW9t6/WkSzJ/FQjzJkyfvW6/9oMqTmwqLQ8B5hcfnpmUnRcQOYAe4\nZ1GXhUIB5g8Dh4SV0dQ5iyeA1ZIukHQGsAnY1dBrTa2Y9VWkOb7MqmikZxERxyXdBDwELAHuiIhn\nm3itaTNX78FBYJPQyKXTkRvhYciCHBBWk5mIuHjcJ7d2gtMWN9fQwqwtDosMOSQsRw6LTDggLHcO\ni5YVQ8IBYTlzWLTAAWFd5D9RnzAHhXWVexYT4pCwrnNYNMwhYX3hsGiIQ8L6xmFRM4eE9ZVPcNbI\nQWF95p5FDRwSNg3cs6jIQWHTwj2LCoZB4ZCwaeCwGIN7EzaNHBYjcm/CppXDoiT3JmzaOSxKcG/C\nrGJYSHoJ+DlwAjgeERdLWg7cA5wPvARcGxGvVWtmexwUZgN1XDr93YhYU/hsv1uAPRGxGtiTHneS\ng8LslCbmWWwEdqb7O4FrGniNxjkozN6qalgE8F+SZlKFMYCVEXE43X8FWFnxNSaqWIPDQWF2StUT\nnB+NiEOSfg3YLel7xZUREfN9zH+O5QsdEmbzq9SziIhD6fsR4JsMqqe/KmkVQPp+ZJ7n7oiIi6vU\nMaiTg8JsYWOHhaRfkXTW8D7we8A+BmUKN6fNNgP3V22kmbWvyjBkJfBNScOf828R8Z+SngDulbSF\nQdXma6s3s1lBuEdhtoipL1/ooLAp4vKF4/A5CrPRTOXnWcxVaNjMFjZ1PQv3KMzGM5U9CweF2eim\nKiw8/DAb39SEhYcfZtVMRVg4KMyq631YOCjM6tHrsHBQmNWn12EBDgqzuvQ2LHzlw6xevQwLDz/M\n6tfLsAAHhVndehcWHn6YNaNXYeHhh1lzehMWDgqzZvUmLMysWb0IC/cqzJrX+bBwUJhNxqJhIekO\nSUck7SssWy5pt6Tn0/ez03JJulXSQUnPSLqoycab2eSU6VncCayftWy+eqZXA6vTbStwWz3NnJt7\nFWaTs2hYRMR/Az+dtXi+eqYbgbti4FFg2bDgUFMcFGaTMe45i/nqmb4H+GFhu5fTsreRtFXSk5Ke\nHKcBnnxlNlmVP7B3oXqmizxvB7ADxq8b4l6F2eSM27OYr57pIeC8wnbnpmW1cq/CbPLGDYv56pnu\nAq5PV0XWAUcLw5VauVdhNlmLDkMkfQ24HHi3pJeBvwK+wNz1TB8ENgAHgTeAG+pusHsVZu3oVK1T\nXyo1q6RSrdPOz+A0s8noXFi4V2HWjs6ERRAOCrMWdSYszKxdnQgLXwExa18nwgJ8rsKsbZ0JCzNr\nV/Zh4SGIWR6yDwvwEMQsB1mHhXsVZvnIOizMLB/ZhoX/DsQsL9mGhZnlJeuwcK/CLB9ZhoVPbJrl\nJ8uwMLP8ZBsWHoKY5SW7sPAQxCxP45Yv/JykQ5L2ptuGwrptqXzhc5I+1lTDzWyyxi1fCLA9Itak\n24MAki4ENgEfTM/5R0lLRm2UhyBm+Rm3fOF8NgJ3R8SxiHiRwad8X1KhfWaWiSrnLG5KldLvGFZR\np2L5Qp+vMMvXuGFxG/A+YA1wGPjSqD8gInZExMVVPprczCZnrLCIiFcj4kREvAl8hVNDjYmULzSz\nyRsrLIZ1TpOPA8MrJbuATZKWSroAWA08PtLP9slNsyyNW77wcklrgABeAj4JEBHPSroX2A8cB26M\niBPNNN3MJimr8oWuDWLWqH6UL/SVELO8ZRMWZpa3rMLCQxCzfGUVFmaWL4eFmZXisDCzUrIIi7Ws\nbbsJZraILMJihhmf3DTLXBZhsZa1nmdhlrkswsLM8pdNWHgYYpa3bMLCzPLmsDCzUhwWZlaKw8LM\nSnFYmFkpDgszK8VhYWallClfeJ6kRyTtl/SspJvT8uWSdkt6Pn0/Oy2XpFtTCcNnJF3U9E6YWfPK\n9CyOA5+JiAuBdcCNqUzhLcCeiFgN7EmPAa5m8Kneq4GtDGqMmFnHlSlfeDginkr3fw4cYFBlbCOw\nM222E7gm3d8I3BUDjwLLZpUOMLMOGumchaTzgQ8DjwErI+JwWvUKsDLdL1XCsFi+8Ef8aMRmm9ko\n6vhDzdJhIemdwNeBT0XEz97SkEE9gZFaUyxfuIIVozzVzFpQKiwknc4gKL4aEd9Ii18dDi/S9yNp\nuUsYmvVQmashAm4HDkTElwurdgGb0/3NwP2F5denqyLrgKOF4YqZddSi5QuBjwB/DHxX0t607M+B\nLwD3StoC/AC4Nq17ENgAHATeAG6otcVm1opFwyIivg3zftjElXNsH8CNFdtlZpnxDE4zK8VhYWal\nOCzMrBSHhZmV4rAws1KyCQvXDTFrRl3HVjZhYWZ5KzMpq3EuX2jWnLqOrSx6Fi5faNacSF9VZREW\nM8y03QSzXqujd5FFWJhZ/hwWZlaKw8LMSskqLHyS06xedR5TWYWFmeUrm7DwPAuzZvRqnoWZ5c9h\nYWalVClf+DlJhyTtTbcNhedsS+ULn5P0sVEa5JOcZvWo+1gq87chw/KFT0k6C5iRtDut2x4Rf1vc\nOJU23AR8EPh14FuSfiMiTtTZcDObrCrlC+ezEbg7Io5FxIsMPuX7kjoaa2btqVK+EOCmVCn9jmEV\ndUqWL5zz5/uKiFmt6jymqpQvvA14H7AGOAx8aZQXLtY6HeV5ZtaOscsXRsSrEXEiIt4EvsKpoUap\n8oXFWqdvW+eTnGaVNHEMjV2+cFjnNPk4sC/d3wVskrRU0gXAauDx+ppsZm2oUr7wOklrGFRPfwn4\nJEBEPCvpXmA/gyspN45yJUTIPQuzGtR9DlCDaoPtkvSWRgzDwic8zUa3wPEzM9ewvyzP4DSzUrIM\nC/cozKpp4hjKMiyGfO7CbDRNHjNZh4WZ5SPbsBh2o9y7MCun6QsD2YaFmeUl67DwiU6z0TR5zGQd\nFkMeipgtbBLHSPZh4d6FWTlNHyvZh4WZ5aEzYeGhiNncJnVsdCIsPBQxW9gkjpFOhIWZta8zYeE/\nXTd7uyAm1vPuTFgMOTDMBiZ9LHQuLMysHZ0KC/+9iNlAGx8Q1amwAF8ZMRua9LFQ5gN7z5T0uKTv\npPKFn0/LL5D0WCpTeI+kM9LypenxwbT+/CYa7t6FTau23vtlehbHgCsi4kMMaoSsl7QO+CKD8oXv\nB14DtqTttwCvpeXb03a1cu/Cpl0bx0CZ8oUREb9ID09PtwCuAO5Ly3cC16T7G9Nj0vorUzmB2rl3\nYdOmzfd82SJDS1IZgCPAbuAF4PWIOJ42KZYoPFm+MK0/CpxTZ6PBJztt+rT9qfelwiJVHlvDoLrY\nJcAHqr6wyxeadctIV0Mi4nXgEeAyYJmkYZGiYonCk+UL0/p3AT+Z42fNW76wLPcubFq03auAcldD\nVkhalu6/A7gKOMAgND6RNtsM3J/u70qPSesfjhwqGZlZJWXKF64CdkpawiBc7o2IByTtB+6W9NfA\n0wzqoZK+/6ukg8BPgU0NtPuk4d+MTHKOvNkk5dCrgEzLF44jl39QszrV/L52+UIza15vwsInO61v\ncust9yYswIFh/ZFbUEDPwgLy+sc1qyK393LvwmLIvQvrqlzfu70MCw9HrKtyHH4M9TIsIM9/bLMy\ncn3v9jYshty7sK7I/b3a67DwcMS6Iufhx1CvwwIcGJa/LgQFTEFYgAPD8tWVoIApCQtwYFh+uhQU\nMEVhAd35pdj06NJ7cqrCYsi9C2tbF9+DUxcWxeFIF39h1m3F912XehUwhWEB3fslWf908T1Y5pOy\neqnYw+jiL866p+vvtansWRQNP5bPrEldDwqoVr7wTkkvStqbbmvSckm6NZUvfEbSRU3vRFXFz/E0\nq1OfPh+2zDBkWL7wF5JOB74t6T/Suj+NiPtmbX81sDrdLgVuS9/NrMOqlC+cz0bgrvS8RxnUF1lV\nvanN8qQtq1tXr3rMZ6zyhRHxWFr1N2mosV3S0rTsZPnCpFjaMGu+rGp16PLl0YWMVb5Q0m8B2xiU\nMfwdYDnw2VFeONfyhUpf4F6Gja4YEn0KChi/fOH6iDichhrHgH9hUAMVCuULk2Jpw+LPqly+sEkO\nDBtVH3sTReOWL/ze8DyEJAHXAPvSU3YB16erIuuAoxFxuJHWN8yBYWX1PSigWvnChyWtAATsBf4k\nbf8gsAE4CLwB3FB/syfH5RFtMdMQFNCj8oVNK/Yu+v6msHI6+J5w+cJJ8IlPK+rzicz5OCxG5Mur\n062vl0XLcFiMwb2M6TSNvYkih0UF7mVMh2nuTRQ5LCoqvnkcGP3TwZOYjZnaz7Oo01w9jGl/Y3Wd\nQ+Lt3LOokXsZ3Td7SOmgOMU9i5q5l9FdDomFOSwa4tDoDodEOQ6Lhjk08uWQGI3DYkIcGvlwSIzH\nJzgnzCdB2+WgGJ97Fi2YLzD85m2G/43r4bBo2XyzQP2mrmZ2r83/ntU5LDIxu7fh4BidA6JZDosM\nzf4jNZ8QXZhDYjIcFhmbLzSK66bVXCeHp/3fpGkOiw6Y64ToNAaHA6JdDouOWSw45tquixa6rNz1\nfeuq0vMsUqGhpyU9kB5fIOmxVNP0HklnpOVL0+ODaf35zTTdNOurKOb4ytVibV1oP21yRpmUdTNw\noPD4i8D2iHg/8BqwJS3fAryWlm9P29kEzD6oZh9c8x2UTYfJXK81VzAs1n5rV9nyhecCvw/8c3os\n4ApgWBR5J4PaITCodboz3b8PuDJtby0pexDOdSDX8TVueywvZc9Z/B3wZ8BZ6fE5wOsRcTw9LtYz\nPVnrNCKOSzqatv9x8QdK2gpsTQ+PcapIUd+8m1n7novhATrmgTrWfnUkFLL9nVX0m1WevGhYSPoD\n4EhEzEi6vMqLFUXEDmBHeo0ncy1jWFVf962v+wX93beqdYXL9Cw+AvyhpA3AmcCvAn8PLJN0Wupd\nFOuZDmudvizpNOBdwE+qNNLM2rfoOYuI2BYR50bE+cAm4OGI+CMGBZI/kTbbDNyf7u9Kj0nrH44c\nyp6ZWSVV/kT9s8CnJR1kcE7i9rT8duCctPzTwC0lftaOCu3IXV/3ra/7Bf3dt0r7lUWtUzPLnz/8\nxsxKaT0sJK2X9Fya8VlmyJIVSXdIOiJpX2HZckm7JT2fvp+dlkvSrWlfn5F0UXstX5ik8yQ9Imm/\npGcl3ZyWd3rfJJ0p6XFJ30n79fm0vBczkpucad1qWEhaAvwDcDVwIXCdpAvbbNMY7gTWz1p2C7An\nIlYDezh13uZqYHW6bQVum1Abx3Ec+ExEXAisA25Mv5uu79sx4IqI+BCwBlgvaR39mZHc3EzriGjt\nBlwGPFR4vA3Y1mabxtyP84F9hcfPAavS/VXAc+n+PwHXzbVd7jcGV7uu6tO+Ab8MPAVcymAS1mlp\n+cn3JfAQcFm6f1raTm23fZ79OZdBgF8BPACozv1qexhycrZnUpwJ2mUrI+Jwuv8KsDLd7+T+pi7q\nh4HH6MG+pa76XuAIsBt4gZIzkoHhjOQcDWdav5kel55pTYn9ajssei8G0d3ZS06S3gl8HfhURPys\nuK6r+xYRJyJiDYP/iS8BPtBykyorzrRu6jXaDovhbM+h4kzQLntV0iqA9P1IWt6p/ZV0OoOg+GpE\nfCMt7sW+AUTE6wwmF15GmpGcVs01I5nMZyQPZ1q/BNzNYChycqZ12qbSfrUdFk8Aq9MZ2zMYzBDd\n1XKb6lCcxTp7duv16crBOuBooUuflfSXwrcDByLiy4VVnd43SSskLUv338HgPMwBOj4jOSYx0zqD\nkzIbgO8zGDf+RdvtGaP9XwMOA//HYEy4hcHYbw/wPPAtYHnaVgyu/rwAfBe4uO32L7BfH2UwxHgG\n2JtuG7q+b8BvA0+n/doH/GVa/l7gceAg8O/A0rT8zPT4YFr/3rb3ocQ+Xg48UPd+eQanmZXS9jDE\nzDrCYWFmpTgszKwUh4WZleKwMLNSHBZmVorDwsxKcViYWSn/DwhFsJWfch4TAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x112de9940>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "img = np.zeros((400, 400, 3)).astype('uint8')\n",
    "\n",
    "def draw_point(img, pt, size, color):\n",
    "    cv2.circle(img, pt, size, color, thickness=-1)\n",
    "    \n",
    "n = 1000\n",
    "r = 200\n",
    "x0, y0 = 200, 200\n",
    "for i in range(n //4*3):\n",
    "    theta = i * 2 * np.pi / n\n",
    "    pt = (x0 + r * np.cos(theta), y0 + r * np.sin(theta))\n",
    "    pt = (int(pt[0]), int(pt[1]))\n",
    "    draw_point(img, pt, 1, (255, 0, 255))\n",
    "\n",
    "plt.imshow(img)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x114972940>"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQsAAAD8CAYAAABgtYFHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF29JREFUeJztnW0MJVV5x39PebXViijdUKEBdVuDTV0JxSX6gWKsy7bp\nYmLMkqYQQ7I2wUSjaQWbtJq0SU2qtKQtKRbq2liR+hI2xNaui0nbDwKiK/IislaMkJWNCqgxoV14\n+uGeuztc79w7d17Py/+X3Nx5u3fOOXPmP885c2b+5u4IIcQ6fm7qBAgh0kBiIYRohMRCCNEIiYUQ\nohESCyFEIyQWQohGDCYWZrbDzB4ys0Nmds1Q+xFCjIMNMc7CzE4Avgm8EXgUuBu43N0f6H1nQohR\nGCqyuBA45O7/4+7/C9wC7BpoX0KIEThxoP99KfDdyvyjwGvrNjYzDSMVYni+7+5ntP3xUGKxFjPb\nA+yZav9CFMh3uvx4KLF4DDi7Mn9WWHYMd78RuBEUWQiRAkP1WdwNbDWzc83sZGA3sG+gfYkavONH\niCqDRBbuftTM3gF8HjgBuNnd7x9iXyUz9Am97v9t4P2LuBjk1unGiVAzpBFTn7xT71905h53v6Dt\njyfr4BTNWXWSjnmCLu5rMV1es53IA4lFpCwTiNhOwjrxWEx7bOkW7ZBYREIK4rCOJuKRWp7EcSQW\nE5ODSNRRzceicOSSx5KQWExAiWH6onAo2kgPicWIlCgSy5jnW9FGWuh9FiOxeCXVibH+7oqIC0UW\nA6NwezXLogyVU5xILAZCTY7NqIqGBDZOJBY9o4reDUUa8aI+ix6RUPRHtV9HD7bFgSKLHpBIDMei\nYKh8p0ORRUckFOOiCGM6JBYdkFCMx2KzRIyPxKIl1YFEEorxUD/GdEgsNqRaSSUS06AoYxokFhug\naCIuJBjjIrFoiKKJOJFgjEenW6dm9gjwY+AZ4Ki7X2BmpwOfBM4BHgHe6u5PdEvmtKgipoFurQ5L\nH5HFb7n7tsq7/a4BDrj7VuBAmE8WNT3iR30Y4zBEM2QXsDdM7wUuG2Afo6CmR1pIMIalq1g48B9m\ndk9wGAPY4u6Hw/T3gC0d9zEJEoo0kWAMR9fh3q9398fM7JeA/Wb2jepKd/e61/zHal+ogVbpY2h4\n+BB0iizc/bHwfQT4LDP39MfN7EyA8H2k5rc3uvsFXXwMhkSVLG0UYfRPa7Ews18wsxfMp4HfBu5j\nZlN4ZdjsSuC2rokcCzU9hKinSzNkC/BZM5v/z7+4+7+b2d3ArWZ2FTPX5rd2T+bwSCjyQ82RfpF9\nIRKKWBm6UhR4vGVf2AcFVpzoGNJLddl/69WHm1G8WChEnY46cRjieMybJMv+W7aLzShWLKZvfJXJ\n1A5sywRDtovNKP5BMlWC8Vh25Z7SBX7VdssEpPQLTJGRRekHfUymjiSW0bTpKa/W51KkWECZB3tM\nYm3/t72dWufVGku+xqA4sVBUMTyptPXbdm6XaohUlFiUeDUYi1gjiWXMo4s+/mdOCdFGcR2cuR7I\nKUnR9LnvZ0dKeKdGMWKR6wGcmpLC8E3I8e5JMWIh+if1t4gNNfhr2V2UHChCLHJvS45Nbu3zIU7o\nHL1asxeLnCp1DOTW7Bg6Dzn1ZWQvFqI/Um92rGKsJ1xTjjKyFgtFFf2hsuxO6n0ZWYuF6JdchWLs\n51NSbZZkKxa6EvZHKY/xj3nypigY2YqF6I5eSTcsqQlG1mKhSt6eVCpwX0xVV1Lq+FwrFmZ2s5kd\nMbP7KstON7P9ZvZw+H5RWG5mdr2ZHTKze83s/CETX0fshR47Od/1WMcUdSeVjs8mkcVHgR0Ly+r8\nTC8FtobPHuCGfpK5OaVVcpE2KXR8rhULd/9P4IcLi+v8THcBH/MZXwJOmxsOjUWsBZ0KJXcMx3Cy\nxpCGOtr2WdT5mb4U+G5lu0fDsp/BzPaY2ZfN7Mst0yAGokShEOvp3MHpM+ORjYVwSPtCVfZ2xHg1\nK5FYo4u2YlHnZ/oYcHZlu7PCslGIrXBTouTmR4zEKBhtxaLOz3QfcEW4K7IdeKrSXBGRElOFFMeJ\nTTDWvlbPzD4BXAy8xMweBf4M+EuW+5l+DtgJHAJ+CrxtgDQvRVfG7qjs4iMmv9ZsvE4lFu3I7ZHz\nPoixLvWUpk5ep1mP4BTNienEEHEisRASigSIof8iC7GIMWxMgekboGITphaMLMRCtEcCK5qSjVio\n0ou+iLkuTRldZCMWYjPUBEmXqQRDYiGEaETyYqErZHtiDrfFaqaILpIXCyFKZWzBkFgUiKIx0QaJ\nhRALpCSmY0YXWYiF2t5CDM/ap07FODS9MvQljBLYfJg/mTo0EouR6XpQ636vk18M/Ri7xGJgVolD\n1wPrNdN9/LdIizGiC4nFACw7aEOcvMv+0xe+h9p3rqT+UOKQ0UXSYhFbr3UMV/dlhjWpnwCiGUO/\nVStpsYiBGASijsXbarGJq0iLtvaF7zezx8zsYPjsrKy7NtgXPmRmbxoq4TEQs1BUiTVdIi3a2hcC\nXOfu28LncwBmdh6wG3hV+M3fm9kJfSU2Jhb7BGI/IY3nplNRRp4MeXzb2hfWsQu4xd2fdvdvM3vL\n94Ud0hcdVUelFERiFSk4d4+JymI1XUZwviM4pd88d1EnY/vCnEQCFGWsIodjOwRtxeIG4OXANuAw\n8KFN/2BI+8K+yfk2ZFUwJBp50ffxbCUW7v64uz/j7s8CH+F4U2NS+8K+yS2aqENRxnFyOcZD5KOV\nWMx9TgNvBuZ3SvYBu83sFDM7F9gK3NUtidOTSwWqUrooLKLyWE9b+8KLzWwbszJ+BHg7gLvfb2a3\nAg8AR4Gr3f2ZYZI+7BDXEgcyxWSVJ/qhz2OZvH1h3xU75/6JOeuEsEShzFEglxxH2Rf2RQn9E7A+\nb6X1X+Saz77rsMQiUOrVtI7SBKOk494WiQVlCsUm5CwYOedtTl95zEIsuhSGhKIelYmokoVYdKXE\nk6JpnnMuG10oNqNYsdCIxRlNy0BllSZ9CmGRYlHC7dEmbBpd5CQYiio2J3mxaHuwc789ugk5iUAT\nShSKPo5x8mKxKaWdGOtoGjXkcmLp+LcnG7HYpBLkUvGnIJeTTXVgc7IRiybkUtH7Jsc+iWWU2Pzo\nkyzEokllV0VZTZMyTLnsdPy7k4VYNEUVZTVNyye1CKR0oegr39mLhcZTbE6O5VWqUPRJ9mIxR5Wl\nGTn1X+jdHP2SjVjkVMmnZlVZplLOsacvRbJ2JCu9rdqFvt6a1eak7cswWse9X7ISi2Wv2VOFaU8b\nwWgiDov/tfibLsdQQjEcTewLzzazL5rZA2Z2v5m9Myw/3cz2m9nD4ftFYbmZ2fXBwvBeMzt/6Ews\nok7N/mhiFeBL1tuKz7J9rNrOWb6PZWmopln0S5M+i6PAe9z9PGA7cHWwKbwGOODuW4EDYR7gUmZv\n9d4K7GHmMTIaVjMt2lNnFbBKIPrcb5141KVBDEMT+8LD7v6VMP1j4EFmLmO7gL1hs73AZWF6F/Ax\nn/El4LQF6wCRKOuigjH2vygadetE/2x0N8TMzgFeA9wJbHH3w2HV94AtYbqRhWFq9oVClE5jsTCz\n5wOfBt7l7j+qrvOZn8BG3QQp2ReKGcsO8Lq+hL73X9fkUD/V8DQSCzM7iZlQfNzdPxMWPz5vXoTv\nI2H5pBaGi21q0Z26zsNVHZF97reub6QuDWIYmtwNMeAm4EF3/3Bl1T7gyjB9JXBbZfkV4a7IduCp\nSnNlFNR+7Y8mXipN7mCsEpJ12627m7KYPgnGMKx1JDOz1wP/BXwdeDYsfh+zfotbgV8BvgO81d1/\nGMTlb4EdwE+Bt7n7yn6JLo5kVRavfrqV1o268tukXDUoKw7CWJlOjmTJ2xdWWVZJVHHasarcUihT\n3U59Ln2IRTbPhqRQgVMhdaGA+NOXItmIxTqmj5/SIBUxaMK8H0PHvh+yFwt1dm5OjuUlwehO9mJR\nRRVmNU3LJzUxKf0uSV/5zkIsmoTOpVeYdTQpw5TLTse/O1mIRVNSuyKORU79FKuQYHQjG7Fo884D\nsTm5CIrqwOZkIxZNyaWy90XTqCKXk0vHvz3Ji0XbSqwHj45T2glUYnOkj2OcvFi0QQ8ezWia9xz7\nNEoUjK4UKRag8RdzmpaByipN+hTDYsWiSolXl02jihxRdLEZWYhFl6ueKkw9KhNRJQux6IoEYzU5\nN0FyztucvvIosQiUKBhNRmuWcDJBWce9LRKLCk08MnJgXd5KE4pc89l3HZZYLFD669lKE4oqJR7v\nTUhaLIY8uCUKRslCkSt9Hssu9oXvN7PHzOxg+Oys/ObaYF/4kJm9qcf0TkKOgiFBeC4qj/U0MUae\n2xd+xcxeANxjZvvDuuvc/a+qGwdrw93Aq4BfBr5gZr/q7s/0mfAxWNaHkWOl0vsqjxPeVZk8Q1zg\nutgX1rELuMXdn3b3bwOHgAv7SOxU5Dw8vMmr/kWa9H08u9gXArwjOKXfPHdRp6F9YWosdnymLhq5\nR0tdyOHYDkEX+8IbgJcD24DDwIc22XGqXqc5iYaiieeislhNa/tCd3/c3Z9x92eBj3C8qdHIvjB1\nr9PUfDYXnb50YuTJkMe3tX3h3Oc08GbgvjC9D9htZqeY2bnAVuCu/pIcD3VOXbERa7pEWjS5G/I6\n4A+Ar5vZwbDsfcDlZraNWV18BHg7gLvfb2a3Ag8wu5NydYp3QpqyKsKY+uq9KBIljh0R/ZG8fWGM\nt7qWZWasNDbZd4xlFgspN9MapL2TfWGTyEJsyLJbrUMIyCqFTbGyx8DcwSxVQR0yzRKLgVnVr9Fn\nTJdixRb9MUb7QGIxMnUnda5uYGI8hq4bEotIGFsEUg2zxc8yVq9j0k+dzpm+i1aI/MlCLITok5Qi\nrjHv3kgsCiSlk0HEg8RCiEQZe0xI8mKhq2R71NeTLlMMHkteLIQQ4yCxKBRFZOky1ZD0bMRCIbXo\ni5jr0pTPrmQjFqIdMZ8YIi6yEAs9et0ONUXSYuonYrMQC9ENiWz8TC0UILEQAQmGWIfEonDUFHku\nMVzBF4klTdmIhfotuqOyi49YhAKavbD3VDO7y8y+FuwLPxCWn2tmdwabwk+a2clh+Slh/lBYf86w\nWRBdiaEiip8lJqGAZpHF08Al7v5qZh4hO8xsO/BBZvaFrwCeAK4K218FPBGWXxe2G4VYCjVFFJnF\nRWxCAc3sC93dfxJmTwofBy4BPhWW7wUuC9O7wjxh/RuCncBoqMK3I6aKWTIxCgU0Nxk6IdgAHAH2\nA98CnnT3o2GTqkXhMfvCsP4p4MV9JnplWsfaUcaUKrYx5DtWoYCGYhGcx7Yxcxe7EHhl1x2nal+Y\nO2qOxHmixsBGd0Pc/Ungi8BFwGlmNn+HZ9Wi8Jh9YVj/QuAHS/5rUPvCkiu7SI8UrCWb3A05w8xO\nC9PPA94IPMhMNN4SNrsSuC1M7wvzhPV3+MhORrEWdirkYvzchinqTkwudqto8nbvM4G9ZnYCM3G5\n1d1vN7MHgFvM7M+BrzLzQyV8/7OZHQJ+COweIN2N0Bus2zM32ymFqfIaezRRJXn7wjpSOgixU4Lo\nTlFfJthnJ/vCbEZwLqKOuv4oJcrIXCg6k61YiP7JVTDGzFcKHZl1ZC0Wii76Q2XZnVQ6MuvIWixE\nv+R8l2Tok7caTaQoFFCAWOiK2C/Vip5DmQ6dh5SbHYtkLxYgweib6tUxhzId4iReFInUhQIKEQsx\nDKk3S4ZI82JZ5CASc4oRi5wOWkzk1izpi1yiiSrFiMUcVej+WRSMFMq4736EnPom6mgy3Dsb5oOL\nShiRODZ1ghFjOfclZov/E2Ne+6QosYByRiNOSbWMYxaNtmnKtU9iHcWJxZyYK3EOLN4tieUEa3vc\nS4sillGkWCi6GI9lHaBTC3XT/S6rIyWKxJziOjgXkWiMx+KJNnZnaNN9LUtXjnc3NqXIyALyGlSU\nEqsijWXbDLn/xXQ02bZkihWLOfMmiSrG+CyLNKrf67bfhHX/3fX/S6B4sZgjwZieOvGom+97f2I1\nEgs0/iJW+jgWU3em5kQX+8KPmtm3zexg+GwLy83Mrg/2hfea2flDZ6IP1IeRHxKKfmkSWcztC39i\nZicB/21m/xbW/ZG7f2ph+0uBreHzWuCG8B09ijCEqKeLfWEdu4CPhd99iZm/yJndkzouijDSRlFF\n/7SyL3T3O8OqvwhNjevM7JSw7Jh9YaBqbRg9ub2roUQkFMPQyr7QzH4duJaZjeFvAqcD791kx7Hb\nF0ow0kRCMRxt7Qt3uPvh0NR4GvgnZh6oULEvDFStDav/Nah9YR9IMNJCQjEsbe0LvzHvhzAzAy4D\n7gs/2QdcEe6KbAeecvfDg6R+BFJ/G1QJxP5IfC50sS+8w8zOYHZ8DgJ/GLb/HLATOAT8FHhb/8ke\nFz14lgYSimHJ1r5wCHT1ig8dk42QfeFYqEkSFxKKcZFYbIhurU6P+iimQWLREkUZ05CbF0dKSCw6\noNfgj4eiiemRWHREgjEuEorp0CPqPbCsSaJK3Q+xvOhXKLLoFUUZ/ZGjV2jqKLLombqOT1X2ZuiV\n+/EisRiIZb4Zqvj1SFjjR2IxMOrPWI0iiXRQn8VIpGgePDQSirRQZDEidZZ+pZ0kEok0kVhMwLoo\nI8eTp4Q85o7EYmKWPWeSS8Qhr9C8kFhEQhMf0NhPNNkA5o3EIlKaeIIubjc2EoeykFgkwKpH4scc\nn7DuDo5EIm8kFgkxthfoOiQOZSGxSJh1J2tX8ZAYiCqNB2UFo6GvmtntYf5cM7szeJp+0sxODstP\nCfOHwvpzhkm6WId1/AhRZZMRnO8EHqzMfxC4zt1fATwBXBWWXwU8EZZfF7YTQiROU/vCs4DfAf4x\nzBtwCTA3Rd7LzDsEZl6ne8P0p4A3hO2FEAnTtM/ir4E/Bl4Q5l8MPOnuR8N81c/0mNepux81s6fC\n9t+v/qGZ7QH2hNmnOW5SlBsvYSHvmZBrviDfvP1alx+vFQsz+13giLvfY2YXd9lZFXe/Ebgx7OPL\nMdsYdiHXvOWaL8g3b119hZtEFq8Dfs/MdgKnAr8I/A1wmpmdGKKLqp/p3Ov0UTM7EXgh8IMuiRRC\nTM/aPgt3v9bdz3L3c4DdwB3u/vvMDJLfEja7ErgtTO8L84T1d3gMtmdCiE50eZ/Fe4F3m9khZn0S\nN4XlNwEvDsvfDVzT4L9u7JCO2Mk1b7nmC/LNW6d8ReF1KoSIH70pSwjRiMnFwsx2mNlDYcRnkyZL\nVJjZzWZ2xMzuqyw73cz2m9nD4ftFYbmZ2fUhr/ea2fnTpXw1Zna2mX3RzB4ws/vN7J1hedJ5M7NT\nzewuM/tayNcHwvIsRiQPOdJ6UrEwsxOAvwMuBc4DLjez86ZMUws+CuxYWHYNcMDdtwIHON5vcymw\nNXz2ADeMlMY2HAXe4+7nAduBq8OxST1vTwOXuPurgW3ADjPbTj4jkocbae3uk32Ai4DPV+avBa6d\nMk0t83EOcF9l/iHgzDB9JvBQmP4H4PJl28X+YXa364055Q34eeArwGuZDcI6MSw/Vi+BzwMXhekT\nw3Y2ddpr8nMWMwG/BLid2SM+veVr6mbIsdGegepI0JTZ4u6Hw/T3gC1hOsn8hhD1NcCdZJC3EKof\nBI4A+4Fv0XBEMjAfkRwj85HWz4b5xiOtaZCvqcUie3wm3cnecjKz5wOfBt7l7j+qrks1b+7+jLtv\nY3YlvhB45cRJ6kx1pPVQ+5haLOajPedUR4KmzONmdiZA+D4SlieVXzM7iZlQfNzdPxMWZ5E3AHd/\nktngwosII5LDqmUjkol8RPJ8pPUjwC3MmiLHRlqHbTrla2qxuBvYGnpsT2Y2QnTfxGnqg+oo1sXR\nrVeEOwfbgacqIX1UhCeFbwIedPcPV1YlnTczO8PMTgvTz2PWD/MgiY9I9jFGWkfQKbMT+CazduOf\nTJ2eFun/BHAY+D9mbcKrmLX9DgAPA18ATg/bGrO7P98Cvg5cMHX6V+Tr9cyaGPcCB8NnZ+p5A34D\n+GrI133An4blLwPuAg4B/wqcEpafGuYPhfUvmzoPDfJ4MXB73/nSCE4hRCOmboYIIRJBYiGEaITE\nQgjRCImFEKIREgshRCMkFkKIRkgshBCNkFgIIRrx/5jkS9pXSZ0yAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x113d56f60>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "img = np.zeros((400, 400, 3)).astype('uint8')\n",
    "def point(img, pt, size, color):\n",
    "    cv2.circle(img, pt, size, color, thickness=-1)\n",
    "    \n",
    "n = 10000\n",
    "r = 100\n",
    "r2 = 100\n",
    "k = 4\n",
    "x0, y0 = 200, 200\n",
    "for i in range(n):\n",
    "    theta = i * 2 * np.pi / n\n",
    "    alpha = theta * k\n",
    "    x1 = x0 + r * np.cos(theta)\n",
    "    y1 = y0 + r * np.sin(theta)\n",
    "    pt = (x1 + r2 * np.cos(alpha), y1 + r2 * np.sin(alpha))\n",
    "    pt = (int(pt[0]), int(pt[1]))\n",
    "    point(img, pt, 1, (255, 0, 0))\n",
    "\n",
    "plt.imshow(img)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "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.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}