p1/izzivi/04 veliki kih/veliki-kih.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Veliki kih\n",
    "\n",
    "Neka visokomogočna in precej velika bitja, so zbrala 65536 zvezd najrazličnejših svetlosti. Razporedila so jih na dvodimenzionalno mrežo velikosti 256x256. Če so stvaritev pogledala od daleč, je bila videti kot črno-bela slika. (No, v resnici črno-oranžno-rumeno-bela, saj gre za zvezde). Zvezde so torej predstavljale piksle neke oranženkaste slike.\n",
    "\n",
    "Slika je bila imenitna in znana in imela celo ime. Rekli so ji ... bogvekako, tudi ime se je izgubilo.\n",
    "\n",
    "Dolga leta je bila razstavljena. Potem je eden od ogledovalcev kihnil. Močno kihnil.\n",
    "\n",
    "Zvezde so se razletele vsaka na svojo stran. Ene počasneje, ene hitreje; ene levo, ene desno, ene gor, ene dol; večina, pač, nekam postrani.\n",
    "\n",
    "Okrog zvezd so se sčasoma naredili planeti. Na planetih je zrasel najprej en sam plevel, kasneje pa tudi solata in žirafe. In neka manj mogočna in manj velika bitja, ki pa so sčasoma odkrila teleskope in spektre in tako naprej. V datoteko [meritve.txt](meritve.txt) so zbrala podatke o meritvah. Vsaka vrstica datoteke predstavlja eno zvedo: prvi stolpec pove njeno svetlost, drugi in tretji njeni koordinati, četrti in peti pa njeno hitrost v smeri x in smeri y (torej: za koliko se spremenijo koordinate zvezde v enem stoletju). Ker so ona, prva bitja, bila in tudi ostala zelo diskretna, so vsa števila cela.\n",
    "\n",
    "Kako je bilo ime sliki?\n",
    "\n",
    "Namig: za začetek bo treba ugotoviti, kaj je bilo na njej ... pred velikim kihom.\n",
    "\n",
    "### Tehnična pomoč\n",
    "\n",
    "Izziva se lahko seveda lotiš v poljubnem jeziku. Ker se na predavanjih pogovarjamo o numpyju in Pythonu, pa je tule pomoč za numpy in Python.\n",
    "\n",
    "Podatke najlažje preberemo z"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[    161,  -73865,   37242,      -6,       3],\n",
       "       [    161,  -61517, -110923,      -5,      -9],\n",
       "       [    157, -123325,   98876,     -10,       8],\n",
       "       ...,\n",
       "       [    104,  -12324,   12370,      -1,       1],\n",
       "       [    104,   98860,   86456,       8,       7],\n",
       "       [    108,   74256,   12459,       6,       1]])"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "np.loadtxt(\"meritve.txt\", dtype=np.int64)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Kako to premetavati naprej, odkrij sam(a).\n",
    "\n",
    "Če želim v jupytru pokazati sliko dimenzij (30, 50), ki ima neko bolj temno zvezdo na (5, 7) in neko svetlejšo na (25, 33), nam je storiti tole:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7fb8de199550>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAADrCAYAAAB5JG1xAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAACz9JREFUeJzt3UGInPd5x/HvL5JDS12o3MRGSEoVgigRoZFBCINzcENTFDdUziEQQ4sOBeUQgw0uxfUlbqCQQ+O0hxDYxsKCpg4GJ7UwpalQXZxDcSM5biyjFBvjJIqFhDHF6qVB9tPDvkFbVauZnZmdkZ/5fmCZmXfenffhj/1l+O+MnapCkvTe975FDyBJmg2DLklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpia3T/HKSg8DfAFuAb1bVV0ac79dSJWnj3qyqD446aeJ36Em2AF8HPg3sBe5NsnfS15Mkresn45w0zZbLAeDVqnqtqn4BfBs4NMXrSZKmME3QdwA/W/P43HDs/0hyJMmpJKemuJYkaYRp9tBzjWP/b4+8qlaAFXAPXZI20zTv0M8Bu9Y83gm8Md04kqRJTRP0HwB7knw4yfuBzwPHZzOWJGmjJt5yqarLSe4DvsfqxxaPVtXLM5tMkrQhmef/scg9dEmayOmq2j/qJL8pKklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNWHQJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCYMuiQ1YdAlqYmtix7gvaIeWf+5XOc5SZqXqYKe5HXgEvAOcLmq9s9iKEnSxs3iHfrvVtWbM3gdSdIU3EOXpCamDXoB/5zkdJIj1zohyZEkp5KcmvJakqTrmHbL5c6qeiPJrcCJJD+uqufWnlBVK8AKQJKa8nqSpHVM9Q69qt4Ybi8C3wUOzGIoSdLGTfwOPcmvAe+rqkvD/d8HvjyzyW4wfjRR0o1umi2X24DvJvnl6/x9Vf3TTKaSJG3YxEGvqteAj89wFknSFPzYoiQ1YdAlqQmDLklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNWHQJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCZGBj3J0SQXk5xZc+yWJCeSvDLcbtvcMSVJo4zzDv1x4OBVxx4CTlbVHuDk8FiStEAjg15VzwFvXXX4EHBsuH8MuGfGc0mSNmjrhL93W1WdB6iq80luXe/EJEeAIxNeR5I0pkmDPraqWgFWAJLUZl9PkpbVpJ9yuZBkO8Bwe3F2I0mSJjFp0I8Dh4f7h4GnZzOOJGlS43xs8Qng34DfTnIuyZ8AXwE+leQV4FPDY0nSAqVqftva7qFL0kROV9X+USf5TVFJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNWHQJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCYMuiQ1YdAlqQmDLklNjAx6kqNJLiY5s+bYI0l+nuTF4efuzR1TkjTKOO/QHwcOXuP416pq3/Dzj7MdS5K0USODXlXPAW/NYRZJ0hSm2UO/L8mPhi2ZbTObSJI0kUmD/g3gI8A+4Dzw1fVOTHIkyakkpya8liRpDBMFvaouVNU7VfUu8LfAgeucu1JV+6tq/6RDSpJGmyjoSbavefhZ4Mx650qS5mPrqBOSPAHcBXwgyTngS8BdSfYBBbwOfGETZ5QkjSFVNb+LJfO7mCT1cXqcbWu/KSpJTRh0SWrCoEtSEwZdkpow6JLUhEGXpCYMuiQ1YdAlqQmDLklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNWHQJakJgy5JTYwMepJdSZ5NcjbJy0nuH47fkuREkleG222bP64kaT3jvEO/DDxYVR8F7gC+mGQv8BBwsqr2ACeHx5KkBRkZ9Ko6X1UvDPcvAWeBHcAh4Nhw2jHgns0aUpI02taNnJxkN3A78DxwW1Wdh9XoJ7l1nd85AhyZbkxJ0ihjBz3JzcBTwANV9XaSsX6vqlaAleE1apIhJUmjjfUplyQ3sRrzb1XVd4bDF5JsH57fDlzcnBElSeMY51MuAR4DzlbVo2ueOg4cHu4fBp6e/XiSpHGl6vq7IEk+AXwfeAl4dzj8MKv76E8CHwJ+Cnyuqt4a8VpuuUjSxp2uqv2jThoZ9Fky6JI0kbGC7jdFJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCYMuiQ1YdAlqQmDLklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDVh0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNTEy6El2JXk2ydkkLye5fzj+SJKfJ3lx+Ll788eVtBFVte6P+tk6xjmXgQer6oUkvw6cTnJieO5rVfVXmzeeJGlcI4NeVeeB88P9S0nOAjs2ezBJ0sZsaA89yW7gduD54dB9SX6U5GiSbTOeTZK0AWMHPcnNwFPAA1X1NvAN4CPAPlbfwX91nd87kuRUklMzmFeStI6M88eRJDcBzwDfq6pHr/H8buCZqvrYiNfxLzHSHF3v3+8kc5xEUzpdVftHnTTOp1wCPAacXRvzJNvXnPZZ4MwkU0qSZmOcT7ncCfwx8FKSF4djDwP3JtkHFPA68IVNmVDSxHwXvlzG2nKZ2cXccpGkScxmy0WS9N5g0CWpCYMuSU0YdElqwqBLUhMGXZKaMOiS1IRBl6QmDLokNWHQJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCYMuiQ1YdAlqQmDLklNGHRJasKgS1ITBl2SmjDoktSEQZekJgy6JDWxdc7XexP4yXD/A8NjuRZruRZXuBZXLPta/NY4J6WqNnuQa184OVVV+xdy8RuMa3GFa3GFa3GFazEet1wkqQmDLklNLDLoKwu89o3GtbjCtbjCtbjCtRjDwvbQJUmz5ZaLJDVh0CWpiYUEPcnBJP+Z5NUkDy1ihkVJcjTJxSRn1hy7JcmJJK8Mt9sWOeO8JNmV5NkkZ5O8nOT+4fjSrUeSX0ny70n+Y1iLvxiOL91aACTZkuSHSZ4ZHi/lOmzU3IOeZAvwdeDTwF7g3iR75z3HAj0OHLzq2EPAyaraA5wcHi+Dy8CDVfVR4A7gi8M/C8u4Hv8DfLKqPg7sAw4muYPlXAuA+4Gzax4v6zpsyCLeoR8AXq2q16rqF8C3gUMLmGMhquo54K2rDh8Cjg33jwH3zHWoBamq81X1wnD/Eqv/Au9gCdejVv338PCm4adYwrVIshP4A+Cbaw4v3TpMYhFB3wH8bM3jc8OxZXZbVZ2H1cgBty54nrlLshu4HXieJV2PYZvhReAicKKqlnUt/hr4M+DdNceWcR02bBFBzzWO+dnJJZbkZuAp4IGqenvR8yxKVb1TVfuAncCBJB9b9EzzluQzwMWqOr3oWd6LFhH0c8CuNY93Am8sYI4byYUk2wGG24sLnmduktzEasy/VVXfGQ4v7XoAVNV/Af/K6t9alm0t7gT+MMnrrG7HfjLJ37F86zCRRQT9B8CeJB9O8n7g88DxBcxxIzkOHB7uHwaeXuAsc5MkwGPA2ap6dM1TS7ceST6Y5DeG+78K/B7wY5ZsLarqz6tqZ1XtZrUN/1JVf8SSrcOkFvJN0SR3s7pPtgU4WlV/OfchFiTJE8BdrP7nQC8AXwL+AXgS+BDwU+BzVXX1H07bSfIJ4PvAS1zZL32Y1X30pVqPJL/D6h/7trD6RuvJqvpykt9kydbil5LcBfxpVX1mmddhI/zqvyQ14TdFJakJgy5JTRh0SWrCoEtSEwZdkpow6JLUhEGXpCb+F5Mq2HMkPYZsAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "\n",
    "\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "slika = np.zeros((30, 50), dtype=np.uint8)\n",
    "slika[5, 7] = 64\n",
    "slika[25, 33] = 128\n",
    "\n",
    "plt.imshow(slika, cmap=\"afmhot\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Z `dtype=np.uint8` v klicu `np.zeros` povemo, da matrika vsebuje cela števila (osembitne nepredznačene `int`-e). V tem primeru `plt.imshow` pričakuje vrednosti med 0 in 255, kar ravno ustreza svetlostim zvezd. (Kakšno naključje.)\n",
    "\n",
    "Argument `cmap=\"afmhot\"` v klicu `plt.imshow` določi barvno paleto. Če bi napisali `cmap=\"gray\"`, bi bila slika črno-bela. Ker gre za sliko iz zvezd, bo ustrezna paleta `\"afmhot\"`.\n",
    "\n",
    "Namig: ni težko, samo simulacijo je potrebno pognati nazaj. Zadnja dva stolpca povesta, za koliko se drugi in tretji stolpci spremenita v vsakem stoletju. Odkriti morate, za koliko stoletij **nazaj** je potrebno iti, da se bodo vse številke v drugem stolpcu razlikovale za največ 256; v drugem pa ravno tako. Ko dobite takšne koordinate, pa izrišite sliko po gornjih navodilih."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}