2d cellular automata-for github.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "f969e83d",
   "metadata": {},
   "source": [
    "# Two-dimensional (2D) cellular automaton\n",
    "This is an extension of the elementary CA. Here, we have cells in a two-dimensional grid. Each cell has\n",
    "- **state:** which can be zero or one.\n",
    "- **neighborhood:** famous neighborhoods are the **Moore** (with eight neighbors) and the **von Neumann** (with four neighbors).\n",
    "\n",
    "Again, we have a **set of rules** that are applied to the current cell and its neighborhood to update the state of the cell.\n",
    "\n",
    "**Hint:** A well-known 2D CA is the **Conway's game of life**.\n",
    "<br> Here, we use the **von Neumann** neighborhood to create some interesting patterns.\n",
    "\n",
    "###### by Hamed Shah-Hosseini\n",
    "Explanation at: https://www.pinterest.com/HamedShahHosseini/\n",
    "<br>https://github.com/ostad-ai/Computer-Graphics"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "9840e83d",
   "metadata": {},
   "outputs": [],
   "source": [
    "from matplotlib import pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "f9dac8c7",
   "metadata": {},
   "outputs": [],
   "source": [
    "# the next generation of cells\n",
    "def update(neighborhood='neumann', no_alive_neighbors=[1]):\n",
    "    global nextCells,cells\n",
    "    if neighborhood=='neumann': #von Neumann\n",
    "        Nset=[(-1,0),(1,0),(0,-1),(0,1)] \n",
    "    else: #extended von Neumann \n",
    "        Nset=[(-2,0),(-1,0),(1,0),(0,-1),(0,1),(0,2),(2,0),(0,-2)]        \n",
    "    for x in range(gridW):\n",
    "        for y in range(gridH):\n",
    "            noNeighbors=0\n",
    "            for x2,y2 in Nset: \n",
    "                    xn=x+x2; yn=y+y2\n",
    "                    if (0<=xn<gridW) and (0<=yn<gridH):\n",
    "                        if cells[yn,xn]==0:\n",
    "                            noNeighbors+=1\n",
    "            if cells[y,x]==0:\n",
    "                # alive stays alive\n",
    "                nextCells[y,x]=0 \n",
    "            else:\n",
    "                #dead becomes alive having no. of alive neighbors \n",
    "                if noNeighbors in no_alive_neighbors: \n",
    "                    nextCells[y,x]=0               \n",
    "    cells=nextCells.copy()           "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "d0200ca9",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxoAAAF0CAYAAACkIU9RAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAATVklEQVR4nO3dy27kOBIF0OqB//+XaxazGLdgK0nq8qlzdgUWRTkfCgQYjPzn79+/f/8AAAAE/Wf2DQAAAOeRaAAAAHESDQAAIE6iAQAAxEk0AACAOIkGAAAQJ9EAAADiJBoAAECcRAMAAIiTaAAAAHESDQAAIE6iAQAAxEk0AACAOIkGAAAQJ9EAAADiJBoAAECcRAMAAIiTaAAAAHESDQAAIE6iAQAAxEk0AACAOIkGAAAQJ9EAAADiJBoAAECcRAMAAIiTaAAAAHFfs28ARvjnn3/+9e+/f/9OuhMA+B+xidPZ0QAAAOIkGgAAQNw/f+3TcajrlvQdXwMARhCbeBM7GgAAQJxEAwAAiJNoAAAAcc5o8ApaCAKwGrGJ09nRAAAA4iQaAABAnNIpjqWFIACrEZt4EzsaAABAnEQDAACIk2gAAABxX7NvAGrctQJsHbuO14wBgNgEP7OjAQAAxEk0AACAOO1tWV5NK8DRfH0A3klsgs/saAAAAHESDQAAIE6iAQAAxDmjwVa0EARgdSPOb4hH7MCOBgAAECfRAAAA4pROcayarWtfAwBSSstxP6kpAYYV2dEAAADiJBoAAECcRAMAAIj7mn0DvFOPWtNe7W0T96OWFmB9qdg0ovV6K7GJkexoAAAAcRINAAAgTntbhilt6VfzkVzp11e10wXYz92zu/VZvVJsuurx98Jv7GgAAABxEg0AACBOogEAAMQ5o8EUrS0EW9sEPmkh2NoKUAtBgL20nrWriWmp83ylcav1+pBgRwMAAIiTaAAAAHFKpximdSu397zr3NHzAJinR4wZXTp1N09sYiY7GgAAQJxEAwAAiJNoAAAAcc5o8EiPdrOjx67jK48B8FkqVtzpMW/0GUWxid7saAAAAHESDQAAIE6iAQAAxDmjQbWaWlSyfF0BfiY2zSM28Rs7GgAAQJxEAwAAiFM6xSMrtanV3haAn8xsE9vLzLa8UMqOBgAAECfRAAAA4iQaAABAnDMaDLNyjaxaV4BzlZ596/H8HxGbUnHz7hriFi3saAAAAHESDQAAIE7pFFO0bsmu1t629F59zQDWN7st7eg2uWITvdnRAAAA4iQaAABAnEQDAACIc0aDYRLt9lqv+UTr+ZHENQHoK9FCffTZiuuaNWcUS68JCXY0AACAOIkGAAAQp3SKKbS3BWB1pb/GXRPTUqVMraVb4hEj2dEAAADiJBoAAECcRAMAAIj7mn0DvFOqhW1rjWprK8CaOlx1sADnmNHCtlSqvS2k2dEAAADiJBoAAECc9rYsp0cL217tbbWwBXiHHr++XRObSq95nSc2MZMdDQAAIE6iAQAAxEk0AACAOGc0WMKJ7fd8tQD2JjbBM3Y0AACAOIkGAAAQ55fBmWJ0C9sZ7W21FwTYy6dY8dv/nfGr4TUx7juxiZHsaAAAAHESDQAAIE6iAQAAxGlvyzClNautta5PamTval17zwNgnrtnd+mzuuYaPWJMTUxJXQdK2NEAAADiJBoAAECc0im6WamF7WrtbbUXBJijJlb89v9q1yhdr6bMqndZl9hEgh0NAAAgTqIBAADESTQAAIA4ZzSIqmnb1+JJm8DEmqPXA+C51md3zXmKxHo1StdPxU2xiRZ2NAAAgDiJBgAAEKd0ikdWalN7Sntb7QUBnpld5pRab9X7FpsoZUcDAACIk2gAAABxEg0AACDua/YNsJ+7utDWlnoz513n7jIPgP87sb16q95t2cUmStnRAAAA4iQaAABAnPa2fLRSK9q3t7fVXhDgf1pLUq9WmtdanrXSPLGJ7+xoAAAAcRINAAAgTqIBAADEOaPBj1Zu28fPfJWB04lN+xGb3s2OBgAAECfRAAAA4iQaAABA3NfsG2ANK/2uhd/ReD720zjAbmqeazXPx9+MnvfEyr//8Z3Y9G52NAAAgDiJBgAAEKe97Uu1bpfWzN1l3nXu6fMAVvXpuVb6LLu7Tk2ZT+m91Ky30hqtr0VNTEldhz3Z0QAAAOIkGgAAQJxEAwAAiNPettBd27xdtLZGPX3sOn7C2HX8hPaCJ/wNkPa22FT6945oN1t6zVRsStzLE4nzMdfrnPBcP+Fv6MmOBgAAECfRAAAA4rS3/cUJrUNnbKX2XnOnX1/dcb3VnPA9hKQTWoU+ea4lym5a139SytpDa+lwy/U/6fFanPJ53uXv6MWOBgAAECfRAAAA4iQaAABAnDMahXZpITi6NeoJY9fxE8au4ye0vt3lPmGk0vr32d+X0hiaOtsw4nxB63M0dfZs9H2XXr9Gzes08zPsvGA7OxoAAECcRAMAAIhTOvWLXbbJWu/zyd/XulW/6rzr3NPn1cyd/XjY5XsIo+zynWhtw9ujrOjT/y1dv7U8qWaN1HVKY0WqxKz0vmrWX7mV8y7fwxXY0QAAAOIkGgAAQJxEAwAAiPuafQMjrdb+tNUurVF3GbuOnzB2HV+p9e1q7yHMlvrct645OjalzkGk7q10/Zr3KdVetscZxdEtbO/WT51R7B2baohN/2ZHAwAAiJNoAAAAcce3t+219Zcwe/typl7b4ytp3R7eRWtbxNUc/ghkUT2+E6PbWM/4Xo8uu1n52dWqx6+N311/9Gs44vM8Ir6fEpvsaAAAAHESDQAAIE6iAQAAxB1/RuO71dqflt7ria1RVxq7jp8wdh1fqfXtau8hrGalcwGt36Wa7/LKep9neOL0cydv+Oy9IR7Z0QAAAOIkGgAAQNyrSqdqtLY5G73e6Hk1c3eZd517+ryauaPbZV6N/h7CCkaU5LSWLP72/z5p/aXsVk/KTu+u85tev3a+0nsx+n369H9/0/pezP6uncqOBgAAECfRAAAA4iQaAABA3JZnNHrUuM1onak16hpj1/ETxq7jK43d6fUettIWlxqzY9NVa33/LucZUq/FqucgUvcz+r3/dLZh1dd7xOfyjbHJjgYAABAn0QAAAOK2KZ3q0VZz9i93nrbeiDVX+/XV09Ybsebs9/COdrrUSpV+lF7zTq9Sl5b1ev0ydY9Sl9b1Zv/69uhfLZ/52l/1aL07+vvz05qlejx3erGjAQAAxEk0AACAOIkGAAAQt80Zje9Wa5152th1/PSx6/gJY9fx08eu4yNaUMNV6xmf2bX/p683Y83Uervet/Xq12w927J6bLKjAQAAxEk0AACAuG1Kp1q3kHrPu849fV7N3F3mXeeePq9m7i7zrnNnvKa8U4/Pdo+ym9PXG7Hm7Pew5jo9nmurrlez5uj1UmvuHJvsaAAAAHESDQAAIE6iAQAAxE09o7FSe8y3tT/VGnWd19t7+K73kPWlPqN3nLfKzrvO3XXeLq/329+n0nk1c3vMWyE22dEAAADiJBoAAEDc8NKpXr+8CLADpVRrEpuAN+sVm+xoAAAAcRINAAAgTqIBAADEaW9bMHYdP33sOn762HX8hLHr+Olj1/GVx1if9rb7zbvO3XXeLq/329+n0nk1c7W3BQAAKCTRAAAA4iQaAABA3NQzGjXUG64xr2buLvOuc0+fVzN3l3nXuTNeU96px2f77rPV47N9wnoj1pz9HtZcp8dzbdX1atYcvV5qzZ1jkx0NAAAgTqIBAADEfc2+gVI17eBa5p3SOlNr1Pqx6/gJY9fx08eu460t/VZoBcheSssbUm1xa/S6rvXWd/prc+J6PcrBVmBHAwAAiJNoAAAAcRINAAAg7vj2tolrPjGiDneV9UasOXq965qnrzdizdnv4Z2dWwgyR2uL09Zr3mltv9pjvVR717s1e7UjLV2v199YqnT9Ea/F3X2lPl931+z9eb77v6nPc40ez51e7GgAAABxEg0AACBum9Kp72q2L1PXTLTOvFtTa9Q9W6OuOnYdX2nsTq/3sJX2ttSYHZuuEqUurevdSZXdpF6L1tfpTuJeUvcz+r3/VEq06us94nP5xthkRwMAAIiTaAAAAHESDQAAIG7LMxojjG5z2bre6Hk1c3eZd517+ryaual5rbSb5Y16tCa92uXMQqtdz53czZ39Xox+nz7939+MOJNSI3F+cWd2NAAAgDiJBgAAEPeq0qnV2p+W3qvWqNrb1o5dx09oYTujBTXMkCrZSGj9Lo0ouxmh9G/sVXZzp/V5vMvr/4bP3hvikR0NAAAgTqIBAADESTQAAIC4489orFaP912iVduuZtSzjrZyXWhCqp3ibIc/AllUj+/E6HbUM77XPVqDn3CeoUbvcyez4/uIz/OI+H5KbLKjAQAAxEk0AACAuONLp75brf1pq9Nbo2pvu8d72Gq19xBmS33ur1pLZFqNLrtpLXma/Yvmqfe3dV6PX9gunXedO6L8rVSqpE5s+jc7GgAAQJxEAwAAiJNoAAAAca86o1GjtW5wtBHt2Ea3uJvRUq9Hzeiq82rmzn487PI9hFF2+U70qu9vvU6Pcwmr/Y2/zetxXmWlv2+GXb6HK7CjAQAAxEk0AACAuK/ZN7CqXVqQ7doadaWx6/gJY9fx0S1sU3a5TxilpnxnZtwa8evQI0qCvo/VPEdHay3rqnmuzv57Vy0FvBKb/s2OBgAAECfRAAAA4iQaAABAnPa2vzihHm9EDeWIOtyZ613XPH291ZzwPYSklVt+lnryXEuc22pd/8mZuR5az+G1XP+THq/FKZ/nXf6OXuxoAAAAcRINAAAgTulUoVXb29bQ3vbnsev4CWPX8RNaw57wN0Da22LTd7N/1bpVa2zqfS9XvUrTxKZ3saMBAADESTQAAIA4iQYAABDnjMZLPWkbWjp3l3nXuafPA1hVqlVo69mOHs/jEbEpFTdLr9F6PubJddiTHQ0AACBOogEAAMQpneLPnz9rtZvV3jbf3hZgRzXPtcSvU88o3W21Uulujza5nMGOBgAAECfRAAAA4iQaAABAnDMa/GhEfSlZvsrA6cSm/YhN72ZHAwAAiJNoAAAAcRINAAAgzhkNPlrpNy/8jsbvYwBvUvPMv7PSvF1+/0NsopQdDQAAIE6iAQAAxCmdotqOW7k1c3eZB8D/3T1L756dPWJMYr0n1+kRY1r/Xt7NjgYAABAn0QAAAOIkGgAAQNzX7BtgP6u2qX1De1sAftZ6LqFUj2v20vtexSZK2dEAAADiJBoAAECc9rZEjd6uHbGV3Xs7/m49AJ5rfXbXlOAm1qtRun4qbopNtLCjAQAAxEk0AACAOIkGAAAQ54wG3azUwnan9rYA9FMTK377f7VrlK5Xc56i9H5qriE2kWZHAwAAiJNoAAAAcUqnGKa0pV5rK74nLfxaW9gm5gEwz8zypE9mlnVBgh0NAAAgTqIBAADESTQAAIA4ZzSY4u3tbQFYz6dY8dv/fXIOo1VNjPtObGIkOxoAAECcRAMAAIhTOsUSRmwzj+arBbA3sQmesaMBAADESTQAAIA4iQYAABDnjAbLOaW9LQDnqIkxd1pjU+k1r/PEJmayowEAAMRJNAAAgDilUywvsXVcM7d1HgBnaf3179Jf5v6k9TqtvxoOaXY0AACAOIkGAAAQJ9EAAADivmbfAO9UUzO6cnvbO9oLApwj1d62hydnFKEnOxoAAECcRAMAAIjT3pZhWtvUJq75ROn9aIsLsJ8ev749Ojb1+NVySLCjAQAAxEk0AACAOIkGAAAQ54wGU9S0ty2dp70tAE+sdLZixP2ITfRmRwMAAIiTaAAAAHFKpxgm0UKwx7zr3NHzAOirtJS1x/N/RGxKxc27a4hbtLCjAQAAxEk0AACAOIkGAAAQ54wGj/RoNzt67Dq+8hgA9VY+I9jK2UJ2YEcDAACIk2gAAABxSqeoNmJLmJ/5ugL8TGyaR2ziN3Y0AACAOIkGAAAQJ9EAAADinNHgkZXa1GpvC8CfP7lYcafHvNHtdcUmerOjAQAAxEk0AACAOIkGAAAQ54wGw8ysNa2ZO3oeAPP0iDF3z/hUrChdX2xiJjsaAABAnEQDAACI+5p9A7xHTUu/lnkz2tu2XhOANZSWGaXa4vaQKjmGNDsaAABAnEQDAACIk2gAAABx2tsyTI+a0RF1saX3o4UgwH5a29S2XjOlx72JTaTZ0QAAAOIkGgAAQJzSKaaoaW+buOaT9raJ+/E1A1if2ARZdjQAAIA4iQYAABAn0QAAAOKc0eBY2s0CMMPdOYjW2NTj/Aj0ZkcDAACIk2gAAABxX7NvAGq0tgl80kJQK0AAWj35lfARvzAOPdnRAAAA4iQaAABAnEQDAACI096W5a1co+rrA/BOYhN8ZkcDAACIk2gAAABxSqfYiva2AKxGbIKf2dEAAADiJBoAAECcRAMAAIhzRoNj1bQe9DUAYASxiTexowEAAMRJNAAAgDilU7zCpxaCADCa2MTp7GgAAABxEg0AACBOogEAAMQ5o8GxtBAEYDViE29iRwMAAIiTaAAAAHFKp3gFLQQBWI3YxOnsaAAAAHESDQAAIE6iAQAAxDmjAQAAxNnRAAAA4iQaAABAnEQDAACIk2gAAABxEg0AACBOogEAAMRJNAAAgDiJBgAAECfRAAAA4iQaAABAnEQDAACIk2gAAABxEg0AACBOogEAAMRJNAAAgDiJBgAAECfRAAAA4iQaAABAnEQDAACIk2gAAABxEg0AACBOogEAAMRJNAAAgLj/As6v0t3eMO3bAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 1000x500 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gridW, gridH = 100,100\n",
    "#holds the cells of cellular automatun\n",
    "cells=255*np.ones((gridH,gridW),dtype=np.uint8) \n",
    "#initialize grid cells with one dot in the middle:\n",
    "cells[gridH//2,gridW//2]=0; nextCells=cells.copy()\n",
    "for _ in range(45): update()\n",
    "cells1=cells.copy()\n",
    "#---another one\n",
    "cells=255*np.ones((gridH,gridW),dtype=np.uint8) \n",
    "cells[gridH//2,gridW//2]=0; nextCells=cells.copy()\n",
    "for _ in range(45): update(no_alive_neighbors=[1,3])\n",
    "plt.figure(figsize=(10,5))\n",
    "plt.subplot(121);plt.axis('off'); plt.imshow(cells1,cmap=plt.cm.gray)\n",
    "plt.subplot(122);plt.axis('off'); plt.imshow(cells,cmap=plt.cm.gray)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e75dec1d",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.11.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}