{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tissue neighborhood quantification\n", "\n", "In this notebook, we will analyse neighborhood-relationships between cells. \n", "We count the number of neighbors for each cell and take a look at this number in a parametric image.\n", "Afterwards, we average this number locally between neighbors.\n", "\n", "## Test data\n", "Let's generate some cells and take a look at the borders between them" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pyclesperanto_prototype as cle\n", "import numpy as np\n", "from skimage.io import imshow\n", "import matplotlib" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAz5ElEQVR4nO2dedxV0/7HP8uDBikNIiVNT2j41S/3RqMr0mBImYq65vwoRVzKcOW6DRIRiQchmqjETRJxy5DbpG7ziIqSSqaItH5/nLNO69nPHta0z977nPV+vc7rOWefvdf5djzP2/f7XWvvTSilsFgsFhkOizoAi8WSPKw4LBaLNFYcFotFGisOi8UijRWHxWKRxorDYrFIE5o4CCEdCSHrCCEbCSEDw/oci8WSfUgY6zgIIQUA1gNoD2AbgEUAelBKVxv/MIvFknXCyjiaA9hIKd1MKf0NwGQAXUL6LIvFkmUOD2nc6gC2cq+3ATjda2dCiHDaU+2UWupReVDwy6/Gxiq/fbexsZwcqEVCG/vrw44BANAfzXwX+yvXNDIOAFTftMbYWE4OljvC+JiVq1Q0Ot7n+37yfO/7beWMfhZ+376LUnps0G5hicPtN7yYHAghvQH0lh34uvEPqMbkScWV642N1eGBl4yN5eTbonASxPtKd0XV9PM/3l9rZMxN1441Mg4A3Hfxn42N5eTnlscZH/Ov119qbKxeSz9GjYB9Zt/e0tjnHfj6H1+K7BdWqbINwInc6xoAvuZ3oJQWUUr/RCn9U0gxWAS4r3TXYq8L2p0SUSSWJBGWOBYBKCSE1CaEHAmgO4A3Q/qsWHH3V89FHYIlR+i19OOoQ/AkFHFQSg8A6AvgHQBrALxKKV2lO+6IOmVRvvtNKN/9Jt2hjDOvdxvM690GQEoecRfIfaW7lsg2GHHLOv73imOiDiHryEij4yOfaH9epc23oNLmW4T3D20dB6V0FqW0PqW0LqV0iM5YI+qUxYg6ZYtti6M8nJiSx6ayjbGpbGMjYwElyxM34iKPt6e1B5Aceay+8R4MLKiPgQX1lcfIdqYhIwxGKOs4pIPwmFVxysKLHybrN+J0GqQs0/BjaPXrlcZ2CqP8WPXETUQYPDqN0pZDNuDlTe8pHctk4cZnE/cqRuSOyebo6hvvcd0+/A+53y1Vcag0SZ3S2Fn62CUifcfYLjkXlQagnn3sebJM5rHp302w6d9NpMcQkYYq2c4yTNFyyAYAQK+652QepjCdeRz1yTfaY6y+8R5PaQDQyj5kkC1ZVDINRiwzDhlp8MhkHnueLOP7ft2/LA8cQ1YaIlmHiCxksg5dYchmHUwaMvCZiV+2wWMq86jyjwuLvf5y9n+kjvcThhORzEO3TBHNOrykIZpxxE4cqtJgiMgjSBpO3CSim2k4JSKbXYjKI5viUJEGzxUj5Ba46crDKQ2GiDxkhOHETyDZEIdfppE4cegKg8dPHrLS4GECMVGeXFf4uPYYQfIwUZ6IiENXGDxhy8NLFk7c5KEjCy94iZhqivrJI6g8SZQ4apQqoP1rlDY+7g+Tx2qJwottE8trHW9CGoC3OMLoZ7gJxKQwAHlpuOElElFhePH2Seb6NE7uLT0ZN334i/FxmUBkehmJb46aIAxpxIkfbmpYYltYTdC4TM8G4dY81ZVGmNxbejIAYGwb87+rJtZ3eJHT4nh75lVRh5BVwp454eVhOtswTZV/XJh5xBUmDUYY8tCZOfEjZ8Uxdm+qBMsXeWRrurWg3SmJkIZpOn2ptibFC6c0GGHIIwxyUhxMGoxclsd7E+fgvYlzcOa4+K+ktaTwkgYjCfLISXG4YUwe84ZpHW6qMeqGlUeKaz58LOoQPAmSBkNXHmN6zsKYnrOw6KgKWuN4Edb1OCLDmW0YJS2NGtUHYdtXegIxwcI2o0ps4+Uxz+A1MZLCOd8uBZCSxwttbo02GAei0mCMbVNGarZlTM9ZsiEpkzcZB6CZdTgyjRrVB2lGEz5hZCDTXv0Ud94YeIEoJRZOvUbreCYNRpwzD1GCMg+WWWRTGkCOreMQzTY6nS9xlS6B0kQ0+zBZprhlG7IEZSTTXv3U9/0Rz3yr9fn1ml/r+V7zS14QHscpDCe6mce5n12UeX7bRd6X8fNDNtvgYVmHjhz+/PP3Qvvl3QIw2RLFVx6SfQwRcZjubZgQB8MpkCBh8KjKw08aDBF5BEmDoSoPXho8MgLRkQYAHDxpmtbxDBF52AVgAbiWLfOGaTc/3Ri5ZBROnXzQ2HgmpQGkShpW1shIA4BS2SIiDSC4dBGVhipe0gCAUTPELhKsKw2TmGyU5kTGodMQ7XT+S0Zk4cw6Ri7x/+Ne013d2abFwehXTr0HJJJ5iArDiVvmISsN2YzDTxo8fpmHKWmYyjgYfplHojKOo2vUQtthL0bz4YYzjJFLRgVKAwBOnXww85AhjtIAgjMPVWkAJTMPlUwjrEapV+YRp0zDiYnMIxbiYLQd9qK0QEKdfpVAVBhuqEokbnjJQ0caDN0ZF0BcHqLZBmPUjHLFBBJnaZgilus43OQxf9DVJbbFRRomYfLQKWVk0c02eO688dhiZYsJaTAWTr0Gd58ZzrkXDFlp8IyaUQ67u8f7ItWMRUdVEJ5pcSOW4nCDl0mPm3J3CTmDzz6YRMIqU0xz543HYvpnF0QdhjQ60kgiOvKIVakiyqSx4d0tLY6cOvlgbHsbScOrXMk3aTBU+x2JFMfIn+vhtJHxvVmNJd445ZGv0mCoyCOR4mDkkzzKDpoedQg5BZNH3KWxclQlrL71htA/p82O1WizY7Xw/okTx8if6xV7fdrIj3NeIK+VPxOAlYdpwpBG5clq989xY+WoSpnnYcpjwKYF0sckThxe6Ahk5qqjDEcTHvkqj0n9J2JS/4lo3LRF1KFkBV4ajDDkoSINIEGzKmHBpBEkj/Mb/lzsGPZ65qqj8MWpm4Hzv8HGmeZuoMRg2QZP2UHTsW9YN+OfFUcm9Z9YYhsvjxXL1H7xw+TFhc8DAK5ufp3S8W7SYPDyaPDYs0rjmyBRGYezTHFDJuv4Y8tpwvvOXHVU5sG//uLUzZl96p2/Qng8XcoOmh7r7KPHylbaY7hJw4lqBvJVuX8rHRfEI3UOxcMEIoOfNJysvvUG5SxkwKYFytkGEJNzVU6tW5e+OGy47z4i0vBiyR2pX2IZUYjCi8OJbgbilm34IZOF/DIn1Um/q9tFUp8RRKmXS/53nNRIroQUEYYffllIQcf7ir3e8tHZWp/F4IXhRlD2ISOMINwyEVFJTKt1ZnJOqw9bHEBKHtkWB0NFILLSYHjJg4nCDVPycJMGQ0QeusJw4hSIUxpOZCUSJAs33ARiUho8TCAymYUVh4OFl/XTOt4NEWk4EZGIqjSc7BvWzVcYTnQE4icNhp88TEuDsWLZgkBh8IjIQ0UYblzd/LrQpAEAz90iv9xfVByJaI7qSiNOsD6Il0BMSQPwzzLceGj6DCl5iMiCp8fKVq7yCEsaQHCW4aRm67me8jAlDEaY0gibRDVHc4l656/IPBgmpaHKQ9NnCO0nKw2Gs2kapjRWH69228aareeiZuu5mdeP1GlhXBpJx4rDUgI/eZR6ebiyNBhMHmFKwwQ1W89NrDBUyhQZYl+q5FKZ4kW981cYXwOy+yy961fwZYuuKNwIWxqq2YZFDK2MgxDyBSFkBSFkGSFkcXpbJULIu4SQDemfFc2EmtvEoUxxYiK78GL2C4tCGdc0j5/wz6hDiCUmSpWzKKVNuU7sQABzKaWFAOamXyuRD9lGPpMUeYTFJfc9EnUIyoTR4+gCgF0w4yUAF6kMEndpdF9aHgMnNI06jERyZMvmmef5Ko/h7VJ/ekmVh644KIA5hJAlhJDe6W3HUUq3A0D6Z1XZQcOQRvUGU7XH6L60fObBsPKQg5cGI9/kwaTBSKI8dMXRilLaDEAnAH0IIW1FDySE9CaELCaELN77ww+aYfjz1eATAOjJg5dFEmjU+VLtMXrWudFAJIdwkwYj3+SRdLRmVSilX6d/7iSEvA6gOYBvCCHVKKXbCSHVAOz0OLYIQBEAlK5fj/asVOrQm5W2AgCabj1ROTYmCydMHl+tviRwjKTJglHt4dStAht1vhQrZ72mNVbPOjfilc3PaMfkJw3G7BcWoeM1f9b+rDjjzDYYLOuY+uDt2p8R9lQsoJFxEEKOIoQczZ4DOBfASgBvAmAXsrwKwBuqn7HsxK3Sx3w1+ARPafB4ZR9u5UgQcSlXqj08KyMNk+hmHiLSYJjIPJI8FXvJfY9kHioUjr8ND/3Z3B3bvNApVY4D8BEhZDmAhQDeopTOBjAcQHtCyAYA7dOvlVGRhyhOeSQxw2CyCEMYPKrykJGGpTiyAikcf1vmedjyiMVJbqXr16MnPvWo0L58+SKSWYjS5hUz6d3wK5dJH3PCJWszz2/p96TvviqCkC1X1h8ntqYkqITRlYZs2TLsx4LM8y6FZ2l9NgC0W7bHdXvjqmK/qzxeJYoszlKGl0UQdy0KvhVCos6OlREHABx73enGY4hKHLw0GG7y0M0oROUhKg0eN4GYyjT85MGLwgtVgXhJw4moREyJg7G8cJvysX4CseKQxJQ4gGB5uMnCCS8PU2WIiDxUxMHgBWKyRHHKQ0QYTmQEIioNJ14SMS0NQE8cPE6JiIrDnuQWAs5m6QmXrC32EOGJ0X1DiCxcWB/EdF+DNUyH/VigJA0AeGPDB0L7qUoDAFbsHKB8bFSo9kKsOEJi4ISmUqJw44nRfY02PYPWduhkGwzTaz8YqsLgCZKHjjQYTnmEkW2YRkUe8f9XJZTxr8fzjuVe8jAhjbD47JIzjI3lJQ8T0mDkQ+ZhxRFz7i5lfprVKY84SyMM3tjwQeYBmJUGI4ny+L3vF8L7xv56HG68XLEOen0nf73PbBHXbIOHyWP6EteFvXlDGNJgJKFMAYABpy+XPiYZ/zKO2Xf0AJCShyU/OLhE+BQoiyQq0gASKA4ek/L4sOc4Y2NZzLGk9p0AkimP98tfHnUIoZEocbBsI86EUaaUKRP+uQemePX4z4yNxaSRRG7vvR1AvOWhmm0ACRHH7Dt6eEojX0qWMmUqJEogurhJI4ys4/BJ5m/kzKTBiJs8Bpy+XEsaQMzF4ScMHl15/PT3r/DT37/C2+sfVB7j+ikX4vopF2L+b+FehJcJxIRETDdGa715Mu4s6q49TrYyjTHPpVbnmpSHUxpxQ1cYjNguOVcpS2RnWn76+1eu2zvVF7uJz/VTLvR8r+2RV0jF4seow+S+i19+CT6ZiWFCHrXePLnEthG95Uo2GVkcdtp8qbHdYNJwcqCH+h3gg6TR7ocpymMzHrutbLHXZ81cL3ysiDTenHlVcs9V0ellBMnDSxZO/OThJwweE/KQlQYPL5AyZSpkXr+2sOStEwoKHlL6DDdpOAmSiGyGoSsOL2kwZOUhk2WoysMpDB4ReYhmGokTx6qyZlI8Jg5RQYgy7TP1C8yoCkRHGk7GzxO/kz3DTSYiogji8mE1tY7XEUeQNHiCBKJTlogIxE8WbngJRKY8ERVHLBaAFW77HKhf2shYpoUB6EkDAOb/NtFo6RIFJoRhioNL2irJQ0YawKHeh1MgJvoY75e/3FMessJgfHB+/RLyMNXTcBILcVjCRSXbAIA//rgLAFD3rRkGo4kGWWnwMIH0/+AfpsIBUFweqrJw8sH59TPP//Wt3vVm/Yj1rEouITvbYrJMiRtTBm3RHkNmalZHGmETt6laUaw4skjYU7VuqGYbPJvOu0g/kIiIszQYprKNbGLFkWVE5JHL2UY2SYI0kortcURALjRLo+aYRl8A+/VmZyzqWHEEsLZGCzSskXq+6l+jjY3b78vdmeejT6qceR7XbGPTeRclrkk6qNcIAMCwl5N7zktcsaWKD2trtAhl3KnbOxR73e/L3ej35W7j0tjb0tyVs+LEMY2+iDoEozx9STgXUvq99LBQxgWsODzJljR4Pv/cTE2+t+UZGWlcOOhrI2P2eagAHVdebGQsyyFufjuc+yZP/7EZgPDkYcXhgpc0Gl7QT2tcP2nowmQRRpbR56FDFwpOojxYyRI3eGmYzDqYNBhhyCOnxNH1lg+1jl9bo0UkmQaPbNYhIgudrIOXBiNKeaiWKXGVB48JeTilwTAtj1iIY3/96M+XEUU265i6vYN0piEqj7B7GG7SYHRcebGyQGoPOR4LD/6GhQd/Uw1NiTjJw6tEefqSM5UEMv3HZp7SYJiUR2xmVTa9/3vmed12R0gf/2upJphU1Bc9eqv1CWQyjYYX9AucYQm7LIkLHVdejNmNpgnvX3vI8cVeu8mj+WFHuh4bRVN04JCTDr345PnM019aXqc8pkhfg8nj/6bO890vSBZOmDyO+HWQ1HFOYnF2bP2mh9En3/d2mJdIfi3VxHW7jDx0ShOnPEzKonZt9zu56UjjzWFyN+n2yziciMjDKQ0ReImYEIfo1GwxYXigIg+VZqiXPGSl4YZTIIk6rT5IHDzVOwX+m1zhZWKyj8HkEUaGwcvDRJYRpjicMJGoyMKLc//HzAyRnzxEhOFERCAmZk+qvvCj9hhu8PIQFUcsehzZYFJRX0wq6htK8zOssoT1OkyVJjJNUh1pMExKwyTOXsfAISdlHiqU+eR5lOHKmKSh0vuITY8jqYTZywCi6WfksjQYg3qNAD1ljNExmTycGUhYazVM8nvpYVJ9j7zJOCyWbMFnH0mQBgCcWaE5Wh43V3j/vBPHZ38/LeoQIiWoXDGRbeQ7jz3ZFc9csQfPXLEHTV4+EHU4gZxZobn0MXknjiQxtup8DH3X/BmgppahZ5Mpc+bhupEbog4jkMee7Bp1CFKoSAMQEAchZBwhZCchZCW3rRIh5F1CyIb0z4rce4MIIRsJIesIIeE2AHKYsVX1bwHgh5s84pptTJlzaDryupEbYiuQpElDB5GM40UAHR3bBgKYSyktBDA3/RqEkAYAugNomD7mKUJI7H4bb634aPBOEeKUxtB3a2YeJmHy6PNQQSylMWXOvGLS4NEViOnGqJ804liunFmhuXK2AQjMqlBK5xNCajk2dwHwl/TzlwD8G8Bd6e2TKaX7AXxOCNkIoDmABcoRcpTqsxS76gBVNvdWHuOlJ64EUFwej303QHm8+3ufjAeK1ikf7yQo0+DlcXd7vWt33l/wC3DvJgD1A/fNNl7CcMLk8fwdhUL7f3n4p6knG09DzXpLlGJzkrRMQ0cYDNUex3GU0u0AkP5ZNb29OoCt3H7b0tu0KdVnaeb5rjpFJobMcGvFR6WzkPWzemQePWoMRo8ag7XjkC1PdDKQ+wt+yTy/9kfxu4GJMvaJt5SPFZUGj0oGsmWjfqM8m9LYec3RWsfrZhk8ppujxGWb69JUQkhvQshiQsji73fLr15VkQfLNrwQFcj6We4X3NERSNg9DR5eGoxrf1xvTCC3nnM3gJQ8dASigp88MtkGh448ZKQRZbliUhgM1QVg3xBCqlFKtxNCqgFgNyDdBuBEbr8aAFxb+JTSIgBFQGrJud+H8dkGz646RVplixd+ZYyXNHic8pi0bbDrfgwdaciWLm7S4GHyGHe0fPnChOGEyeOmW84LHEMl2xDBTRqMLQpli0qm0eTlA1jeS+1PbvS+9E2g0ndT+OeUckrjmELoXJV0j2MmpbRR+vXDAHZTSocTQgYCqEQpvZMQ0hDARKT6Gicg1TgtpJT+4Td+0LkqXuJgiMojKOPw47HvBghJww8vgZjONrwEEiQNN0QE4iUMN/zkYVIaz99R6CsLL0QEolOeyIojIwwf/CQim2kMfH68mZPcCCGTkGqEVgHwDYD7AcwA8CqAmgC2ALiUUronvf89AK4FcADArZTSt4OC8BNHkDQYQfLQkQYA3DxB7zaQPLxAwixR7m6/RUkWTvzkISMNhps8TGcafx+4O3gnD/zkodvTkBGHiDR4eIGolibGxJENvMQhKg2GnzziJA5G29/aGR+TZ39Hc7M9QEmBqEiDwcsjTtJguMnDVCM0SB6ywuD5cLbe75SoOHLqJDevnoeuNCwpWP+jfNcXtcdifY9KJ0dbq3vhbJpOn21uDY1br0NHFjxtOr6vLQ8RYrvkXDbbYJieqrVYwsaUNLJJLMWhKg0GLw+bbVjiCJueTaI0gBwrVXh21SnCW7eFM7VnsZggqdIAYphx6GYbFoslfGIljnySxmXlNuH4Ss9GHYYlB1k9NPxrzsRCHFt/a4ABm5dHHUbWuKzcpszzfJfHnnU/RR1CTvFM278BCF8esRAHo8/tBZlHnDC5hoOXBiPf5WGa1te731pCh24d9c5EzgZMGozVQ0/LPEwT2+Yok8eYR3xXq5dg6piprtvLbDxWO6Yw4eWxY88NEUaSff5nYZnM8/8211vpOmdj6rSA1tf3xUfPmbmJN2Nnw8dRdVV/o2OawikNJ6uHnoYGd5u5jAAQk5WjZRo0onXHv+a7j6hAvMQByMlj3pTNJbZNOaC3ItEt2/BDVyDzytyGM878P60x3DCxAIzxl+/9/5vIiIRJg8eUPJ7e+nrgPjJSqdG7PDrslvt98CJIGk78BJKoJeci4mD4CcRPGiU+00UibrJwoioPWWkwZOUxr8xtxV6bFocpaQQJww0/ibhJg0dHICLScOImkRq9y5fYZkIesuJguAkkZ8XBYAKRkYUXC4fI3SFLRR6q4mCICMQpDR4TAinzcvEbGR0xvarHnsGoiIPhFEiQNHhkBaIiDTeatbrP8z1deaiKAygpj5wXBwCcVWeK9mfLSoMhIw9dafB4CcRPGk5kJeIUBo+KPHSkwRhZ6SLtMQB/kZiSBuAvDoaKQHSkAaiLI1azKjKYkAYANL9H7XJslx9eGZcfXtlIDDK4zcDISAMAPp33ND6d93TgfmVeHuErDQD4vdtO3/edmJCGSVpf3zfz4DEpDVHeqVwX71SuK7y/rjQA9WnbxIojLkQhEF4estLg8ZNHkDB4fu+2U0ggJqVxx54ZxsZiMHlEIQ0eEXmYkAZDRR6JLVVMZRyqpYoXfAljskQp8Tm/rDY63llbzN+qsP0LDY2PyWOqXHHS8/6zjY4nUqb4wZcwJoXhxvjj+uff9TjiwF9rnnzoxZ7wxGFJBrrSAA5lIF807KY9lilsqWKQUidUKfa6dKXOEUUizwc1S04VWuJDvZrNcM6PX0QdRgYrDkM4pcFIkjySxNON30e56v2iDiPrxEUeVhyWxPF04/czz+MsDxNlCpDKNnjO+fGLyAVixWEAr2yDEUbW0bib2eZokomzPMIkSnlYcWgSJI0wWN0ptVo2H+XBZxs8uSoPZ7bhxKQ8xtTsILyvnVXRINvSYMKwuFOuej/89NXoqMMwRpA0TCEjDIYVhyKy0mDlyq97Zkl/lp8wWNaxYnoD6XF5Vi4eiCcA3FJV/V4pYeOVbfDoyOOwZqkFcRP/9TSuuED+FAiTyEiDZR3vHV1LaH8VUTixpYoCOplGWLMsOmXLysUDM8+f2DnURDjGEZEGQ7ZsOazZ0xlpMCb+61KpMZwUPPKb8rGqmYZf2TKmZofMwwSJWTna+NmXSmyrMqy01ueqrBo1UZ7IZB2y5Yls5sFLg8dE5mFq5aiMNHjcMg+nIIKQzTycwmgy/UGp4wH9EoXPPGRFcfT+U5JzdmylY6vRDt2uUjo2m/Iw1dMQEYdOP0NEHl7C4NGVhwlxqEqDwctDVhoMN3mIZhSy4jDV1+ja8Eal40TFkfhSZdegX6MOQZqgckW3CRpUtohII1dgZYiqNIDiZUvBI79JlSHLu5lZyyGDqjRkSLw4AHV5VD3hLpw/5p9C+5qeQfGSh6mZEzd5rFw8UEoaOv2ORZ9uwdBT3848VNDNNkwiKwweUXlkaxbFBDkzq7Jr0K+BZUvVE+5y3c7kMbPPva7vhzXtWrpS52Jli+np1sbdVmPF9AZaGcYTO4dKlSyLPnW/GjiTx91rOinHYokPie9xOGHy8JKELO8OeczIOH6QxRRL+lwXytiT7h5vbCynQLwkIYKfQExmGvuO05umZvRaP0F7DL9+h8lsQ6dUyZseh5OqJ9xlTBrZgCxOifu0Mc9HHEkwfOmiIw0AWiVMrpGkEoWRc+IwTft7bg1tbCYNRhLksejTLdrS4HEKJE59DdO49TqSKA0gh3ocScMpDcZpY54PrWyJM/tnbgUAXDOsEC8M2hBxNOHQ5fKbAMeF6w7gaxz+t/OjCUgDm3HEkLhmHmds7hnKuLc/3CiUcU2xr+k+PHNZV60xulx+k+d7Bx6eiQMPz9QanzFoaC+ccuU+I2P5ESgOQsg4QshOQshKbttgQshXhJBl6Udn7r1BhJCNhJB1hBAz61sjxnS54pVt8MRRHnM/Gm58TDdpXDOs0PjnqLKvafh/hAwdeQwa2guDhvbKvA5bHiIZx4sAOrpsH0UpbZp+zAIAQkgDAN0BNEwf8xQhJF53kFbk/qWNcP9S/f8zikiDoSuP8aW+xO+PnKk1BuPn8dsAmJWHX6ZxzbDCSAWyr+m+EtLQzTpEUJEHLwyeMOURKA5K6XwAewTH6wJgMqV0P6X0cwAbATTXiC926AhERhoMFXmML/Ulxpf6MvNaVx5MGgwT8hAtT3Tk0WTXU2ixSv7O9X5ZRtzk4SUNRljyEFrHQQipBWAmpbRR+vVgAFcD+AHAYgC3U0q/I4Q8CeBTSukr6f2eB/A2pdR3ZZOJdRzvbOySeX7bx021xnLS8uL3lI99oFmqwlORBo9Xw5QXhAhH3D5PaD+nLNw4u7X6wjKVvoZf07TJrqeExljQ0PvObbJlyY2vBt9/xa+3IYpb8zRIGE7WTigrtF/Y6zjGAqgLoCmA7QAeSW8nLvu6/sUQQnoTQhYTQhbv/1XPirw0AGBUq2Va45nk/qWNtKUBFM88WEYhKw1ALPsQkUYUuGUfTXY9JSwNAGixqm+JLMStLBEhG9mHE2cvQxTTmYdSxuH1HiFkEABQSoel33sHwGBK6QK/8QsKKtFyZc9B5yvqCAfulIUbJjIPnWyD5x8HVxgZp/+Ac4yM48w8VGShmnHozqIsvUH97nVO5l4pfsc6N/yyDhPZBuO+383cgCwo8wg14yCEVONedgXAZlzeBNCdEFKKEFIbQCGAhaLjzpq4udjDCxFpAPHKPOIGyzx+Hr8tthlGEogi69DBVOYhMh07CcACACcTQrYRQq4DMIIQsoIQ8l8AZwG4DQAopasAvApgNYDZAPpQSv9QDc5NIqLSYFh5eKPbNFVpksZ9zYYKbvIwmW2Y5pQr92kLJHDlKKW0h8tmz1Y/pXQIgCE6QbnB5FGQU3M0llzhmcu6ZsqWOEuD55Qr9wk3TZ3kxcpR1azDVH/Dkj8kRRq6JEocBc3Vm2K2ZAkHmXIlF8sUniRKQ7VsSZQ4dJGRh802xBGRR65LI+nIyiOvxGGxWLyRkUfenVY/qtWywPUdNtuQZ8KcHf47LM9OHLKcPeFO7bUcjM50B2aR442MFQWHt3gfGCe2b15mHEnudzz+aPyktmPf1YH7NGsSv7hNE5Y0TC3+8uPwFnIXUMq7jAMA7my3D2j3ief7i1AWf16dvdOpZTC1chQA3tu7BW2P1RtDRBqMZk3ew9Ll5uKPG0nNOGSlAeRpxiHCogZq89th0uuY7dgz7uXMQ4X39m7JPABg/oY2yvHISIMRx8zj7Al3Ghmn3nOT0O/ZUUbGyhYq0gDyMOO4s514JrGoQTwyj17HbHfdvmfcy6h0rdgJT0wUbszf0AZtCz9Uik2FuGUeWx/8J+oDWL/2SKXj6z03qdhrXh6jDZ5XYxJVYTBicXsEdpKbHzOblIyzS6kBUp8jIw0nTCDvLVe/uOwnjUve/9YLL1mIsLSb+lJyUYGoZBs8suJo9UDxa1Q8sW2u1udvfVDsRlxO1q89soQoZJARSeMa3rd2uOjz+4XHkZHElzfckZx7xwaJw00aPKIC0REHAHw/6RSt4wF/eejIgqEjDUaQPHSlwQiSh1MWTlTloSoNADh4xzTlY3mCBOInDZ4ggchmFjkjjiBpuOEmEl1pAGbEAZSUhwlh8JiQB+AuEFPSYJTppnfvX1l56EgDMCcOHqdERKXhxSXHq/etRMUR6x6HijQA4I39j5bY9jH+TysWU9JIEs7eR9ykAQC31DhbWB660giLfs+OwugbbtMWRjaxsyoR0HLFocskms42TMNmXUxLwyS31Djb9/2tD/4zttJgJEkaQIzFoZptJIWWK66KvTQYcZYG45YaZ7sKJO7CCIOpO8KfIYulOHJdGklj4VU3GB/TRJniBi+PfJRGtohFj6Nyxe9xeZfZmdevpn9eti4n7ueUWPbefeisVyaP5i89G1U4Qtzw1C+4AS3Tr2ah4vLOvvvnKlN3fKjVJA0ilhkH49WT34k6BIuDMLIPU9zw1C8ltn3XZFYEkeQ+sRaHJTr4bMPJwqtuiJ1A3KTBsPIwT+zFYbOO7OMnjTjiJw2GlYdZYi+OOBDGGo5Fp76PybuSfeVl1awjrMZoEFYe5kiEOHIt61h06qFlwHGTh2y2EXXJIpJtxJ1p5UZiz8PfRx2GFIkQB6Anjy3jPsSJl15pMBp1eGnEDdUSRVYeZ55SU+lzcpFp5UZmnpuWR5jrORIjDkBeHlvGfYgt4w59eSrymDx6ON7+9mrp49zwksbkXc2NZB57jxuGvccNQ52PO6LOxx3ljtXsa4g2TDveXRVAMuVh+jwVXhqMpGQesTjJrWqVAnp5l6OE9nVb28HLQZStr01w3T55tPgfUKdjXxTaTyXL6F5F+M6Z2HvcMOnxN7c6tG4mjGYov96DyUKFeWu9ryOiU6bIrO/YOfvaEtuqvHee0ue6ycKNSn+roDS+GzLrORJ1dqyMOICUPFRk4cXW1yZICcMNL4nolCYi8lCRBs/St/6idbwXTBw60vBi3totRnobXvJwE4UbMvIQFQbDikMAWXGccZfZy/o93KaxkXGc8jDRz/CSh64wGGGJAwD+Xv2NUMat32+dsbF4eYgKg0dUHrLiAMzI4/XXXy+xbcJNJc8eZ+TEafVJg/VCOh37orEm6ORdzUvIw5Q0LClUhMHYdc5bmedeElGRBpDqd6jIw00WPFeOPXS9Gj+J+JGo5mhSMD1zwjdPrTTMoiMNEVSlwZBtlgZJwwkvERlsxpEgrDTiza5z3kKV987TloUssrJwwstjKO4QOsZmHBaLQZImDVVsxmEJhS4rnsbyFdvRpGO1qEPJGgXfzEa72qnn739+krFxxw3bCQC4dlBqhioMWcx9/4PUk2PE9rcZR57T7Lx/hzr+8tnJuMqZadrV/tLIODMO7M88HzdsZ0YiJslIQ4LEicP0VCwAFE3tYnzMfKXLiqfRZcXTxbbFWR439X7MyDgF38wusU1XHrw0wkJFGkACxTHzrX5Gxyv85l4AVh66uAmDJ47yuP671gD05eEmDYaqPPyk8d3aVkpjmiSwx0EIORHAeADHAzgIoIhS+jghpBKAKQBqAfgCwGWU0u/SxwwCcB2APwD0o5Rqnd66q9HAYq9nvtUP5583WmdIAIekke/c89Che84Mkcjo/EThxvLZ0fc8mCyc3NT7MYwtujW7wXiQjUwDUM82AIGVo4SQagCqUUqXEkKOBrAEwEUArgawh1I6nBAyEEBFSuldhJAGACYBaA7gBADvAahPKf3D6zPcVo46ZeGGqjz8hNH7Ev3VjlXa/lV7DDcqfDwSLXp/bHTMTm39T67zE4msOBg68phYZkbm+eAbTpU+3kscDFl5+GUbPKLNUhlpVDxF/XfBSxr/OeZsoZWjgaUKpXQ7pXRp+vmPANYAqA6gCwB2S7KXkJIJ0tsnU0r3U0o/B7ARKYkIsavRQCFpAPJlS+E39wZmGXEtWSp8nJrmW1BkLk0NkgaQykbYg0dVGoBa2TKxzIxi0gCAwc+uweBn1wiPESQNWUSlIcKMA/sTkWkwpM5VIYTUAjAfQCMAWyilx3DvfUcprUgIeRLAp5TSV9LbnwfwNqV0qte4VasU0LPveVDtXwCxzEO2LFHNPKoUrU09eWWo0vFuMGk40ck+RKThxsrO45U/k0c063DKwg+3DERFFiJZh6w0gjIOFWmoZhx+4hDNOITFQQgpB2AegCGU0umEkL0e4hgDYIFDHLMopdMc4/UG0BsAylascNqF990pFIcfvEBM9C9k5JERhheSIvGShRsyAlEVBsOUOBhOgciIwottl5i7KplTIjpZhlMeuhmGrDhEMg1jpQoAEEKOADANwARK6fT05m/S/Q/WB2ETzNsAnMgdXgPA184xKaVFlNI/UUr/VOoo8TNj/WCli6mmp9GypefdxR8+yEgDECtfOrVtri2NMGBli1spokqNqebu/cLPuJguTbKJifKEJ1AchBAC4HkAayil/Kl0bwJgN0G9CsAb3PbuhJBShJDaAAoBiF+VRhPTMyUi8gjMNtzwkIisNBh+8oijMHhMCSMsTK31YFOzSZcGIJZxtALQC0A7Qsiy9KMzgOEA2hNCNgBon34NSukqpG7GthrAbAB9/GZUkoCXPKoUrVWThpO0PFSlwVhQ1KqEQOIujaRgKtswKY0o13MEruOglH4EgHi87XqbcErpEABDNOKKHUVTuxTreRgRBoeuNHiYPExP3VqSRxjZBhCTlaM1D1aMOgQpTEsjLGy2kfv4ZR1hSQOIiTgA4PFfLoo6hECKpnZJjDQsZpl28wVRhyBFmNIAYnZaPZNH/5g3yyz5Q4fdf8s8Z/K4+Kl/RRVOIGELgxGbjMOSv/Trau7aFabosPtvxaTBM+3mC2KZgWRLGkBMxaFTtrS+ohRaX1HKXDAcdSedjinvTA/eMSZcuu7nqEMI5OzxOwCYl4fOWg4vYTiJizy+W9sqq9IAYioOQL/nEZY8ACRCHmvu2xt1CIEwaTDikHmISoMhK483tps9F+qDJ27GB0/cjO+nfY3vp5VYZxkasbivStPqjeicPq+6vifS75j26LGe7300UX/evO6k0123X96hm/bYjAo3VTI2llMar51sZmUuABQ9PhMt39mjPY5TGjyjXzdz9SyVpeey4uAJ6n2cNP/GEtu6VFM7J+qDJ272fb/CxScojWv8XJUwURWHnzB4dOThJQ2GKXmYEkdQpqEqkaLHZxZ7rSoPP2EwTIij+yeNir0e+WgL3/11hMHjJQ83aYjgJpYgaTiRkYjRc1WixK1kmfboscLSANTLliBpAPEqW0TKE5m+R9HjMzMPJ590kBediDQA/ZLFKQ0AuGPAAs/9TUkDcC9dVKUBHCpt2E9ZaYRF7DMORv8yM6Rk4YVo9iEiDR6dzGPHjkKc/MBu5eNV+hnOzMNNDkHIZB2i0uCRyTzcZOGGM/MwKQ2eAd3VSgU/yk+trnysaNaRM6UKo8ZYs8unvQQiKww3giSyY0dh4BiiIlFtgrafbbYL75SIiii8cApEVBJ+rLgnvEwxDGkAeuJww00mOVOqhIVb+WJCGoB7+bJjR2HmIcK6+ysH7qMzc/Jux7OUj3WDL11MSgNIlS7dP2mUeZig8RBzje2kojMTk7fiAIrLw5Q0GFPemZ4RiKgsnKy7v7KnQOI43fpJh0rGpWFJYTrb4FGRR16LA0jJw7Q0eFSlweOURxylYUk2svLIe3EkBSaPOEvjyk/GYXOfT6MOw6KIjDwSIQ7TjVGeMx4aia6te4Y2vkniLg2GlUfukwhxhMUZDx26eE5S5BFHeGmEQd3DB6DFKHMXILboE6vT6t0II9vghcHTtXVPvP7RK8Y/L1fxEwbLOuqMOUN5/LqHDyj2msljwW1mLkY8YvV63NmgvpGxAOC1e/sCAFqkL3u7YGN8FgcGMaJsQwDAhYL7513G4SUNhs084oFTGqbp1XU4gJQ8RqxeH8pntKhnZsq36LA/MPKyLUbGcoNJQ4ZYiyPM3oYfVh7xRrdsYdLg0RUIyzacmJIHgFDkoSINIObiME1QtsETR3mMmnNF1CFkEO1rqDRKRbKNsHoeTCAyIvGSBkNVHkWH/YGiw4rfICDMzEOGvBJHknlkdmqtiUl5qK4elW2GhjXLko2GaZBEgqTBaFGvm7BA3IQRBqrZBpBH4pDJNhhxzDqAlDxMCETlfBXVGRRRecj2NmTl4VamiMLL47V7+wpLgydIHiLCMJF16EgDiKk4Gh65CQ2P3IQK/Y9Hhf7Ha4+nIg2GrjzOnaz3Hwg4lG040RWI6fNVggiSR9gNUROMWL1eSRg8bvKQzTKiLllidXZswyM3Be77/eNy50LoSINHdpp2yJriN2ee032V8md7icPJbedOFNqvfWXH9zxB7JfQxHoNr+lZXWmITNHqZBs8tSeODt5Jgms2v6Z87B2v1pQ+xi/buLBHs+ScVn/EiYX0mP6PBu8oyPeP7zAmDB4/eThFEYSISESF4cWsHg/KHeAhkDAWeDGBmM4y3ARiShiMOIkDEJeHSHkiKo5Yliq6mChvZJCVBpAqYfzKmKxLAwCurJl68JtCXBUaRmliV5i6o9vTcJKT4ggLvt8xZM34zEMHN4HoSgMAOk+6T/3gtECSJg03TGcbcSSo32FaGoAVhzRdW/fUloUbTB4mpAEoZhw5Ass68kEajGw3S3NWHGu+Nt/jCBtT0rDklzQY2ZRHzoojn8nnbCPfccojjDIFsOKwWHIOJo+wpAFYcVgsOUmY0gCsOCw+TBgzOOoQlKg3v27UIeQ8VhwWiySfX9Ev6hAiJ6fFEcbMyp+ePdf4mJZkYXrlaBKJxZJzQsi3AH4GsCvqWCSoAhtvmNh4w8Ur3pMopYH3Wo2FOACAELJYZI18XLDxhouNN1x0483pUsVisYSDFYfFYpEmTuIoijoASWy84WLjDReteGPT47BYLMkhThmHxWJJCJGLgxDSkRCyjhCykRAyMOp43CCEfEEIWUEIWUYIWZzeVokQ8i4hZEP6Z8UI4xtHCNlJCFnJbfOMjxAyKP19ryOEdIhJvIMJIV+lv+NlhJDOMYr3RELIB4SQNYSQVYSQ/untsfyOfeI19x1TSiN7ACgAsAlAHQBHAlgOoEGUMXnE+QWAKo5tIwAMTD8fCOChCONrC6AZgJVB8QFokP6eSwGonf7+C2IQ72AAd7jsG4d4qwFoln5+NID16bhi+R37xGvsO44642gOYCOldDOl9DcAkwF0iTgmUboAeCn9/CUAF0UVCKV0PoA9js1e8XUBMJlSup9S+jmAjUj9d8gaHvF6EYd4t1NKl6af/whgDYDqiOl37BOvF9LxRi2O6gC2cq+3wf8fGBUUwBxCyBJCSO/0tuMopduB1H8oAFUji84dr/ji/J33JYT8N13KsLQ/VvESQmoB+F8A/0ECvmNHvICh7zhqcRCXbXGc5mlFKW0GoBOAPoSQtlEHpEFcv/OxAOoCaApgO4BH0ttjEy8hpByAaQBupZT+4Lery7asx+wSr7HvOGpxbANwIve6BoCvI4rFE0rp1+mfOwG8jlQa9w0hpBoApH/ujC5CV7zii+V3Tin9hlL6B6X0IIBncShVjkW8hJAjkPojnEApnZ7eHNvv2C1ek99x1OJYBKCQEFKbEHIkgO4A3ow4pmIQQo4ihBzNngM4F8BKpOK8Kr3bVQDeiCZCT7ziexNAd0JIKUJIbQCFABZGEF8x2B9gmq5IfcdADOIlhBAAzwNYQynlbwAUy+/YK16j33E2u9MeHeDOSHV9NwG4J+p4XOKrg1THeTmAVSxGAJUBzAWwIf2zUoQxTkIq9fwdqf97XOcXH4B70t/3OgCdYhLvywBWAPhv+he5WozibY1U6v5fAMvSj85x/Y594jX2HduVoxaLRZqoSxWLxZJArDgsFos0VhwWi0UaKw6LxSKNFYfFYpHGisNisUhjxWGxWKSx4rBYLNL8PwxK431w/FSAAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Generate artificial cells as test data\n", "tissue = cle.artificial_tissue_2d()\n", "\n", "cle.imshow(tissue, labels=True)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAmjklEQVR4nO1dX8ymRXX/HRHY4J+URSALrAUNa4IGqN3ghY3RNrjKzeJFE7xoMJrsjSbatKlLvChpQ0ptam9sm2Ag0mohBrWSVl2BaIgJdgWzLiDdZRWr627YWjSSkmDF04vvfWR2mJln/pzzPOd53/klm+/73u+ZM+eZOfObM+ec2Y+YGR0dHR0leNncCnR0dCwPnTg6OjqK0Ymjo6OjGJ04Ojo6itGJo6OjoxidODo6OoqhRhxE9C4iOkJEx4hov1Y/HR0d04M06jiI6AwARwFcC+A4gG8DeC8zf0+8s46Ojsmh5XFcA+AYM/+AmX8J4G4Ae5X66ujomBgvV5J7MYAfOz8fB/CW2MNn0dm8Da/IErzryudw9PA5bdp58gCIyJSUNYfsAa19SM7R0sZUWmZqLDX0fxY/+ykznz/2nBZxUOCz085ERLQPwD4A2IZz8Bb6gyzBBw4cwp6Lrm7V7zR5AERkSsqaSvaBE4cAnIc9F129+r69D8k5WtqYitqTMzdBPLr1jKT+9/M9/5XznNZR5TiAnc7PlwA44T7AzLcx825m3n0mzlZSoyMHg+ENXwcC6bAPl/CnhBZxfBvA5UR0GRGdBeAGAPcq9WUOS1l40rtVhxykvL+aPnOgclRh5l8R0YcAHABwBoA7mPnxWnnuC0m61JLwF6GkjtJElNJtGF9rpGJNnzG06FtjO3L9Hctqo1bHwcxfZuZdzPx6Zr6lRsZgwANCLrWF3T00aVI6Sr6fq0vKyOYmu0FPX1cLc52Cr1/N3E+9Kdb2Z7Zy1CeM2ML0n63to3aBp9o0BxmdSW1dPCECTkHi7OwTQA2WEH/xSc63V22da8emhaRUCsBK8Wrazm5WpfSFal27FHJklfSb60r6erUef0oJw29b02boK9eQ/XdM9Sm9I6fGu7R9rG2pndToIdXP/XzPI8y8e6y9OeKoHTipyQkZe+q5UnKLydTuV9sQx54v2Q3H+hNLG0d0KiV4qedb3yu3feq5xRHHz0+e95ufNQeudHBd+C5hy+4Ug+Su1BpQrBmrFoIqaS81BzECb9lUcvv2209BHGPPLIo4dl+1jQ8e2Cr7kIwLhD6v6aPVYDTkjRmAxK5cYoitfeX2J/H8gNKxk3rXkB7SXlSt7osiDj/GIYHWs+uY3ClTbSVyNKLyqb4sjGuLF9EiVwKSxBuSO5BSruxc4jCbVZFCKCOzDpCIf5T25WehrIxrLJORm34ek6sFzXSzK1PjHdaeOCym7zSgnf+3VjsTQyibY4XgXPjzZTndHMJaEocmk1tD665aAle+xcXoIuSFtMiSRGy+lmSva0ccS2fyFmzCOy4dYyQvZa/atrB2xBGCVrFQbXvJHczfWSWqNZeM4b0tvn/u/NeSx5Rzv5bEEZqYkkrGEKwdf3JSunPrODWszZGL0k0j9x1iZOFuKBrjoPUf+ZhFTZQ5dPwZJsvaWT+UWUg9k4OQDK13r5W7pDnKRWzRt1S8SmFt6jgkym39ZwZIVqJKTK5ULcmAMU8lVSiloUPJ+2kXwvn1EANKK0alC8Zq2ue026g6jpKBSbmAMZevFVNlPXIRq3kYvo89F0KNG+yOR6iPEjfdfT6ElmNLLNBeIk+6NqO2HkVal7XwOGoWZmwnaZET+jyGuStPQzIHSLz/WD+t3oSkV1LTpsRzlagSnmoDW5THsevK56rZsHZyQixcw+ghwogFqkLPWQjgtdRnlHoHJX3EZGt7cBKejDUvE5D1PEwQx4C5FpJEoZCve8gNj7nmuSSiZYytcksWkRQxlcjRyrL4KfABFkljgNRYmMiqHD18zmlHhtxMgMUJKt1NASzG6FKIZTEk3seVXSsrNy1ZkzaNxYfWGaY8jtSOZJ3VW7yWkCcylfclOY7+bqYhewrUeEUtx72pIeF1mPA4fEwRcLKK2tSoFcQ8DwsYdFuK96oJlzxq3tmUx5EDrfOqVUyxk2nItbwApw62WkfNWloccYSOKx0dudA8TmlA28ZjAd4xLI44BtS+8JKxKe+pDZ88lkAa1uZ+kcQRq+RrrQWxjE07ovnQem9J0pCco1AKW2MMalPlJoOjJZAa1JJKTzeg5Lq8GjuXK9dy4FELoXlZQtqzZY5CnpCfkvZ/34K1D46mJmNgzBoCyS1bdicu5OFI7wipehbrRzTJXRdIF89ZRY1uqeNTTpVyST8t5LaYuyq559GSOwQ58mJtQjKkz8wpeRJpWk0vaYDkrjvWV25/U1XglvZT+3xum5zCyrX68wg1Ex0yqKkWttSiLmk79nxqZ9Icj5L3kCSckIzYwq7tM7f92HtJ2OWYd5orf+OJQ7J9rdwchq+RO9bO7StGGKHFI2G4tQtD8uye4x2Gns3tu2ZeY200bLPW7oBOHKLtW2WWTKTUu+b257dp9erGnp2yarOUFEo3ghZPslVOST8l8teGOCTduLmII9TOhUaMpOV4V0OGtZ6UZi1F6xFX6kiT6kcrK1Qrf1H/H8cmIZYv1zakHOSktt1IfizTkSvfGmn4z4fSn1bTv1Nj8XUcS0Vo17WAUK2IpHvt1yNYXIhaWbKpMIXepj2OpU7c3Ggdt1StiMSuO7S35qZ35KOJOIjoh0T0KBEdIqKHV59tJ6L7iOjJ1ddzZVRdf1gydHdxlx5JcmHJ0woh5+i2qZDwON7BzFc7AZX9AB5g5ssBPLD6uRh911hfjB2DNglLfX+No8peAHeuvr8TwPWlAqySxuC6u5O91ImfCzlZi3XHWCHaEtBKHAzga0T0CBHtW312ITOfBIDV1wtqBGvntlvbuq77Eid+DoQ2hE0bQ38MrG2OuWjNqryVmU8Q0QUA7iOi/8xtuCKafQDw2otfjgMPHzrt9633KGJVlNq1CpYgcRdF6j5LatxDmZxNglYFqSaaPA5mPrH6egrAFwFcA+BpItoBAKuvpyJtb2Pm3cy8+/zzzgjKr9mFYpkA92tKdug44re1Cq00Z6s3kKPTJngeuaX3IfuT7EMC1R4HEb0CwMuY+dnV9+8E8BcA7gVwI4BbV1+/NCZr+PMIPlp2oVgbP9WYIpOcflP3QaaAVgmzH8CsnYcSQ5YcyyUQvYsQcZaSgN/Wah3HhQC+SUTfBXAQwL8z81exRRjXEtGTAK5d/VyFkl3INe5cI3Vl11ZDDm2mJg+tGosQUnUdKdTsfnMTcQxT6VSb+vbtd/hMS2/zd1UGjA2AhbssrXcjBqRK0lPPSOiU606P6dI6pjV6hCDZf4t31+IBtL5nie2szSU3H6F0XiumJo7URI4ZrOZxocU19ttJBWZD+sQWUuuxc6zfkv7956VtdazPnLah9mtPHBrBP63sQWrHyNmtNRdh6TOpdkNbjYUSOsKUZsZKnpcgzxp5JX1KeTCunFzi6JfchOGnfWPPlMiwHujTPFOHxrM0ZuLqJk2c/gaxhJSyhM6mL7ktFf4kuAGv3AmSNr6xQLMESWkvmJbg75TvHwq4W0VtGrwThwJCEW4LiBmJZQOX1M33PjQW+BLrUWp07sRhGNqpVferRdLQwBQxiE0Yy8XGOKyeJZewECWDrkuET56aY2B9fF86Bsey2i3S41iiO9jRhj7Xeqght0URh7tDSpLHUnaFTUMoLb00WNa5xeNcBHFols5KQtP1XcIYaNfDLAXWCU/ClkwTx1j+Xmpi/EBhiwwNhPLu1oxREiHS0FyEkjJ93a0Rn1TltdnK0ZIXlKh4dFFbNVgjI6eP3CrUkr41vINSHVIyNO6/xGRJLaaUflKBaM1SemDhJec1BlJ7T8RvU1Ka3SIjFymDGyOR2N0RiYUe0iEmv3ROxsZeo/y+hTxy2kqXibuQ3CgWSRytzB8yqpoBjxlyqYGPPTMGaU8qBzXeVso7iMmv3TlbiWNsTFvuqmh4xmPjVLLR5fS7KOLYfdU2Pnhg5+yu3JjMWjKraTtXoDG2EKTGVYpQpTaW1LMDUpuRhk3UyI/JLB2rxV1y07qb0YpWQ/fP0NYxFk+xEuwrXRClc+DOW2ssLCZb8vjoypwCJrIqRw+fIyrPinEPmJM8pDyXPRddbW5cc1E7BrF3lhiLMcKokR/TVQMmiGNTkEseVmoYNPqXkFmSmpUYSze1Kjkm7nsMspdC8J04JsKcxUASi2ZJx60BVgh4HWEmxrEJCJ1tfXRjL0NJSrpDDt3jmBhL3sGtoI/d/OgeRwJ+daH07iUdrQ/Jl4opjHlKc2DqCtmOF7FWHoekYYfuHEjsdK7cUCBLqm5D4v6NZXQiyIPW/K8VcUhBa5eKyXUJpHaiQ5ffpI5Fbvt1JaI5ELsvIynX/VkSnTgcuAuv9n5ISnZKbi1CRT+SKTnfQ3I/mwOlXqUFnUMYuy8jJVfr/U0Qx64rn2uWITkwkos7h4z8fkvrPYa2qaNPzfiE9HbltYx5Tfva/qySh4vWcU3Zmcb7m7qrMmDKOyFue8lgW8vdDum7DC33ZSTvU/jtYtBKVdfEwDQC2Kn3KLWbkudz5mxRl9xeTdv55yfPi/6+9NhQa8gSF4padBmT7/fRqm9Ou1zDrSWw3DqMqS/HtRBcjtwSQpAk7zGiWRxxhP4/jhLE3LFc0qndiUr7rZE99nmLzLE2JYSQkl9DMCG0HiNLyDDVX+141rRJbVAt3leo/aKJIwYrrmZItmYGRmq3LZFRm9qOZQokFryWHIlj5Vi7llIByU0pJnuQmUscJoKjmtCqldCGpIsekpdCSxBtGO9QNsYC/CBk66L0CSg0di2k4fchPZa1gdO1J44BmsarvTCmXHiSaWNLhOEiVkdTq6/fdml1LzXk0UvOOzYakuQWS6laJdABNeS2MR4HsIx8vibG3l+rSG2TEDuqWEXtkXKjiGNJ0DK2JZKn7/pb1d1fhEsh4Bo9R4mDiO4golNE9Jjz2XYiuo+Inlx9Pdf53U1EdIyIjhDRnmKNOl5igNILJUQeVr2N2I5olTysjV8KLUHbnBjHpwF8EsA/OZ/tB/AAM99KRPtXP3+UiK4AcAOANwK4CMD9RLSLmV+o0i6A1gi1K2eApcmOlXlLFZf5MpdAGC5i5FFTACX1zqkx1CD+VkjM+ajHwcwPAnjG+3gvgDtX398J4Hrn87uZ+XlmfgrAMQDXVGvnQCJSnSqsKpWp4T6PGWBOpWppXxaRa9g5qdCcflqQo6sV8pDcKGpjHBcy80kAWH29YPX5xQB+7Dx3fPVZE0LuautEhM6hucaXCi5Kk5oPifNzbDyljHuQ0yIz9/1KjjCDtyplR1N6ay36amR5pNOxFPgsWJpKRPsA7AOAbSj78wium91SsON/5hp87PmUsYTax57NlZlCTXXm2FGoZTH4utTIrD2K5maM/Da1dlQzTlLVo61jJIGsknMiuhTAvzHzm1Y/HwHwdmY+SUQ7AHyDmd9ARDcBADP/1eq5AwBuZuaHUvJTJec5C1W6jDplZDn9xYw31k6yHHmsn9TvQ/Jqxit1pMrpu9XAh/HMfRcJ8s3Vq7RdjocxNqYl/YneVQkQx98A+B8nOLqdmf+MiN4I4F+wFde4CMADAC4fC47WEkfO793nWnfRnH5KZUgFe2P91OpfuthzZEtuAjn9uCghey1yq9noYjqNvWPt0USMOIjoLgBvB/AaAE8D+HMA/wrgcwBeC+BHAP6QmZ9ZPf8xAO8H8CsAH2Hmr4wpESOO0h2+lVxq5ZfK8SF9Rpbsp/bYNiYvZORTZDla2k5lR7UeUAilui7+dmwtO8854aXyJGXG+tDawQdIeGC1clLypYLHknLHZPuoPQrVtB2wuD863Qo3OKa5m0nASnouF9Lj6QcyLc0NEA+0SugZCshKHVentCuTJedSHkJHRwta0945WKqtmiQOoJ40QnUJ1na0js2GdY8rB+aIQ8Jt82sUOjqsYqn2aYo4luq2dXRsGkwQx64rn1u021YKvwqwo0Ma2nZlgjg2CaF7IpsKifsiHS/C33w1x9VEOvbo4XPEUmB+e0s1HGP3RNzPNg3SRXaSFbktd1qmgj9+ofUkqbsJ4nDhv7D0/YWS51val8Anj006tvloIVF33LQWunQxnARS9qKVZTRbOTqgxJBShpIjp/VCUQolk9ayS2h7Lxrl97E7Fj5K75BI7bY+IeUg176kPNma0vRQm8WXnPuIvXQpseRCunS9ZvdrvbPgGrpG6blGWfdYn6E2OdcNYm3HUEv4Y5CaG2m7WjviGCDlKtayf+3FLqkq2NKLfNJ3IFplDfKk7rjk6lFz8U/istxYX61eUattbQxxDJBi6xoZJZOt7drn9lFDIGOkVUuipW1K9SppG5Mj7amVzE1unxoey1oTh3amo6TdgFJPoBaxlFvt4smJCeUa/dQkqnX80pgz6diW1lh24phI1thka2RIpGINvu6tO1+orVawVntcpeVKemWam9HGXauvgcQEjGUENFN0rbL9VF2tvKkIUxslWZNc1MiLka7/eym4/Z2xI69Nrxzt6FBCK3nuuejql7S1shF14hCCayTuhEvvYBpYkq6bDitz1IlDELHdwcpkrwtCVbabgFAV6FzoxCEAjdLmjjSse3XSMZ5BTm6VrTY6cTRibPI0CcXigtGGT9KWyUMCoU1JizRLZHXiaMDUmQPXRd1ED2csy7Bu5DHVplRz9OnEUYmSBSyR5gvVGGikD5eKdSOPXPtqsQGfMEqIaKPrOGpRs+u3VCSmYigSlY5+sNGiN5Pzji55tFTStta1tKJ0PnNsQLrGaDEexzCZc+8oLYu0ZlfMMeCW3dZPI9fK0URrWb3/uf9vgB+ArB3PlvGrta+YziF9QvUhpTBRcr77qm188MDOqratu+xcF7RKLphpXHpKyZ/yPbXkhNqVVmKm+s6JP8xhY63VpmtXch6bvFqXstYlncJ9rT135r7TmHwrsZNW7y6229bIyPVC5iZbf+607NWEx1FzrX5Ay2CX7tBTXNDS3O1LDarlnVsWrC9jTh1Cckpk1XiBc5JPrsexmBhHDCXn0djZduxcKr37ps6j7u+l5LdEz0vfPTaWpWd/ybSzxJm+VpYGwVjA4okDSJNHKqjqG0GIWLQmNLS4JfsJjUmt0Y8t+Fig0f3nPrtU1M6NlaOfJBYT4xhDy5m2JjgmAV9nDXJqPWL55BEjpFw5OW2tpoRbEXuvJb7vWngcIUi6p1pILUZrcHUMpXBrZPnyhp+tz5skLM95CmtHHNLnWUBncn3SWAJ5uJCIw6QIZN3gE6+mpzkF1o44loCYwSyFPCQNfYmLpgZjwXfrc+5jbWIc64La+hJtaO2OsUItK+8uuaDHYmmS6W/tMRz1OIjoDiI6RUSPOZ/dTEQ/IaJDq3/XOb+7iYiOEdERItqjpfhUCJ3tW5CzAKU8D80jlqSsVHHfXAjFcrT0cY9tNX3EMmea45dzVPk0gHcFPv87Zr569e/LAEBEVwC4AcAbV23+gYjOkFLWAloMumTXbjUkP+DYAo0gbmos5k7fprI+tfOhNeexQLU2eYweVZj5QSK6NFPeXgB3M/PzAJ4iomMArgHwUL2KduDuOqkJkaoIzd3lYs/43lKp6xpaQJJHqVyvq7SvWGHdWJ/uszFC0z5GlqTox57R9JKySs5XxPFvzPym1c83A3gfgF8AeBjAnzDzz4jokwC+xcyfWT13O4CvMPM9KfktJefA6QMofRavkRcrNpPSoaZWJSSn5bnW83gtieV4AKWLrbQmJaVP6TM5/Uyh6wDRP8gUII4LAfwUAAP4SwA7mPn9RPT3AB7yiOPLzPz5gMx9APYBwGsvfvnvPvXwpU2LMzS4EuTRKktKlzGvolWfUjKamjiGdjGULqIWGb4s7QJCiY0oVxfVuyrM/DQzv8DMvwbwKWwdRwDgOAD3fvwlAE5EZNzGzLuZeff5522FQUI57hDGjFw6uChRedkK//zaWnzlHrlqPJhaN7hlTmLHh9LjX8hepsgWSUBq3ltRRRxE5P69p/cAGDIu9wK4gYjOJqLLAFwO4OCYvKOHz0lG1kOBvrHFs5SaiDkQGpsWMpoaErpKvOfUNiYVU5LIWOWkY+/CVnDzDUR0nIg+AODjRPQoER0G8A4AfwwAzPw4gM8B+B6ArwL4IDO/kKuMaxAhV7qUxTt5xOEGOyXiLyXPLoGcchGyMcvvKZWxysmqvDfw8e2J528BcEu1Rg604hYdHZIIZdss26nEhrqYkvPaiagdpCUYQIcdaGT0tNGSrl0McbQG1TRz2puKElJe2qKqxdLez68bycViiEMK3chlkUMefTxto8Yr75fcOsTQPbrlotQj3yiPo++OHS40sm5LJc9SvTfO46g901nAFHclSlFSxm5Jby1o3mHRwulzeCyrzUZ5HEA+s1okllCtQI0Mlzhb3rPEO9uUmpqlvV+th71xHseAJe2QMWMs0S9236H1PUvLvS2NKyCjkz+2NRfRUrI1xqpVx40ijtIdstagtM7MNdW0Ofd6at6z1qAtpcVTCz6E3IuB0lWkGpc2B9TKXMRfcnMHruWmoMatzBBi925iv0v103rrs5QMctpILYbcflxIFFrVtm9ZdKVtQzrW2Edpv6LX6rWRIo7UJJcsllaXT3PnaCGMHPk1MmJypHa/VjJuXfwWxndAzbuVjF+JvmtBHKU74Bi0d8hcOb4uGq6oBnlIp6pLjwax9qVe2dSeSkrWAN+bbvU0S+UMWPxfq6+J2PttpXXRgKVA4QD3rD7op0UaLfJy4zOSwUop+IRcq6OUJ1aKtUzHDndThn+W4C9Ky/CDfdbGEhgfT4uk4WPQy6K9xmCSOCwbqgTc91rKO1rWMxZAdO3Isv6SXudU9TImjiq7rnwOBw4cOu0zyxO9idA4UknKi3ke3Y50YII4Qv914Lp7HUuATxbW58TXz1qx2VSYotDO5FEF2JwSZasIGZ3l+EyI1LoN6cEscXTMh9RCk0pFSu6EKU+ok4cOTBNHn/TpkXMcGYKNFuYlV1/32Y52mCYOC9A414ci/0tFqf5zva/1OMfS7MA8cazbbmE5RVhKknPPjfVgbQ6WuomYJw6g3kD9IN7cExML4EkFG1vetXYR1gZM55gLK8erAaFNREo/bVJfBHEAeQY9GG8ol186kCHSaZmE3ItLtX3EqiZz5LXu3GM3OH1dfLKxmKXRxtKzQKYvuYXg351wUXsxKCQvJis2qWPPlyzKkja5hDQG6RhOTbl67oJpvTsz1nZsjlvm1JdVo18JSjNYa3E7NoQWwhiTUyprzBimMu7W99eIEUjJ10oNh8a3xCZK51Z6M8hF6XpZe+KwdEvTlTPIkpAbkyG56LWqCzUDl1I6t24eJV5qqb4Sm04OfPm5xLGYGId1aJRmhwKPlrMyS0fpmEoev2KyW9Pdwzu5/yT0M3FXZV0g5WnEZK5D+tEiWsbT3zBK4mU5skuC27n9+Tq77c/Ykadb9zgWgFDkvcMmpEjDbZ8ij9b+avXrHkdHhzCmIHfJOFfNsagTR4cqtIKvViF9XJU+/sT6eFFm/0tuHZmQLDgazswa1ZCW4RKk1Htrp8zdPkrlL5I4tOsOOsoRq/5cAnlILvABre+dqiyVRo3cRRCHX5YsZYQ9S/Eiasq+x8r73c+G5+eGr7PGAh/QmlKVlBnrp9b2R4mDiHYS0deJ6AkiepyIPrz6fDsR3UdET66+nuu0uYmIjhHRESLaU6NY7A6DxMCFDMfaBai5UHt/ZKwGQtp9l5gr6aOFBKbazFrfdbRylIh2ANjBzN8holcBeATA9QDeB+AZZr6ViPYDOJeZP0pEVwC4C8A1AC4CcD+AXcz8QqwPt3I0t0RW476AZrVnKzTOuyVl16HnNErpx9qN6TUmY+yOiMadohy5pfOrdZ1B7A8yMfNJACdX3z9LRE8AuBjAXgBvXz12J4BvAPjo6vO7mfl5AE8R0TFskchDOS8zIKeIxQ3EjSFnYkplTgV3okPeUqtMF6nKQqldvkROyi5yCuNKvaZSe5K2k6k8jdZ+iu6qENGlAB4E8CYAP2Lm33J+9zNmPpeIPgngW8z8mdXntwP4CjPfE5O7+6ptfPDATgD1O1HqnsAAjQtJsb5K+yzRpdX7kHg/zTsypeMYIpfaucgZm1pvV8pz1mwn/icgieiVAD4P4CPM/Asiij4a+Owl7ERE+wDsA4BteOmfR8iFfz71jaZGbo3nMWYYoT5SslLPhUhE2hVOQWJX9OsdWvQLeUlS9hT6rFa2q9fUgXnp/rKyKkR0JrZI47PM/IXVx0+v4h9DHOTU6vPjAHY6zS8BcMKXycy3MfNuZt59Js6u1f83CE1IyyCVBMxShhrTIxaALNU/R88p6gFqkFqkrWg9VsXqUCT0WzppAHlZFQJwO4AnmPkTzq/uBXDj6vsbAXzJ+fwGIjqbiC4DcDmAg2IaJ+BHyaXkxYywJPMwGGKIEHLO6rl6thLRXLCqnzu2khvR1O8q3V+Ox/FWAH8E4PeJ6NDq33UAbgVwLRE9CeDa1c9g5scBfA7A9wB8FcAHUxkV64iRh8RxKLSrtbrYgyyrXkaH3FEPGPestAL9OVmVbyIctwCA4P++w8y3ALglV4ldVz4HPJr79PTwz+BWI+q+MVnLDnVMC826FDOVoxaKb1KQTIVqw7JuHXKQzvyUwARxHD18DgD75DGgL8wOy5hiczNzrV4qCNXRIYFUOtb/fE64RWtTesMmPA4XpZWFISzFc+l4EZbmLBZY9jM/FnWeiszMEccAy+RhyWCWDN/Ypca1RV5O+jpW4zEH5urfJHGUTopbTKUxqZqyNRCryLQEnzQsjWvurp2qn0khVvxXglj7Vrm5MP13VUruoQyoKcdOIea2SrqGmrKkzuUacmIVtZLjkOor9HxN3zk1M7VkPnaNIVZMOCYjhrX4g0yxySwxCgmDyDGGFmjLqV30MSOs1Te3KK2VPErvB0kVy42N/xgJxFATAK19p7UgDiDN1DULoLSN5k4lJaPk/XLJxUfo+RJ9S96xZjzGbKNmUdeglLRq5NUQdm6btSEOQMZNllhcErJj7VzUyqghxRRh5IyT1s5Xc0wd6yd2jJPOREinbVs9ohJ91pI4pHcFjd2nZEdvPQLVuvRSrnlMlpZ8F3N6eHPI1oiDDXBl9r8dm0DIwKUmJhRpD01Uqp+caH1L5FzSqP1MU046s1R+zmc1Mi1kcOaARD3KRhIHEF6cUgsqln0pWUypidXcMWvh6qKllxQZLQ1a891SXrCxxAHoGruE7NDEWiSNjmWjhjw2mjiWAu1gngRCx76O5aDUpswTxxSBLP97y7BMGpt6lFgXlMydeeLQQjd2GYSIvbQ6Mke+/32HPErG18y1+hA0vI2YzJaKyE1FKh3cOp6ptPYU9tAqa0m2dLqu/a/VvwRjk7npaTorSM2ThDcTkp9KfefKdL1XDVvSsMtagjNLHHMxdiePPEjv/DX9t7RzF3iIQEqIZKxIrUVXn5A07LJmHs0ShxZyBskyeVjQKVeHlnHMLZCTQizWNUYiJV5sCRGlKo0l40e15L8xxDH3DtmC0B2LVuOpHYtST7DU2HPnSYvc3UK91LWBksLBkipNv1xfy2Zbx80kcUhH0lt2PAs7vAupc3lNu9rjY+5YlupUO0cl+vsk4o95CYGm9M0t15ewS4kwgKlLbjmDUfqyrYNU2z51kahGh5gbXdNHazuNsWwd55x2El6nxKJL2fkUc5Fqu6jbse5fqwfCKbgcxNpNZSyxc2kLieQYe+o8PPZcrk6SC0Z6nnLaT20LuXIGSBNAbZtFEcfYtfpcSO7yIdkxeaUTmOPqSi0mFzkLy39Oq5ZGo9YhlmbV7mNOebleVG4/ucRhugCsFNopK3eSWkjKXTQaxl7TPqST+ztJaBVISdR4LA2+XYagMdYmg6OWEQpI1kxIqH7A/doyyS1BQ79fzUyUhuwxUl5XjM219Bh04iiAO/ixlF2NTH+hz30O10wDToVNIo05MoBrSRxWU6k52ARDnwqbNJYxm9eqX1pL4thkbNJO23E6fPLQ3Dg7cXR0rBGkj70xdOLo6FhTaHqdnTg6klhinAjoRzVtdOLo6CjAkgPvklhb4tCY4B547Og2sAUTJedE9N8A/hfAT+fWpQCvQddXE11fXcT0/W1mPn+ssQniAAAiejinRt4Kur666PrqolXftT2qdHR06KETR0dHRzEsEcdtcytQiK6vLrq+umjS10yMo6OjYzmw5HF0dHQsBLMTBxG9i4iOENExIto/tz4hENEPiehRIjpERA+vPttORPcR0ZOrr+fOqN8dRHSKiB5zPovqR0Q3rcb7CBHtMaLvzUT0k9UYHyKi6wzpu5OIvk5ETxDR40T04dXnJsc4oa/cGDPzbP8AnAHg+wBeB+AsAN8FcMWcOkX0/CGA13iffRzA/tX3+wH89Yz6vQ3AmwE8NqYfgCtW43w2gMtW43+GAX1vBvCngWct6LsDwJtX378KwNGVXibHOKGv2BjP7XFcA+AYM/+AmX8J4G4Ae2fWKRd7Ady5+v5OANfPpQgzPwjgGe/jmH57AdzNzM8z81PY+mOh10yh54CIvjFY0PckM39n9f2zAJ4AcDGMjnFC3xiK9Z2bOC4G8GPn5+NIv+BcYABfI6JHiGjf6rMLmfkksDVRAC6YTbswYvpZHvMPEdHh1VFmcPtN6UtElwL4HQD/gQWMsacvIDTGcxMHBT6zmOZ5KzO/GcC7AXyQiN42t0INsDrm/wjg9QCuBnASwN+uPjejLxG9EsDnAXyEmX+RejTw2eQ6B/QVG+O5ieM4gJ3Oz5cAODGTLlEw84nV11MAvogtN+5pItoBAKuvp+bTMIiYfibHnJmfZuYXmPnXAD6FF11lE/oS0ZnYWoSfZeYvrD42O8YhfSXHeG7i+DaAy4noMiI6C8ANAO6dWafTQESvIKJXDd8DeCeAx7Cl542rx24E8KV5NIwipt+9AG4gorOJ6DIAlwM4OIN+p2FYgCu8B1tjDBjQl4gIwO0AnmDmTzi/MjnGMX1Fx3jK6HQkAnwdtqK+3wfwsbn1Cej3OmxFnL8L4PFBRwDnAXgAwJOrr9tn1PEubLme/4et3eMDKf0AfGw13kcAvNuIvv8M4FEAh1eGvMOQvr+HLdf9MIBDq3/XWR3jhL5iY9wrRzs6Ooox91Glo6NjgejE0dHRUYxOHB0dHcXoxNHR0VGMThwdHR3F6MTR0dFRjE4cHR0dxejE0dHRUYz/B/uvQrCOSluYAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "membranes = cle.detect_label_edges(tissue)\n", "cle.imshow(membranes)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Analysis and visualization of neighbor count" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAdzUlEQVR4nO2da6wdV3WAv9VAeCOS5oHJozHgIps+wLHSSrdFUBQIkei15RQFySilkfyHqqUCqbFx1EhJMEWCtlIFUqpEpDUlIBvbCeWpqBXiCgi2FUJiN8RAWkysGEooCKTQhNUfdyYZn8w5Z2b23jNrz1mfdHXvnXPOnDV79nyz9p69Z0RVcRzHacOvDR2A4zj54eJwHKc1Lg7HcVrj4nAcpzUuDsdxWuPicBynNcnEISJXiMiDInJcRK5L9T2O4/SPpBjHISJnAN8GLgdOAN8A3q6qR6N/meM4vZMq47gMOK6q31XVXwJ3AMuJvstxnJ55VqL1XgB8v/L/CeD3pr1ZzjxHef4ljVb83Ff+Iigwp56zeAyAkz+5IMr6nvsS30+xKPdNHSePx9lfT/G/h3+kqufOe1sqcUjNstPaRCKyHdgOwPMuhj881GjFa+86EhqbM8FW9j71900H3x9lnWuXfT/FoLpvpnHTW+PsMwA+I//V5G2pmiongIsq/18IPFJ9g6reoqqbVHUTZ84VnJOIyYq5a3nnQJE4OZFKHN8A1onIWhE5E7gauDPRd5liD9uGDsEZCU2yjaFIclUFQESuBP4eOAO4TVVvnvrel2zSJk2VXXc9fTbcx1XhQUakThjb2DNAJM2YVylDmyzrIzZV9rDNdFmmoK00Qpsr64sugGNy6WFV3TTv/cnGcajqZ1X1N1X1FbOk0YRdd+08TRpg28YlsbKPreyNur1N1mWlyVKWYS6ZXLmvQvZX33V7fYd+w2QZR6sgpmQck7KYxtDZR5NK3fWMOVmJQra19VksIOs4vLyh8zbPKk/Lmce08m27z7qKo0vWMSmNphlHqqsqwTSVBqwWdJcDKkbzIuWZsO8sIxaHlzcAp5dNrAPeYrNlXtl2rZ9t2XXXzlby6JJplJjMONpIo0qbnTPvgG9SOdtKo8k6mxzgbbYzVBhts45SGm2olkvTMo0popD1tinfJvsteH81FMc0aTTNOMyJo6s0SprsnBgHfGimMbnOthWmqTz6FEcXaYQQKo9p+zCW4Kcxa9/1IY5ZmUZ24tj1kzdFW9+sHRNywJcVKkbzZD+bg9cxTx4xmidNxNG3MKqkalrWrTdFc6+6D2Otf5Y85jVPshLHyzat0WsPvTP6evdxlcne+BjSgOniSFHB6wQypDCmMU0kofUg1j6r48aV3Vy/tCP6ekuBtOnLGPxyrAUsSiMmdYJI1Qlq5fLsPOr2ueV6cOPK7tN+xyS02T+LUYtj0Uh95aQqD4vZRpU9bHvqxyqTskghj5ArJ7MYrTi2cGDoEHqlr8utu5Z3ZiGN2MSuT9MkkUIeKRilOBZJGkdXLuXoyqXZVDhnvhxy2JejFIdlUneyObZpuo9i7ctUmejoxLFI2UZd5bpxZfdTP44t2u4Ty/twdOJwniZFxVskMVvAqjzMzlXpgvVKHbOZEiPlnTd2wHp59s0WDnTeh0MKIMVcmdFkHNYrecq+ja5Ma9Js4YD58hyKLmVjIWuI3dcxGnFYJ+aBGLsiVgXiwmhG03KyII2SmPIYhThyqezl2cpqvFbjssq88rIkjZJY8jAhjsc4a/Cb8fRNV4mkqoxHltYnWe/Ymbb/LEqjJIY8THWOlvJos2FjOEtWt8FiX4gzm3L/lfvOsjRiYUocJXXZR51MxiCNSSYrYR94thGHLRxg48qxocNoROiVFpPiqKO6kZYnLsWiLgtZhDOZ0x8h8jDRx9EWa/ecTM0WDnjfhpOErv0dWYpjK3u9L8BxItFFHlmKo2SR5OGZwWJyZGl9L315bZ8Fk504JjduP5sXRiAuD/vE3EfVdVm7EJCdOKaxKAJxeSwGdfvZkjyyuariPM2RpfXZXPZbRMqDvus+mnVysDLmJ6uMo0kbbBGyDlitXJ592KbL/mnzmSGnL2STcbTpuJmUh6UULzbVitbmDOfS6Yem2UfI/pis332cPLMRRwj72TxqeZTMasK4KIZllkBi75s+Rh8vhDgWickMxIVhi7Hsnyz6OPp80vqYyLVSLgo5758sxOE4ji1cHI4zQlJ3kJoXhzdTHMceQZ2jIvIw8DPgSeAJVd0kImcDnwQuAR4G3qaqj4WF6TiOJWJkHG9Q1deo6qbi/+uAu1V1HXB38X8nPNtwHJukaKosA7cXf98O3RpbLg3HsUuoOBT4oogcFpHtxbLzVfUkQPH7vLYrdWk4jm1CB4AtqeojInIe8CUR+c+mHyxEsx3gWRe/NDCM2SzCqFHH6ZMgcajqI8XvUyKyH7gMeFRE1qjqSRFZA5ya8tlbgFsAnrdpg1Zf63K380lcFs6i0sdclc5NFRF5gYi8qPwbeBNwP3AncE3xtmuAg12/o8uNVC0/8MhxUtPX7PCQPo7zga+IyDeBe4B/U9XPAx8ALheRh4DLi/87s2gPanKcrvR5S4nOTRVV/S7wuzXL/wd4Y0hQk1TlUW2+jCWzqM6YzHn+wtAsYjkOdf+Z7GbH7uOqUT1XZXKadc4zJvtm1j0uytfGWpZD37AqO3GMgXk3dXF5zKbNDYvGLpChcHH0QJd7T7o8nknIfVa9PONifpJbzmxcORZc2Z1VYpSFl2c8XByJ8Eoaj5hl6fslDi4O43hFj0+ZCXrZdidLcQzdozwPr5DO2MlOHOU4DuvycOJx/dKOoUNwJshOHFVcHuOn3Mc5ymPMmWdW4shhun2KyrKogsx5u8t6MFZ5ZCGOreydKo2cK1cbFuWh2iV125oi60hRpnWjgceGaXHMEkYVCwdUeWCnTqnL77Gwzanoa9vKfRXz+8YoiTpEVee/KzHP27RB1x7ac9qyLs2Svie9zapwN67sjvY9bWWU4+S/NgdvjLKdVqYhZZfy+bDTiC3Zm2T34cr9g6diMuPo2pfR15mqyRk/VubRZT118VX/r45jsDCeoe/saVaZdomlaRmmKOehThJmMo73HtocZV2pCjKkcnc9Q8Zs9nSJIdXcjlBRhGQcbcp0Xl0KEYHV7COrjOMsbD92JXSH5Hgp0TJdy7Pt56ZlljGytNyzD58duwB0PUOPaUp6iLxLecTst4I0M3ar8kjZBDSRcSwCbSvumLOUGGfGNuVjuSyH7l/qioujR4aowDHOkrlWbrAtjZxxcfRMk4rslT0OXo7pcHEMgFfocGL3NzjtcHE0oM9ntViVSo7NlRtXdrtgEuHiGIBpldmqNKwxNhnkOH3AxdEzsyp9igMi5gjWHCu4dVKdLFLvKxdHC0KbK7meKauVO0d5WC33VOVq+tmxFrE8uatp5bWWddR9dkh5dC0fq/KoEqNcc3h2bDQe46yhQ2hMWzl16aCzUsnnTQYLqaRDPBzcSrnC9LLtWq5932rBzJDzac+HbUqflXALB+buJEuVNCX72RxU9nWfnVa2Q5dp9WBPPdGuLIN5ZevPjq0w+YT6aSIZsmlSJ4+hK/Ysrl/akSy+UHlM0td8ixC6lmeXiXbTynbIsjEzrX7yRj7TsPbA6VQToCB+j3vbGEO+P4XUY40lGaJcc7nUntW0+txJdSaPvd6xTAwLJUW5jrm86nBxBGK5edKVRTgIUmWIdWU3xvJ0cThOZKqiyEUa+7jqGX2Ls3BxLBjzKnIuFd0y1Wwmh4y0jTBKXByGObK0PkknY65yyOFOZDmIokoXaUADcYjIbSJySkTuryw7W0S+JCIPFb/Pqry2Q0SOi8iDIvLmTlE5yQ+SXNviR5bWmxVIbtIIoUnG8THgioll1wF3q+o64O7if0RkA3A18OriMx8RkTOiRbsgTB4Y5SjL2NlHKYocrwqECqTPS7IWhdK2T2OSueJQ1S8DP55YvAzcXvx9Ozw1EmUZuENVH1fV7wHHgcs6RzfBVvaaGxAU++w3b30xJWKxQrelrUDK98eUcG7lGCKMkq59HOer6kmA4vd5xfILgO9X3neiWBZMdfSoBXmkeJhRWwmFVP7cOvDm0SUDWTT5hmYZVWJ3jkrNstqhqSKyXUQOicihJ37Y/rkqQ8ojxd2w+my311X2mHfLGnIqwKxyrHstlnxjvjc2MYVR0lUcj4rIGoDi96li+Qngosr7LgQeqVuBqt6iqptUddOzzp09O9bSk+qtSaNt02VeBQ4RyBAzXpsyq4y7xNyljELkUd3PFsq46yS3O4FrgA8Uvw9Wlv+riHwYeBmwDrgnNMhZxJ5gNQvr990syyHG7NLyvU06ES1U5CpHlta3ftjRvLKr0mf2MK1sm04EjJ1plMyd5CYinwBeD5wDPAr8DXAA+BRwMfDfwJ+o6o+L978P+DPgCeDdqvq5eUHMmuTWdIp96sqbShopmyj72Zz0ye5gTxoxmHUghpZnzGfXTlKNu6swjsmljSa5mZ4d2/a+HCkrsbUmShNixzxZ6ccojZI6ecTKNObJI6Rct9Fslvk0FlIckK4yuziexuoArJTkJOEQeTQVh9kh512kATYu1TpOG3LM3EyKo6s0SlwejnXKZk+O0gCj4oiBy8OxTq7SAIPiCM02HMdJjylxLJI0FrGD0emHPo4jU+JYFFwaTirKJnpqeZh6PELos1VSMeRENscG5WhUy0z261WPodgjSE2Jo0q5oZYGgTlpqB6QlsVqWR7zLgZsZW9UeZgeAFal76HnKUYODnFQWB+4Ni8+yyKZpG1Zx9q2tlcQZwlkFCNH65glkFBppJzgNNQBEFscsbajS1xjlEiMbeo69KBOIKMVR0kpkBgZRtuC7yKPoSt9DIH0mWn0FUdfzNre0O0JGbM0KY/sh5zPYx9XDSINaH+/SgsVvbxDVpdYYt8gOFRiVvsZuhJyJ7mhBjpmK45Yz5DtKp8cb/Bb0lQEKe4oPraDPiZtBRJDGl2vXmYrDivkLpAur3UlpjTGLKAm2xYz0+giD7OXY/si1g4o5VHt/7DQRJlHyhjHfHC3oUs5TLtEnappUsrjpobvX3hxxCYHWTj9EUOe5TosZbbeVImIDz5zUmFJGuDiiIZLo1+uX9rht04YEBeHkx3Vs69lecTq47GWbYCLIwpDZBuWD5i+8bLoHxdHIENKYxEPmGln37GWRZ/Zxk0H39/4vX5VJYC+pTHWgyMWfT6cqw/6kkYbYZS4ODrSZwWdJYzytTEdMNNociCNRR4ppdFFFJO4ODpgsWKO5YCZRpsDyUJZ7GczG+nWOZpCGjFkUSUbcUwOi41RObqk/pabJxYOmBR0OZCGKovq/rp+acegT6mH+MIoMSGOs3is03j50MqxhQOtDkwrzZN5nxuTPELOvinLIlV/U6xsI5UwSrK/qjLGDsPQbRpjmVhhP5tble8QYzBSSwNGIA7ofqBsXDnWOJXs6ywe66C3II+Qe4CArYFPbYVRpel2WNreeZhoqsSgSVo6bSRfKY9pOy43aVTXZ6XZUsrDZ8yOg9GIA04/ULpU0Lrso4/ZrhtXjmV1tgmhiUDGWBbzOkpz2+ZRNFWqhNyGbQjKWIfufe+bFHcXy5XcpAEjFEdsUkpoct2LJg94pkByPIiaUrdtuW6vi2Mg5vW3LBqlQKz0yaSgbt/mur9dHAbJtTKNmVChzdqnN67sjrrPdy3vjLauacwVh4jcJiKnROT+yrIbROQHInJv8XNl5bUdInJcRB4UkTenCrxPYjdXmqzPojz6usQ75qxjFjnJo0nG8THgiprlf6eqryl+PgsgIhuAq4FXF5/5iIicESvYIQm5jl+ljYRCK1LMA7DvqfxbOGBOIH3Ek4s85opDVb8M/Ljh+paBO1T1cVX9HnAcuCwgPnOECCTWJeJ5xD7oJre3z8FlIdtx48ru6Nmiy2OVRo+AFJFLgM+o6m8V/98A/CnwU+AQ8B5VfUxE/hH4mqruKd53K/A5VZ05EeVlm9botYfeGbAZpxM7zQ/p+Q4ZV9IkhiHnYlia8Nd0n8e8BJzyWcNVYlx5aTwMfbMkfQTkR4FXAK8BTgIfKpZLzXtrzSQi20XkkIgc+vkPf9ExjHosXeLaz+YoZ71qBSwzitwmcIVSt71tOxZjjvOx1pSaRezMo1PGMe01EdkBoKq7i9e+ANygql+dtf7yodNdH0c3DSu2h3hZUKpBUzncYqAkZvMj5QOfY2a+vc2aTZlxiMiayr9bgPKKy53A1SLyHBFZC6wD7mm63n1cddpPKJYyD6vE6vRdVHLKOiBe5tHkcuwngK8CrxKREyJyLfBBEfmWiNwHvAH4KwBVfQD4FHAU+DzwLlV9smtwMSTi8kiHC2eVaU0oq+xa3hkskCZXVd6uqmtU9dmqeqGq3qqq71DV31bV31HVP1bVk5X336yqr1DVV6nq54KiqxArC3GcFFTlYVkaVULksRAjR7tmHZ6tOG3JRRqhZCWOkM5Tl0AavLnyNDlKo2uzJStxhNJGHi6a5rg88qetPBZKHI7jTKeNPEZ1B7AmNLllvWcb7Zk3tsLqTXs2rhyLFtuRpfVZ3URqktUxHs2aWwuZceQsBosVs0nb3mLcsUm1jX3U17Z3Rl9Iccyr6JY7uaydudsO9x4z1vZNU7o8TmEhxdEEi/JIUTFDZ5+2xaI8LA1h75uuz2BZuD6ONpX9xpXdJpo1qStj2yfahRKzXyEG5T6ONXy8um0WRQnhD21qNMktNeUkt1nEmAAXkkXEqFxtKtFQB1ZTgYRmZG23bzKu0O/vekIIlUsskbSJv5UkGk5yy0IcsWbN9l3Z65hVcaychefJo6+ZvqniiHF/lVD6EkjrzGIs4rAiDYh3YE9WGivCmKTuwLV0kyRoH0/o96WYDRu7OSM/CjimE9/Ipxdi358jBKsHd0omDxJr0mi7Dgv9VXXkWLdMi2OsVCuK9UpTysPiVaaSeUK4fmmHWWnkitmrKpayjRRYF0YVy9IoKcUwGesiCkPPkbDmSgNMZhxjl0ZupJBcqgO6ut5FlEZfmBSH37THHrk8JLqaceSQKaVCz6m7b3g8TIqjxOVhD8vyGNOzWa1jWhyOTSxmH/OezerExbw4POtw5tFEDC6PuJgXhwWsnV2t0LVchuq0dHnEIwtxeNZhl6GlOgYZ7GdzdleAshAHhMljCwcGr+Bjpm3Z5vYQo5RUh/XHlkfKKyvZiAPay2Py+apd5OHCaUbTDtPyPTnKI3bMdXOBcsk8zE9ym6RucFiXHTptYpGLIg7V8g0p01TPZm1zgMYURl+3LajSahRpw0luZoecT2MfV7GVvcE7c/JmKy6MNISWa91+3s/m4ANr1k2aUmRDY3uERHYZB8AetiWMxrFOzGnoVXmkbD51EUeMrGPD0uFnLDt2cOP0D4w143Cc2KTub+mabTR5lEcddbKosn75yFN/z5TIDLLqHHWc2FiVRknbztJ50pikKpE2eMbhOJHpuz+jrSwmOS0DafgZzzgcJyK5SaMrnnE4TkTKpk9KgaSQRdu+Ds84HCcBsfpO+hhS36WDNDtx+KVYJxdC5WFVGpChOMY2kMYZjj4G/XWVxyxpHF25tGM08ZjbxyEiFwH/DLwU+BVwi6r+g4icDXwSuAR4GHibqj5WfGYHcC3wJPAXqvqFkCAnh5nvZ3OWcx2c4clhhHBfM367ZhvQLON4AniPqq4Hfh94l4hsAK4D7lbVdcDdxf8Ur10NvBq4AviIiJzRNrCt7H3qpw7PPBYXi3cgm0Wbk1xTaYRmHSHSgAYZh6qeBE4Wf/9MRI4BFwDLwOuLt90O/Afw18XyO1T1ceB7InIcuAz4apOA2tzh3DOPxaaUR9Mh6JZl0+d9RUKlAS0vx4rIJcBrga8D5xdSQVVPish5xdsuAL5W+diJYtlMuj4SweXhzBKIZVnkTGNxiMgLgX3Au1X1pyJTbxJS98IzZtKJyHZgO8CLL35x0zBqKZstLpDFxrIktnDgGc3rvu9eFiPTKGl0VUVEns2qND6uqp8uFj8qImuK19cAp4rlJ4CLKh+/EHhkcp2qeouqblLVTS849/ld4z8N7/dwciFnaUADcchqanErcExVP1x56U7gmuLva4CDleVXi8hzRGQtsA64J17IjpMnQz2HN7Y0oFnGsQS8A/gjEbm3+LkS+ABwuYg8BFxe/I+qPgB8CjgKfB54l6o+GT1yx8mQmNIYcjxHk6sqX6G+3wLgjVM+czNwc0BcjuNEIEW2AUZGjp78ydyLLo7j1DAr60glDTAiDoCbDr5/6BAcZzSklAYYm1ZfymPX8s6BI3GcPEktjBIzGYfjOGH0JQ0wKg5vtjhOc46uXNqrNMDI4xHklZuUDx16xvIuTRYfPepYJ9XjHSDCSbfh4xFMiwOaycNl4eRKV4nMu/t5Z4EswnNVXBhO7sybX1MnliaPTKiecFM0/U32cVSp2+jJh0k7zlgpxVL+tvJQavNNlZJdyztdFs5CEzKJs3HWMbamikvDcbpT11cY0oQx31RxHCcNu5Z3dh5s6eJwnAxIea+ZLvJwcTiO01oeLg4nGtvYwzb2DB2G05E28shCHIeXNyRbt1f2OFTL0Mtz/GQhjlR4ZY9D6rJzudvD/DiOFNnGrEroz6ZtTpODOaQ8p60/1j7axp6o+3sy3pjrTn0j7pveWlya/UyzcRwLl3HMq+x+ZrNBH1lM+TvVd8Va7zb2sI+roqyrjqek0QLT4kjZtzELl4dtQvdP3edDBTLtszHrUgp5dJEGGBdHbNrsRIvysBRT01i6xNzkMymzhOpP08+EvD4vliopM482LJQ4cqaaWg9N2xisNwXmfccsibSRS4r3htA124AFEkeqM98QDHmVIeTsmWL9fUps8ipcijrVZJ0xso4QaYBxcVx68CiXHjwavJ4Ubde+Pj9rHbldphxDx3SMMo/RxzJ0k8X85dhJ2naYxqqMbS+txbw013Qbmn5H19hilOW07wpdd5NtyEFMbdnK3tafmZltNLwcm504mnB4eUOSSjKrcrb9PosVPdVBPeu7Yq+7bhvGKIwqTeXRqHmyyOM4YjRv2tC1rTvrc0NU9rqYUsaRYt1jl0RXQvs0JhllxgGwfvlI1PWV9HGm9MofxiKW5ayso5U0FjnjSEmqDklLl1tzZxHLsu/O0tGKo+8H1MRgkSp6ahaxLPuUx2jF4TiLyKQ8YvdtlLg4HGdklPJIJQ1wcTjOKEkpDXBxOI7TAReH4zitGbU4cryy4jg5YGMAmMgPgZ8DPxo6lhacg8ebEo83LdPi/Q1VPXfeh02IA0BEDjUZsWYFjzctHm9aQuMddVPFcZw0uDgcx2mNJXHcMnQALfF40+LxpiUoXjN9HI7j5IOljMNxnEwYXBwicoWIPCgix0XkuqHjqUNEHhaRb4nIvSJyqFh2toh8SUQeKn6fNWB8t4nIKRG5v7JsanwisqMo7wdF5M1G4r1BRH5QlPG9InKloXgvEpF/F5FjIvKAiPxlsdxkGc+IN14Zq+pgP8AZwHeAlwNnAt8ENgwZ05Q4HwbOmVj2QeC64u/rgL8dML7XARuB++fFB2woyvk5wNqi/M8wEO8NwHtr3msh3jXAxuLvFwHfLuIyWcYz4o1WxkNnHJcBx1X1u6r6S+AOYHngmJqyDNxe/H07JH645wxU9cvAjycWT4tvGbhDVR9X1e8Bx1ndD70xJd5pWIj3pKoeKf7+GXAMuACjZTwj3mm0jndocVwAfL/y/wlmb+BQKPBFETksItuLZeer6klY3VHAeYNFV8+0+CyX+Z+LyH1FU6ZM+03FKyKXAK8Fvk4GZTwRL0Qq46HFITXLLF7mWVLVjcBbgHeJyOuGDigAq2X+UeAVwGuAk8CHiuVm4hWRFwL7gHer6k9nvbVmWe8x18QbrYyHFscJ4KLK/xcCjwwUy1RU9ZHi9ylgP6tp3KMisgag+H1quAhrmRafyTJX1UdV9UlV/RXwTzydKpuIV0SezepB+HFV/XSx2GwZ18Ubs4yHFsc3gHUislZEzgSuBu4cOKbTEJEXiMiLyr+BNwH3sxrnNcXbrgEODhPhVKbFdydwtYg8R0TWAuuAewaI7zTKA7BgC6tlDAbiFREBbgWOqeqHKy+ZLONp8UYt4z57p6f0AF/Jaq/vd4D3DR1PTXwvZ7XH+ZvAA2WMwK8DdwMPFb/PHjDGT7Caev4fq2ePa2fFB7yvKO8HgbcYifdfgG8B9xUVeY2heP+A1dT9PuDe4udKq2U8I95oZewjRx3Hac3QTRXHcTLExeE4TmtcHI7jtMbF4ThOa1wcjuO0xsXhOE5rXByO47TGxeE4Tmv+H9fJr+MDnfl8AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "touch_matrix = cle.generate_touch_matrix(tissue)\n", "neighbor_count = cle.count_touching_neighbors(touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Average the measurement between cells to reduce noise\n", "Mean of touching neighbors" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAmmUlEQVR4nO2dX6xd1X2gvx+egGzjCLC5Vza44xsHUByRMYxlj5whIoOpCS+Qh5HISJSHRvSBSI3UPhD1odFISJmR2nmaVKJJVMp0gpDaNEiD6oKVKLSesWMSB4Mtgx27g3ute2MyVhxskWLWPNyzrtfd3n/W37PXPmd90tU9d99z1l5nn7W/81v/RSlFoVAouHBd3xkoFArDo4ijUCg4U8RRKBScKeIoFArOFHEUCgVnijgKhYIzycQhIg+JyAkROSkiT6c6T6FQGD+SYhyHiKwC3gYeBM4CPwa+pJQ6Fv1khUJh7KSKOHYCJ5VSP1dK/QZ4AXgk0bkKhcKY+VeJ0r0NeNf4+yywq+nJsmG9uu63NlslvPq6y2E5q2ENl6KldYVV0dKqcgMfJEv7Oq4A8D43RkrvoyjpAKwa5S0FKT6v2J/TZVY3/u9XH3086rk++unPziulbu16XipxSM2xFXUiEXkSeBJANt/OmtdesUr47rVHgzNXZTtHoqV1gZuipVVlKyeTpLuOi8uPDzX73Ynr+U2UdABu4kK0tKqk+Lxifk7H2db5nH3v7412vl/fOPNPNs9LVVU5C5ghxO3AvPkEpdSzSqkdSqkdsmF9omwUujClAbCTgz3lpDAkUonjx8AdIjInItcDjwEvJTpXVuzmQN9ZKEwINtFGXyQRh1LqQ+ArwD7gOPCiUuqt0HT3rt3HJubZtDJ4yYJtHGMbS51GuzmQvUDWcfGaaEOTW9Sxg8N9Z2HsuEhj79p9wee7e+1Rp2aAVG0cKKVeBl6OkVbdhdnEPPNsipF8MnZzgAPsDk7nPn4EEO39NgnDZCcHo7V3hLCHV4EleRxmR8+56UZ/VuD/eY070vBpN8x65OjetftabRor+jjC9qDX60ijjtDowyyIoe+1LcqIzd34N2Lv4dXlH5PcIw/zs4J45dMG36jDt7MhWcQRisuF8I0+/hP/85pjxxxt3yaNUKoFMYRxCQPg3/I6sFIeR7k7StqxI4+buBDcs9L1OY0rOt67dp9TD0tID2WSkaOurLp3uzK7Y33t6fLh1EnDxEYgrtKwqbbYyMLlfYYKw7W6oqXhgimVapTRRCx5VD9D16qli9xtPrfQaoqtOJqk8b/lgdeVUp0XN7uII6Shx9bsXdKAlQWqTiI+kYZZZakWUJcCmGv7jo804GpkMsuC9WtiRB51n6Ftu5RPNKirLX1/djHGQmUTcXzxJ38YLb22D8ZGGk1ogcSonqziw+A0ugpgjOqJTcThK4w6XOQB7pGH7WdXJ4+YVUeN+RnGahRtizq6pGEbcWQhjpt3fEI9cPg/R093nk1BokhFDGlAszhStGfUCSSmMMBdGnU0iSRU9rE+szrW8Wv+lkejp6sF4hJh2Ioj616VUHKURkzqWuxTNYLmNrajibqel5QN2KGs49cAPMrfRk87xviOJiZaHDEnrw2B1D0npjxiRxux0QPyhiANTQp5pJjbBRMsjodHY8+mRR7j6m7dycFBSCM2VyL3I1SloUkhjxRMpDgergxYnWR57OB1dvA6d/F231kpWNIkDc0Q5DGR4qgjljxC56CkbGQr8lgi53lCXdLQhMrjAV7lAV7lcf4yKJ0mJk4c1WgjJrpA5lIwN3HummN38fbyzzSi11bJ5TMysZWGxlUeWhYPWA6iC2HixNFGSNRRLYg5FswqKeSxh/3sJU1rfeiCStXXD+Ez6qJLHuOUhclEjeOwjTYuscY6TZvCZztMOWY1pS7acOUEd7b+fw/7W/+/j7CVp9oE5DLxsEs4oTOUzTJwkJ1eabhGGyZ6jEeIHJ7nd6yeN3UDwFyrKG3ycP2msimYsds2YohDUxVIlzBMfOVhE7XYyMM2SvGVR1NZcBFIiDQALkZaB9ZGHmUAWAd11ZZUC/Ds5gC7OBQtvZjSAFa0ibhIA+wE4PuaLinEXCu2jrayYPt5hkojJjEbSici4ghpEL3EmiiyqH6jdaXpG/JCfHFoNq9YmN4Nm8jDt22kLvJwlYZrxGFbJto+x1jSiBVxaNoij0FFHB/jX5hhsZdzx44wbKOWXRxa/nEhR2lAtxRCGlSrkvCJNFI1lDZ9fjlFGlViRB5ZiEMzw6KzQFJ2v7oQUs3xlUhuNMkhRi9MjGqJ7efj+jlWP7ucpRGLrMSh0QIxf+rIRRox6UMgodGGSVUSMbtut3Ok13aNLnZxaDDSCI06slvIpwlTHrmvPRkDUx66Hp2qmhKbVOM8UjMJ4z5ceJy/tO6mrZJlxNHFEFa7jskuDmXbtjE0muQwbdLQ+EYegxTHFk5znrL7W8GPIY4CTomPPAYpDs00yWMaqmfjJLd5R01s4DxznEl+Htdh64MTxxZOr/j7POsnXiB6Sb0ij7ikkEbMtUA2cH75cUp5+CyePDhxNBEikIusi5ybdEyrPPRmUuPcH6ZPTGloUsjDd8X1wfSqpEJLo0seZoG9yLrlvy+ybnntzxQCqlvAdyjbIcagThTVzyI3do2WWDzouYVmnTQ0pjxOs8Ur/RgMKuKoVlPqcIk6tnLK+rkXje88829zweBxfhvu4PDERx821zO3CMTMzy6PBZ7bpFFljjPeUcj8aINKXwYTcdhIQ1OVxwbeA9xE4UvMb0Ob7QJMebhEIbounuM3tqsMqs9ve08uzw3Jg8Y2+nARRpWqPOoikdibQA1GHCGcZ/2yPMaJWZ1xxWePkbYqTFOj3TouZiOPWNFD03WPXe1xyW+bQEKkUYcWiRZIip3jpkIckCbasN2J3LVwhmxMVI1AbFr5QwQXixRVDi1F27RtJRqSV7P6cpBd0aVhknKryUG0cbhUU3Knq2cgxm5mGteuwb7aC1KeN7TaU/1fzLymlEZqpibiyI26KCSmNHwZd9Ult8ZNuDYCyzGPfTOIiKMwXsZ1o+R+Qw553Mg/Bq6z2kX24pikakoTKQqny4LMdaS+YYZ6QxaWCBKHiJwRkaMickREDo+O3SIir4jIO6PfN8fJ6mSTQzVlnOTSk1PwI0bE8Xml1HZjncKngf1KqTuA/aO/vZiGaGOamXZ59LVcZgxSVFUeAZ4bPX4ORptCOJK7NHIa/zBkpvUa6vc9VHmEikMBfy8ir4vIk6Njs0qpcwCj3zOuiaaQxim2BqdRN9FqWgt+TKbtGlbf7xDlESqOzyql7gW+ADwlIp+zfaGIPCkih0Xk8OVfvB+YjXb0qNEQeUxjY9443/O0yWPoBIlDKTU/+r0IfA/YCSyIyEaA0e9anSqlnlVK7VBK7Vh969oV/zszmr4TwgbeW/4xOcVWa4EMdSp3aI+KSZFHXJreo88K/02k7oqFAHGIyFoRWacfA78NvAm8BDwxetoTwPd9z+EjjzpZ1NEkDx9Z5FLgL7EmqjQ0RR7jo2tl/y7OsyHJZuNVQiKOWeAfRORnwCHgfyml/g74BvCgiLwDPDj625vQyKONqjyGFlnAVVmkEIbJEK/N0HEVyHk2LD9OLQ/vIedKqZ8D/6bm+HvAAyGZqmLKw2w4jTHjVcsjdL+OkOHJ5jofXeM5UguijVRT0qv4XkszP7mJLuRamfJYrPQ1mLKoUpVHdXPxEAY3V+UMcxO1gE21QC0wWyuPPoXRROpZtV0zW9vOO6nzTLREjrHN+bVaJDEEMjhx5I7NNO6uG60qjxylYZJSIHXX0+U8OQgkt3YbMxLxlUj2c1WGSLWgmMsO2haiBWZTZC0pqW7O6nKNIWkUVuLbFlLEkQhXUdSxwGz20UaV1PLoO40hnNMVH3kUcRSmiiHcyH3gKo8ijswpBT0+MaLBaWeQjaN9LT5cKLgwFDH5NJAOLuLQ4zgmfdvHwlUmrUs1J6ayVyWmPI6wPVpahXhoaQxRHpdY3XcWkjEoceS+Rkchbng+RFloFkYjPHOWR8hAsEG0cbQJo7R3TCZNmycNod1goTIs/BKrWcPlnnJzLRM/ctQ2wgiVh95YaZEZ71mJZXGfqxsAhUYKQ15lvSqN3Ig1XyVbcbhWS3zkUbcTm4s8mgreUL4ZY1HdMcxl9zRN7N3WYp3PhTZpxIo69KTMG/gAgA+4wfq1Ez/Jzbctw0YeNts2dsnDdhf1IRXwWc9Iq2mbQZeZquNuy4h9PtsoI0QeTevH3MAHVvKIKQ3ISByxGz5t93Vtojp92Wfv2dDJXzELeFfhtpWJz36k1fcf+tmEMM5rWoduLLURiO1KdV3RR2xpQCbiuH70xmOQolCGblg9CVWXlBsYu+J7PWNJI0Y7Rlv04bs2bl30kUIakIk4CmnxLej6dVcmoJjEkEbshk9THjFW4Yer0QfAG9wdJc06BjWOY8ikaiwcIjGiF5frk/O1zHmcRxtFHGOkjwIc41tyFR9GyEk/5CwNTaxoY5wUcYwZ2x6ZQjjlOqZj+JXXATIJjaV9YxNJFXGko0QcHYT2qLiSa2EfYnVlgZnsR3IOlSKOFlJJoynayFUauTFpMkhVzlKWpyKOBsYtja7/+eI7IrQunT4Hbk0qqSa/nWYLkE4eRRw19CGNnDHlM0R55BqhmNKIWea0NDQp5DFR4vgMR4Nev3W0JXUKbKWRW9RR99o+5RE6mC1nYpS9qjQ0seWRhTh+4zDDb2j4LIqbS2TSJpxNzAcJZH6UwjjJSR5NVRTfL6/TbGmUhiamPLLpjm3aH9aWDbzHPJu8C3PsSCOXmz81m5gPEkDda5s+w75vfFOkIXmxadfQ5bFrcFiXLKrE2nUvG3GYVHeobxJJ3RR6V3nEFEbOsphlMdmNFyqPKmZaubap+F5P18bQrZxqlIerNExCBZKlOKqYIrHZcLrrWyxFO0ZsafgshtMn1Rs8lkjMdGKMJVlgJmpPk06zi5DeE11eX+UB7zSa8B2MmEUbxzjQdeohSCNVui43TKybK0diR16zLA76evl8QU2NOFKRc/XElyHfBLakqLY1CSSnhYqbcJVHEUehEBlTHkOQBrh/AU6dOF7jvr6z0Ctd0cQ0RBupucIqNvAeG3iPNVzqOzudlDaOCWOWhSQFb4hy+CQnmeNM39no5Aqr+s6CE75V7U5xiMh3RGRRRN40jt0iIq+IyDuj3zcb//uaiJwUkRMistcrVwVmWUic/rXyyFUon+Tk8uM5zmQrkKFJIwSbiOMvgIcqx54G9iul7gD2j/5GRLYBjwGfHr3mmyKS3dXMaeHdOqrSWMOl5Z+451lc/p2jND7JyRXSMAkVSOz32yaNHKsrPiOaTTrFoZT6EfDLyuFHgOdGj58DHjWOv6CU+kApdRo4Cez0zl2FLZwO3mhad8fOLw+aDpNI7LEWXZFGTIlcZF22jXdNwqjiKpAUkhxapBGjJ9C3jWNWKXUOYPRb923dBrxrPO/s6Fgw5ujRmLvUg9+8iSusWv6JdSO7Vk9CzmkWnhRdyiEjPm2lYdJXFWac0tjD/qDXh0YZJrEbR6XmmKp9osiTInJYRA5f/sX7zifykUfX4C9bgTQVlhCBpG7TMKkrPDELlSZ0IpwPbfLoM9Los7qS4rP1HXK+ICIblVLnRGQjLH8iZ4HNxvNuh/qSo5R6FngWYGbH5lq5aJrmqqTaqb5troRNYakWkkusaX1+iDTMc3WdB7qjC/3/mFUwfQ1tpOwTbdiQQ/VkDZesPqM69HwVLcaQeSox8I04XgKeGD1+Avi+cfwxEblBROaAO4BDYVlsJ3a1pYoZhfiGpSkaNn3O4/KtkyoCaSOWNPTNpdszcpCGL6dGE+2r6KpZV/Us1cjmzohDRL4L3A9sEJGzwB8D3wBeFJHfBf4v8B8BlFJviciLwDHgQ+AppdSVkAzaTLG3iTxC56jEKCz6pja/dVJUUczzhBac2JPtmmbSxo40UvUSjVsaNtRFIamnQnSKQyn1pYZ/1U7VU0o9AzwTkimNy7ocqaotKdA3durZr7EKT+zqS1UesaXhWx3oItaXR1f+fDdomuNM0m0fTQYxrd6WJnmMe4uDSSWmQHS1JccxDnXEjDTq5BFrN7fPcHQs8sh2yLnPKmCQvs2jUIhN2QIyEr7S0JjyKNFGIUd0pDVEaUCm4ojBedYXaRSyZqjSgAzFERptFAqF9GQljmmSxiquJGv9LxRSk5U4poVVXB3aMu3ymPb3H5up7I4N3VslFTEHaZnS0FxizWC6JYdAihXiV3El+1mw45IGZCYOEy0RV4GY4zjMxqfcG0rNb95pk8gCs8uPQyWt05o2eYxTGpCxODS+Aqlyiq1O8qibVxFaaOqijTq0RKZBINXRraZEwE0k1demkkcTLuUjZhVt3NKAAYhDYyOQriHnXRFI1ySskG8cW2mY+FRhqjdPbGLdiLbD4W2jkab3nWK2bxPmZ1xXTqqyaNulLXcGIw5NVSC+81PMD+w+XrN+3bjDVZfoI7U0qjfuuBs29fur5sPmfY9TIHDtF0WTKIcqj8GJQ3OGOavtILtwkYbGVR4+0UaVLoGklEbTN73v+hKhk+9C3mv13Dlss2m7wXQdfVRTYMDdsTGkAf77rOiFA8dN3Y2aShqzLFitgepCbjvf6XVHcsjX1tHqG7b0JQ0YsDhyoQ+BmPJIKQ1bbBcqyuHmbCOX/NnIo09pwICrKrHwqarUoeVhVmFSCiW2MGKMVanKI3WbSwpS9MT4CMmUh1mF6VsYmqkXR2z6qL4U8iVGFKMl8r3lXUj6p1RVCsAwo4Np4iA7x75ifBtFHIVBops0p41c5FHEURgcpjBylkesxtaDlc0Q+9ivpkoRR2Hw5CyPlPQpjyKOAZJLt2EfNEliUuVRjTaqxJTHcbZZP7f0qgyIaRaGDeu4OFHXqEsasXARhqaIYwCM42ZYZIZFZribo8nP5YtNVBEiD/26i6zrvQ3BRRouW2yCnyiqlKpKgUVmlh8fzWSAURWXqohrtaVuyLntTZgC30ijTXbH2bb8E4PBRBzV6fQxdm6LNWo0JamjDVMamqPcnVXk4dN+0RR5uFzPeTYFRR59VJ2qu+TFEkWVLMRxPR94LdQTKo/XuC9befQhDJNc5BHS6Fm9cX2uaag8XInRrrGJefazJ0Jumhl8VWUSd27rWxqTRIzZryHVlj56e1JLAyZAHOAvj828yxljh+8cSCkN3QBqS0h7xyXWGLes380zKV2stu9jXL0oMciiqhIDm2rLZt6tPa7lsYUzkXPlRipphEQYrlWWpoV99M0zSd2l08zEiANWyqNJEm3URR/jkMkCM6zhcpK0Y1RLdORRFYjP6l82ApmUSMOkq6F0SNEGTEhVxWQz73pJoy8WRjf2JVb3nJNuzKpL6Hqj0zpJrY6hSQMmUByxSdkGslCJBoYgj0uj9b5iURXIJMuk7r0NURpQxNEbVWlohiCPFExDBFJXVdnFoR5yEk4RR4bkKo/cd8NLRYyBXG2v38WhaAKZZxOf4liUtNroFIeIfEdEFkXkTePY10Xkn0XkyOjnYeN/XxORkyJyQkT2psr4OIldXWmKNkxylMc09oiMMwoKkcf8aJUOTWp52EQcfwE8VHP8vymlto9+XgYQkW3AY8CnR6/5pojkudmmIzdxgZu4EJyOjTQ0ofK4yDpWR+qtmWFxOc1poK7qNI737iOPpgFqKeXRKQ6l1I+AX1qm9wjwglLqA6XUaeAkDLT1p4EQgbhIQ+Mjj+pIyVB5aGmY6Q+BBWZYxYfOr2uLMnKTR9eo1lTyCBnH8RUR+R3gMPAHSqn/B9wG/B/jOWdHx5KziXPA0qK7MZb6N9nOkWuO2crjAjcBftLQXGJ14zgP24Ks5XHZUkRVWdSdN6fGzLbra8rjSkuRt30/LttJ+orGlEddz4vLMPhPcSz6ZDffxtE/A7YC24FzwJ+MjkvNc1VdAiLypIgcFpHD7/8i7BtRS0OT04rdN3EhSBoaM/IImX9hE310SSMnFphxur6r+PCaKMS3R6ePyKvalmFL7MjDSxxKqQWl1BWl1EfAn3O1OnIW2Gw89Xaon1qolHpWKbVDKbXjulvXc4R7nPKwiXPLP3XEkkddtOHKXbwdnhGW5BGjsNbJY2Z5Jkv+0tCyCBGyFkjKqCmWWHZxyFsYJjHl4SUOEdlo/PlFQPe4vAQ8JiI3iMgccAfYV9iOcM+KnyaaZFElp8gjN7Q8hiKLXBlKe48mljw62zhE5LvA/cAGETkL/DFwv4hsZ6kacgb4PQCl1Fsi8iJwDPgQeEop5b21mSmP7fwUsJeGJkWbx6QQ2miaWztHX9Rdh5yFouUR0u7RKQ6l1JdqDn+75fnPAM9456gBLRFXcRQK48CUR87SMAlpNJ2KkaO+VZYY7RuF6WIo0ghlUOJ4mJe9X1vaO9IwLTeKDUO8Fp/imFe7x6DEEYqLPEq0Yc8Qb5jCSlzlMVXiKBQKzbjIY6JWALPBppelRBvudI2pmM20y7dtVO60segQkU9lxDHk9g6fuRepsZlPE2P07LTyPR5Nfg4XacAURhxw7eZOVU5wZ7TRnrFpm2vhSowxGC6T8BaYyTbymGZcpQFTKg4bcpRH9abzGYIcc8CWz8zdHOUxzdUVH2nAFFZVXHaMO8GdCXNizyyLtTebyw5jOS3Nl1u1ZZbFbK7NuFgclSpfBhNxbOenzLNxxTHXUaQ+20xqeejoQ0+Td2WWRacbxvZbuU4e4+geHfcKZeYyBjdxIVjqddc35TR5H46NRnXWRb4u7z9EEE0MIuLQ81SqzLPxGpmk4gR3ektD0yWD2eXvAf9QfijSsJGoXjSpbu2TkGpkblWlOo51DAW/i7etrkEKacAAIo4maZjYRCI+0UYK6iKPmAU59Q7pMSONBWaCBHAXbztHHpMgDZO66/ca98XMTi1Zi8NGGnXURSGh4hhCgUtN7OpJjMZnF3nk/hm6CKNvBlFVmTTMApx7YdbkuOq6pktAodW/cTAkaUDG4vCNNobCEAqzJmdpaJrq/EO5xjG5j9eSnyNLcUy6NIZGipsv1RgZM91plMa4yEIcq7iyogX9DHOcYa7vbE095tiPoURI77KZNVxmDZenetZu6qgjC3E0UeSRHznL490V62QvMc3ySEnW4ij0R9tgqByjjzppaIo84pO9OErUMX6GNvy6TRqaIo+4ZC+OHEjx7eq7oVJO+F6XviYPDv1658QgxDFpUYdZgHMrzK7RRt9VFptoI3fWcZFdHOw7G04MQhwQJo85znAucBesWOQmChPfKoqrPELn/EwS5jWPLY+UPSuDEQe4y2Nu1LGr8ZFHzG/UJmnEqraEpBParmHbYKqfM0R5xJZ+3TUfSuSR9VwVW0w5dKHlsbFhLYtU7Rmuz3O5kUMLdMzGUH39zIl8TdfURh51M2M1PtUUl53mNXULJvleM5vX7eIgB9nllX6V+3gtyaS3wYnjDHNs4bSTLOowo4+NzCerq/ve1LbbK+YkjTpCr2udXG7iQnDbRtv1tVlZzXX7y6H1VHUxqKqKJlQaVXKThs3rc++VSTnWI1aDaPX6ue4In/L6x6qyHGA3q7iy4icGg4s4hkKsQjW0DY2HiM/arRqb6qVvtOFbZTnA7tb/m/K4wirn9GGgEUfuxL6xzeiiSCMuIdKwIbSK4hp5dEmjim8EUiKOAVGkkTc6Ohx3e4arLKr4yKNEHIVCRIYmDV9KxFFIgk+359Ax13uN+b4vjxZSWj3a+yWFLFzbOkrEMeWkrv5Ma/Uq1vv+IfcvP77M6mWJxMSngXRw4vg8P4ieZuoGsmmirps4Z3m4bGrVRl2EEfq+TWmkYmp6Vd7g7qjpzYzGGhR5hNE1riRHeeibPVQebdUS3/fdJo17Mlhas1McIrJZRH4gIsdF5C0R+f3R8VtE5BUReWf0+2bjNV8TkZMickJE9oZmcj3nl38gnjxmMluMpi8uLvcFuBVy19flII91Rq5NYkUeMRhHpAH+0QbYRRwfAn+glPoU8O+Ap0RkG/A0sF8pdQewf/Q3o/89BnwaeAj4pog457Aqiyoh8phhsVYaJepYKYMUN3qMNFN1efrII/a0AFtphEYdIdIAC3Eopc4ppX4yenwROA7cBjwCPDd62nPAo6PHjwAvKKU+UEqdBk4CO20z1CaLKq7yaBKGSZHHSpokEiKAWEJyFUhs2cRM74fcP4hIQ+PUxiEiW4B7gIPArFLqHCzJBZanQ94GvGu87OzoWCsuwjCxlYdLtSTWEORJI2UkEkKbQJqqJk3YRh3T1M1ch/U4DhG5Efhr4KtKqV+JSONTa46pmvSeBJ4E+Phvfdw2G7VoeXyGo8vHQtsv5tnkFLrmdjPlTorxDrHSMj/36peI7znq5hyNK8LQxIg0NFYRh4h8jCVp/JVS6m9GhxdEZOPo/xth+U49CyumL94O196BSqlnlVI7lFI71t66xjf/K9ACidXoWaot6clduKZEYldNxklMaYBdr4oA3waOK6X+1PjXS8ATo8dPAN83jj8mIjeIyBxwB3AoXpbbid1TYrs2Q2FyidXjosvJ0KUBdhHHZ4HHgf8gIkdGPw8D3wAeFJF3gAdHf6OUegt4ETgG/B3wlFIqziIAPdEkjxzr+4U0xIo2Ykqjz/EcnW0cSql/oL7dAuCBhtc8AzwTkK/sqLZ5FGEUhkCKaAMyGTk6tJtwaPktTC5tUUcqaUAm4gD4IZ/vOwudzLOpSGNKOcL2vrPgREppQGbT6rU87k8wka1Q8MFcIVzLYztH+smMBamFockm4ihMLymW7w/lNe5rzNcRtmcZgYxLGpCpOEKqLc/zOM/zeMTcXOUUW7MsME0MIa86jznJwzYvuVzfe/jpWKUBmYoDwts8UskD8ikwbQwxjznIwzUPrtc51vYEmn3sZR97eY/1vMf6qGm3ka04wF4e59i0/GMSSx6n2Moptq44lvONaeYtdj5jpdeUTg7ycMW26qKfE2OPEy2MKuMSSFaNo67Y7AX7PI/zOM97n6MqDJMjbM+uoayuAJvHfPNblVGMdFJSFVDXBswxhNV2XZred5c86qogdcKoYspjPe91Pt+VrCMOqI866qKLNnwjjzZpaHKKPFy+9WzTa/o29Xnftq8JvYnrXt+WZswoJ9a10mix6N820hgHotQ1E1fHzrodd6rth/9763Pu5wdeu81XsY0+bKRhEhJ5hEYuPgWzer4YabThk35XlGBie/NX0xxS1ehlHvZ+rW3U8ax89XWl1I6u5w1GHHdxIuo5mwTiKow6um4om5vI9qbMJeKJIaImhnyzxyREHHXUycRWHINu4wihru0jhjSgPoJwvZFsopBcpAEr8xs7X9MqitTodhCfNpDs2zhSYrZ9xJKGxmwb8L2R2lrrc5KGJteBUZNA7GjDxKcXZqrFAUvyiC0Nkxg3UjWNcnMWYuMqj6kXx1AIjV7GwWF2cJjO6nEhU1zkMQhxxG4YNbnEGl5lT7L0Y5K7NOoeFyaTQYgjFZe4utbpUOSRI6lFUSKZ/Mi+VyVFtGEKw+RV9rCHV6Ofb1Jpu5n1/3ZwOFr6MdKsphcrLZ2eScy0U7PvfbeBZVMXcTRJQ1MijzwYRxSjf6c6V6x0v8WXk6647yoNyFwcKds22ijyyJvQG7Lu9aECaXptTCmlkIePNCBzccSmK9owyVEeOdXzbfPik2eb16SMEswf29eE/L+Jb/FlvsWXVxzLZa+fqRLHkDFD675xzUPuVYGuc7RJxEUuts+tE0YKfKMNmCJxuEQbmhyjDui3l8H3vKkilHFKrNrlnCKashFGjKgjRBqQqTj0YiQH2M0Bdgen5yMNTag8Ytzgbd92OUQgtqQK6cdJjGte93rXKKPvKktWs2NtRq7t5oBT2iHSMHHtpo3ZNWdbUG3P4Zu3GDd207lC07Z5D7mKKWRgn8/2lG3Rxq9vnBnOtPrVO7apucP/I1p6uzkQTRgmbfJwLZQ5FvRUN3XbuWKnXfcechWGJnREsK08bKontuLIsqoSSozqjQu+dV2bAVTjpC5PKfORIu3cJdEXoW0aVbIfOZoT5sjSWAW0bjRk34V/XL0V46DvazkOqvsaV4ktDZjQiCMlr7In6TflNBT01EzjtRx3Y+nEiiPnmaRNTFNBT800XstxymNixVEoTCNVeaSopkARR6EwcWh5pJIGFHEUChNJSmlAEUehUPCgiKNQKDgz0eJI0bNyExeip1koDI0shpyLyC+A94HzfefFgQ2U/Kak5DctTfn910qpW7tenIU4AETksM0Y+Vwo+U1LyW9aQvM70VWVQqGQhiKOQqHgTE7ieLbvDDhS8puWkt+0BOU3mzaOQqEwHHKKOAqFwkDoXRwi8pCInBCRkyLydN/5qUNEzojIURE5IiKHR8duEZFXROSd0e+be8zfd0RkUUTeNI415k9Evja63idEJO3YZPv8fl1E/nl0jY+IyMMZ5XeziPxARI6LyFsi8vuj41le45b8xrvGSqnefoBVwCngE8D1wM+AbX3mqSGfZ4ANlWP/FXh69Php4L/0mL/PAfcCb3blD9g2us43AHOj678qg/x+HfjDmufmkN+NwL2jx+uAt0f5yvIat+Q32jXuO+LYCZxUSv1cKfUb4AXgkZ7zZMsjwHOjx88Bj/aVEaXUj4BfVg435e8R4AWl1AdKqdPASZY+h7HRkN8mcsjvOaXUT0aPLwLHgdvI9Bq35LcJ5/z2LY7bgHeNv8/S/gb7QgF/LyKvi8iTo2OzSqlzsPRBATO95a6epvzlfM2/IiJvjKoyOuzPKr8isgW4BzjIAK5xJb8Q6Rr3LQ6pOZZjN89nlVL3Al8AnhKRz/WdoQByveZ/BmwFtgPngD8ZHc8mvyJyI/DXwFeVUr9qe2rNsbHnuSa/0a5x3+I4C2w2/r4dPDaKSIxSan70exH4Hkth3IKIbAQY/V7sL4e1NOUvy2uulFpQSl1RSn0E/DlXQ+Us8isiH2PpJvwrpdTfjA5ne43r8hvzGvctjh8Dd4jInIhcDzwGvNRznlYgImtFZJ1+DPw28CZL+Xxi9LQngO/3k8NGmvL3EvCYiNwgInPAHcChHvK3An0DjvgiS9cYMsiviAjwbeC4UupPjX9leY2b8hv1Go+zdbqhBfhhllp9TwF/1Hd+avL3CZZanH8GvKXzCKwH9gPvjH7f0mMev8tS6PkvLH17/G5b/oA/Gl3vE8AXMsnv88BR4I1RQd6YUX7/PUuh+xvAkdHPw7le45b8RrvGZeRooVBwpu+qSqFQGCBFHIVCwZkijkKh4EwRR6FQcKaIo1AoOFPEUSgUnCniKBQKzhRxFAoFZ/4/CleIevQ7pOEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "local_mean_neighbor_count = cle.mean_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_mean_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Median of touching neighbors" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAVBElEQVR4nO3dX6wc5XnH8e9TQ1BCiDDFgDGmOMiRMBclzpF74SqijQiEG9vCjRzJrYssuRegNlIixSgX4QIErZRUvUkkR6C4PREOwjFYbZSEWq1QrCZwbAGxOXVwwC0ntmxSQEFEJcU8vTizeL3eP/PnnZn33fl9pKOzO2d39jnvzvz2mT+7a+6OiEgRv9d2ASKSHgWHiBSm4BCRwhQcIlKYgkNEClNwiEhhtQWHmd1hZsfM7LiZ7azrcUSkeVbHeRxmtgT4BXAbsAA8B3zB3V8K/mAi0ri6Oo51wHF3f8XdfwfsATbU9Fgi0rCLaprvCuC1vusLwB+NLOLKy/3iG66tqZTpsZQ3a5v35bwFwFtcHmR+b7I0yHxSVOfzNCj0OP/voflfu/uySberKzhsyLTztonMbAewA+Ci669h1dxsTaVMj7t4opb5buJJyBbAfWwMMs+9bA4ynxTV9TyNEnKs5+1T/5XndnVtqiwAK/uuXwec7L+Bu+9y9xl3n7loWXdfndq2GBqjr4sMU1dwPAesNrNVZvYhYAuwv6bHisosW9suQaR2tWyquPt7ZnYv8CNgCfCoux+tOt/+FjC2Vrg/MHqXtxLv5te4zmITTwbbZAlhlq1Rj2Xb7uKJyutDb916IOftazuPw91/4O6fcPcb3f3BKvO6iycu2G5sejuyjFDdxyxbg3YyeTZHYtlk6f3fqXRys2xlE09GM355lFmXajmPo6gPz6zxYTtH8/5DbXcfeRbqsq+Yg/Ou0gkUXZirPlao/7lfzJ3HqLqb6t7KrAeD69gD9tAhd5+ZdL+6jqpUViQFy7Zqw57oogtmna+ETXcZoR+rv/5QK3yMmy2TnqemNv2KrgdVuvYoO46y/1CRQZv0ZOdZOIuu2KHmWWQhrBoYRRf4Mo/XPy55xzRkEFWZb5FloInwyLsOjFrH8nYc0QVH1X0XeQYuxApftRsYnGfR+eVdCJsMjqa366uGx6gxr+NFo1+dAZJn+R+3jiUXHF+e2xhsfuMGr8oT3lugYtlRN2kBDLEi51nI29wRWFeHUMeLxTB1hMi45X/SC3NSwXHtzHLfPnd38PnuZXM0K3kdRi10dazIwx4rxiMHo4Ik9uWgzgAp0sXnDY6p/jyO2BeWqoatuHWtzDGGxDDDnvMUloM6xrfOUxaiPaoixdW9cvcfHYg9SFIIi0F1HH2pKzymtuOIfcEOran/N4WTm1IMjZ7Yx7ZnKoMjlcEPIYUVWYpJ4fmcyuAQSV3s4TF1wRH7gIvkFfOyPHXBITJNYg2PqQqOWAdZZNpMTXAoNESaMzXBISLNmYrgULch0qypCA4RaVbywaFuQ6R5SQeHQkOkHUm+yU2BIdKuJDuOmD66X6SLkgwOdRwi7UoyOES6JMYXyuSCI8ZBFKlbbMt9csEh0lUxhUeSR1VEuqo/PNo8SJBUxxFT4oq0rc1Pf0um41BoiAw3uG400Ykk1XGIyGRNdCIKDhEpLIng0GaKSFySCA4RiYuCQ2QK1b2DNPrg0GaKSHwqHY41sxPA28BZ4D13nzGzK4DvATcAJ4DPu/ub1coUkZiE6Dj+xN1vcfeZ7PpO4IC7rwYOZNdLUbchEqc6NlU2ALuzy7uh3MaWQkMkXlWDw4Efm9khM9uRTbva3U8BZL+vKjpThYZI3Kqecr7e3U+a2VXA02b2n3nvmAXNDoCPXf+ximWISJMqdRzufjL7fQbYB6wDTpvZcoDs95kR993l7jPuPnPpso+c97d9bNTHA4qUFPV7VczsUjO7rHcZ+CxwBNgPbMtutg14quxjKDxEimtiU79Kx3E18BMzewF4FvgXd/8h8DBwm5m9DNyWXS9N4SFSXN3hUXofh7u/AvzhkOn/A3ymSlGD+sNDO05F8qnz7fbRnzk6SB2ISDkh326fzAf5iEgYIT5+MLmOQ0TCKduBKDhEOq5MeCg4RKRweCg4RKQw7RwV6bgyO0iT6zh0HodIODqqIiKNSSo4utptrD0433YJMoWqnEyZxD6OrgZGv154HF5/U8uVTJe1B+c7N6Yhzr6OOjgUGBfq7z66tsCH1hvLLoVHqLdsRBscCo3JRm3CdGUlqGJw7FLp6AZX/CLrydS/yU2hcU6Z/RtrD86PXDFk/FjEOk6jPtwqbxiEfnOouXvQGZZx7cxy/+e5pW2XEZ2QC3Hsr6RNKTKmbY9Z0ZV91Atukfk8YA8d6vvGgpGi6Dgu5622S5AOKBrEwzq3JpT96MwqHUlR0e7jkHDafuWMQZUAaGL/R6gVvKnPq4mi45ALxbqtnSKNZXgKjinX9W4jhdBI8VPtFBwRSmFhT4HGsT4KDhEpTMERGb1KhnN4/U2d31Sri4IjIgqNbkrxhEcFR0TqeHVUGEkdFBwy9VLYXEmt61BwREZdRz0UHmFFERxvcXnbJUQlhYU8RSmMa8hvW6tTNKec6/thJXZNfhZKbx2I9eSwKDqOQb03+ZR9s49cSJsrYTU1nrG+iEYZHIO6GCIptNUpCjmuTb17NsbNlySCo6tCh4e6jkV1jGvXxlbBIZ1U19GrrgSIgkMksC6Eh4KjY7qwULetC/unFByR0wlhaelCaECO4DCzR83sjJkd6Zt2hZk9bWYvZ7+X9v3tPjM7bmbHzOz2ugqXahQe4XUlNCBfx/Ed4I6BaTuBA+6+GjiQXcfM1gBbgJuz+3zTzJYEq7ajem8P11GWcEL/710KDcgRHO7+DPDGwOQNwO7s8m744ASLDcAed3/X3V8FjgPrwpQa78kwTQoZIl1b2PuF/N+7OI5l93Fc7e6nALLfV2XTVwCv9d1uIZtWmULjQlUW2C4u7INSC9+Y1oHQO0dtyLSh3/hkZjvMbM7M5t55/beBy5BxFBrnKHzLKfsmt9NmttzdT5nZcuBMNn0BWNl3u+uAk8Nm4O67gF2w+E1u4x4spqSNTf/Cm2e7vcsL+yiH199UeJ9H0+MY29stynYc+4Ft2eVtwFN907eY2SVmtgpYDTxbrUTJa9K+D4XGaEX2GzU5jrG+R2tix2FmjwG3Alea2QLwNeBh4HEz2w78N/BnAO5+1MweB14C3gPucfezVQpUt1Fcb8Fee3BeYVHQpO6j6dCIVTRfOr197u4Lpis0pC3DwqOp0GgzMPJ+6XQ0H+QjEpOu78OYJNpTztVtSFekFhoQaXAoNKQrUgwNiDQ4RLog1dCACIND3YZI/KIKDoWGSBqiCg4RSUNUh2P13SoiaYgqOPr1QqRogIza4aQgEgkn+k2VUOfqp7wHWyQ20QdHT54AyfP3WN80JN2Tchcc7abKKIObMGVDQPtTJAabeDLJF7LkgqMn1GArNKRtsX/B9DDJbKrUJaUnS6ZbjN8RO0rng0MkNimER7KbKqGk8CRJ9/QvlzF2xeo4RCIX4yaMgkMkAbF1HQoOESlMwSESudi6DVBwiEgJCg6RiMXYbYCCQyRasYYGKDhEohRzaICCQyQ6sYcGdDw4YjupRiQVnQ6OFJJduiWVZbLTwSEi5XQ6OLSpIjFJpduAjr07VkEhEsbUBodCQlKSUrcBUxgcCgxJTWqhAR3fxyHSthRDAxQcIq1KtUNWcIi0LMZP+JpkYnCY2aNmdsbMjvRNu9/MfmVmz2c/d/b97T4zO25mx8zs9roKF5k2KYVHno7jO8AdQ6b/vbvfkv38AMDM1gBbgJuz+3zTzJaEKrZNaw/Os/bgfNtlyJRLJTwmBoe7PwO8kXN+G4A97v6uu78KHAfWVagvOgoQqVsK4VHlcOy9ZvYXwBzwJXd/E1gB/LTvNgvZtMbsY2PwgR8WFHnD4/D6m4LWIt0wrV+P8C3gRuAW4BTw9Wy6DbmtD5uBme0wszkzm3vn9d+WLGO4mAZa3YlMo1LB4e6n3f2su78PfJtzmyMLwMq+m14HnBwxj13uPuPuM5cu+0iZMsYKFR4hVnyFh1QR46ZLqeAws+V9VzcBvSMu+4EtZnaJma0CVgPPViuxvJg6D5FpMnEfh5k9BtwKXGlmC8DXgFvN7BYWN0NOAH8F4O5Hzexx4CXgPeAedz9bS+U51bHPQ6TrJgaHu39hyORHxtz+QeDBKkWJSNw6ceZo2U0W7ZsQGa4TwQHa3yESUmeCA4qFh7oNkdE6FRwiEkbngiNP16FuQ2S8zgUHaH+HpCXG5bWTwTHpvA69v0RkvE4GRx4KD5HROhccRc4iVXiIDDd1n3I+TplTz3vh0dthqjBpx+AO6y49D8OW27b3e3QqOKro0oIai3FHt9YenO/0c9ILk7YCpDPBoTe6xanKoe+uhwe01410bh+HxEOfdZIuBYckT+HRPAWHTAV9iPQ5TWyWKzikFXWt5AqPZig4ZKp0fWdpT91dh4JDpoZCozkKDpkKCo1mKTgkeQqN5ik4pHHagZk+BYckTd1GOzoRHDrdPC5a2ZtR53I/1e9VUWDE6/D6m7TJkrCpCA4FRJrydB7jwkWdy2SbeLKWN70lGxwKi24YFg56V2z7kgwOhUa3KTTy28rsBdPu4onK800yOERktGFh0W8vmz+4XDZEOnFURaQrJoXGoP4QKUIdh0jiiobFoPPD46Fc91HHIZKwqqFRljoOkcTUERa9fR0P5Ly9Og6RjiuzgzS54NChWJFwdFRFpCNm2dp2CZP3cZjZSuAfgWuA94Fd7v4PZnYF8D3gBuAE8Hl3fzO7z33AduAs8Nfu/qMqRarLEAmvyolgeXaOvgd8yd0Pm9llwCEzexr4S+CAuz9sZjuBncBXzGwNsAW4GbgW+Fcz+4S7ny1SmMJCZLRZtlbaSVr17NGJweHup4BT2eW3zWweWAFsAG7NbrYb+HfgK9n0Pe7+LvCqmR0H1gH/kacgBYZIvRo/5dzMbgA+CfwMuDoLFdz9lJldld1sBfDTvrstZNPGUmCIpCN3cJjZR4G9wBfd/TdmNvKmQ6b5kPntAHYAXHP9RcDSvKWISAkhOo2eXEdVzOxiFkPju+7+/WzyaTNbnv19OXAmm74ArOy7+3XAycF5uvsud59x95mly3QemkidQoYG5AgOW2wtHgHm3f0bfX/aD2zLLm8DnuqbvsXMLjGzVcBq4NlwJYtIEaFDA/J1HOuBPwf+1Myez37uBB4GbjOzl4Hbsuu4+1HgceAl4IfAPUWPqIjIZG2ez5HnqMpPGL7fAuAzI+7zIPBghbpEJIA6ug2I5MzRt7i87RJEkjSu66grNCCS4ABq+UBVka6qMzQgsrfV98JD53SIlFN3YPRE03GISDVNhQZEGhxVNlv2sVGbPdIps2xtNDQg0uCA6vs8FB7SBW29UEa1j2PQPjbm2t8xauDy3l8kNW2/MEbbceSRJ23bHmCR0GJYpqMPjmGDVLQ9i2GgRUKIZVmOPjjg3GBV2Z7TTlNJXUzLbxLBAeEGTQEiUl0ywRGawkOkvM4GByg8JB2xLaudDg6I7wkRSUHng0PC2cpsa99lKs3qfHDsZbMW9gD6x1DjOf06HRx72fzBZS3s5dU9dupk4hP1Ked16Q+MfluZjeLr9VIxbmXu/a3KeA7OP8Q8B+cX8vkerDelZencOvFQrtt3ruMYFRo9emWLQxNdTO93XY8Var5bmZ243FZRZt6dC448FB5xq/r8DLt/1QAZdd+Qy1Id4VF2np0KjiKDFGN4xFRT3lrK1JznPnV2Cf0/ee9T5e+TaulXZ+dRRKeCI2X9rXXbitYQ+6bApMcYFyJFwqWO21ZRJYQ6ExxlBimGlXSYNo8yVHn1rGP+TYbY4CHnOrqpPPMM0XVUnUeUwbGXzef9hJhfWXVsT4eaR2qHKetq6ZsUYsxD7GNpe5PF3C/4PujGXTuz3LfP3Z1rMIp+tmKoAS56aC3kobm8C1TexyhbW4gVe9RjVZ13nv8hhWAqqsxnjY5bJ+btU4fcfWbSPKIIjg/PrPFVc+Ge1Lt4opZEHrdwFl0oY1zQ61qpxz1W6HkP+x+mMTD65Q2PPOtE3uCIclOlqqbbuLLbunlOoGrSsJrqrKOOeU97SJQVep2Yyo6jTk28Umrhr6aLYzmu6ygSGp3uOOpU1w7JmA63pq6LY9l0l63giEiXFvS6dXEsmwwPBYfIFBkMj7rCRMEhMmV6YVFnB6LgEJlCdW+2KDhEpDAFh4gUpuAQkcKiOAHMzF4H3gF+3XYtBVyJ6q2T6q3XqHr/wN2XTbpzFMEBYGZzec5Yi4XqrZfqrVfVerWpIiKFKThEpLCYgmNX2wUUpHrrpXrrVaneaPZxiEg6Yuo4RCQRrQeHmd1hZsfM7LiZ7Wy7nmHM7ISZ/dzMnjezuWzaFWb2tJm9nP1e2mJ9j5rZGTM70jdtZH1mdl823sfM7PZI6r3fzH6VjfHzZnZnRPWuNLN/M7N5MztqZn+TTY9yjMfUG26M3b21H2AJ8Evg48CHgBeANW3WNKLOE8CVA9P+DtiZXd4J/G2L9X0aWAscmVQfsCYb50uAVdn4L4mg3vuBLw+5bQz1LgfWZpcvA36R1RXlGI+pN9gYt91xrAOOu/sr7v47YA+woeWa8toA7M4u7wY2tlWIuz8DvDEweVR9G4A97v6uu78KHGfxeWjMiHpHiaHeU+5+OLv8NjAPrCDSMR5T7yiF6207OFYAr/VdX2D8P9gWB35sZofMbEc27Wp3PwWLTxRwVWvVDTeqvpjH/F4zezHblOm1/VHVa2Y3AJ8EfkYCYzxQLwQa47aDw4ZMi/Ewz3p3Xwt8DrjHzD7ddkEVxDrm3wJuBG4BTgFfz6ZHU6+ZfRTYC3zR3X8z7qZDpjVe85B6g41x28GxAKzsu34dcLKlWkZy95PZ7zPAPhbbuNNmthwg+32mvQqHGlVflGPu7qfd/ay7vw98m3OtchT1mtnFLK6E33X372eTox3jYfWGHOO2g+M5YLWZrTKzDwFbgP0t13QeM7vUzC7rXQY+Cxxhsc5t2c22AU+1U+FIo+rbD2wxs0vMbBWwGni2hfrO01sBM5tYHGOIoF4zM+ARYN7dv9H3pyjHeFS9Qce4yb3TI/YA38niXt9fAl9tu54h9X2cxT3OLwBHezUCvw8cAF7Ofl/RYo2Psdh6/h+Lrx7bx9UHfDUb72PA5yKp95+AnwMvZgvy8ojq/WMWW/cXgeeznztjHeMx9QYbY505KiKFtb2pIiIJUnCISGEKDhEpTMEhIoUpOESkMAWHiBSm4BCRwhQcIlLY/wOIjQwrriWFGQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "local_median_neighbor_count = cle.median_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_median_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAW7ElEQVR4nO2db6wcV3mHn7eBhACpSJo/ch3TGHAlmy/BXKUfaFFaBAmRomsrpjKSW1eicj8ECaRWwgmR4EMS0UrQbyCZEmHVNGkU17kGIQqNqFClQnJtBYhzZWJICsZWHJS0oCJCHd5+uDPOZDO7O3/O7Jwz83ukq92dOzv77plznnnPmTOz5u4IIUQdfqvvAIQQ6SFxCCFqI3EIIWojcQghaiNxCCFqI3EIIWrTmTjM7GYzO2lmp8xsf1efI4RYPNbFPA4zuwj4AfBe4DTwGPBBd38y+IcJIRZOVxnHDcApd/+Ru/8aeABY7uizhBAL5jUdbXcj8JPC69PAH0xb2S6+0nn9dR2FMp/Xve2XvX12LPzqv1/fdwhTed2bxr1/fnVqgfvmf479zN2vmrdaV+KwkmWv6BOZ2T5gHwCXvhn+aLWjUOaz+cvHe/vsGFhb2d53CDPZvDze/bN263bYsMAP/Ir9V5XVuuqqnAY2FV5fC5wpruDuB9x9yd2XuHiu4ERHxC4NESddieMxYIuZbTazi4HdwNGOPkuIQbJ2a7xS76Sr4u7nzezDwL8CFwH3ufuJLj5LNEOZRtzELA3obowDd/8q8NWuti+ak5I0to5wfCN2aUCH4kiJtVu3s3UEA6SpCGOMskiNUYtjDLLISUUaIg1GK46xSCM1YYw920ihmwIjvchtLNJIjbFLIyVGJ44xSSO1bCO1eEOTSrYBI+iqjEkUOak2wDFnHClJAwaecUgaIhVSq6uDFsfYkDTSJiV5DFYcKe2EEEgawyCVejtIcaRS+CEZ8/jA0Eih/g5SHF1yiD19hzAVyWOdmPdRVdrKY+uXj1/464LBiaNLW+cVMuaKuXX5+IW/MZLCPqpK3brctSyKDE4cXTFZEVOomF3I467lO7lr+c7g24X2ZZriPprHPAksUhZFOrlZce0g3rTkIe4A1kXhVal8ezgU/HMXwbwB1XmCuHvl3laff2x529T/1SnTefuo7f65jYcuPD/MrlbbakI+x6NN/a48T+Qrdszdl+atNhhxhJRG3SNVquLImRRInYyiqTxmSSOnSrlW3VdN91FRGkX6EEhbKsmjojjUVSlwiD2dpLe38dDUChgDxTGRut2QJt2WKtKA+VLouisya5/FvD+nEfLgOoiMo22BhKiAk0e0eRUrxiNWm8ZQJfOoKoxJyrKFrrPCqmUR436cx8zMI6WM43Vv+2Vv565DH7WqZhf5erEcudrGMS/zaCoNCDPo2VV2Esv+q0OIthaFOHKajA7HMlmmjQRik0hTpsmjjTRyQjT8qtuoux+GsO/qEkVX5dKlbb55dXYqWZZehZBGqCPREXYE2U7OIlPg0JW+2G0JIY2QzOuytC2LlLoupV2Wil2VZC6rjyWzWBR9nwJsw13Ld7KTh/sOozZjyxq2fvl448v5o+qqiHK6rNBjayzTMsyxlUNO0wOyxCFGx6Q8xiqNnCbykDgSYQjTp2MiL8/YpbGogde6JyYkjgTIxwskj7B00SBDbrO4rdgEJ3EQ/oxIl0ge46BMFDHJI5mzKl2RS2OePIpnCY6w48Lr/PlOHu5EQGVnJw6xJ/nrY4ZMLvfQ18dM/q/Ps23JzOPoghAGn2zYoeUx77RmyCs/QxL76diusszJ793VVPciIQWyZu9M5+rYRYmji0Yyq4G0rZx1G1+TS9EXLbrYCPX92wo+ZN1sIxKJo4RFiyOnSeVs2gCnVdBZYyOLajyxU7ccmnzfsv3TVdbXRCASxwR9SWOSKpUzVAPcw6Fag6ltBJK6NHIWvX+6HPDsUhyjHxxdNMVB1Vn/D0HdMzB1B3iHIosis8og9PeN6SxJXSSOnpg8SzO5rC+qyiOGWLtiUu5D/q5NkTjEq1jkUTdmUv6uXZ+q1QSwCIixghZjyuepxBin6IdWGYeZPQP8AngJOO/uS2Z2BfDPwHXAM8CfuvsL7cIcPjE2yi5jKk6iE+kRIuP4Y3e/vjASux94xN23AI9kr4V4FSlN9RevpIuuyjJwMHt+EIZbO1Txm1Est7GWYerfu604HPi6mR0zs33Zsmvc/SxA9nh1y88IQleDRalXgEVTVl5jK8MhfN+24niXu28H3g/cbmbvrvpGM9tnZqtmtnr+ucUMgaR2C742xFg5Z8UUY7xiOq3E4e5nssdzwBHgBuBZM9sAkD2em/LeA+6+5O5Lr7nq8jZh1OIwuwYvkKpX/C6SKrHEFG9XLOI7LqJ+NxaHmb3BzC7LnwPvA54AjgJ7s9X2Aittg+yCkIUbS4U/wo5oYilSJ6YY40+NRcxIbXM69hrgiJnl2/knd/+amT0GPGhmHwJ+DHygfZjdcJhdSU/7hfgbWuzxDZXbeKjTzGM0F7nNI5RA2l74Nu/9i7zSti2Lvq1AnXKsu70iXV3c2DVVRKKrY2vSlzjKKlTZNhbdCNuyiEv1u7ympmr8VbcfgzhyZglE4qhJyC5LiKyh7CK4tixSHiEbSoi7rNX57k1jn/YZMUljkkmJVBWHrlXpgMmKkg9a1hm8jLmyzSN07MWzRE23vYhyT3GfNT1gKuPISH2QtCqLyDpibkBd3uqx7DNiLosieeahjEOUkkpF7opp3z9kuaRYxnUPnBLHCEmxYodkstvYRXkMvYx1I5+RortbrdNlA09FHk3meyjjEGLENJ0kJnFkDP36FSFCInGIoKSSnot2B0uNcQgxMkJk1xIH45nDMUnogVFlG/ETqks+anGMVRg5oW4YLGGkQchxvNGJY+yymKTtaVlJIw1CD/4PVhwSRD0WdZWuWDxdnDEcnDgkjG6QMNKkq2kGgxOHCIdkkTZdzk3SPA4hRG0kDiFEbSQOIURtJA4hRG0kDiFEbSQOIURtJA4hBkjXt4mQOIQYGFH/dmyMaNaoGDuLuiGVZo4KMQAWfQe7QWUcuv2fGBuH2dVLvR+UOEDyEOOhz7o+yK5KWYFq/EOkTkwHxcFlHNPoK6UTYoiMRhxCiHBIHEKI2oxOHOquCNGe0YlDpMlOHh7979zGxFxxmNl9ZnbOzJ4oLLvCzL5hZk9lj5cX/neHmZ0ys5NmdlNXgYvxUBSGBBIHVTKOLwI3TyzbDzzi7luAR7LXmNk2YDfw9uw9nzWzi4JFK0bFLElIIP0yVxzu/i3g+YnFy8DB7PlBuHBX22XgAXd/0d2fBk4BN4QJNQyaz5EGVaUggfRD0zGOa9z9LED2eHW2fCPwk8J6p7NlQlSmiQjGIJCYDnqhB0etZJmXrmi2z8xWzWz1/HMvBA6jnJgKXnTD0OURC02nnD9rZhvc/ayZbQDOZctPA5sK610LnCnbgLsfAA4AXLq0rVQuYnyo4ZcT2zSCphnHUWBv9nwvsFJYvtvMLjGzzcAW4NF2IYqxEEoaQ5JPrJdKzM04zOx+4EbgSjM7DXwC+BTwoJl9CPgx8AEAdz9hZg8CTwLngdvd/aWOYq+FuilxM6TGHooYhZEzVxzu/sEp/3rPlPXvAe5pE5QYF6GlkfpPV8YsjJxBXlY/ibKNeFGmsU4KsiiiKedC9Exq0oARiEPZhoiZFKUBAxeHpCFiJlVpwMDFIYToBolDiIGxtrK988+QOESvpH7qNDZyaXQtj1GcjhViDEzKovh66/LxoJ81aHEUB580UBovxayj7byOsWYw8zKMtZXtQeUxaHEUOcwuySMBJht+HZFIGtXWCyEQc+//wtRLl7b55tVDC/1MSSQ9ZkkkVWmEOCXbdDyjTCBr9s5j7r40772jFUcRSaQakw23r8YaSxyhaCuPNgOhk/KoKo7RdFWmIWnMZ9Z9P/totKmLYpK8DjYRyCJOvZYx+tOxKc/e65oqt+PTRWrhuI2Hah3IQkij6TZGLw5RTh0hjOF+n4ukijxCZhpNtjX6MQ51VdbpouEPrUvRF8WsuPOuyQ7TGIcQQyA/uN29cm/PkbyMuioCUHYQOzFJAyQOIUQDJA4hIie2bAMkDiFEAyQOISImxmwDJA4hoiVWaYDEIQrozEo8xCwNkDjEBJJH/8QuDRi5ODRrtBzJQ8xj1OLQBW7TkTz6IYVsA0YuDiFEM0YtDnVVZtMm61DGUp9Usg0Y2UVuEkV9jrBD9/0Ur2Kw4pAkwpHLYGi37IuJlLINGKA4JIzuKGYfkkY4UpMGDFAcolskjLCkKA0Y+eCoEH1z1/KdfYfQCIlDiJ65a/nO5AQyVxxmdp+ZnTOzJwrLPmlmPzWzx7O/Wwr/u8PMTpnZSTO7qavAhRgaKcmjSsbxReDmkuV/7+7XZ39fBTCzbcBu4O3Zez5rZheFClaIoZOKPOaKw92/BTxfcXvLwAPu/qK7Pw2cAm5oEZ8QoyMFebQ5q/JhM/tzYBX4a3d/AdgIfLuwzuls2cLQj0uLIVCUR4xnXpoOjn4OeCtwPXAW+HS23ErWLf3hFjPbZ2arZrZ6/rkXGoZRji5eE6JbGonD3Z9195fc/TfA53m5O3Ia2FRY9VrgzJRtHHD3JXdfes1VlzcJYyaShxgKMXZdGonDzDYUXu4E8jMuR4HdZnaJmW0GtgCPtguxOZKHEN0wd4zDzO4HbgSuNLPTwCeAG83seta7Ic8AfwXg7ifM7EHgSeA8cLu7v9RJ5BXRmIcQ4ZkrDnf/YMniL8xY/x7gnjZBCSHiZhQzR9VlESIsoxAHSB5ChGQ04gDJQ4hQjEocQogwjE4cyjqEaM/oxAGSh0iLIU05TxrN6xCiHbp1YELoZsEiFkaXcaSYbezk4dKfKKjzswVChGRUGUdK0qgqhbL1lIkMgxjHNnJGJY7YCZFBSBrDIGZpwIi6KjFnG9O6Ik23JdImdmmAMg4hoiEFYeSMJuOIFWUIAtKSBkgcvSJpiFSROMRcJDgxicY4xFSKwtCPTYsiyjhELZR9CJA4xBRmCSLk6WORJhKHeBWSgpiHxNETQ2icQ/gOohkSRw/E3ODqxhbzd0mJGH90aRajEUcsN++JuaE1jS3m75QSKcljNOKAZvJI8TqSJjG3jU0DpmFIRR7mXvqb0Avl0qVtvnn10EI/c9pFb3Uqf9U5DTE2qGLsXcSn+R7t6G0K+g475u5L81Yb7QSwYvZxGw81ajzF90xrKDFKA+KNS6TBqLoq0wjRiHSHrlcy5u8egti7LBJHQIr9fDUc0ZaY5SFxdICkIUIRqzwkDiFEbSQOIURtJA7RGeqyDReJQ3SK5NGOWG8pKHGgyUpdE7M8Yo4tVmlAguJYW9kedHv5DFLJo1tibKAxnzqPWRpQQRxmtsnMvmlma2Z2wsw+ki2/wsy+YWZPZY+XF95zh5mdMrOTZnZT2yDXVrZf+MtfhyDm31oZIjE00HyuzWQsMcSWE7s0oMK1Kma2Adjg7sfN7DLgGLAD+AvgeXf/lJntBy5394+Z2TbgfuAG4HeBfwN+391fmvYZZdeqVJHD1uXjc9cpY5YwYqpAQ6VNdldlmn/V95fRd+bZuzRCXavi7meBs9nzX5jZGrARWAZuzFY7CPw78LFs+QPu/iLwtJmdYl0i/1kl7jrZxNrK9lryqJJhHGGH5NExO3m4dgOdNaU/5YsNc3oXRk1qXeRmZtcB7wC+A1yTSQV3P2tmV2erbQS+XXjb6WzZTJp2P6rKo063RPKIhyr7YZZAmtxeoO+sIwUqi8PM3ggcBj7q7j83s6mrlix7VX/IzPYB+wC46s1Vwygll05RIG3HLySPbpnW2NuUeRf3TelSIqllGUUqnVUxs9eyLo0vufu/ZIufzcY/8nGQc9ny08CmwtuvBc5MbtPdD7j7krsv8dtXNY3/FeQCCTXoqSNP9xTPbMQo6q5iSlkaUO2sigFfANbc/TOFfx0F9mbP9wIrheW7zewSM9sMbAEeDRfybEKfKZE8uidGYRQJHV/q0oBqGce7gD8D/sTMHs/+bgE+BbzXzJ4C3pu9xt1PAA8CTwJfA26fdUYlBSQPEYohSAOqnVX5D8rHLQDeM+U99wD3tIgrOjTmIcTLJDdzVAjRPxJHDdRlGS+x/LxGLIz2ZsVCVGEPL89ozuWhSxWUcYgIKDbOWNjDoalxHWbX6DOQwYnj7pV7Oxu5PsyuKCv5NFKINY8xplirxjJmeQxOHDldnvaKqZJPI8UYY4i5bgx15dH0wsyp/MPE34KI4pfc7G1LzqdXO9l2iLtET6sch9jTettdMFn5Q8a5h0NBtjergfZZrm3kNW/so6weNb5FxDxJ/GWzzVa9Onbw4oB28ph3RIlNHvMqftN4Q8moSsPsQkzzthkq25kmj6bdmlKx1M0s6kikojgG21Up0rTbUmVnx5Be51SJpU68+QBh2XuafO+q72lbpnXjDbkPy+pMm7GQvGtzoYuzwO7ILEaRcRSpmn3U3dltjpJt0/8mFX/y80JsYxaxbD/E9+6LtVtb3PmuatahrspspgkkxEh5iLS4aqOJpeJ32SBTbuwhaSWOMspkInHMZ1IeIU+vLeKIHlsDyuONLa6hEFwcOUWBSBzVyOXR1Tn5Q+xp3ZDKBKLGOS46k0ZOLg8Njlbj7pV7O53IE6KBxzjfQQyMmoOuoxdHKsQ4w1IMjBrykDgI/yNPXSFpiFiQODJSkYcQMSBxFJA8hKiGxDGB5CFio/MzKg2QOEqQPISYjcQxBclD9M3ardujzDZAtw4UIjpilUURZRwzUNYhFk0K0gCJYy6Sh1gEMXdLypA4KiB5iC5JSRg5EkdFJA8hXkbiEELURuKogbIOIdaROGoieQghcQghGiBxCCFqI3EIIWojcQghaiNxCCFqI3EIIWojcdQk+K+NC5EgcfyuitlzwP8CP+s7lhpcieLtEsXbLdPi/T13v2rem6MQB4CZrVb5IZhYULzdoni7pW286qoIIWojcQghahOTOA70HUBNFG+3KN5uaRVvNGMcQoh0iCnjEEIkQu/iMLObzeykmZ0ys/19x1OGmT1jZt83s8fNbDVbdoWZfcPMnsoeL+8xvvvM7JyZPVFYNjU+M7sjK++TZnZTJPF+0sx+mpXx42Z2S0TxbjKzb5rZmpmdMLOPZMujLOMZ8YYrY3fv7Q+4CPgh8BbgYuC7wLY+Y5oS5zPAlRPL/g7Ynz3fD/xtj/G9G9gOPDEvPmBbVs6XAJuz8r8ogng/CfxNyboxxLsB2J49vwz4QRZXlGU8I95gZdx3xnEDcMrdf+TuvwYeAJZ7jqkqy8DB7PlBYEdfgbj7t4DnJxZPi28ZeMDdX3T3p4FTrO+HhTEl3mnEEO9Zdz+ePf8FsAZsJNIynhHvNGrH27c4NgI/Kbw+zewv2BcOfN3MjpnZvmzZNe5+FtZ3FHB1b9GVMy2+mMv8w2b2vawrk6f9UcVrZtcB7wC+QwJlPBEvBCrjvsVhJctiPM3zLnffDrwfuN3M3t13QC2Itcw/B7wVuB44C3w6Wx5NvGb2RuAw8FF3//msVUuWLTzmkniDlXHf4jgNbCq8vhY401MsU3H3M9njOeAI62ncs2a2ASB7PNdfhKVMiy/KMnf3Z939JXf/DfB5Xk6Vo4jXzF7LeiP8krv/S7Y42jIuizdkGfctjseALWa22cwuBnYDR3uO6RWY2RvM7LL8OfA+4AnW49ybrbYXWOknwqlMi+8osNvMLjGzzcAW4NEe4nsFeQPM2Ml6GUME8ZqZAV8A1tz9M4V/RVnG0+INWsaLHJ2eMgJ8C+ujvj8EPt53PCXxvYX1EefvAifyGIHfAR4Bnsoer+gxxvtZTz3/j/Wjx4dmxQd8PCvvk8D7I4n3H4HvA9/LKvKGiOL9Q9ZT9+8Bj2d/t8RaxjPiDVbGmjkqhKhN310VIUSCSBxCiNpIHEKI2kgcQojaSBxCiNpIHEKI2kgcQojaSBxCiNr8P6EJRtuxRzkBAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "local_minimum_neighbor_count = cle.minimum_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_minimum_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAARvElEQVR4nO3dXawc9X3G8e9TN9DykmLKm2M7xYl8YXNRYh25lU6FqBDh5cZwkcq5qKwIyb0wUiK1F3YjFCSDTCuF3hHVUVCsKsW1lAQsFbVBViQUXwSMxYuN5eAADSd2bFJAIa1CYufXizNbr092z+7svP1n5vlIR2d3zuzsb/8785z/vCsiMDPL4/eaLsDM2sfBYWa5OTjMLDcHh5nl5uAws9wcHGaWW2XBIeluSSclnZK0s6r3MbP6qYrjOCStAH4E3AksAC8Cn4+I10t/MzOrXVU9js3AqYh4MyJ+DewHtlT0XmZWs9+vaLqrgXeGni8AfzZu5CuuuyKuufmPpprwNXxQqLA2u+KXv6r8Pf73qj+o/D26ro7vaeD0VTeVOr0zL/3s5xFx/aTxqgoOjRh2yTqRpO3AdoCPf/LjPHDkC1NN+H6eLlpba206fKLy9zg6v67y9+i6Or6ni37GQ/O7SpvaI9rzX9OMV9WqygKwduj5GuD08AgRsTci5iJi7srrr6ioDMur3pne2qqq4HgRWC9pnaTLgK3AwYrey8xqVklwRMR54EHgP4ETwIGIOF50un1eTamTex3tsvvwntrfs6ptHETEs8CzZUyr74HRxIK86fAJjs5vqP19264voZv0kaP387RDoyczYhc0+V3V3etINjj6Hhjg0GiLTYdPJPFd1Rkela2qFNHn0EhhBhzw6spkKX1fdUouOPoaGn2dAduq799XMsHRt8Do+4zXRm34znYf3lPqAWHjJLGNo2+HkbdhBhxIZf29aW1qg92H91S+vSOJ4LD0tWnBsUVVhoeDw6bm8LABB4fl0tfw8N6lSzk4zCw3B0fN+vof25pR1XYOB4eZ5ebgsNzcazIHR426tMB16bNM0vZjWapYXXFw1KTNM944XfxMS3XlM5YdHg6OGnRl5hvFn609ygwPB0dFBt3brs18o3TxM3bxM0F54eHgKFlfwmKpPn7mtiojPBwcJelrYAzryufvyueokoPDStX2ha7t9U+raK/DwWGla+vC19a6Z1UkPBwcZvQvNAZmDQ8Hh/VeX0NjYJbwSObSgWZN6HtoDAwu3fnIlOO7x2G95dBYNMu1RtzjsFYYXsh9UZ3yzNqWDo4S+D9XdUa1rUOkeV5VKcihUZ1p2nbW9u/793Z0fkOh0HWPw0pVRg8g70K9dPzlauhrYJTdM3NwWGmKzpxlLdSD6Sytx6FRHgdHAX2dEUdJJTSWTvPo/Ibefk9Vbv9xcNjMypoxq1yw+xoaVfPGUZtJG0LDquMeh+VSZvfXoVGdqndTOzhsoipmQodGuxUKDklvAx8CF4DzETEn6Vrg34CbgbeBv4qI94uVaU3wwVU2ThnbOP4yIm6NiLns+U7gUESsBw5lz82sQ6rYOLoF2Jc93gfcV8F7NM5d7dm57dqv6DaOAL4nKYB/joi9wI0RcQYgIs5IuqFokanwDF+c27AbivY45iNiE3APsEPSbdO+UNJ2SUckHXn/3fMFy6heH2f4sj9zH9uwqwr1OCLidPb7nKTvApuBs5JWZb2NVcC5Ma/dC+wF2Dj3h1Gkjqp4Rr949GUZ07F61LFRe+Yeh6QrJV09eAx8FjgGHAS2ZaNtA54pWmSd+nQjpWkVbQu3ZX3q2hNWZFXlRuAHkl4BXgD+PSL+A3gMuFPSG8Cd2fNW8Aw+nk9fT1+du89nXlWJiDeBPx0x/L+BO4oU1QTP4JPlOX191PhWvqaOtfGRozazcaevD//NqtH0wXk+yc0K8zah/nFwWGkG4eEQ6T4Hh5XKodEPDg4zy83BYWa5OTjMLDcHh5nl5uDINL1f3KxNHBxmlpuDw8xyc3AM8eqK2XQcHEs4PKwNmj7M38ExgsPD2qKp8PDZsWP4hsXWFsudpVwV9zim4NCwNqhzPnWPw6xDhsOjyh6IexxmlpuDw8xyc3CYWW4ODjPLzcExgfeomP0uB8cyHBpmozk4xnBoWJtVfTCYg2MEh4a1WR1HkPoAsCEODGu7ug47d3CYdUDdJ2Y6ODLubVgb+d6xDXFgWFs1efmHXgaHw8LaKKXrxPRur4pDw6y43gWHmRXn4DCz3HoXHCmtJ5q1Ve+Cw8yKmxgckp6UdE7SsaFh10p6TtIb2e+VQ3/bJemUpJOS7qqqcDNrzjQ9jm8Cdy8ZthM4FBHrgUPZcyRtBLYCt2SveULSitKqNbMkTAyOiHgeeG/J4C3AvuzxPuC+oeH7I+KjiHgLOAVsLqfUcnh3rFlxs27juDEizgBkv2/Ihq8G3hkabyEbZmYFpfRPr+yNoxoxLEaOKG2XdETSkfffPV9yGaOl1PBmbTZrcJyVtAog+30uG74ArB0abw1wetQEImJvRMxFxNzK63t55LtZbk3fM3Zg1uA4CGzLHm8DnhkavlXS5ZLWAeuBF4qVaGZLDQIk2XvHSnoKuB24TtIC8BXgMeCApAeAnwCfA4iI45IOAK8D54EdEXGhotpzSSGlzarQxL1jJwZHRHx+zJ/uGDP+o8CjRYoys/w2HT5RW3j04shR9zasL+qa13sRHGZWrs4Hh3sbZuXrdHA4NMyq0engMLNqODjMOqaOnraDw6yDqg4PH+tt1lHD4VH28R2dCw5vEDX7XWUfHOZVFbOeKPPclk4Fh3sbZpOVESCtX1VxWJjNpsjqS6t7HA4Ns2a0NjgcGmbFzboctTY4zKwcs4RH67dxmFlxecPDPQ4zy83BYWa5OTjMLDcHh5nl5uAws9wcHGYGwEPzu6Ye17tjzXouT2AMODgMuHi9Bh+R232zBMVSDg675ESno/MbHB4dVEZYDGtNcHhmrsaosyMdHt1RdmAMJBEcV/zyV55RazbpdGqHR7tVFRgD3qvSQ3XenNjqV3VogIOjV47Ob8gVGg4YG8fB0QN5A2Ppa82WcnB0XBkLfpHgsW5ycNjUHB424OAws9wcHGaWm4PDpubjOtph9+E9lb/HxOCQ9KSkc5KODQ17WNJPJb2c/dw79Lddkk5JOinprqoKN7Pxqg6PaXoc3wTuHjH8nyLi1uznWQBJG4GtwC3Za56QtKKsYi2/snoJ7m20T5XhMTE4IuJ54L0pp7cF2B8RH0XEW8ApYHOB+qwERRd6h0Z7VRUeRbZxPCjp1WxVZmU2bDXwztA4C9kwa1je+4UOxndotF8V4TFrcHwN+DRwK3AG+Go2XCPGjVETkLRd0hFJR979YMYqLLdpgsBh0T1lh8dMwRERZyPiQkT8Fvg6F1dHFoC1Q6OuAU6PmcbeiJiLiLnfrLmplhNzbNGoYHAPo/vKDI+ZTquXtCoizmRP7wcGe1wOAv8q6XHgE8B64IVpp7s0POrYrdRXgzuVOyj6ZffhPfVcAUzSU8DtwHWSFoCvALdLupXF1ZC3gb8BiIjjkg4ArwPngR0RcWHW4oY/oEOkfA6NfhosS0UCRBEjN0HU6hNzq+KBI1+YalwHiFl5lobHI9rzUkTMTXqdjxw1s9xaFRzubZiVa/fhPTMtV60KDjOrRt7wcHCYGZAvPBwcZpabg8PMcnNwmFluDg4zy83BYWa5tSY4fAyHWTpaERwODbO0JB8cDg2z9CQdHA4NszQlHRxmlqZkg8O9DbN0JRkcDg2ztM106cCyXcMH3M/T//98cHNjX6HKLE1J9jgGfHd0szQlHRxmlqbkg8O9DrP0JB8cXXV0foND0VqrFcHRtQVs+PN07bNZP7QiOKA7C1hXPof1WxK7Y6fV9juPjQuNsnY/L51+m9vK0taq4GiraXsZw+PlWegnBdIoDhUrojWrKgNt6+rPWu8sYZN3+m1rS0tH64KjTYoumMu93gu+NcnBUZGyFupR03FgWNMcHBUoe8Ee7l04NCwFDo4WcWhYKhwcZpabg8Mq4Y233ebg6LmqF26HRze1LjiqOHDJM3d5RvU03L7d07rgqGKPRRXT7ZtJqyZu326ZGByS1kr6vqQTko5L+mI2/FpJz0l6I/u9cug1uySdknRS0l1lF13lMRJ9NFjo87ZH3te5vbtjmh7HeeBvI2ID8OfADkkbgZ3AoYhYDxzKnpP9bStwC3A38ISkFWUXXmQmHDeze8a+NAyqaA+3cTdMDI6IOBMRR7PHHwIngNXAFmBfNto+4L7s8RZgf0R8FBFvAaeAzSXXDeSfCadZGDxjX2pciBQNbmu3XNs4JN0MfAb4IXBjRJyBxXABbshGWw28M/SyhWxYJaroJnuhGK3Knoi1y9Sn1Uu6Cvg28KWI+IWksaOOGBYjprcd2A5w0yeLnd0/6noWZZxgVsap7Taab4HRblP1OCR9jMXQ+FZEfCcbfFbSquzvq4Bz2fAFYO3Qy9cAp5dOMyL2RsRcRMytvL6cy4KUvYfEYVA9t3E7TbNXRcA3gBMR8fjQnw4C27LH24BnhoZvlXS5pHXAeuCF8kpeXlW7a+t8T7PUTfOvfh74a+A1SS9nw/4eeAw4IOkB4CfA5wAi4rikA8DrLO6R2RERF8ouvE7jVlscGNZXE4MjIn7A6O0WAHeMec2jwKMF6krO0vBwaFifte7I0RQ4NKzvHBw5eFdkf/l7v5Svcm62jFE3z/IuZPc4LAEp/jdfrnfpnqeDwxKR0oLok/Ymc3BYMlJYEGc5/6mPHByWlDYuiH1cdfHGUeukvPfRLWPBb/u9jfNwj8OSU8Ud8Oq6Ollfeh4ODktSFVcjq+viTX0ID6+qWLLGHTeR+vVS+rC64uCw5PXhP3jbeFXFzHJzcJhZbg4OM8vNwWFmuTk4zCw371UxK1kbT79/aH5X9mjPVOO7x2FWkTKvtl/lLumLoTE9B4dZhVI/nH2W0AAHh1nlZl3gy771Zpm8jcOsBnm2e9QVDrP2NsA9DrNalXHT8zKCpUhogIPDrHbjztKt68boZfCqilkDmrxwUNHeBjg4zFotT3iUERgDXlUx64EyQwMcHGatN2m1p+zQAAeHWSfUvbHUwWHWEXWGh4PDrEOWhkcVqyng4DDrnEF4VBUa4OAw66QqQwMcHGY2AweHmeXm4DCz3BQRTdeApHeB/wF+3nQtOVyH662S663WuHr/JCKun/TiJIIDQNKRiJhruo5pud5qud5qFa3XqypmlpuDw8xySyk49jZdQE6ut1qut1qF6k1mG4eZtUdKPQ4za4nGg0PS3ZJOSjolaWfT9Ywi6W1Jr0l6WdKRbNi1kp6T9Eb2e2WD9T0p6ZykY0PDxtYnaVfW3icl3ZVIvQ9L+mnWxi9LujehetdK+r6kE5KOS/piNjzJNl6m3vLaOCIa+wFWAD8GPgVcBrwCbGyypjF1vg1ct2TYPwI7s8c7gX9osL7bgE3AsUn1ARuzdr4cWJe1/4oE6n0Y+LsR46ZQ7ypgU/b4auBHWV1JtvEy9ZbWxk33ODYDpyLizYj4NbAf2NJwTdPaAuzLHu8D7muqkIh4HnhvyeBx9W0B9kfERxHxFnCKxe+hNmPqHSeFes9ExNHs8YfACWA1ibbxMvWOk7vepoNjNfDO0PMFlv+ATQnge5JekrQ9G3ZjRJyBxS8KuKGx6kYbV1/Kbf6gpFezVZlBtz+peiXdDHwG+CEtaOMl9UJJbdx0cGjEsBR388xHxCbgHmCHpNuaLqiAVNv8a8CngVuBM8BXs+HJ1CvpKuDbwJci4hfLjTpiWO01j6i3tDZuOjgWgLVDz9cApxuqZayIOJ39Pgd8l8Vu3FlJqwCy3+eaq3CkcfUl2eYRcTYiLkTEb4Gvc7GrnES9kj7G4kL4rYj4TjY42TYeVW+Zbdx0cLwIrJe0TtJlwFbgYMM1XULSlZKuHjwGPgscY7HObdlo24BnmqlwrHH1HQS2Srpc0jpgPfBCA/VdYrAAZu5nsY0hgXolCfgGcCIiHh/6U5JtPK7eUtu4zq3TY7YA38viVt8fA19uup4R9X2KxS3OrwDHBzUCfwwcAt7Ifl/bYI1Psdj1/A2L/z0eWK4+4MtZe58E7kmk3n8BXgNezWbkVQnV+xcsdt1fBV7Ofu5NtY2Xqbe0NvaRo2aWW9OrKmbWQg4OM8vNwWFmuTk4zCw3B4eZ5ebgMLPcHBxmlpuDw8xy+z+RKjtXLKLQWgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "local_maximum_neighbor_count = cle.maximum_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_maximum_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAkzklEQVR4nO2dXcwdR3nH/w+2k9hgV3kdYrm2lTeAhRJEGiwrIFHRVKgkzo0Bicq9qHwRyb0IEkitVKdcwA0VrQS9AzWIiKiipJEAJRekJbIoCNTmU8H5cIMNeSGvbdnFruqIOMFOphdnpx6vZ3bn45ndZ3bnJx2d9z3n7Jw5s7O/fWZ2doaUUqhUKpUQ3jZ2BiqVSnlUcVQqlWCqOCqVSjBVHJVKJZgqjkqlEkwVR6VSCSabOIjoTiJ6iYiOEdHBXN9TqVSGh3KM4yCiNQB+DuBPAKwCeBLAnymlXmT/skqlMji5Io7bABxTSv1SKfU7AA8C2JvpuyqVysCszZTuNgCvGP+vAvig68NEmxWwwzNpSsmXg98xpvUmY1pt1mVMW1eFi0zpce6nnGWaI23u/XR1x3tc+0vzs98opd7Z96lc4rDVmsvaRER0AMCBxX/bARzyTDrHwbPCmNY5xrTabBkg3TNMaXLup7OMabXJsb8499N7PD7DWT7X/crnU7maKqu4PITYDuCE+QGl1H1Kqd1Kqd3A5kzZqPTTruR1X1T6ySWOJwHsJKIbiegqAPsAPJLpu4SRszlRmRc+0cY4ZGmqKKUuEtGnAfwbgDUA7ldKvZCe8nrjb+62XSrnjb+1PC6MkRFPusLpzeBrsnBwHpfv+zkQIo0lpDdXzvd/xCBXHweUUt8H8H2e1GyVZi3kyaPNOvDIg7tvwyc9KfI43jyXIg+OfTV0pBEmDSCjOHjoqyhcVwKWkdZB2lXwqdEHpzRyda7auID4Zttxx+vS5TFk+baJjTrCpQGIFkdIBYmNPo4kfi8QW/B+lCoNjSlLrr4fbnlsQvqVlTGFYRIqj/i6K1QcMRUjVB42aQCXCtMnDyEF79tsGaNZ4oKzudL1202puKINE055nG99f2hkyL2/hmqmpJ3wBIojpUL4ysMlDROzYG15iin4rgoq5aw1BrosTgdswyEP2z4cS/BDkh4lCxIH1xmkTx4+0mjTjkI4midDVLzSKvf1yC8Pn33nkkfu8uSINvqaKzxNayHi4B5OYnaaxojCBVd/xhJTOi44K7geEDbEFZYQaWja+8QlktB91+6TybnPzmJxexcXOq9aIPz9cBOfj4NTGqWQ66xYyohS20GSswM7FX1wP5Eh7Xyym7g4ct43IpHcoXQp8gAWstAPqbSbFDnkkef3T1gcP22e5yKPofozSpBHjoOF+0Y7V3o55MHPRMXx09b/U5bH2eYxx2ZZqfRJSL48JioOG1zyeDRx+5ydbFUeC1L3UU58I5dUeTzTPE4lpmNnguJoRxucPNp6HhtbJTxiPOaItH1kEtrcCZXHM8YjLxMURxcpUUe7IkqsmG1yyOMx5KuYqWVa4j7qo08ew8nCJMtkxcGZoF0K+AlDSr7RxqaANH0q3x7PtDibKRyddTf1vP9Yz/u7Er//gY73fMsU6N9HIWnZWDX+viUyjZT9pcd4pMjBt/N8x9OLybW6mZA4QpsoXfIIPVP5VEzuvg3OXv62QPqEYRIrjy5paHzK1Xdfxcpj1fF6iEBS9xXXOE0fefiJY2ZNFRNbs+VR5AlvVwEcZkyP+9Kg2ScSIg0g7izoIw2gf1/kboq4pAH478+c86WGwtdROpGII6VDdBN4KmD7jNZV6YD4kBfIVxn78tyFT+ThK4w2tmghR1Ro4lsWXfuRaz9x3xnSFXkUFXH8FuNdu+Y+a63Cr9IdNh4hSJQG0B95xEoD4On0zBWduPafpEijTXrkIUQcmicQLpCcl19D8BWGjViJSMMljxRpaDgOfN80Qvdje99JlgYPQu6ObWOTh+3uQSnS4ERXwJSmTCip0YbJM7i82cIhDc0Ql1dTyuIwFiuBlMAppNymIFQcNkyZSJ49nAvzDKYlUsqZ7BkAz42diQg4BVoC8fIQ1lTxZW5rlxyG3L6N0nBFLXMrB01cf0eh4liPsEFclYpJWx5zlYYmXB6FikMzJ3kcHTsDE0PLQ7o0zmGYJurjzcOPAsXRnhpuE6YvED2/RJUHLzmkwZmmOUgxpzxWgrcoUBwuUgSyjTMjmZmrPLYZjzlgG9mcQx4rUVsVdFUlF9tazy7M9T62Gf/r7c4iz1wbttmsjgLYmeG7JGLbL+ZrPuuwDE3qJfWuu7hNeeSe9NpNYeLwmQo/ZGWulYDvbldgW4XOJQ8bOvKYskB8ogtT4iG8P2IbH8wD+zDC5REy9YP+rpg6txKxzSUKEkfI+hntJoveGSs8WemE84zgM3em2XQJkYg+K3KLjmP8RmhzpP35LpGYESLA9/tdzQjf6CNlrpj2d9t+00pC+ldSkDhS4FgfNIaUyhkz4W5XE6brngqugydVGlz9FzqdtkBcUaImtBxC+hy6BMJdN9v1boU5/dmIAxgm2nARWjlTZuluRyA+979wnH2lSKOd5vGAtH0lmtJJ2R4RnPOEtpIt5ULEwbk6+dj0HaScU/vH3HkbIg+uYeU5r5SEpt1VBtxXNcqdfb8QcUwRWxQiYfEgX3mUII1Y2nIv5R6h4ajiqFjokgfnzWsSpWFSsjDyRjMFiGNKzRQXOS7jpkYvZp5y3OkqXRqVLpLEQUQrAF4F8CaAi0qp3US0BOBfACxj0Tvzp0qp/0nL5hyQ0Expk/PW+CWUfUafNxxDzv9YKXWrMU/hQQCHlFI7ARxq/o9kDtHGnBlv5KMMyu0czXGvyl5cmvbpAQAfj0tGujSWUCt+LOaozbmWof7dZcojVRwKwA+I6GkiOtC8tkUpdRIAmufrw5PNIY1lhjSWcKUw5lrxY7EN9Z5bGbZ/b3nySO0c/bBS6gQRXQ/gMSL6L98NG9E0somf+9APvWOWET8oprTKvR3pt3jH3gfiouv+kNrnURJJEYdS6kTzfBrA97CYUfgUEW0FgOb5tGPb+5RSuxd9I7/Xevc80jsLzxkPk2X4Rx+2CKMENjfPHBPncl398LmprLRyjsH1G211NZb8EUy0OIjo7US0Uf8N4GMAngfwCID9zcf2A3g4Pnsx8vDdAcuO12NkIaXCb8YlaXCSKo+QO1GllOVYuE54vpwFcJEvOw5SmipbAHyPiHQ6/6yU+lciehLAQ0R0N4BfA/hUWhbPI19H6TIub7qUWGlziMKGtNvX54CWh+8EVWZT7yJyDtOKTlkp9UsAf2B5/QyAj6Zk6krMyMOUCEdItsyUVsrw5GXj776JY0NFwdHXoQm5fR2Il0ZsWS4bf69EfreJq6kXU54pJyWzbrYl0lVG7ciDTyQFjBxtcx7TWldlufX/FtjlkRJZcMrDxHX7OsATafR1mC57vLcS+d1d/UPme0NPdqwlEtMc0dukH/YFikM6PlcHlnveb8uDozmSSx5At0BSsZXncsD2+rMrAduEdCr7SERaE9iUTpwCJjRZsSTaFWW59fAh9yXqHGiBcPdr6PJcRvx4HN/tUq5ElbL8o0lcR2oVRzaWkFbRgYU8ODs/h6jYuW5eWx4gDY7yaachLdqwES6PKo5sSF3rpcSzImdZLjte5yyXEss4TB5VHOLJ0Q1VYsXmZBmXR4M5ymPaZVyoOKSezTXS8wcsKva0K7cfOcughGYKsDg5hZ2gChSHHsdRwsFZ4SH/SMj5MsurKpzyqCKSyYbmuUR5DDWqd3gKE4f0OTqAPAIqJeQFePO6of8jYtF1VbI84vvPChHHerilMZdIocS7dFOwSSNH1JGjTNt1VZo8wvs0bCkIxjfCSF2pzRylGRvVmCF1zmI1K3rq/BXcI0l1R+NriekMFWnokw7nXCDSo2KeuilYHKE7IEYetntCQu7GdVXw3PLQpC4cxDUMvX1lYgPC5eErC66ybUeqHGuodNWbzQDOJKStCbnJrc3kb3KLtbaPPPruPgX65eFTyXNVcB9sB4F5VjV/m15r1lw6MgTX5UyzjPokMnRfRleZxkQfvvU1RR6uPPvml/dQFyQO7hDPRxBdtCcRimmnpt6NmNp/Y0YkfTK0LVZtk0nMuAdJnZw+ZeobfcTUWV2PfATiu//78st/mAsRB2cfbaowbKR2bg3VdMmJpMFiseUZKmLXAclxkuuKPmJPGLboI0+9K702V7yInb9VRyESF4sKJSV6y7W2rykPrquDZpSZb+7RQi7HToHQS4lzucwcS0h5Si5LaZdq/ajiGJQxRj9ynCWlX2LsQrI0NCXk8XKqOAbHRx7lVSSZ1HLMRe3jGIUpdJaOTclRUPnU2ttLrjaoeTbsmsVaCutRXiepzm+VDDdVHJ3kkkZ7xKdUWUhlaiLIJeWYEbx+1D4OJ0NJIzdcB1lp0UYp5Dp36760PIPvqjisTEUaXJQuDan5N6XBGUW1O+D55TExccTeb6HJtfYqMK40UiqltIMu9rdI+x02OOThumrHKw8hfRxvjZ2BjJQaZQBlHGwh5FyHOBTXoafzF1r2Ppf5+fo8hIgDcK8P68s5uJdP9IE70ihZGFIZ+6A3vz9Fqj6Hna9AQgcV6sgjTSCCxGHSLixXhbGNxQ+VB6cwJMuixMupkoktz9BDrut7UkYipwlEqDjamAXns+C0bflE7rVY20iWRuUSnM2VkGZFyqGmv+fVhDRcxDVfChEHB1ompcxbmYOQs+SUoxPuvo7YfgkphMtjYldVxqAUaYRQ6gEQQo7f6JpUu4Tzc9hVlyqOSoUdUx4lSEPjL48ZiqOkHZmDvhB9DtFGbswyLHEhqX5mKI6SOIc8sziNfVkzhjPgX84hB/MQb684iOh+IjpNRM8bry0R0WNEdLR5vtZ4714iOkZELxHRHbkyPn3yTfu2wCYPqZXenJtzFXIFIrX8+PGJOL4J4M7WawcBHFJK7QRwqPkfRHQzgH0A3tds81UiWsOW29nQlsY55Ik+zKsBEiv9Gbgn9E0VCHfU1VV+02uu9IpDKfVjXDl18l4ADzR/PwDg48brDyql3lBKvQzgGIDbeLIKLHZ26i3o3DuRawUwTZ8cOCVyDn7jYsbAd/2RUIFoSXLuN4nSzUtsH8cWpdRJAGier29e3wbgFeNzq81rDJhnCAnzV1wwHqfAsyxDqAxS5GFuK63DOGbRopgIhEMeQ0pj44Df1Q135yhZXlPWDxIdIKKniOgp4H8jvipGHlzRhussnSKQ3H0afd+VvhCxDLrkYTvIU+QRIo1pNVdia8opItqqlDpJRFsBnG5eXwWww/jcdgAnbAkope4DcB8AEL3XKpdLdK1UP+QBB/iF9m152IbAm3BFDj4y7fsuXSXGqugc66va6DrIzyJ8IF9MpJEy1+xS6/lXkenwEBtxPAJgf/P3fgAPG6/vI6KriehGLFb0eSIti30M2WyJ7Q/gasr00df3ESKoMSIQLmnoqOM8/Dt+z8I/+hiyebIEu9RuMB7D01sziOjbAG4HcB0RrQL4PIAvAXiIiO4G8GsAnwIApdQLRPQQgBex0Os9Sqk307Lo0/vtE3mknkE5OhG1PMwIJEfEpNPkiMjWYpjogzvSiD24+6KPoaXhg5bHcFEIKdXTShgiE/ReBfyj5Z3QS2ZdB4kEcbQpbSBWLoFwS4Pj7mfbQcsljb7zdcr9T6nyuOtppdTuvk9NoTfMwHWGnVbH1Hjk6P/I1aeRSo6FpjW2vg6umyVvwBCRh+Ah57E7SsKl2kolhPLusBYqjlS7m/Ko0UZFIrpelicNQKw4ONiEKo2KbMqUBiBSHKV1GFYq80OYOOYkjXWoEVElD4y3hzkQJo65sM74e+7yyLUA1lzRY4TyykPY5djUtVVywTmGY53ltZShyJUreQ38yx5KWszJRfvWBlMevAO4BUccHHNELME9ZFcaF43H3FhnPFLRB3eOVdol3z7fdz8UbwRSwGlO76xU2y8h/k5IXaFTIw/fA0PLI3X3nIP8cS3tAXvtMgop83YdyRV5+H5/F5zRS580NFoe6dFHAeLQcAjEjDxiJLIO8fKIOZvGNGFy3y3MFRH55tMst66yd9ULHXnwr9h+JX1N7fZrr2KcOTbSBVKQODScEYgmRCIp8oghJPqwHYzmDW+ptO/wHbpj0xX5+dSFIQUCXBmZuJrLY8kDWAgkTh4FikPD1Vk1ROTB0XbvE4jPlIOaUIm4pgQ4gzh5pEZF6xBfddt9H0OJpAu9tGOMQHybKbwI7hztg6uNGNtxytWZF4qtqRAz5aDPNj7ziITepDb0xEt9vGY8xuZVhK0PyyGNuE7TgsUhhTEEYsqDa/awNiETD3XNRu77faHkuPokQR6Anzw4I41weRTcVOGCa7ZrW/s7p1C4DkKdDselxrY8xojIUslxJSYmqjXlYTZhcjVNwuRRxcGOWUlCws7KNOEYQ6Tr0XsY0uKhNlVYaXduyZnOvh/poyLnzk2QFMFVcbDhkkRJ8iiJTZBxRWRoZMijiqNSIOblZMny4LrV4abW/2Nd0btEFQcLfVFFjTryIlkeORlPHlUcyYwhBSmXDcfANXhtqvJoRxttOOVxzPuT9apKEkNLY87C8GEDplVGfdLgwl8YmiqOaEKloT8fc4l2iINB8i3jGp+h8inykBS1hEgj9O7tcFG0qeKIIiXS2Ah54zumIg1NqDwkCQOIjzS67qFKl4VJQeLYZnktddRnzPZTa55MTRoalzyGlETMYtaptOXBKwyNEHFcBbsY+kiZnCdm+yGlUYWxIGU6gLY8pEUWNjj6NdYBOMKQjpsJXFUpYVrANn0CqtLgY4PxGAOue6FCyCsNYBLiAOLlcQH+UQR3tOFKL6c0OOZx9WUzFtGCfsQgfdpDX3zlMdRVlHSENFU48Gl2uDqO+q545GqitDtKc0ljyAjDNbGPloC0+TgqMUxIHMDl8oiZ3m+MwVwbschzrukIc0ojZvYvH4FMJdIw6esoLSfaACbTVDG5gGHnBE1Fi07GzUv+pM43mtKEmRplSQOYpDhKot20KkEem8E7SXFbIFOWia0pXZ40gCqOEXH1x5Qgjxxsax5TviHQ1lTJfwUkB1UcIpEqj1zLIUiXxXrwLsfR5gj4BLIJwAeZ0nLTKw4iup+IThPR88ZrXyCi40T0bPO4y3jvXiI6RkQvEdEduTJeNj6X5yTKI8dVH5s0JIlkyJnRUuTRbvLllYdPxPFNAHdaXv8HpdStzeP7AEBENwPYB+B9zTZfJaI1XJkdF651XUMXf0phA/iihBxrsnYJYmPP+7mxRRlDSCRGHq5+oXzy6BWHUurH8K/tewE8qJR6Qyn1MhYD5XlXux2dFIHELv4USnukZKo8bGuypjLEwLvYJkbXNtLk0deZnEceKX0cnyaiw01T5trmtW0AXjE+s4q4m1AiWNs8cnydTRQXPR+alKHHXfLYYHnYCL0ash7dB96Q8174RB/rLQ/bez5p9BEipCXEjWw+gu6+j5DL2fzyiBXH1wC8G8CtAE4C+HLzOlk+q2wJENEBInqKiJ5KH8/fHsc2kKu8uAie+xVMeaTcf+EjD6kzntvkERpVuJog3JFJLmLHv/DKI0ocSqlTSqk3lVJvAfg6LjVHVgHsMD66HcAJRxr3KaV2K6V2A9cAeCYwF2uNhw0ueXD0a3CNTVgHnpu1bPLwPStzkDqfCUdeh/i9XDdgHgHPgDk+eUSJg4i2Gv9+AoC+4vIIgH1EdDUR3QhgJ4KWw36m9XDhO1JeUuQhDS2PoWQxVUorOx559B6BRPRtALcDuI6IVgF8HsDtRHQrFs2QFQB/AQBKqReI6CEAL2Jxqr5HKfVmfPZMeezyzXKLbQCOx2dh0mxGWl9FzHKJki61crEeV94TJHm6By2Px6NTIKWsXRCDQvRuBfyt56djLtKkiINrcWOuu0K5x3ekdnKOIY4cC05zoOWRQxq5Lq225fGhpxfdB93MZORobJNFagWtyEVypMFHYeJIGRJS+zvyEBKxTLGZYlKiND6ImGimMHGkEiKPGm344yOPqUujdMLkMTNxVCoVN/7ymKE4fKKOGm2Ec6HnMQdKnwTa/0LJDMUBlN3fIfEg9LnSM8Zs30OTa0xH/tvkQ6QBTG7OUV9+NXYGEuC8HMtxv0lIfsZYoGhIzqO8AWFAqDSA2YqjVE61/t8ekQbnzWkxEqvykEXcOK4ZNlVKjDZO4UppAItbg3x5DXJWcpfWbOFeb+Y8hl3DJgaFWGkARUUcu3Blp2Vo9jmkEVsZuhYEtmEThQ2bPIY4ow89Qxl3h7VtP/rs29Boop1myPbvb55T18DlHx1eiDh2OV7XlWmon5F6BumTh68suihFGj5NllxXt8aMBPR39wnk/T3va5n0CSTPLSUFiMMlDROfSERKE8UmDw5haHL3IXBGGmcx/HIIUpoPZj7aEumThoktGsnfxyJcHD7SsJHjTCWlwo0Jd/NkrtJooztUQ4QxLsLFMVXMqIMz2siJxFnXfZEqDJNypAGIvqoSG22UwjpUaQxBCdLgJv9vFiqOqUujNHIIbohmyhylMQxCmiprAWwx/tcT75Q8NHwKmOMttDy22D4oCLPjeS3me99R3oFoQiMOTZ3yTx6Sm1e2S91Czo0TQ7g4KuPRNbrTNZJ1TLrGx1R5cFOAOGrUMTzShoT34TMit8qDkwLEIYEcnWyPo3wpxkYdQ4/f0FR5cFGIOEo/wNqYM0tL+22h0cbYTRaJ85OEch7Aj8bORBCFiANIO8DOQs6tzvFrWeQntokSKo96mfQSZllwyyNfORckDiBcHmdx+cEQI4+jCLt9vQuXNI6DJ/LQ6TzXPEJI7dfw7TDVg8lKlAf3ycdWBmVEHhNp9IVUer3zXRX3qON1mzx8J9LxjTJMeYSMYXFJp0se5hBnzs5Q23gP18jT1NvYY5opusrHjO/oujEtJg0XPwLwR5Hp276PP9ouUBzHsTioUiu7WZjn4RZGF6ZMXBKJbZro3+nzuRi0VHIPsksdrm470NYjvW8jdXBY6AFZYoTlprCmiob7cmGMNNrYIpLU/owuKXA1b3JxCvnuceHqEE09b+aUAVOT5ZoNwDV0+YOBQsUhlVVcEghXJ6hNDpKFMTd8pgmMFUykPK7ZcOlhfT9dIlUcWeC+cmJGF1UaZZEalQTKwyUL5+fj5FFgH8ecqdKQje73GLg/I1QWV2xvyON1v01qxFGpsFKYNCKpEUclE7p5NcQqZFIw5//kPKBbExznkMXrYZMa14hj9uRu/kgeKZsTrjVsbjP+zrRWS6A0gCLFIWXoeMWO7TKxZHlwjQq2iSJVHrf1fySVCGkARYrjeub0DjOnN1f6xpVIlMdK85wqjy5BxMqjQxqvn4lMk49ecRDRDiL6IREdIaIXiOgzzetLRPQYER1tnq81trmXiI4R0UtEdEd6Nm8wHgCfPKo0FhxH3ICy0O0kyGPFeJhwRR4cDBBpANHRBuAXcVwE8JdKqZsAfAjAPUR0M4CDAA4ppXYCONT8j+a9fQDeB+BOAF8lojXhWWvLok2KPA6jSsPFccSLxIdUeazCfuBzECMPn4giJOrwlEZq1JEgDcDjqopS6iSAk83frxLRESxucNgL4PbmYw8A+HcAf928/qBS6g0ALxPRMSxK4z/8suQShY3rAZwO+HyVRTiuG+9SpPI4wq+22A7qleZ52TONld5PhMG5iPdAUQaQLA0g8HIsES0D+AAWe35LIxUopU4SkQ4BtgH4T2OzVXjdSRUiDBNfeQwljSkP0hrjt/lEASvN83LHeyHf53PXM6c0ysNbHET0DgDfAfBZpdQ5IudQVdsbVyiOiA4AOLD47/d9s+FAO8sUyNDRxZSFkQPdZGlHHil9DSsJ27ry0JZIrDBew5VjOwaMMgCWSEPjdVWFiNZhIY1vKaW+27x8ioi2Nu9vxaWjdhXADmPz7QBOtNNUSt2nlNqtlNrNt0iyFkhtkpSDFoh5g6AkzDwV2jQBWKUB+F1VIQDfAHBEKfUV461HAOxv/t4P4GHj9X1EdDUR3QhgJ4An+LLcxxjSqNFGGhKFYcI91qNsaQB+TZUPA/hzAM8R0bPNa38D4EsAHiKiuwH8GsCnAEAp9QIRPQTgRSyuyNyjlHqTO+MyqMKYDzlGgiby+hngms186QXgc1XlJ7D3WwDARx3bfBHAFxPyVQBVGpUCyBBtAGJGjpZ2r12VRkUIXeM5MkkDECMOABgn5AqnSmOe+E5MLYSM0gDEneq1PMYfi1+pLLjF+FvLQ3BnbmZhaARFHJX5ckv/RwbnFrjztR0iI5CBpAGIFUdKs+WTzSMHtwDYkyntHJSQV51HSfLwzYsQebx+ZlBpAGLFAaT3eeSSB1DWASmZdh4lyCM0D6Hy4F5rd3PzILgvfvIjWByAvzw2GQ8TLnnYwlbJB+Yex9/caedIR4I8QvFtutzUPJ+C/5KZLrQw2gwjEFJq2BDHmgn6gHJPA9/VUdoWhYvv9n/ESV9FfjQh7Rz0Hdix+W2ny5WODY7Rv+391pcml7BcHac3OV7vY4vltdBoPOQYf9vTi9tAej4VmIMRsBWSLbroIjby8KlMkiIPn7yE5HeP8UhJJ3Sb1IPYtn1XmpxRji3yiJUGcCkq0c8yhi0UEHFoziBMFi58o4/QypQSeexh2D6U9vdxpNFFTPohkYfv/mqnmatptDFDmjcnbOt7nPtFHAWJg2u9UI1LIBwVqe+A8jmIfA9KKREPh4hcDHWwc5FDGkCaOGzYjv0qDg/a8uCskEOc0aVIQ6PzKy1fQ1OKODSmA6o4PNHyyHUWexTpB5JNIHM/OCUjrZnig/bAZDpHc/NJ5A19OQ7wdhpVGhVuwi7hVnEUw57Ws0TWoy6YVTL+8ihEHLmaKcBiST2PuZRFIF0atr8rU6QQceTCXIezFHlIJLcoaiQjDWG31dvIEW24Fu7dhjrfRghdB7N+L2WR5Hb6HGm20+NcxFnn72LzXMDh9f+ETWUxw4ijr6LUyEMGQ0Qx+jnXd13s/4gX2wGcY0rLRvj8N8LFkbNvo4sqD9mkHui27VMF4tqWSx5AHnnETZolXBzchISlEuUhqZ3vm5eYPPtsk6ss1rcevtt0ESsP2123OSMPf0pqhM0cM7TmbJfHEHrQ5srzEGVh/lbbd/mWRUi/x1ATBMVP0TmjiCOmgkmMOoBxrzLEfm+uCCX3513bxu6DvujDRxocUUfavL5CxXG+eVwETxsx5ayUKg+OA9yVRmmXKfvyWsJv4ShzW50Oncd03CaLMHFoYbS5iHiJcISyMfIwK5i0cQ4x7Xi9XS5ydExyfgc3Zp2ObZrEyiN9FQEhfRxvwf8A95HH2oD0fOkb4xF6NvXJ31BheN/4CM4Drv1dXGl3/QZJwuDmHPznqeFbdkRYxMEF5yUwH2KvHPgMoBoSW56kRhlDpjkFeNcqmqg4cmE2WWJCfRtDH6w+cP22vu8YgrHLcgj6miz8C5xVcQSzDXnPlHOo6LmZY1kO21k6YXGsGzsDEcypoudmjmU5nDwmLI5KZY605ZFnHeYqjkplcmh55Fu8vYqjUpkk+aQBVHFUKpUIqjgqlUowExdHjisrZzOkWamUhZB1Vei/AfwWwG/GzksA16HmNyc1v3lx5fcGpdQ7+zYWIQ4AIKKnfBaCkULNb15qfvOSmt+JN1UqlUoOqjgqlUowksRx39gZCKTmNy81v3lJyq+YPo5KpVIOkiKOSqVSCKOLg4juJKKXiOgYER0cOz82iGiFiJ4jomeJ6KnmtSUieoyIjjbP146Yv/uJ6DQRPW+85swfEd3blPdLRHSHkPx+gYiON2X8LBHdJSi/O4joh0R0hIheIKLPNK+LLOOO/PKVsVJqtAeANQB+AeBdAK4C8DMAN4+ZJ0c+VwBc13rt7wEcbP4+CODvRszfRwDsAvB8X/4A3NyU89UAbmzKf42A/H4BwF9ZPishv1sB7Gr+3gjg502+RJZxR37ZynjsiOM2AMeUUr9USv0OwIMA9o6cJ1/2Anig+fsBAB8fKyNKqR/jyiGtrvztBfCgUuoNpdTLAI5hsR8Gw5FfFxLye1Ip9Uzz96sAjmAxo5PIMu7Ir4vg/I4tjm0AXjH+X4XMxUwUgB8Q0dNEdKB5bYtS6iSw2FEArh8td3Zc+ZNc5p8mosNNU0aH/aLyS0TLAD4A4HEUUMat/AJMZTy2OMjymsTLPB9WSu0CsAfAPUT0kbEzlIDUMv8agHcDuBXASQBfbl4Xk18iegeA7wD4rFKqa7otEXm25JetjMcWxyqAHcb/2wGcGCkvTpRSJ5rn0wC+h0UYd4qItgJA83x6vBxaceVPZJkrpU4ppd5USr0F4Ou4FCqLyC8RrcPiIPyWUuq7zctiy9iWX84yHlscTwLYSUQ3EtFVAPYBeGTkPF0GEb2diDbqvwF8DMDzWORzf/Ox/QAeHieHTlz5ewTAPiK6mohuBLATwBMj5O8y9AHY8AksyhgQkF8iIgDfAHBEKfUV4y2RZezKL2sZD9k77egBvguLXt9fAPjc2Pmx5O9dWPQ4/wzACzqPADYDOATgaPO8NGIev41F6HkBi7PH3V35A/C5prxfArBHSH7/CcBzAA43FXmroPz+IRah+2EAzzaPu6SWcUd+2cq4jhytVCrBjN1UqVQqBVLFUalUgqniqFQqwVRxVCqVYKo4KpVKMFUclUolmCqOSqUSTBVHpVIJ5v8AjG9xcrAxvYYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "local_standard_deviation_neighbor_count = cle.standard_deviation_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_standard_deviation_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ4AAAD8CAYAAACGnEoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAZTUlEQVR4nO2dbcwcV3WAn1MHUj4Vp/nAcUJjwJXs/GhwXqWVXKG0KCREit5YSSsjubVQJPdHkECiUp0QBFKcmFaC/gPJiAirpqRRTOKEIkpkUUVYLcG2AsR+a2JI2hhbMRRHIBChSU9/vDN4vJ7dnY97Z86dOY/0andnZ2fPe/fOM+d+zIyoKo7jOHX4nb4DcBwnPVwcjuPUxsXhOE5tXByO49TGxeE4Tm1cHI7j1CaaOETkZhE5JiLHRWR7rO9xHKd7JMY8DhFZAfwAuBE4AXwH+ICqHg3+ZY7jdE6sjON64Liq/khVfwM8BCxG+i7HcTrmgkjbXQ28WHh9AvijqUFccpG+7uorKm3418ff2C6yEn73Xb8Kvs0YrORMtG1fxMsAPM+aINtbw/NBtgPhYuqK0L/TGVYG3d4sfn1o6aeqeum89WKJQ0qWndMmEpFtwDaAC97+NtYc3FNpw0u3bmgd3CRrnjgcfJsxuJ1Homx3E49BVjm3UO13mMcetgTZDoSLqStC/k57uYOLgm1tPkty3X9VWS9WU+UEcFXh9ZXAyeIKqrpLVRdUdeGCS7szqnMuy9I4S8gd3hkuscTxHWCtiKwRkdcDm4HHI32XKXzHc0Kxlzv6DmEqUZoqqvqqiHwI+FdgBfCgqh4Jse11WbMiRpOlDUVh5M8tp9iTmUaRPWwxFbu1eLqga2nkzasdFdePNo9DVb+mqn+gqu9U1ftDb39dAv0SVrOPWdLIsRa7tXhmsZc7Wu34fUmjDknPHF33xOEgAmmbvcyq1HvYEqzSt61Qm3iskjSskpI8oL1AuqBpR26sUZVOWffE4UY7/6En1p+3rG5KnEpl7lIYMcskxWbLXu6INiLWhjYxRZk5Wpc3LKzXqsOxs6gjjzJpFKlSOevuIKEqfJ0fvK0wrIo01rBxTClV+d26ylCmxbJDdh5S1YV5n0+6qTJJ1WbLPGnA2SbGtB2hyQ4yb5tVsZr+dpl9hfiusm3E/B+sNF1CZD+DyjhyZmUeVaQxjfxoFKJyhTiyzasAIZonMTKvkMTKiLpqDhV/wy6kMq/OVM04BimOnKVbN7QSRSxiN1li9GeUxWyxf2da2baNNXYTJqY06mQYo2yqTGJRGiEpq2yxOkEtSqKMrpsfbcl3aoudp7MYxKiKs0zskZPiiIblnRHsxwfnyyJG5hFLSIPNOFKoOCHparg15LyUWMSIL/Q2p+3QqWQegxSH9YodktQndY2ReXJIQR6DFIdlUpu85ISlqhTayuN2HokqoMH1cYwp2xjT/zoE6u7Idfs8usxUBicOxxkS8+TRV7NmUOKwfgQO2Uyx/r8OkTbnyYTYwS31fQxGHNZ3JO/bGAZNrrUSor/CGt452hEhxWZdkmOg6m9gcacPwSDEkcqOFOokN8cG837HoUoDjIhjJWcGXchlNJWIS8cW036PoddnU30ceWHXGYIawo5U/B+8LyQ9Jvs9hi4NMCaOnLKCL5PJEKQxSQoXOnbK2cMWHuW2vsPoBJPiKKMokzFMsS7LQoYoSidNkhFHkUe5bRTyyHFhONYw0TlalzFJw3EskqQ4HMfpl+TE4dmGY5mx1M/kxOE4Tv+4OBwnMGO4uFJS4hj6j+EMiyHX12TEMeQfwRkuQ80+khGH46TM0ASS5AQwx0mVojxSnp6eRMYxJFM7zhBIQhyO49jCxeE4Tm3Mi8ObKY5jj1adoyLyAvAL4DXgVVVdEJGLgX8GrgZeAP5CVc+0C9NxHEuEyDj+VFWvVdWF7PV2YL+qrgX2Z68b4dmGM2RSrt8xmiqLwO7s+W5oNuaUcqE6TlVSredtxaHAN0TkkIhsy5ZdrqqnALLHy+puNNXCdJwmpFjf204A26iqJ0XkMuBJEfnPqh/MRLMN4K1vf2vLMBzH6ZJWGYeqnsweTwOPAtcDL4nIKoDs8fSUz+5S1QVVXXjTpW88571HuS3pWXWOU5fUpqQ3FoeIvElE3pI/B94HPAs8DmzNVtsK7Gv6HS4PZ2zkAmkqka4E1CbjuBz4loh8F3ga+BdV/TrwKeBGEXkOuDF73RiXhzNW6kqguG5seTTu41DVHwF/WLL8f4D3tglqkqI8UkrnHCcEs06Mm7U/TL4X8iCc3NmxY7s1guMUaVP388+GEEhy4nAcpx0hTu03f66K4zjxaJrBuDgcZ+Q0kYeLw3Gc2vJwcTiOUxvvHHWckdOkgzS5jMOHYh0nHD6q4jjGGPKs56TEMdZsY8OBpb5DcGqSS8OyPNrElkQfx1iFUSSXx+GN63qOZFhsOLAUvEwnd0hrs50HP3PUUmFboZh9uETakZdlSHlYzjAgXHxmxeHSmM+0JowLZT6TZRcio5u1U4bKOuqc5Dbvs20wKQ6Xxlma9G+U7QTFo2pZBRpTmc8q0ybZR9Udso08pn1H1W2GzoREVYNusAlXLKzSrx5c2XcY5gjZKfrxjXfX/swQZVKnTOcJpM3OGGNnn7bNOtvZITsPFe5YMBUToyoX8XLfITgjoK6INxxYKv1MiEtbzmvWNNl+2Wdi9bmYbKo4YWmSbcDZSjeEzKNN9pZ/tmk5TqPYzAi1g3fVOWsi43DOx+duhMNyWVofhZmGi2PghDhKplq5wbY0clIsXxeHQVKo7Cng5RgP7+NwkqTKUXoDLo5YeMbRE3nP+eQOYPUomWI6/fGNdwfv0HSWcXH0QNm5DI9ym1lpWCNFic0ixVErF0fHzKr0MY6O9x3YGWQ7KVbuFEi1XF0cHZLqkTLVyp1jtbnS5Z3XQuPi6Iiq0rCWdVir0E3la1UeRayV9SxMiONlLuo7hGg0mT5spZKnVJGrYKVcYXrZpnLXejPDsUO7P2yqzRLL9F2mxcytjYSq1O/QU9FDYyLjmKQ4VGm14MqwHHOoTlJnmablWfegaPUgalIck1jeIXNCx2YprR4SIcv1vgM7KwukTRPEYvMlCXFYJ5bQQsujzlHSWkUNSYxyHVtG5+JoieUsqClDlkZOrNGrMoEMsTxdHI4TmKI8higNcHGMjnkp9VArepcc3rjOZL9ESFwchonVbk61PZ7CjjiWK8zPFYeIPCgip0Xk2cKyi0XkSRF5LntcWXjvbhE5LiLHROSmWIEPndg7d6ptcctH8rFIA6plHF8Ebp5Yth3Yr6prgf3Za0RkPbAZuCb7zGdFZEWwaEdCVxlB/j2Wd8ZptI05dBmPSRpQQRyq+hTws4nFi8Du7Plu+O3QwiLwkKq+oqrPA8eB68OEavOIGDqmLpsReVs8ZeoKJF8/5I4+NmlA8z6Oy1X1FED2eFm2fDXwYmG9E9my1liu4KGO2F1LY0g0+Q1ClEGX5WhpHwjdOSoly0rv+CQi20TkoIgc/OVPfhU4jH5oIxCXRhhmlX/Ze23KYsjlOI+mJ7m9JCKrVPWUiKwCTmfLTwBXFda7EjhZtgFV3QXsguU7uc36MkumrcJkvPMmibk04jOrDh3euK721de6LkdrEw2biuNxYCvwqexxX2H5P4nIZ4ArgLXA022DTB0rZzqORRqbeKz2fVrzsqkikC7Lse86M4254hCRLwM3AJeIyAngEywL42ERuRP4b+DPAVT1iIg8DBwFXgXuUtXX2gSYWrYxizKBxM42xiKLSZrWm3nZh0tjGTM3nb7z4AfPWz4kaZQR++LEY5VGCMp+m67Ks09hVL3ptJkL+TiOJcbehzEPs1POh55tOE5OatIAo+JwaThjIUVpgFFxOM4YSFUaYFAcnm04jn1MiWNs0vBRDycGO/Y9EP07TIljjLg8nJDk0ogtD1PDsUO7t0pVmkx5dpxJJmVRfH3v4j1Bv8uUOIrkEqkrkGkdTtZFVMw8XCJOXeZlGDv2PRBUHuabKqHup5JSD/bhjeu8CTNwQtbHqs2SHfseCNaEMZtxTFIlA5n3Y6TWFErhrE2nOlYOXrk82mQgyYgjZ1IgTX+MVCRi9axNpx5WpFGkTfMlOXHkhPohLEujyDyBuDTGSRdDr2WY7+OIjcUjwSzKBOHSGCchpNF0G6MXR4oUReHSGCchM40m2zJ9PY4uSKWp4qRLHyMojblN/HocfVDsg/BswElKGjXwpkpAJjsufSKXE5LQsz/b4OIIxDRJuDyckFiRh4vDcSIRa8Tu3sV7eheIiyMA87IKzzqcGPQpDxdHS1wKTp+ElEedbfmoSgtcGs4QaCIfF0dD6kojX9+HaJ3Q5Dt+1eHaEFmKi6MBbTKNDQeWXB5OFO5dvGeqPEZzIZ9J9nLHectu55FW22wya9SbJ05V8nvYdsmkPGJ1oJoQxxlWlophHnu5o5U86t6Y2KXhpEAXoy3Jj6o0EU7fuIDGw1DPhUpeHNBOHlVTydA7u8tjPAxRHoMQB7SXxyyBxNrJXR5Oqpg4rf4NC+t1zcE9QbbVtsN0kvsO7Ay6vTJ8lGUcpHDRqB2ys9Jp9YPJOFLFpeGkiItjDh/feHe0bbs0xsWQ+jpMDMeOEZfGuYzhAkgpNFWq4hmH0zvWO4lD3BSsK2mEvOnSLOaKQ0QeFJHTIvJsYdknReTHIvJM9ndL4b27ReS4iBwTkZtiBd4loZsrQz2iNqFMGpZEkkqWMCmM2PKoknF8Ebi5ZPk/qOq12d/XAERkPbAZuCb7zGdFZEWoYIeAS+MsswSx4cBSrwIpyzKsSmSaJGLKY644VPUp4GcVt7cIPKSqr6jq88Bx4PoW8Q0Kl8ZZqkqhjTyaNjFmfcaaPKrcbDoGbTpHPyQifwUcBD6qqmeA1cB/FNY5kS2LTuj5G0VCzOVwaTSnyiUJqu7sbe49PLlelVGSWKKpI4TQd6qH5p2jnwPeCVwLnAI+nS2XknVLZ5iJyDYROSgiB1/9yZmGYSwTUxohcGmEoSz7qJtVTGuChM5MYtG08zN05tEo41DVl/LnIvJ54KvZyxPAVYVVrwROTtnGLmAXLN+Q6XYeqTVtvCtZhMg2/Boc59L2eiYhOqu72Omt3VMlZObRKOMQkVWFl5uAfMTlcWCziFwoImuAtcDTVbd7O4+c8zdrPcexgLU+j3mEyjzmZhwi8mXgBuASETkBfAK4QUSuZbkZ8gLw1wCqekREHgaOAq8Cd6nqa02DKwoiz0ZcGmljaag1FGXXdbEslFwebbKPueJQ1Q+ULP7CjPXvB+5vHNEUXBiOZYrysCyNIm2aLj5zdAZdnBnrDItUpNEWF4fTGUNsphRJURpNR2lcHFPwbCMsQ5dG6tSVh4vDcRygnjxcHCV4tuGMkXWLhyuv6+JwOsEnwDWni9Pk60gDXBylxLzq15hxedikrjTAxTEVl0ccXB62aCINcHHMxOURB2vy2MRjg7oeaBXWLR5uLA0wcnuEKxZW6Z0HPzhznRj3jq1D3mFqrdKnStvh2bZSbyoKK3M16vR71BHEklxX6fYISYhj3lmzXQlkbEel2PQljza/oxVx5MwTSN2sYjDiaHKHtlgicXGEp2t5tP0NrYmjjDZ3NawqDtO3R2haADGaNS6N8ISYTXrfgZ2V5eG/YTi8c9RJnnkT9sbY+Rkbs+Jok2454+O+AztLBTJGYXTR52dSHC6N4RPrpLeiPMYoja4wKY55lw50nGkc3rjOmybEzzpMiiPH5eHUwefYdIdpcThOVVwa3WJeHJ51OPNwaXSPeXFYYOzt5dD41cDSJwlxeNbhTGMI2UaKnblJiAPaycPFY4sh7OyhKAojtDxi1vtkxAH1CyLEsG5qR4JUSFEeoc9TKatbqdQ30+eqVCWkWVP54YZAFXnM6g/pWj7FutFUIlXq1yYeCyapuvdkrkpy4sgLInQa5sKwSZkcLNzEu+7OPbT6lZw4IHzbbWg/6tDpWxo5ITODWNsuu15Hmyt/5SQpDsexQpXmS9MDU1N5zLu4z9K+Db997tccdRyDtM1m636+7q0UihKpg2ccjhOIPEPouunb9r4rTeThGYfjBCQ1aTTFMw7HSYwYssj7OqqeDOAZh+MYpouLI/ud3BxnQFiVBiQojtBtSJ/D4VhkljTuXbynu0CmMFccInKViHxTRJZE5IiIfDhbfrGIPCkiz2WPKwufuVtEjovIMRG5qW2Q+dmD+U4eamd3aTgW6ereLW0mglXJOF4FPqqq64A/Bu4SkfXAdmC/qq4F9mevyd7bDFwD3Ax8VkRW1A1sUhZl7zclxdOYnXFQVRpts462s0fnjqqo6ingVPb8FyKyBKwGFoEbstV2A/8G/G22/CFVfQV4XkSOA9cD/14loDo79NjPF3CGQ5d3iAsx5bxWH4eIXA28G/g2cHkmlVwul2WrrQZeLHzsRLZsJk2zgKqfcWk4Tjgqz+MQkTcDe4GPqOrPRWTqqiXLzrtBrYhsA7YBvO3tFwArJ1epTC6ForVdFI51ur4PbYhMI6dSxiEir2NZGl9S1a9ki18SkVXZ+6uA09nyE8BVhY9fCZyc3Kaq7lLVBVVdWHlpmHlooTtPHScWKUsDqo2qCPAFYElVP1N463Fga/Z8K7CvsHyziFwoImuAtcDT4UKejUvDsU7q0oBqGcdG4C+BPxORZ7K/W4BPATeKyHPAjdlrVPUI8DBwFPg6cJeqvhY8csdJkJDS6HM+R5VRlW9R3m8B8N4pn7kfuL9FXI7jBCBGtgFGZo4uvXxN3yE4TpLMyjpiSQOMiAPgun1H+w7BcQZDTGmAsdPqc3kcWlzfcySOkyaxhZFjJuNwHKcdXUkDjIqjTbPlun1HvdnjjIp7F+/pVBpgVBzQvs/D5eFYJPTNkXbc+gA7bn2ApVs3sHRrswsPN0FUz5sN3jnyrgXl0wdL36vS3zFLEt5f4lhhC3vOW9b0HkE7bp19+cB1TzTLQJbkukOqujBvPVOdo3WpklVct++oy8PpnTJpwPwMpEws86QBnJN9NJXILMw2VXLK5FC3H8ObLU6fTJNGFXKx5I9VpNEF5psqOYcW1wcRgGcfThe0kcU02vRhVM06BtdUCZU1+FwRJzYxpNGWMum0acKYb6rEwpsvzthpMxIzWnGAy8NJh5hDrU22PWpxgMvDcaC+PEYvDiccW9hjsn3vVKOOPEYvDu8kDUNRGC6P4TNqcbg0whBbFJ7J2COZ4diQuDDCMGtnzt/bw5Zg2w+xzcnthdpWvr1UyYdml6bevOBcRpdxuDTSoIssJn+0vsPHPnmtyXyO0YnDSZ+2O3rZ59sKxLp8ptF0EtioxJF6tmGpclaNpUnMVT4TqyxygdQRSaxYuj5Vvg6jEkfKFFPrvqkbQ8ydPDbzJBLrnJQuhOFTziuQerZRpM92edPvjXX07lJik0PO1k5kq0PbU+1HIY4+pRGics062lnIQKoyL9YU/peYwrDaLClj8MOxfUlj8ugUcthv2ndV/Y7Jit/0cyFpmwlU+R+siqlrYYS4sE8y1+OoS9fCqFspLVb0aTHFiCP/rtDbLvsfrAojx1LzpOr1OEbRVIlN05GDKhOouqQsJqtZRpfbHAKhLx84+KZKLEJV0LJmRt+Vv6vRii7ouywtMMprjlok5pHSK3p7vCzjM1hxpHidDa/o4fCyjMtgxeE4TpxmCrg4HGewxJIGuDgcZ5DElAa4OBzHaYCLw3Gc2gxaHDFGVry33nGsTDkX+QnwS+CnfcdSg0vweGPi8cZlWry/r6qXzvuwCXEAiMjBKnPkreDxxsXjjUvbeAfdVHEcJw4uDsdxamNJHLv6DqAmHm9cPN64tIrXTB+H4zjpYCnjcBwnEXoXh4jcLCLHROS4iGzvO54yROQFEfm+iDwjIgezZReLyJMi8lz2uLLH+B4UkdMi8mxh2dT4ROTurLyPichNRuL9pIj8OCvjZ0TkFkPxXiUi3xSRJRE5IiIfzpabLOMZ8YYrY1Xt7Q9YAfwQeAfweuC7wPo+Y5oS5wvAJRPL/h7Ynj3fDvxdj/G9B9gAPDsvPmB9Vs4XAmuy8l9hIN5PAn9Tsq6FeFcBG7LnbwF+kMVlsoxnxBusjPvOOK4Hjqvqj1T1N8BDwGLPMVVlEdidPd8N3NZXIKr6FPCzicXT4lsEHlLVV1T1eeA4y79DZ0yJdxoW4j2lqoez578AloDVGC3jGfFOo3a8fYtjNfBi4fUJZv+DfaHAN0TkkIhsy5ZdrqqnYPmHAi7rLbpypsVnucw/JCLfy5oyedpvKl4RuRp4N/BtEijjiXghUBn3LY6ye2NbHObZqKobgPcDd4nIe/oOqAVWy/xzwDuBa4FTwKez5WbiFZE3A3uBj6jqz2etWrKs85hL4g1Wxn2L4wRwVeH1lcDJnmKZiqqezB5PA4+ynMa9JCKrALLH0/1FWMq0+EyWuaq+pKqvqer/AZ/nbKpsIl4ReR3LO+GXVPUr2WKzZVwWb8gy7lsc3wHWisgaEXk9sBl4vOeYzkFE3iQib8mfA+8DnmU5zq3ZaluBff1EOJVp8T0ObBaRC0VkDbAWeLqH+M4h3wEzNrFcxmAgXhER4AvAkqp+pvCWyTKeFm/QMu6yd3pKD/AtLPf6/hD4WN/xlMT3DpZ7nL8LHMljBH4P2A88lz1e3GOMX2Y59fxflo8ed86KD/hYVt7HgPcbifcfge8D38sq8ipD8f4Jy6n794Bnsr9brJbxjHiDlbHPHHUcpzZ9N1Ucx0kQF4fjOLVxcTiOUxsXh+M4tXFxOI5TGxeH4zi1cXE4jlMbF4fjOLX5fy/dRoucgZu2AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# most popular number of neighbors locally\n", "local_mode_neighbor_count = cle.mode_of_touching_neighbors(neighbor_count, touch_matrix)\n", "\n", "parametric_image = cle.replace_intensities(tissue, local_mode_neighbor_count)\n", "cle.imshow(parametric_image, min_display_intensity=0, max_display_intensity=10, color_map='jet')" ] }, { "cell_type": "code", "execution_count": null, "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.9.9" } }, "nbformat": 4, "nbformat_minor": 4 }