{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Count touching neighbors\n", "\n", "In this notebook we visualize which cells in a tissue are neighbord by drawing a mesh. Furthermore, we can change the intensity of the mesh so that it corresponds to the distance between the centroids of the cells." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import pyclesperanto_prototype as cle\n", "import numpy as np\n", "import pandas as pd" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Test data\n", "Let's generate some tissue-like structure consisting of cells which typically have approximately 6 neighbors." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAccUlEQVR4nO2de7BV9XXHvwuqTIvEQASLPIRYogI1kd4QWxUTHVRQg3lVqI34GGkyxKpJRnEcElsGY2zFdJz4wJHIpFFrRlDbBBpH09CgUTFR5FEeRo08BBNNFJL4gNU/7t6Xfc/d+5z9+D3Wb+/1mTlzz9ln37MXR3+f/V2/3z7nEjNDUZTm0s93AYqi+EUloCgNRyWgKA1HJaAoDUcloCgNRyWgKA3HmgSI6Ewi2kREW4lonq3jKIpSDbJxnQAR9QewGcBUANsAPA1gFjNvMH4wRVEqYSsJTAawlZl/yczvALgPwAxLx1IUpQJ/Yul1RwB4JfF4G4CPZe1MRMbiyLHve7+plzLCH0cN911CR/a//Ubmc2+9+a7DSvIxZsAA3yW05Y39bZ7cu9dZHa28+Nvf/ZqZh7ZutyUBStnWa6AT0RwAc0wf+J4TTzP9kpVYf/M1vktoyztbvt9xn0ce3emgkvwsHXuU7xLa8sDe9sNq/5qnHFXSm/OXPfRy2nZb7cA2AKMSj0cC2JHcgZkXM3MXM3dZqkExxNTT5KcZKXQSAAD065rsoJL82JLA0wDGEdFYIjoYwEwAD1s6lmjOmTih5yaNPCkgRkXQmTwCiOnXNVmMDKxIgJnfA/AlAP8NYCOA+5l5vY1jKe5QEZhHggisXSfAzD9k5g8x81HMvNDWcSTzwUf+vddjSWmgSApIoiIwj28R6BWDjpEkgrKoCPpSpBVIw6cIVAKWaE0BSeoiApWBWXzNE6gEPOFzsrBsK5CGiqB6CmglloErIdROAgv/8BPfJRRC6spBEVQE9nAhAiufHShcRIUrBj/78cNSt1/7p6eUrqcq7VqBIvznOvMLKiZTQBouLiySdrGQ6STQibIXG52/7KFn0q7LcVu9IbIGvgRMCQDoO3dQVQq2BQB0pwJpVxjWjWQ6MHH1YVDtwGc/flhuAYTWFuShSuvgQgAxNicNpaUA35hoF8QmARNn+4V/+InXtsAWsQhstAsmiUWgyUA24pJAkbN9HlwmApOtQB7yJgOXKSANU6ngkp9/2sjr1I2qaUCMBEwP/iYRwgpDVRHEAlj1wIdNlFM7qohARDsw+BC7ZdS1LWjlnIkTRLcIZScNWxNAlgimfOa5UnXVhX5dk0tNFIpYIhwy6CA+rev9To5lSwauW4FOxDLw3QpkkUcGJuK/CzG4XiLMQ5oMspYIGycBwLwIpAkg5oHl1/kuoS3tROCy/68iCokCiGkVQZYExMwJKIov6jrPkHeeoJESqOM1BCEi6XLjVQ98uOeWF8kpICbPZxDk/yssYWqyUGorEAppk4W+lwI7iSDECch+XZOBZQ+lP+e4FlFoIpBBMhH4FkAeiiYG6TQ2CZhAU4A5pp42HKP/5a99l9FIGi+BOA3ErcExy79R4Le7/7TCO4eM6rCfe8ZPfAyYPwULFqzyXYoinEa3A0mOWf6NggI4wMF7Xum8kyfmz5/iuwRFOCoBQ0gSwfiJj/kuQQkIlQCA76583MjrSBBBmgA0DSivXrE58zmVQI3QBKCk0U4AgErAOL7SgApASaOTAABdHTDWCiSJReBi1SDv4J+vKwWNIc/AT9J4CdjEhgyqnPHjuQGVgRmGrpyA1870/9HtooO+FZWAAw7e80olETQh6uuFQsWpOvhjVAKOKCMCW4NfWwNzuE4DpgZ+kkZLwMZ8gAlcnPmliOC1gd2fcHvtun34q+v6e65GFjYGfBqNXR3wIYB2KwfjJz7Wc3OF7+sHYgHEPHPdPk+VyOLVKzY7EwDQUAn4TABpIvDZ87sWwWsDJ/fc0ghRBENXmvmSV9eDP0ZEOzBmwHtYMvp1AMDFvxriuRr7xCL4ixO2eK7EPlmDvR3PBNgaFJ0b8DHYsxCXBJaMfr1HCDaQOg/gExtpoN3ZPg8hJYKzNi7AWRsX4MKbl3Xc19fZvh0ikkAaRUTQhPQQAlUGfRrSE8FZGxf02Xbhzctw95XpX4wibfDHiEsCZcgrDE0B2VRJA1XP+u0IKRHEpCUCqQIAKiYBInoJwFsA9gF4j5m7iGgIgP8AMAbASwD+lpnfqFZmZ5Ii0GRQjqLLhrYGfiuSEkHa2b8TkgUAmGkHPsHMv048ngfgUWa+gYjmRY+vNnCc3NicU2g6rgZ+K75FUHTwx22BdAEAdtqBGQCWRveXAjjXwjEUS7RrC3wJIMZHaxBP+pUhBAEA1ZMAA/gRETGAO5h5MYDDmXknADDzTiIalvaLRDQHwBwAGDWwYhWKUebPn4J/XPRH32WkYjMRlB3sWcw7ch1ueHmi0de0QVUJnMjMO6KB/ggR/V/eX4yEsRgAJh1G/v8WmhIMVURgeqB3IgQRVJIAM++Ifu4mouUAJgPYRUTDoxQwHMBuA3UqSi+KisD14A+J0nMCRDSQiAbF9wGcDmAdgIcBzI52mw0g/c+eKEpF8swRVOnpm0KVJHA4gOVEFL/OPcy8koieBnA/EV0C4FcAPle9TEVJJy0RSBv00luC0hJg5l8C6PO3mJj5NwBOq1KUohQhFoG0wR8KtbhiUFGkC2Deket8l5CJSkAJnq+dd73vEoJGJaD0YfuKF3D1sdt9l6E4QiXgCanfJbB9xQs991UEzUAl4ImtPxvnu4RcqAjqT6Mk8O6cY3yXIJbtK17olQKSqAjqjdgvFak7EtqBrEGfxtXHbsc3N46wWI3ii0YlAUBOGvDZDrQ76yvNo3ESkMItd03G5Zee7/y4VQb/1cdu19aghjRSAj7TwC13TcYtdx34XP7ll57vTAamzv4qgnqhcwKOSA78NGIR/Nud37NyfNPxX+cI6kNjJfDunGNw0OLcX39QmE6DPou0VFBVDLb6fxVBPWisBGxRdvC3o0pK0AlApRMqAUPYGPytFJWBCkAWUj9S3GgJVG0JXAz8NDrJwOXg15agGNN2D86974ph1r+pH0CDJfCL1Vd135mwpGfb5PUXt/0dX4M+i+T8wVWf/mePlSh5eO7pkwvtX0QYWeQRSeMk0DP4U3gqIYQkT375dlvl1AJNA3JJiuTKjH0acZ3AL1Zf1XOrKzcu+5rX4+u1A+0pmgJcUnsJmBj4H1v0BQOV1B8VQZjUsh2o8xlfOq5bg+PuOBeAfrNQFWqXBGwJQNNAflwkguPuODcSAPDgqU9ZP14VJLcCQM2SgCaA+hMP/FYePPUpnPuYrNWbUKiNBFwI4GOLvqArBTkx1RZkDfo00hKBiqEzwUtAz/5yqSqCIgLIIhaDLxlIbwWAgOcEfC356dyAG0wIQMlHUElAz/rhUSYN2BCA6zmD2d9aDgBYtBr48omLnB23DMEkARVAuORZLYhn+20mAFerCLEAYhat/rKT45ZFdBLQgV8fshKB69hvKxG0DvyQIGb2XQMmHUb8k7MOhJIQBr/EVYIQPkS04oq5vkvoQxkplBn0vtuCY8Z86xlm7mrdLi4JhCAApV4k24Tln2zepc/BzAkoim1sC0Dq3IAoCYSUAnSpUKkLYiQQkgAUpU6IkMDv9/y57xKUhtPEuYCYjhIgoiVEtJuI1iW2DSGiR4hoS/RzcOK5a4hoKxFtIqIz8hTRXx2gKN7IkwTuBnBmy7Z5AB5l5nEAHo0eg4jGA5gJYEL0O7cSUf88hRx6bfctBM4ecgHOHnIBFtz9OBbc/bjvcpSKuEwBEicHO0qAmVcBeL1l8wwAS6P7SwGcm9h+HzO/zcwvAtgKoPYf41IZKEWQJoKycwKHM/NOAIh+Dou2jwDwSmK/bdG23EhPA2cPuSDzOd8y8P09g0qYmJ4YpJRtqZckEtEcIlpDRGve2PP7Xs9JF0EnfMtAKhKvFlTKS2AXEQ0HgOjn7mj7NgCjEvuNBLAj7QWYeTEzdzFz1+BD/qzP86GLADggA5dC0DRQDF+rApJagrISeBjA7Oj+bAAPJbbPJKIBRDQWwDgApT+6JU0E7VqBTrgSQQifH1C6kSKCPEuE9wJ4AsDRRLSNiC4BcAOAqUS0BcDU6DGYeT2A+wFsALASwFxm3lelQGkiqIILEUhNAtoK+GNVvxlY1W9G5vMiPkU4YcxwvufrF3bc73cL7deSRZUU0In5F/6NsdeSmARuWrim5/74177jsZK+SLhIyPSnC7MG/JzRnwjjU4RNZMHdjxsTwY3LviZKBEkBAMCGoRcBkCeDEGl3di+CiMuGpWMzBcQ0bTUhloFSfG6gU7wvSlAS8DE/4EIAMSZEIDkFKOWJB77JwR8TlASaQFURSJkYzCMATQP5sDHwkwQnAZdpwGUKSFJFBMtOOAMn7FhtsJriFEkAKoJuWlsCm2f+VnRiUChFJguXndD3w5qxCH52xIlG6+pEmRZAJwsP4GLQtxKkBOI0YHPJ0FcKSBInglYZpA36LMqmgiLyMNX7t6YC21KQsDyYxIcAgEAlEHPotX6vHfBBEQFUIW+SsDn5l5SCaSFIE4BPgpYA0AwRxK2BKwEkSSaJpBBcz/y7TglNIriJwTTqdGlxFj4E0EosBAlLfxuGXtRzU6pRCwmYRsJ8gFR8rzykUVQG2gr0pjYSMJUGVADhosmgHLWRAFBdBFIFMOH2f/JdQi++f9HbvktoSzsZaAroS/ATgyaQOviVathcXagTtUoCQPE0oAIoh/Q0oOSndhJQlDQ2DL1IW4EMaimBvGkghBQgbT4gSWhp4KsnT+y5SaPr5Be8HbuWEgA6iyAEASjmOOeqB3s9ligCX9RWAkC2CEIRgOQUUAekiMBnCgAasDqQvKxYB795Tr3qf/DYjR/3XUYmrQmgla+ePBH/+r/r2u5Td2qdBGIuXncKLl53Coatetl3KbXitqFPAOgWgUQ6CSDGZyLwnQKAhkggybBVL/fcJBJKCogFECNNBHkFECOlNfBB7SXwmVmnZD4nVQTSaRVAjBQRFBVAjGsRSEgBQM0l0E4AEgkhBWQJQAplBRDjSgRrrz8dS6Z90cmxOlFbCeQVgKaB/OQRgM80UFUAMTavJ1h7/elYe/3pPY8liKCWEiiaAFQEShamRNA6+JMsmfZFrzII6s+QdcJE/N895cjKr1EGya1AmRbAxbKhqTN/ES6YdoST41y84jbjr5n1Z8hqkwRM9f+aCnojdQ6gzgIA3LYJtZGASVQE3VQRgM25AR8C8IErEQQvgc/MOiW4VYAQkJoAfOEyBSRxIYIgJRAPfJuD32UaGEkbnB0rD6YEYCMNNCUFJLEtgqA+O+D6jB+LoNNkobRBXAWpCaCJg98VQSYB12RdZjySNhgTwO+++DkjryMNE2nAtwB8tQJJbKaBIJKAlJ5/2KqXcfApe629fiyCQ2/7vrVj+CAWQdFlQ9+DH5AhgJhYBKaXD8UnASkCcImvVGC7FSiSCiQIQCrxxUWm0kFHCRDREiLaTUTrEtuuI6LtRPRsdJueeO4aItpKRJuIqNSfzXEx8Scd1yKQOhegtMeECPIkgbsBnJmy/WZm/kh0+yEAENF4ADMBTIh+51Yi6l+kIMkD32YrkIYrEbgUQJ40oCmgGFVF0FECzLwKwOs5X28GgPuY+W1mfhHAVgCT8/xi08/8vpCWAFQA5agigipzAl8iorVRuzA42jYCwCuJfbZF2/pARHOIaA0Rrdm8f3+FMtzgOgU0ERVANcqKoKwEbgNwFICPANgJ4KZoO6Xsm/oJJWZezMxdzNzVf+ihJcuoP3VcOkxrCVQAZigjglISYOZdzLyPmfcDuBMHIv82AKMSu44EsKPMMSThOwXYEoHPViAWwTlXPagCMEzRlYNSEiCi4YmHnwIQrxw8DGAmEQ0gorEAxgF4qswxlN6YFoGEuQAd/HbJK4I8S4T3AngCwNFEtI2ILgFwIxE9T0RrAXwCwJUAwMzrAdwPYAOAlQDmMvO+cv8EGfhOAUlMiUCCABQ55FkdmMXMw5n5IGYeycx3MfPnmfkvmfk4Zv4kM+9M7L+QmY9i5qOZeYXd8ptHHecIFHvkSQPirxhU+lJFBJoCmkcnEagEAqWMCFQAzaWdCFQCDUEFoGShEgiYPGngtqFPqABKIukThDZRCdQYHfzlqZsA+LKrMp9TCQROVhoIQQB78X7fJTSCdgIAVAJtmT3neMw6+iTfZXRElw2VLDoJAFAJZDJ7zvE990MTQQgpIEZiGqhbK9AJlUBOQhFBSAJQ7JInBQCBfMegS5IJoJVYBPdu+qmrcgoxd/nXsRqDAAAn7nvLczX5uP3SW/ts+8qdf+ehkm7qkALyDv4YlQDaD/xQmLv8670er+4/SLQI0gZ/zE2X3gPArwxCpagAAG0HSglAWmvQKoCY1f0H9dykcPult7YVQJJYBq4IPQWUEQAgJAkMW/u27xIKM+vok0S0BVkCaCUWgY90kHfQp3HTpfdoIshBWQEAgpLAZaNfwmWjX3J6zKptgO9EkFcASVyngioCiHGRCK6cJCctFaWKAAAhSSBJLIJbfjXG6nFMzQNISQSSMDHwW7E5TxAL4AO73sJvDg9XBmURJ4GYZCowLYQ6TgQWwdakoY3B34rt9uADuw68L5KFUPXsn4SYU78H1CmjDzqSvzL4mo77VZWBzcHvKg1UGfxpVJWBi4Gfh6JiKBP/fUrBxKB/8gsXP8PMXa3bxSaBNMq2Ci7O/C7aAtMCAMqnAimDPyZPu1C1749TguSEUIagkkAreWTgOvrbEoENAaTRTgjSBn47WmVgY+LPlQxMRf9aJIFWOs0b1KH3B9wJAMhOBiEJAOhOBtue+Qerx3CRDEz2/lkELYEkl41+qZcIfAmgDqsF8TLi8xd/03MlYZCcTEyjrCRcCACokQSAA8ngzTM/5bUOkyJwmQIUO5RZcXAlAEDQxUJKX1QA5bHdCpSlU2oA3AoAqFkSkESVNKCDvxpSBRCTJoLfHD7I+eCPUQlYpKgIdPA3F18CALQdsE7ezxeoABRfqAQc4PuDRorSDpWAI9qJQFOA4hOVgGdUAIpvVAIOaU0DKgAlZu768d6OrasDjpl19EkYcsNpvsuoNT+YtgpnrZjiu4xcnDfrgZ77c9ePx7cnbHBegyYBpVbcMexYAN0ikE5SADE+EoFKQKkNsQBiJIsgTQAxrkWgEvDA9LM3+y6hMUgUQTsB+KCjBIhoFBH9mIg2EtF6Iro82j6EiB4hoi3Rz8GJ37mGiLYS0SYiOsPmPyA0/n7dKAAqAtO0poAkkkSQVwAu00CeJPAegK8w87EATgAwl4jGA5gH4FFmHgfg0egxoudmApgA4EwAtxJRfxvFh0YsgBgVgRnaCUDpTEcJMPNOZv55dP8tABsBjAAwA8DSaLelAM6N7s8AcB8zv83MLwLYCmCy4bprg4qgGnkFICENSGsDYgotERLRGADHA3gSwOHMvBPoFgURDYt2GwHgZ4lf2xZtaxytZ/4spp+9GT/8rw9ZrqZelDn7xyJwvXxYdvC7WjLMPTFIRIcAeADAFcz8ZrtdU7b1+SJDIppDRGuIaM2e/XvylpGL961cbvT1ypBXADGaCPJTNf67TAVVz/4u5gZySYCIDkK3AL7HzMuizbuIaHj0/HAAu6Pt2wAkR8BIADtaX5OZFzNzFzN3HdLvkLL1i6SoAGKmn71ZZeAICe2BFPKsDhCAuwBsZOZFiaceBjA7uj8bwEOJ7TOJaAARjQUwDsBT5krOx3EfHYDjPjrA9WFLCyBJLAMVguKCPEngRACfB3AqET0b3aYDuAHAVCLaAmBq9BjMvB7A/QA2AFgJYC4z77NSfQZj5s/sue9LBqZQEfTG5EqA7TRgaiLQdkvQcWKQmX+K9D4fAFIvgmfmhQAWVqirFMnB30osgrVP2/sLyCZSgJJNSEuBUlcC0qjFFYNj5s9sK4AkcTIwnQ5sCsBHGpD2deO2BKBzAwF+ijDvYM/DcR8dUDoZuD7rxyJo2lJiSGf/mJBSABBYEjApgJgyqcBn7G/ShKErAZhMA7YEYHNeICgJ2CS0CURbIph03h5MOm8PZu+da+X18+I6ATS5LQhGAjZSQBqdZCBp8s+kCOLBn8S3CFxTVQShtQExQUjAlQCSpMlAkgBiqoogbfAnURHkw4UAbLUEoicGfQz+VmIRHPedYR329EfRzx60G/RpzN47F0sHfrtoWcFS9OvJQk0AMaIkIGHQZ7H2ot2iRVB0YBfFlQhCWg3wMfhtfKhIRDtw8PAhogUgnVcHHuHkOLZbA0kCaNJEoQgJKOVxJYAYWyKQJIA8+GwBTM8NqASUwpgWgVQBZKUBCXMAJkWgElBKYUoEUgWQhQQBxJgSgUogYFy3Aq1UFUEIAkimAUkCiDEhApWAUokmXEfwg2mrRArAFKKWCJUwiUWQdwkxhAQQEnEaKLt0qElAMUaeVKACsEfZ1kAlECi+5wOyaCcCFYB9yohAJaAYJ00EoQvggt8O9V1CboqKQCUQIFJTQJKkCEIXQN1RCQRGCAKImb13bq0EUNc0oBIIiJAEoISDSiAQVAAyqGMaUAkoSsNRCRRg7UW7O++k1J6Q0kAeVAIBoK2AUpY8LYFKoCCaBpS6oRJQlJrTKQ2oBBSl4einCAWjcwGKKeauH48nM57TJFACF/MCKgDFFSoBgagAFJeoBIRRNwE8/uKpvktQOqASEETdBKCEgUpACCoAxRcqAcU62hLIRiXgmVcHHtGIFKAikItKwCNNGPyKfFQCnlABKFIgZvZdA4joNQB7Afzady05OQzh1AqEVW9ItQJh1XskM/f5HLQICQAAEa1h5i7fdeQhpFqBsOoNqVYgvHrT0HZAURqOSkBRGo4kCSz2XUABQqoVCKvekGoFwqu3D2LmBBRF8YOkJKAoige8S4CIziSiTUS0lYjm+a4nDSJ6iYieJ6JniWhNtG0IET1CRFuin4M91baEiHYT0brEtszaiOia6L3eRERnCKn3OiLaHr2/zxLRdAn1EtEoIvoxEW0kovVEdHm0Xez7Wwpm9nYD0B/ACwA+COBgAM8BGO+zpow6XwJwWMu2GwHMi+7PA/BNT7VNATAJwLpOtQEYH73HAwCMjd77/gLqvQ7AV1P29VovgOEAJkX3BwHYHNUk9v0tc/OdBCYD2MrMv2TmdwDcB2CG55ryMgPA0uj+UgDn+iiCmVcBeL1lc1ZtMwDcx8xvM/OLALai+7+BMzLqzcJrvcy8k5l/Ht1/C8BGACMg+P0tg28JjADwSuLxtmibNBjAj4joGSKaE207nJl3At3/swAY5q26vmTVJvn9/hIRrY3ahThei6mXiMYAOB7Akwjz/c3EtwQoZZvE5YoTmXkSgGkA5hLRFN8FlUTq+30bgKMAfATATgA3RdtF1EtEhwB4AMAVzPxmu11Ttkl4f9viWwLbAIxKPB4JYIenWjJh5h3Rz90AlqM74u0iouEAEP2U9FdJsmoT+X4z8y5m3sfM+wHciQMR2nu9RHQQugXwPWZeFm0O6v3thG8JPA1gHBGNJaKDAcwE8LDnmnpBRAOJaFB8H8DpANahu87Z0W6zATzkp8JUsmp7GMBMIhpARGMBjAPwlIf6ehEPqIhPofv9BTzXS0QE4C4AG5l5UeKpoN7fjviemQQwHd2zri8AuNZ3PSn1fRDdM77PAVgf1wjgAwAeBbAl+jnEU333ojtCv4vuM9El7WoDcG30Xm8CME1Ivd8F8DyAtegeSMMl1AvgJHTH+bUAno1u0yW/v2VuesWgojQc3+2AoiieUQkoSsNRCShKw1EJKErDUQkoSsNRCShKw1EJKErDUQkoSsP5f6hEM6CgVVXdAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "cells = cle.artificial_tissue_2d(\n", " delta_x=48, \n", " delta_y=32, \n", " random_sigma_x=7, \n", " random_sigma_y=7, \n", " width=250, \n", " height=250)\n", "\n", "cle.imshow(cells, labels=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Mesh between neighboring cells\n", "\n", "Before counting neighbors, we should visualize neighbor-relationships. We can do this by drawing a mesh between centroids of touching neighbor cells." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAfgklEQVR4nO2dW8hmV3nHf49pMuKJOvXQyYEmylQaSzvaISkoYisxMRRGLyrxouQikF5ErNCCES8qhUAqaOmN0hFDY9HE4IGEkppqaJFCazLKeEjSmNGkzThDppoUbS9ijE8vvr2bnT37sNbe6/Cs/a4fDN/77e9993reZz3rv5512iOqSqVS2V1ekNuASqWSlyoClcqOU0WgUtlxqghUKjtOFYFKZcepIlCp7DjRREBErhKRh0XkhIjcGKucSqWyDomxT0BEzgG+B1wBnATuB96jqg8GL6xSqawiViZwGXBCVX+gqj8DbgeORCqrUqms4Jci3fcC4PHO7yeBy8fefJ7s0xfy4kimVCoVgJ/y1I9U9ZX967FEQAauPW/cISLXA9cDvJAXcbm8LZIplUoF4Kv6+f8Yuh5rOHASuKjz+4XAqe4bVPWoqh5W1cPnsi+SGZVKZY5YInA/cFBELhGR84BrgLsilVWpVFYQZTigqj8XkfcC9wDnALeo6gMxyqpUKuuINSeAqt4N3B3iXvecOg7AlecfCnG7Iml90GWX/VEJRzQRCM2V5x8abAjt30pn7Lt1ab/nPaeO/78/tvDdK8sI1TkWIQJzAV+COMw1cldbu/fpCqOl7+pCKH/sKiHrvQgRmMOCOPj05GvL6N6nFYJSxKBr5xJhX4N137gSuq43IQJjhBKHFA18jqmK7w4TLA8RfII39HfoCmWK8mIRQ+yjnB3w5WWyXy1sFur3Ui25A2RJw7aUGczZksvWIb9aqvcuIXz0Vf38N1T1cP/6pjMBH/pO7jrbamBMYWWY4CJgXVtz+3eo3nPb1BLLjqIygVhB4lPZqQMjxHe21NPOvR/S2Lmkzl3fH5pQcT+WCez8Q0V8A68//o5JqPtfef6hszKD2Cwpp+vbFPjUecp675KirJ0eDixV2JRDhZD37O8ziMWaHn1qP0huUg8VUg2PihOBUI6ZCzTXcmIFhuWGGvuesTdChRCZoawlhi9TUNScAKxvHK6BurScEEGRqgcILVqh7tW9Z6w5oBhLkC0h5nBifO9NrQ5YXjIbur9rgKSewBtaHlsqfDGHQjHEKjRr6r37/u78QyqKywRg2azzEtFI3VtYWCLzbXClZS25fGxhuXEsE9i8CKyp9FTj0thjYF9cAjZHUIcoM7efxzqCFHZtTgTALUjn3jdXTqrtwGNYCNixXXWlrZnnFoAu/XpPYdem5gTmlpEspF6udL/LUGMLuVdg6WeGtlPn8q21rGkJ1pZAixSBKXIHqS+pDrWECDwrvi1NCMZ6fStiULQIdAPBSoD60O9hYwZ2yJn1HKlsH18hSN3g5nw0lF3losg5gZY2CGJP4KVcq7bWw/XF1cW+oaDOvfEp9mrPXPlDn7MyMVh0JgB2Uiof5ircihBMzVVM2ec6txHqO7r0pqHKCrUpyEL9thQvAqXh0oAsCNuYACyxbyzjCSkMsTLC0EOfMRtzikLxIuCTpubG5zxCzu8zl16HsM9VGHJs8op1HiDG/UJQvAh0KUEIfMjxfXJOsK4dRoxNsLoMv1zuvxarw8DNiECsVYJQvbLvPXINC1xtTJGtjN17ShhcM8OYvX1pbEYEWlIst/mypjGn/B5WzljM4ZIx9EU097Kmawzk8OfmRADSL72MsTYrSZENdG30tdPKJKbrPEsufOeCUrP5x4vlDtS1wZd7fX0OK9lWn7bh5bZvyUpK6njdvAjA84cIS/H5bOjhSIygCDl3YiEbaBnyfamZYCp2QgRg3X7tJZUYe706xL2sB6cvFr/XUltSCmuxIrC0MYfICqaIMQcRI6gti9QafFYUYrImDlKLWLEiAOsm3MBO4Lqy1t5YqyYWel5Lq0GlscnVARdC7yuInYqu2RmZIk3OtcvR9bultC9EOSntLToTCIHrEGFq1jb1WNR3khLS2JZ6Ztv3u6Vcci2JnReBlqVDhNQC4FNOromyFA3ByhLgEKF3q8amikCHblbg+7nUzNmYSwCsp9sxG1ap8xKrREBEHhOR74jIcRE51lzbLyJfEZFHmp8vD2NqOlwCJefWZNfdZ7kCchcbWsweO3Y2ECIT+D1VPdR5YsmNwL2qehC4t/m9SFzmCXLhcypuK8TYKxGSEpZyh4gxHDgC3Nq8vhV4Z4QyotPtSfsHUaz0Rv1AtmRb6Gwg5PAmxh4JCz5fyqpnDIrIo8BTgAJ/o6pHReS/VfWXO+95SlXPGhKIyPXA9QAv5EW/82a52rncXE73Oduegtwn41IQc24jxxHxnMR6xuCbVPWUiLwK+IqI/LvrB1X1KHAU9h40utKOqIw1tpg7D7fCmkaSooGttS8VMX2xSgRU9VTz84yIfAm4DHhCRA6o6mkROQCcCWBnFvoPnrByMKWLpTPzIUkhACGGLCX7uGWxCIjIi4EXqOpPm9dvB/4CuAu4Fri5+XlnCENTMpeCWkkBWzuG7B3LUnINo3z3N1jw7xwl2OjCmkzg1cCXRKS9z2dV9csicj9wh4hcB/wn8IfrzYzPXC9qOTCHerQhW4eEwVpvm3p4ZXmokqrMxSKgqj8Afnvg+o8B//9JJBP9lN+FXPvkh+zoMmfT0PtTiYKLv3LucFyyEzM2Q8PRGOzsASKfgJs6M5BbCFpC7XSMIQqum69ClZeCmGI5VI7ZicFSGOrtfR3af3+KypliqNyx+QEfXHtrn8+07xmy2ULjd83sQto61thz+GHTIhCi0lyDI2XlTZUX+9DJmuyh39gsZVLgZk+IfQWh7hWK4kTANa1ssTbxtUXGevepv1s7AThXj2vq2GLD71KcCMB8bxTK0a4Vn3Ki0LW3ytnLzolCe81igwg1ZClpz0aRItAl9phyi+vbqenvXxh6PfWZVExlA672WJjj8KVIEUiRXi2ZaY89LCh90023gUzZlvv0ps+8RUk9/hhFigCkcfbaFYSctA3NihD0e0if/Qxj91ry2Tlcxdz6ON+HVacIQ/Ey2a+Xi/v+IiuBPUQs25beN3d6OlV+jlWVKcaGLUN/T00IX8U6RZgNq0IQY1Juzf0srF7MpdOp6nHODmsNPxVFPmMwxZjQwj1CkUsIfLcx56RdsnRZ2dgaRYpAi9UKCh3cqZc8U5dlqR59TmRuhWJFwMqqwBgW0vAuOcbertuJrdG1qZsdbFUMihWBlhiVEno8n/Pzse83dn9fH1poXHPbsbcqBkWLQOheJHTFWuvlUs2llLi06jopaG27cwiKXCLsE2KGOeZS2pre0dJy49T9IMzhmhwNbG25KZZhYy4RFp0JhMbC1uMWq0ugY4SwNcc8Sgg/t9lBqcOETYhACY3Fko2hG5ul7+aD787AObpzBiWxCRFYS8oTgK7vs2TP1OdjzKOkbEQxlnNLywo2IwKlONwKoYK/5Cwgpu0lrSRsRgSWkrrXcXkoiqX99GOfiWlnCTtCXShlWXFTIrDUyal7MyvBsOZ7l9pQQzwgxJe+GFhjE0uELblPzLkyZmfOJbIhe4bel8q+0HW55n4xvneOWN2JJULrjb9lqFfIuSToUm5q+2KUZSk+LE0gbkoEWnxm4XNhKSBbcgfjECHPcliwo09XDHKxORGw2LjG6FZ+brutnPnvYuX0ZIrvnjMr2JwIuGKp17Nii4Wn//RZ20vmtt+HXEOETYqAa+BYCo624nOPEa3MU/RZupRZIqmXFTe1OtBnLIhzr8X3H1jRvxazbB+sCAD4+SmWT0s94NSyuWcMLiXVefoucxXYZi4WZsT7mcDa+4XCNbuzMscSivb7xBSgzYvAkPNiPocg5SaUWIxNWI6lp6mzKhdRLZ2+/2PGyqZFINbSy1R670v72ZjZwBK6E1RdG/sMCUPM7cRT9ZnCf6nKgHRitmkR6LO0AmMEudUsoMuQEAy9p09MUZi6V+xGk2I9P8dwZvMiYKXRDzE0TMmdDYzZtOYeKYYQuf0WglzzGbMiICK3AH8AnFHV32yu7Qc+B1wMPAa8W1Wfav72QeA64Fngfap6TxTLPXA5udcl9+EYa7PQa8XJJVsYe98QfVtKyKrmyClis0uEIvIW4H+AT3dE4CPAk6p6s4jcCLxcVT8gIpcCtwGXAecDXwV+XVWfnSoj1hJhS+pGPsdchefqEVzsSn2AqE93oixX5hSyTN+lzzXlLl4iVNWvicjFvctHgLc2r28F/hn4QHP9dlV9GnhURE6wJwj/utjywFgXAEgz9uzjc94i1TbasfJT2xIDS0uZS3cMvlpVTwM0P1/VXL8AeLzzvpPNtexYPs89RA4hcF16y33wamjZspR6BVsCAOG3DcvAtcHxhohcLyLHROTYMzwd2IznGFrDzxUwJfdcLbl92C27FQQLNrnSxoClOFgqAk+IyAGA5ueZ5vpJ4KLO+y4ETg3dQFWPquphVT18LvsWmjHN2OSTlXPcc6TKBnzFyUKjG1rF6NZrzPpdWi9WO4GlInAXcG3z+lrgzs71a0Rkn4hcAhwE7ltn4jp81rcr7lj131B2YEHsQ9kQ47vMioCI3MbexN7rROSkiFwH3AxcISKPAFc0v6OqDwB3AA8CXwZumFsZiIXFCbglDcdqY4M88xZtua7vs5D5hZoDyHZ2QFXfM/KnwTU9Vb0JuGmNUWvwdXg/tY15mm/pvWOmkWvvnXqZbmlZU3sVYtjen7ewzGafJ7D0MynOGuS6R8x7ppq7CMnQPILPZ0O8xwKbEoGpHW8uxBSCEIeMYhDi3iknCmP4YmgOYel3mTtrYZHNiECoAIyxp93ivUITWwhSNKy1k4qp1v9D+3gzIgDhnB96IilkbxuKGI3KwtJhKFzEYGj3YgqhCs2mHy+2lrXKHqNnCNV4Y082QrjvbSW9nhK3lPYt9cem//ORmJtC1t7fQvD2id1ThxY9K/SHC93rJVP08wRSjMGGjqy6nviKZU+I750ibV3bg1vYYz+363TofaWJQtEiAGkd7rs5JpZtuTbp+BJCCHI1qLlG3T3K3H9PbPEKvTejWBHINU50qQDrDTS135bUVer67dfZmvodmli2nB0UOSdgoZG5BEJM1k5WpsLyikF/T8DYmL//mfa9/WtdrJ5hGKI4EfBNtWJu/JlbNoqN1aDqs0QIYq5c+Db8Kdtcdw4OzSNYoajhgIWJopapMwcW7LNGN0We80/ooYBPqj93n7V29YXAQqwUlwlYcFqX3Onukh1tuUjpqxA9fv9+a/7eJeQ25RAUs1loqQqnPuHWYnHziJVNNzBui4+NoXp4l3JSZS8+MeRbZv2/CCPT7XEszgznzgJccLExVcPvlucjSiGGC+29Ug1/ixABSz3YEP3K6ttqZd3Ykg/HbB7LDubek5vQk34pxaAIEbBMv4KmAjt2MFsXyyG6m26mlt1yfa/c/kwhBsWIgMUA962UfpCHHDJM9UQWfQdnD50sDqGsbHKa25ewhiJEwOL66lpVHhsyrL2n1QbfZ2yDjRWW+DFVnPazg7V+K0IEYFmAxGoQMe47dj9fcegLQQpRmAr8sbK7DSalrS6stSOVz0OVUYwIgI0gmbMhdG88tu00hx/GGvsaW/qbrHLX8dryU2QDOzkcsELOIcnYfMJQwK7tWWM09qEyxmzPJQTWhpxThPRPFQFHLG5ZhvFZ9CVr7kP3z0FOAcj93eeIIZBFiUCuXiJ3ijrF1IpD9/rcZ1OS4/yAC6GHcKGJ5ZOiRCAHS1LEXGI1lCFYEzDfsw4pbLfmoyFiDlWKO0AE6cZubXBYD5AhLNrsk3L7DGssEsPuumOwobvBJHagW2xIU4xtX7Yw3l267t5+tvt7SErJAmLaWGQmkHPdO8d9fMqaWy3IQYhlt/Y+ISlFAGJTpAi0TDlo6ZbPUIGRMrhcespcQmDVn7H9EHK/QOxYKlYEYgWF9Z6hj+84u7Tv1yX0Rhzrvkgl2MWKQEh2QQCGPheblNurfShpGJDCzqJFIETPENvZKdLOJZ8pbWtrqHunHA6t9XMqoSpaBEIRy9kxK7GEPe6xv//Sz5SSBaSiuCXCPks35pQSEH1C2h3j9F4Kvy5ZOsw1BPCNzxx27mQmkFoAQi85xhhnh7AxtV9zL32GJpdQzYqAiNwiImdE5Ludax8WkR+KyPHm39Wdv31QRE6IyMMicmUsw7v4OC5XoK4l9u7FtUKQK7NyFYLSMr6UuGQCfwtcNXD9r1T1UPPvbgARuRS4Bnh985mPi8g5oYydwufUXGkBkaqHWFtGLr+6HEbKiavA5vLfrAio6teAJx3vdwS4XVWfVtVHgRPAZSvsW8TUSbrSBADSbzzybTQWltzGyrdgmws5hWrNnMB7ReTbzXDh5c21C4DHO+852Vw7CxG5XkSOicixZ3h6hRnPMfbkHQsCsHQ2O3VwLF1ytEDfV7kzAFdy27lUBD4BvBY4BJwGPtpcl4H3Dv4XR6p6VFUPq+rhc9m30IznmAvEnIG6pOycwtU9pDVH7gDuMpR2WxEomPZVTjsXLRGq6hPtaxH5JPD3za8ngYs6b70QOLXYOk/6QWAhA/DFis1zy3AW7BxrVEMxAPk7grFMJXddLxIBETmgqqebX98FtCsHdwGfFZGPAecDB4H7Vlu5ACsO9sGizVNHty1sBhoTp/7f5jKb1BOvlup6VgRE5DbgrcArROQk8OfAW0XkEHup/mPAHwOo6gMicgfwIPBz4AZVfTaK5SN0FdeCg1tcZ7At2dwylGGtGeK4lrmUoU1Qa1cQYtSLlbqeFQFVfc/A5U9NvP8m4KY1RoViqGKtOH4I67a5pN8u94nFkDi5ClZKkbA0jwIb2Dbcp630scq3KAyWl7FcGr5V21Od6fcZaljM+DYnAl1cx7E5hSFnr+C7gcViALdMCWlskfXNIqz5b1Mi0O+dlk4u9e/l8n4X+sEYs1HF2KE2NUlolRSnJV1saLHoP1EdXMZPystkv14ub1t9nzVjwrXlDjFmy5xN1raY9sXKUkZgyZY5uvWeQwi+qp//hqoe7l/fjAhYU9cSHibhwtQ+AQt2utqR214LQjomApsaDoyRIwCmhhgWGo8PVvfl+5af295u2RaGKS07+TyBXLRHga1U/hxTjSa3kPn60KJYWYmFzYjAVCVbcHZph1ty95pTrMmoUvu9hPMXmxEBy/QbVH9caA2fVZVc36Gk045zHVRuNjEx6DM5BOn3ALisBFgIBrCdAUAY+1KuGLmWk8KmsYnBncoEci6rjWG5wVkj5O6/FBmMT93mjIPiRWDpwzpi4tsLWJizgDQHg9ayZdHMFQPFiwDYUtySJq2slO1C6HQ5tvCWJFabEAFfLAaAhcnCpYEb25+l3XvpPXNlhEWLwNreIbTD19qTSwisZwFQznn+EL5MXR9Fi8AaQgdA6COrqQIh1OpErF4s9wnA1PfMMYwodonQ0lJWLFtiLx/GsDvUPVMunYZadoRw9sb4/nWJcIS1vVdMMYrZACyJ6Bg5lnTXEHriMhVFikCMsfzSz8WurBhpdqol0jWkbAShsoDQpJooLFIEIFyQhHpASGxCB4Pl7CXXRKXV49+x/VGsCITGx9GpgzT0ZGGqcXbKz63FameQoq6KE4EYQVLCYZRQQmBtu2yMz68p13JnEIviRADiLeu4HvvMfRqtlF7Wt7wSJiv7pLA39txAUUuEqU5/lfBUIB9f5LLd6unOOazZA2FsqkuEHow9AMRSUPjaksN2S2c6fPARrlTE3ERWRaBHPwAsCkDLPaeOzwZF7nGri30WfQvjtpc2eTlHUSKQKlj6qms1SOfssm5/boGawvJzIEL7rSgRyIHVBtTSThpZHsJMBa0F+3zInbnEKLuYR47Hdv6S1M9KAI/NHluwb2pm24J9U4wNDXPT+jSU/4oRgTW4VJ7vXIDLeNyFkDsfu4FhrYH1gzZ3j+pK305LNofyYRFLhEsa8Rq6zrUcrCF7plRPXEpVXkgsDa36+MTnJv4HolzLXFaFYCpd9bE3VFYzVW53WLDElxZS8SU2WIybPkWIQMrHLk3NDViu0G4D8/VVyDPwrrsul5DD/229LxWw2IeSQvikCBFITd+xKUVoCf0Azb2teYiShwOwPAbWHEzyKW+NP2dFQEQuAj4N/CrwC+Coqv61iOwHPgdcDDwGvFtVn2o+80HgOuBZ4H2qes9iCxtyn3yzPCwA242qL1KWBbXL0KRgqhjwHc6twWWfwM+BP1XV3wB+F7hBRC4FbgTuVdWDwL3N7zR/uwZ4PXAV8HEROWeVlQlx3YBjBcvC1GWoMZWKNdvX1v+sCKjqaVX9ZvP6p8BDwAXAEeDW5m23Au9sXh8BblfVp1X1UeAEcNkqKxtiOj/1CkQIpgTASqCWIlJ95vY2WPLvWrzmBETkYuANwNeBV6vqadgTChF5VfO2C4B/63zsZHPNPK4TMRYC24INc0zZuHa1IBYuNlmxN1QMOG8bFpGXAF8A3q+qP5l668C1szYjiMj1InJMRI49w9NONsRKI5c4M2dPUIoAzGH1O/ge0c5ByLKdREBEzmVPAD6jql9sLj8hIgeavx8AzjTXTwIXdT5+IXCqf09VPaqqh1X18LnsW2p/Fko4PJJz3O3Tw1uaH1hiR07bQ8XhrAiIiACfAh5S1Y91/nQXcG3z+lrgzs71a0Rkn4hcAhwE7gtibUNIxy/tVXMEr8X0uU+O9fQQLLE71/xA6EzQJRN4E/BHwO+LyPHm39XAzcAVIvIIcEXzO6r6AHAH8CDwZeAGVX02lMHWGkDqTUzWvv8QSx94UkLm0id1fWR5qIiq/ouqiqr+lqoeav7drao/VtW3qerB5ueTnc/cpKqvVdXXqeo/BLc6IGsqMVXwlrLddk0PlVvc1sZBSj+H9lWRzxOwNI5M/aCT2J9ZSqgTlTnqNeS26ZjEun8RpwiHCJEehxxbxUjXQ33HtfdwKSPWKc6YWK2zsfuuvefmHjS61iGhVTVGpYe4b4ojwqHLSJkRxLC9NIoVgZbYp7RyUNIkYKnEzDZi3DdmLBQtAmtOaMUgZA9WggDEbkix6ilFlhF6GTsmRYvAGiw2svb4aAzbQgdSioYUQwhSZlkhbE/h502IgI+jYjt1beCWkEqmHq5Yn7+JXUZse4t/qIhPo0sVvEsPGVnMTrrkmKvo78VYOwRMbfvaPR4p7N1EJrBkq2dslohTTEKl1hafWpTyHkvL9PV9ykNimxABF3JvSV36d0vkzlTWCFlO232FIHVMbEYELB5dnav8Vu1zNy4XShKrPhZstzw03IwITGHxuKfVzCTUZ2Lhmw1Ys32OHPZuSgSmHJg7Hezalmsz0JoDSNayFZfGYtF2ixnrZkTAUkVPYTEwx7Bqq8WJYB/mhoepKfYA0RCWgtai4vtQ0iPM+nZaioMxcvh37ADRpkQAznZubGcvHYL0P2cpYEsQgJZ+gy9BAGDY7tg2b+L/InQlhkNDzzeM9V5r7hkCSxNpLrTzLd06ty4A8Hy7c7O5TACeE4ElvYKFycXYgjAlkiVlAV1KyQD6pLR7pzKBLkuW6HIHUq4soXQB6L+uuLFJEehmAZYbuyv9OY5Se70luDbqXfBFLDYpAl22Fhx9QRj7my85soDSV1C2wibnBCrzAhFr3mSqzEpednZOYFeZmleYWk6rvfPuUUVgR5hq6EMCUdkdqgjsILWhV7ps5uxApVJZRhWBSmXHqSJQqew4VQQqlR2nikClsuNUEahUdpwqApXKjlNFoFLZcaoIVCo7ThWBSmXHqSJQqew4VQQqlR2nikClsuOYeKiIiPwX8L/Aj3Lb4sgrKMdWKMvekmyFsuz9NVV9Zf+iCREAEJFjQ089sUhJtkJZ9pZkK5Rn7xB1OFCp7DhVBCqVHceSCBzNbYAHJdkKZdlbkq1Qnr1nYWZOoFKp5MFSJlCpVDKQXQRE5CoReVhETojIjbntGUJEHhOR74jIcRE51lzbLyJfEZFHmp8vz2TbLSJyRkS+27k2apuIfLDx9cMicqURez8sIj9s/HtcRK62YK+IXCQi/yQiD4nIAyLyJ811s/5dhKpm+wecA3wfeA1wHvAt4NKcNo3Y+Rjwit61jwA3Nq9vBP4yk21vAd4IfHfONuDSxsf7gEsa359jwN4PA3828N6s9gIHgDc2r18KfK+xyax/l/zLnQlcBpxQ1R+o6s+A24EjmW1y5Qhwa/P6VuCdOYxQ1a8BT/Yuj9l2BLhdVZ9W1UeBE+zVQTJG7B0jq72qelpVv9m8/inwEHABhv27hNwicAHweOf3k801ayjwjyLyDRG5vrn2alU9DXvBArwqm3VnM2abZX+/V0S+3QwX2vTajL0icjHwBuDrlOnfUXKLgAxcs7hc8SZVfSPwDuAGEXlLboMWYtXfnwBeCxwCTgMfba6bsFdEXgJ8AXi/qv5k6q0D1yz4d5LcInASuKjz+4XAqUy2jKKqp5qfZ4AvsZfiPSEiBwCan2fyWXgWY7aZ9LeqPqGqz6rqL4BP8lwKnd1eETmXPQH4jKp+sblclH/nyC0C9wMHReQSETkPuAa4K7NNz0NEXiwiL21fA28Hvsuendc2b7sWuDOPhYOM2XYXcI2I7BORS4CDwH0Z7HsebYNqeBd7/oXM9oqIAJ8CHlLVj3X+VJR/Z8k9Mwlczd6s6/eBD+W2Z8C+17A34/st4IHWRuBXgHuBR5qf+zPZdxt7KfQz7PVE103ZBnyo8fXDwDuM2Pt3wHeAb7PXkA5YsBd4M3vp/LeB482/qy37d8m/umOwUtlxcg8HKpVKZqoIVCo7ThWBSmXHqSJQqew4VQQqlR2nikClsuNUEahUdpwqApXKjvN/o8P8nZ+fSq0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mesh = cle.draw_mesh_between_touching_labels(cells)\n", "\n", "cle.imshow(mesh)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can also combine both visualizations in one image. Note, these images should not be used any further for quantitative analysis. It just serves visualization purposes.\n", "\n", "## Centroid connections and cell borders\n", "A common way for visualizing tissues in this context is by drawing cell-borders and the centroid mesh in different colours." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAulUlEQVR4nO19b6htx3Xfb0VJDE5CchXH5r0jUSlGKXm+UMV9KAEX4zYodkTgJR8a1EIRxKB+sGkMLVjGH+45BYMbiPMtoQoxVUtiRZAYi5I6dUWCKbSxX4LiHD1F0Yul1u++h15S3xK3ASdWVz+cvc+dM2f+rJlZs/fsc/cPLveec/eftWfW/GbN+jObmBkzZsy4uPi2sQWYMWPGuJhJYMaMC46ZBGbMuOCYSWDGjAuOmQRmzLjgmElgxowLjmokQEQfIKJXiOgmET1V6z4zZswoA9XIEyCiewD8GYBHAdwC8GUA/4SZb6jfbMaMGUWoZQk8AuAmM3+Vmf8GwLMArlW614wZMwrw7ZWuuwDwNePzLQA/6juY6K0MfF/5Xa9c3vy+cbv8WhpoTR4XQjK2Kv+Vy+3J1CPWZqPKfucvmfkH7G9rkQA5vttZdxDRkwCe3Hz63vM/S/DccvP7eFl+LQ20Jo+N9XLz2yffDcExY+C5ZVvymIj1+ag6sfofrm9rLQduAbjf+HwfgB36Y+anmfkqM18F3lpJjBnF6JW1J4MZfkgIs8H2rGUJfBnAQ0T0IIBTAI8D+KeV7tU+WpxNU2QyFbelZ2gJue0pPaciqlgCzPwtAB8G8LsAXgbwHDO/VONek8DxcvPTEPtn4xCeoRU0YhVUyxNg5t9h5h9i5ncy8yeKrnW6Ap+utEQbBj6Wz+zwvg3MH3XZYmhEaQ8KDbRpreWAOmhx4lV8WpwMLE0memvAYVZLBnX/nHy62rbH4M9uPkP/+aKjtC0y27TXmVIdmAQJxBS+b4zlovt9en7eKHB05vkgX2G5OD80VVaTLExiHPRZzaVNpp8gRnqTIXYt2G1qfu+AZr9XyRhMFoIucyxEKJr1rME3pOVg3qsf5P0Ad97XNv8EA8nX8TYxRKHp4CuYvSTE7sJysd+2ElQhlhoWUUQ38glg9YebaNwuDocEEgZVKjmkmOo7skgUQyi3pOPFyqHt5Rc+ry1f9nImQ/5qlkftZZFlFZRZABMngR24nChmJ2R2jD1LbeVLafBSpXA828kxYbXmpGs6lWUohU2cuZIUW5HAXCSUZVUB1X0jfLpZRq7W1nhNuq+bBCbhE9giNvgLYCuiqQDJilEik4vM7M+C6/dEtn2uMy6XLQaHg0sy25uyju0LcPX72DL1WOFkP9SsQIrtWgKlM2LCgEnp7Oixyqb2zsDI8CMAwAk2Mq8wkDJ3cqZaMKK2BVTaN6fPnccPYAkEyTFJ39yWQHubipgziYXlKeQPLIy/prK9uZ7dW2dWIIAd9ElH5v0E8eXVmrFa83D5FsdL8JFRPiKMgZttOwRS+jzY7xURvZdCElo7lsDaqqxyDKYsczEwMDXMz/PwJJLX7ZJrFztDPeHKmibuHrFmRBC8M27CdZLvkXg+UHeZJZZR3L4tOwbfdZXx3PXNB4FXX4MIJM6qVIfgyTFheaqzhky6f2iABHMWdNe7wTYtIQLz3NCsN3CSDXC+1NLq9x7JMorat3USuPFTomOzrQEgKcSSQwIAdszgXKUoCp0B4kiJ5mAQXSuHCFwzrcT8jUwmagNWsd97qPb/DlomgYQQYa7i7oRYhM5CLRaWzrxqgzLTgZirfLmEGZQrg7id51s4Od4M1p1QW4kpLyBZQO6AVCEnr0xTcQxGUKKkW++4ZrGGQKFTHEsqSmDKkqDgdlhRgizFTXTuZjkMbSeqgb1Ye64+RPq+73dpu6pZJ4lFSZOzBIA0xds7NjF0mJqqLIWrBkDVWVcQulLNTgzB57T1fK9xz6jTEZC3WUYb+6yDKs7aPV/RgSwHgDQldR6XkOZaiwTMe4QwZjprrFZBbdnicmZWiOhE4+1SKGyyYvd7tWjNDhHQ4ZAAEO5QLSeVOERXGB6SFAalwi5iUgmH1ahatNtQo19Kz8shhEQ5bFTPSlwvvSQwrbRhAcRK2shOP9X8A92zrRbL6H2kqJJGa6ca998FUGX5ZMsUQqLe+Gb9IZOOQs80aRIwFSFLQXsFzJnJlc3t2rnzmmFAdVPWJGRhe6YSgeqAE+hNrI18xWpjYLLLAeCcBIoHT4DZt/nvNhLMVxd8MquHiQBVv4BEvixzdwBHprbjdSc1OnJ/Wxb7+GGKpw6hitABdRa1lHB5em5Se4/vzVmleLwKqZkklWnthHwVIflcx9dcB0tmU617naeJr0AFBVn1B7wckycBVeTOmAlEIBlARcSm5OfwEUCOfD6Lx7xOvzVcbqVjLT+B16wvIFaXjMM4Bt39NnkSSDFTndAqBRU4GqUyZiu061kyLZVeDlX5rGvsoKu9KPU3aBBBrdqKGtcTIaKXkycBE8OsqwIocTQ6oJm6nHJPYBxlXa15Z+mVuozwOVglyy/J9UtRfRnogkAvDoYEiqIEketq5ManXkNsdiuX2EplrB6mC8gSIgapZVg82ysTvioSMyAPhgR6ZIXbNDvSXBZ0SlKyxvc+R0pHC5cEOYN6DOtL4ni0SXSwDD0PpDpQ1J4us1+g2wdHAkBC6KUWk3fX5NMVlug2nsi4z1aRMzvXebz1zKbllKp8LcS4ezk0jsmCgjO4R3J75tY8WJhcFWEqxlbUosrF9fK8nr5HoDpOBMNS0Vg6qeY0VNiVaexQXE4kJXhOb9HZy8CCtjtIS8CGRs57Tnba9vjUFGXrWD4idWXud8TRyiRUNWELZQF2n6tqenFgmVVMsrG2USLNg7cEeuzlaycoX04n7p0TI4IIw2taM32m256VMTTs9igkhDEjGz7kyLJcnJP0HnqdULSaJmsJ5AyKvkNOsAKOV1W24BZVHvr2Aexh/b/GkmZnmTLkxp2hZy3IbuwRiigMSQ4leR4rbHImhoo8tGEJXLmcdVpuctBqzVieDlzFBeyn8pozv4DhS+Xdi5ok7kDjg3jwxxxZBfKMliOi0YYOCzDnXYu5aMMSuHFbzdMphXZeQVIJM5BVNdffJyeMJ5KtAHt5A7n9mbhhh/TZhshrMEPCwYiUDcdzDiJvhzZIANgfHBFmXS5WyQPJdbyUDEKdMvRatMhJaSMjrdgHOmOcgHQ38RRYRqMlODl01LuWt49vKMmojeWAiRynh4KHec9xKMTQBJByn8Fks8zZ7Xq2NJTZX9uBVkKAAPaWck6Hq73sE2Co8HY7loCNSEPtlPhKdqYRWA25JvcYiihNixXJlmsNOEz+FTaKu1rIL5MqV8lsXmwNBPQsmt3Z0OxvoogEiOh1AN8A8CaAbzHzVSK6F8BvAngAwOsAfpaZz8rEjMAOvyk4ulQHmTIkM0SyXNJcBsGaVtXstohgNAegoEJ0C7OuQGkHqprPrLEc+IfM/LCxY8lTAF5g5ocAvNB9ro8K8VNfFVuPMU3RlKo4MSS5DPbxQ8xu3T2C6+1EiNvIl7vhcebtndufk4mpOgavAXhf9/czAH4fwEcr3MeP0HJACHNpYDJxSwkprpJZQEE2e2ciE4nrWa124iM6r8OATAYfRGvthOcOPmejSwATRXsMEtFrAM4AMIB/y8xPE9H/ZubvM445Y+Yjx7lPYrux4Pf+feAj4vuqbbiRiJTa9jHkqeL1HlmJncRWEsmwHJjO/SOBpA1YouHAZoigzh6D72Hm20T0dgBfIKI/lZ7IzE8DeBroNxptF77BFlsuhNBvpzVkUkgyhswkTD231InZYYcIEp832veK4deafoEiEmDm293vu0T0WQCPAHiDiC4x8x0iugTgroKco8DeeMLuiKJO6VNEQ5uYCtBSzbwmoia2dHAFzPrzSEbkGgFMuY17ZJMAEX0XgG9j5m90f/8EgH8N4HkATwD4ZPf7cxqCDonY2no0D7VHDpe8PiulxXBm1vEhItDYf0GAFnRAAyWWwDsAfJaI+uv8BjN/noi+DOA5IvoggP8J4B+Xi1kfsVm0lYHvgsvR5ctsHNpaSE14SVpe2USQMfirl5ZPYEmQTQLM/FUAf8/x/f8C8OMlQg2JnL3mhszrjslhIiaT6/ihSEHSXlnRDXOQmd8py2YfPwRcy9EaaDdjsDJSFM7X+C0QQY8cJfFZC6bTUuP5JLKphDcH9MJnE5VgKzLXfZp1DE4FpmL3dfSpDWofP0TnhOC6r88/kAJanASdlrl+Bp/1pLb7TmFCTjVLJXItUwata6fioElgp9MylUWqHEN2Xuh+xWajnX7tSAl2yRM7xpRtbAL1QSJPqbw5y8/amBwJSM3KHq05vgA0mERiwd78JNG5JnWytjIIgHg/lhCraYmWvL+wFiZHAkB8NlJLVRV2fLKjsMBjLJ2tkmdZFzFlyhkjhf67lkigh9aSxX7m3ISkITBJEjAh7iAlhY7J0qJieyGpp1AIcbnqG+y/Q+cMhZA1IJXHq4/rItGqYpIkMMS6KsfTnlSZBgwSynIenxpPl5YZC+QB4pbK2NWbKX6LQ8jSnCQJAJmNnTirlUYQtOTIQT/Q1CrczBr5DNgzZEo+g+9awH4NRukglJJ5iw6+XBRVEaoJQZd5W1AoQJHZXdkpF5XNN6MK4sc5z6wRIdlB4jVCyzWV5VOCPLHB7Vu2uP6fjMK211lquqsIJ0sCwLix5RCcHeYrZBHUrZcqAJ+uzve903jmxIEXMvuBcQeXLYsNtVm+YRKY5HKgZgoloNPgJ1i5a9VtJRDtslz2rHTGm5JZrfCUcDkjSWMeKgU3hpDzcurmfgzt7TacgGwFiuxkW4T10r3bbGwrrsD/6YzzZe3OU3/ZimA3YClaIQLAbZ3EHJVTxyQtAaDeLJJ9TWswrHCS9yLR0NZoBe9ZIAw366aY+S1ZAz1Muc0ogWbacEuYLAn0qGGuia8nCLWprOf77bulZJCzLbYScgdKC2Z3LB27P0at2lFyzgCY9HJAu+IsKc4v2HlXXante9hyREpqq8hjIJcAxh78gDwaMGi6cyppZGKS0QEbSbNIYK85AP6NJ01k7kVXbXbMCNtpKvIJViplx1lyNbCvv6h/S+RcL8ObokpxTIcTHUiGwGxfnsLt0AucI0XOurcF8zgFmluBt1KRKYW5VDA/q+F4ufvGrRQILIlJLwd6JCXnhLz05osl7OMK12dVFTsx2qHtjNuGHgcwXTWRmhkYg00GW4xZNSp4Kc9BkMAe+rWU2fiBhhCl1irmzkuOG2I2LCWCndCZUjsNHS2o4VQ2U7angIMhgW2DZ74x56JBS/l3rlMz/8JGodOsNtGaVkHvb2oVB0MCAMQzvwkRWw84y42xJs45xyunQltFl3euCEjC/YaaoWlxsvUz8RE1axkcDAksF52HP2PmH9oB14oylDz3EOW8O3AN/r6v7Z2QBNct2SAkFas1+/0FDWD6IUJj9q+e0aVcrOLawWaMoihpu2WFYgsKZpaLzB15HETARxuTPLcSU7tfxsk+dBcQTdcScLD+VEJqrllhzJCg5L7J8pVGU8yqx9RreZaDg7dvgAhbciBOiwRsrz/g7OwUL3wSFJ2MLRKWujLm7kikWfTU6Uc0BySCZDkSwrVj109MhwSE8f4WB5cPZuePLXdsCy2Vmv/YMVZYV6tNvJEjIbKzGH2fPfcYyypoO2OwYrhvbBPMRCvZgVV2/4ntPRDJ5izNIjw/NyJHKUK6mrh1+9ATQ7uWQAEBSM2rrEZW2nTTRj8LjL1GrOKncIXxfMu6QEJXKvbOsaMIGv3oKyYrgKtisSbaiw4oplj6lFjNvE3cZ6+HaxebWnvtlSiR+kwUqHCMIaWdkto0gQicRTwZezukQs9SnNL2YhWz/Gozq7Qk1f5/rSVBTuWi6+/c6zmRmcuRkusvljPXshsoE7V/nppLxrZIoMJ6zdV42o1pKueQSSi14HNY+szT1qr+snZzEszU2ZV89j0EsNu/pq60QwIV1tnVtiA7Iiy73X5K6+j7c8coow3B9FabMtpwEUOtZ4j15xDtN9Q9gOHItQ0SuHJ587uyiZXbgU4l70grd2Zo1Qow4SIC1zE2apJC6Fq1B03xpNJbA5Edmvt7DYU2SACoRgBqg74CXMuUsa0Bn0wl1+jPN98WVGNJNkbKtSbGyhmJkgARfRrATwG4y8zH3Xf3AvhNAA8AeB3AzzLzWfe/jwH4IIA3AfwLZv5dkSS14reIz7pDDXrf/Vz/H4MIfPctJafteYb1VOJbsGWZglUVw6hp47EQIRG9F8D/AfDvDRL4BQBfZ+ZPEtFTAI6Y+aNEdAXAZwA8AuAygP8C4IeY+c3wPboQYSUiUB/kppwZYa9Yh2fNCAPstTdkaNU3sE1HmZrllKh3mlZHauiz7FkzQ4TM/EUiesD6+hqA93V/PwPg9wF8tPv+WWb+JoDXiOgmNoTw30QyCtZMpVBXGNPrK1ByqXd76Nktpd5iiBnLd49Wiq6y0emL1uasGsjNGHwHM98BgO7327vvFwC+Zhx3q/suDZUiBUCh6eiTK7YVeAaSiaD2Zh5QasNC9AU3PYqzLCtPOi70Ow21QACAftqwax8l53qDiJ4koutEdB346/N/KG9R5YrhV1NiHxl0zzLJmctCdhsqZ4L2spikMAXfAJ9u3lFZvH24InJJ4A0iugQA3e+73fe3ANxvHHcfgNuuCzDz08x8dbNGeevuP5WIwOd8ql6x5cqBz3iWoayBVHJqYdC5ohhmv4r6t6CqMLeWgRYnw+7FKEAuCTwP4Inu7ycAfM74/nEiegsRPQjgIQBfyrqD4r5+Kd97kTOTWUUxy/SFURoGDHW1atG4rIMWLIRgMVPJdSSILE+jJEBEn8HGsfd3iegWEX0QwCcBPEpErwJ4tPsMZn4JwHMAbgD4PIAPxSIDQRQQQXMOuOPlxgRMrGBLGmwDzzDi9lMmJ2mbmIQwJhloxf+zzpfsZdBcFaEPQgeOL+c91oDBjqoVfhOGF8XmeoacGn6K4DVCzlTteyXiBJs+3748RVEO8cDPqEi18y6C2Lnu1PcYTCjAyFGSmutc7zVN34Fwl9wgEi0BzWcNXsu1V0CmdaeJ3kGXE2GQ6Fi1JVPMkkzYch+YEgkIEMp4k8BJBEqmbFCGyOCtuf7WuLaXQF3WW87GoY57acLlQ8glnFitxR6kxN0N/B2/kjnYEwe+iemRgKextGaJnc5TWgaIIFAGsbd7hBz4PSKIKXVGae0QBUIlTsVquf9WO/FR3vs1fJgeCQSg1fi9I0nr9VEas4LoGgnLgRqDqr9ev9YWKWojYTIbEjJwZS+qtalj41UcL7P8FzFMiwRKB0kCzNdH5SLLOok8o/OaroiD1h56iUhqN3vPPw9GLa5xRBhMUkg2/12IZZtWtuymRQIe1Ar99FldJdfPUo7EdeLeuRGlqhYq66sEU2arQIYl0FYWoL1cML9XQcLGq5poZz8BKYwBcnJMdYswjpcgQ6FTK75K7w1gz7lGZ4yTLjt7dSw4v4ehYMvFqstZSJDDB0+YM6nCz3WPvn+3ryKLyJEKof8klnXqOi5bH0faz2B6JNCja7CN+bkc5JapyUUq5GRWVhozLZ+usErJQDQUTLRXnl0ZaV1DkuNQVOrbXXObY69pIgurPXu45DdLme1jqjkIDXk0l0iTJYG9pAlNFg1cT9IBVUxYRc9/tpPRIiOpPDkKG82zz9jHIXY/EyX968pQbDXNGpgoCex0QmKoSVOGUMdW6XTL3M5B8oD0+ReEA88cEGrOs14W++9EMkgZ+PY5sYxUV7mz9B5DY3KOQWdjBoggu+Q1gFB24ZCOrMGdZpkzbk42ZnSw+LIQI/13gk3o13xVeWpozzxWmjno8iO0gkmRQJRNSy2ChBklVJjSItuPjRQiSB4oETKww3qrNWfF9DXW4XaosQVMigSAwAArraDLNClr1hzswCNXTkabihyZ+yOoyeGCQQYnx7SZ9WFsQHLWbeaRYdHEZE61cuxEJK3EtBxMxidQLWFEweFmK7eqrBEnZer9xrZSYjKnPIvXArMy61ZnS5REkCT+nxxi20tRP16KnYmaEYLJWQJB5M5Sylte9Z9bMvlU5ai0b4F0qWA72qKmfYHzuISUclBSu5CLSVgCVa2ACgRgy9pK3HhsK8CET2afdRA7JgkVNhfVdvq5iKBW/x2WJQCkpdwqwGUB2AjlnmvKOKjVobwPpE0IWTN+CiruTZGMgCxDWAaTsASARGugN//WSwCORlNKvEllaF/sOOUaoWv7FGTMApwQ7OhK9eQaVwRBY0enwnMkcGUkamESlkBJ2une5p4jEYANc2Zzho0y5Gw1Du2CL7SqOtv38G1u0v/PI1+qHFlyZ/azZqRlMpZATgPzEW0KgCqkmGor6k4K9BpbOVNnR3utPYQV0L9NxyeP73vXUkpV1tgA8xRplcoxRJtf2NoBlQcvJARpuEirg3zrwTHMe/P+ywWA4+5VWmvG6gxZ7Wln36m1XeoMa4ToSu6fZI0ppTuXYlIkkAXbFHQVxAhR3dQOKEVKLnrpzOp7Tpe1sqlkXIrW2PY9fLIXE0Hq4Op8SCcu/1Gj0JwEDpoElgvs18xn7njbUgFIyEmkkZ4bfEbfADOdsTU3ZY0hd3Y9IiwxYGl6gRWgrYOTIoGkWWK9BI5X4UY2a/UDGMS7Xpi2DMj2wvOdK4KiUzV276w2L5Rvm2VYuGQUkaj5udB6KsUkogPJ6Bra56zyHe9CbjroGF56l3e9VxzXTxZipAqotWdWG5b6eUIlyyWwKxzNoidB9WNNfZqUJdAjyIjmbCAhgYAZ22p8XQJ1IlIw81OWVGbYVG75lcm3hWQPhVjJ8jG53zzsW0aZ14zUimhjOq8hM+BUJkcDJg9ixd17su5daZAV+TNy2kSjLwwE+7tHZtslyZWwTCiaQKz205uMpv4aMgPeBikdvH0ZqpKXOClUVHGWzU4sGdAHEMKe/KZc9l4CteSy71kTxjJhiGXlJEmgB5+uVE15Pl2BjzxmXCKGXEZIZvpkIighAEuJNdpi7xpKpFkLxcsxI+O1ti5NkwTWS9AZ76cEF2BnII20b2EOUtfZ4vU1UH/GSwQtTjZWmpJcWVbAgNi+wKWyLk7LMehqjAqlwHv3i93D00nLxapqB+au9cWO1RIcb5JvSOu1WZ1cKyxHX54Mgb5va7x2zMY0LYFuXdbvv18C70CSbGAp2dhSsm1UgeWRo8xRU1WJWJen0CFByyFX0udDhm5LlwQ7YcuKk8l0LIEKG0H0iDoa7RCO7zgL5y8J8ZxnY8DEEe0aB989sEZZ33ksk1yfz1auwvvXhJc4Ko2B9kOEgk5wKURMSYZSCO0wZY305a2MNaMUOQU9keNT2iKb8JTaRHrv6LFFodEphQh7M7sSCw9dB5BkEgay7mrJvXW4FSIoX85ORJF+zw59NgoRWVTIaIySABF9mojuEtHa+G5JRKdE9GL385jxv48R0U0ieoWI3p8sUQbTpQyKoQkg6z6OAWOm/6pjvcRqzSqvYVcrAxYSv5QIWncEJsHOjygkAokl8O8AfMDx/S8x88Pdz+8AABFdAfA4gHd15/wyEd0jksSe+ROTMlKq5gbLnCtBhVi7E8Zzld5DPIuFlDZDoSXFSFlQWoNLHYRaE0YqoiTAzF8E8HXh9a4BeJaZv8nMrwG4CeCR6FlXLp//rRjyc303xRlhCALo4atEDGG0jUAM+O4/hXAgULikKSSCEp/Ah4noK91y4aj7bgHga8Yxt7rv9kBETxLRdSK6jrO/UEnH9O28U0wAClZATiefHNPgb6bJDTkWo0IbT8VXoCJnARHkksCvAHgngIcB3AHwi933Lo11hh+Y+WlmvsrMV3HnbzPFOEdMEcecDXJ391meQiWFORVmBV8MyQrsUlalpZbL7G7JCgi1lYqcmUSQlSfAzG/0fxPRrwL4j93HWwDuNw69D8DtnHvkwFYClSXAwL6AfZmXdWQQlq1WqUw0y7cDMsTgG1QuHQDGnwh8loqqXDYRCNo2iwSI6BIz3+k+/gzON/F6HsBvENGnAFwG8BCAL+XcoxRT9AEEsxeVtu4CkKQgobr+oraN7OoktTB85GT/L2bZDKEnqkvUGBKyDKPJQkT0GQDvA/A2AG8AOOk+P4yNqf86gH/ekwIRfRzAzwH4FoCPMPN/igqRuJ9ADGoNPKAVIJJZS54MMjFJINfZZg/C5cK/+1PJ3gM5croIwpRPe7AO5rDc6Wt3slD7GYMZUGP8gUlgkAq/Amsi5niVYG9jEMW2tdtQxSFsbOwRQrO5Kjv64iaB6dQOCNErgm9QuTpz7CVD0qxQsjTI8BxL1t1jt58PxTX91rV8SFlqDL5M3dEXt4wHRwImpOvYMYkhS0lzqso8OQExjKrACQgR6SCFUgGMHrGI6MtBkYA9O+U6l/prLReb8+11q8Z6tf+cfb0Ex5r5tiBz89XU+4achK1C0xookaHHaO0XsBoPyieQU00YhHD97VMynyzSCscQtlV/IfkU9wjc3tPxORuV91YsQgV/kIZzNRvrJXBMh+0T8DXq0KEfU55Q+nL1cFUlAuj/Hnt2NaFeSlyJAIpkKkFk6XgwJBDCGOZXyPegJkvMSVh5Lz6Vdi2MVqR65cdOGDL/HoRIBYTW5n4CLcCV1lqIvhRYtfN9abhKhVi+QTO2TyC1DceUN2SlDkIEEV04GBIIdXJ2YyvWbFctbjFlVCKsQWbNzLYtsaii7a68FKhSfyGFcDI4GBJQhd14hURgDyjbwaYCl7yZSImqqD1DYtuOVu2YgdgEVQUJ7XkQ0QHprOWdQVJ2M0rcDNR5v8RjRFDajn0UZ5WPZB3ya8jnvYaiFZDar2rtHnyGKe0xWAnBju8RUwBhoo5UCdQ73l7CKC4RVGDLY1sw9lbvBjSz/4ZYi6f07ZgRoTZIwNxZKBE5nbl3TsZ2ZqGBlToLFCmlr+NdG1JGyECrMGgPLjIKtXfAF1PdfB5wa3EbRcRUIH8bJAAUzVjZjJvrRa+0I3CyEsQ63kVunjau6pxyySRBxX0Wa1sDgycCAdkE1gYJ3DD2HRnAfKWFwhbbASIocVolK2bKgBJu4Z2KvQFlz/w5llaP4yVOjmmT+lxBL/h0pX7dXHIpIqYCC6Y9x2CCgyt7dlh3inWqu+eAltMKEMhVkgtgnavhnOTTFejMoUta+QrmtbVTec/Y7evJdKwCZfUlovN9fpUgpuIYdM0YmkzdXUv7RY/aJavB62m0R3cNregEnfH5pqglM7+FLbFWcnhuCQBQ07uSthQNfmXLpT1LwAeL+bJywbtztXGC1WZDUOVr7w1QZedVkdyeWVNr7Z41IxpyiI7tjg/KLLRMtVPC6/T91DcVMc219XJb5htFhABKlXbnfOWdcnbWiBUy2fYULHZ9pVwEKUT9YvtmfO0UIQuvHricqp57VHNcKuvV3r0mYwkY6J16O9txuxpJMHBKvfkq0YYItC0N7yzTI9SWoWNc1x4aIVM5ogdZliWwfX1bjWfW7fupWwIGnA49ewAKZ85cj6xTaSIbfZSAjwik6MfYkd1uN99AEjyXRuityDpLyD7Mvq91j+Wim5DWkeNTsV4Cxyv1vrfRnmMwBbEMuYTGT1FccRGKBtbLrRNT2/noRCjnIAG5sqrH7oUOytzdnfiI4i+IyXFkVnJguzA5EvAqSSxjLgDVYhSlqkP7GlpFR6LzQ6m8AhRHGkYs9MmaDMzJyPVjQqIXA2cvTo4EAMEAzLymtOxz0AIb45lKiWDo3YBS7zdMEY0uxMsHaYjT0+9Vsxun5BgUKUmhAvjukezwypEjIZSZMmCynHUKA6m4ujMHSnKryRPDoA7XqSQLNQDfBiBJHZC6LEjMZUhVhjFM7MGr6JSQQlzFsJdeHlTZg6LDTAIWbAUoYuBUIkhcf/s2M7WPGRMS+VoiABM+2cdq01rtNCkSGEpZbNYtuq+ECDLDijG5Ro/ZRzA2QYUwyD4QmdBut0mRwBhQ6ezULLwE9E4jlSVMJYSUtgX5UjC25VLj3pNJFqrd+DmmX7KPIDOhKQaf97iFARbybLcgXwi+peHY6NtUq/0mQwIlkHReqi9Ash63Yb7WrH81mFYyiK0YrQ0wW2nHnlGlsOVsSWatNpxEiDBpECuHttSVVTGGrTkzOdOwAfV6Be/9tFChYrTa0spu48xMV7lcB1A7MFaYq9VZK2Supsjrsmp8L2NNlcn8vvbWa+aLV5eL7rtE+UtlsBEtUbZDhH0S0UDZgsBESKB2xpSJkG+gRSLoYQ6w1LZyPlefu75YJskgzbrMgXS3pdXi/PPqDCpWYS6BuZ7V95boDVZbApMsFTV0chIkMDTshq1CQopsbyvomLn3PlRfDrjaU3FmzdUB5zJrjV2ystFtf5dyv5L2jIYIieh+Ivo9InqZiF4iop/vvr+XiL5ARK92v4+Mcz5GRDeJ6BUien+2dKYcAyh2LBJQTAS+NaACpmKlVJEzcw+BGFxOQZXJICbT8RKrNYPOeNtmoZ9SSPIEvgXgXzLzDwP4MQAfIqIrAJ4C8AIzPwTghe4zuv89DuBdAD4A4JeJ6J5iSQeCNAEnG/Z+doVE0PoypUeVwQTInGmuTUQLkC17pbV+af9HSYCZ7zDzH3V/fwPAywAWAK4BeKY77BkAP939fQ3As8z8TWZ+DcBNAI8USdnLUtEvkBNGFMGsGPOVO2cqaIgAWolpVyWpVG96QjvHchuS27fCdulZcjiQ5BMgogcA/AiAPwDwDma+A2yIgoje3h22APDfjdNudd81D4myiqMFUgU1iaBaaGgchGQsihbkhC+FYTiJTINUYwr8GVo6IE4bJqLvBvBbAD7CzH8VOtTx3V4yAhE9SUTXieg68NcyGSRmZMbsmtOYQTlyO164+8xUCCAGlZqM3PMC7Zxaou1Faa5FQFZNS09EAkT0HdgQwK8z8293X79BRJe6/18CcLf7/haA+43T7wNw274mMz/NzFc3yQtvzZU/jEomWFBJSpjfvIaHEKSdP2RY1UbKDD+KnJ7BlSPHmMsurYlAEh0gAL8G4GVm/pTxr+cBPNH9/QSAzxnfP05EbyGiBwE8BOBLKtJ2SGr4yOyaO6sGlbeE+QPbUbVUFOSDZjzdiVobuWbIXbPGPwRtS1BiCbwHwD8D8I+I6MXu5zEAnwTwKBG9CuDR7jOY+SUAzwG4AeDzAD7EzG9qCRydhc2ZWLqLbgF2FEBLQR37002BAHqkyigeTJr9Z1gDJW2blBGYA+v8GoQTdQwy83+Fe50PAD/uOecTAD5RIFcaQls02aaflVteMqjMMl7ne/g0cLzcvmchN912KOIouVd0WVB530ANPdi7RoOyujDJ/QS2CmPP7qGdeXJ3fo3JcsbbNM+a+d45JJOtLJW3bPfBSwS1CKC7ntag2speecmijUlUEW5hDFqNtwrvzOCFG4LWMNf3rplRIZclV+Kg07Y2tterPPtv28Z8KWkhtm/HUn5fwAk0ys4PZaPRbkYvbZA9Vk31GViWhbaJ5hy8rg1JIsiOwQtRY7lBi5PtYKpdTeck2FysN+m+mpWL9vVrYDok4DGxSkwkWpykOxAHKPMMzt72/vVacD13JaUTo2I775CXcnvW8MFE33JUgGmQgKdzchvbSRwSMogoiWbMO/psWorr86sIrl/N6bjeWHq11sDe/u/unQVr6aIpe+0QZNskYIf8FOFVXh8ZVF6fAuf1+OKBlaC4e6FMiVN1oIw1H2okEkWtLCCdCDy6oSH79ho1rL8O7ZGAT0EDSGls8bGhQRFAqeImz6ypabCucGqsAs86b+icBW0iCMqdOtg8BKDZNrXbuY3owLuuMp67vvmQMdNKZ88i5R3ZYx6EQGFPjqn4PfcnWBVHZIIItHEp8SSdL+1rgX+odEs156YkEtmcmEJ0IFNBc1I9ayPFIiia6ewMSc/PCifgI1/OlxzVkqIiyq3Rb6rLLKGDODe1eMhJpA0SuHF7MI/70IjdV40ABsA2LKu9LhU+R8lSS3VAJT5/KhEEj6vgG2iDBBRQvXQ1A7HO79l+CvUA1RxUAzlck+GTJ1Pe3HoKJ5T74GBIIIQWyz1btUxE52iGJ83rCZBqDRS1c7+xR49CwkrxW0WhSAQHRQKxjULHgq24qt71xAGUCq+sWkqYOaAkg0WtnRVDxKoWq5L1dDAkMAWTGphWSXBU1hIiKMi8HNQR7KpILURseZgE21rJQBshQmkBUQRVB1hGiDCGIjkrpy8nK2RKWE1ynAC+/lbXgwp+C3Xvv0jGA3gNWQyu9WLtUEvuEqTPDpQcOzSKZqQQOSkPpp39HCwnbEvt6YMtd/Gei7H29+CgSKBHjYF/cryJsbsqxDR2pLHJRLSPYQWoONJcilgx/dsmgikQgCm3GjKJ4KCWAz16hciZFXydslygOOMuVwbdrLH4+wpUHGnAvs9ggDDgGNmMJagit1fWC7AccCEnROc8Z8Cy2pCVsFwAKyyr7VyjFrFIrP8ogSvyoo3ctzRPAQdJAqYV0GrYMAXmxhcnoE36b80cfg2YXusCApAO6kHaYl3+luNBkLgsOEgSMNH0QElBN6BWONm+zTbJjxBBNQdqQAlbzPK8iDhIn0AVDJynX3rvGEHU8JvY0Ny7rxmMqQep2LMELqhP4KIi5FcIhdNUZ+e17LAZldAvCyKENZOAFAVx2BYQGuit5is0iYnrgQszCUwFigo3D/QLhghhHUztwCA4EOafMcPETAKt44DMzhkjQFDkNZPAjBkXHLNPIBVCj2sxphSKmtE2Ijo0WwK5qJlGfCgEoFDrPqM+ZhJoDYdCADMmg5kEWsKhEsBsDTSNmQRawaESwIzmMZNAS5gJYMYImEmgBRxyLkDFF2nO0MFMAmPiwHLQvZiJoGnMJDAWZh/AjEYwk8AYmAlgRkNoZFMR+gsA/xfAX44tixBvw3RkBaYl75RkBaYl799h5h+wv2yCBACAiK67dj1pEVOSFZiWvFOSFZievC7My4EZMy44ZhKYMeOCoyUSeHpsARIwJVmBack7JVmB6cm7h2Z8AjNmzBgHLVkCM2bMGAGjkwARfYCIXiGim0T01NjyuEBErxPRnxDRi0R0vfvuXiL6AhG92v0+Gkm2TxPRXSJaG995ZSOij3Vt/QoRvb8ReZdEdNq174tE9FgL8hLR/UT0e0T0MhG9REQ/333fbPtmgZlH+wFwD4A/B/CDAL4TwB8DuDKmTB45XwfwNuu7XwDwVPf3UwD+zUiyvRfAuwGsY7IBuNK18VsAPNi1/T0NyLsE8K8cx44qL4BLAN7d/f09AP6sk6nZ9s35GdsSeATATWb+KjP/DYBnAVwbWSYprgF4pvv7GQA/PYYQzPxFAF+3vvbJdg3As8z8TWZ+DcBNbPpgMHjk9WFUeZn5DjP/Uff3NwC8DGCBhts3B2OTwALA14zPt7rvWgMD+M9E9IdE1L8v7R3MfAfYKAuAt48m3T58srXc3h8moq90y4XevG5GXiJ6AMCPAPgDTLN9vRibBMjxXYvhivcw87sB/CSADxHRe8cWKBOttvevAHgngIcB3AHwi933TchLRN8N4LcAfISZ/yp0qOO7Fto3iLFJ4BaA+43P9wG4PZIsXjDz7e73XQCfxcbEe4OILgFA9/vueBLuwSdbk+3NzG8w85vM/P8A/CrOTejR5SWi78CGAH6dmX+7+3pS7RvD2CTwZQAPEdGDRPSdAB4H8PzIMu2AiL6LiL6n/xvAT2Dzqs3nATzRHfYEgM+NI6ETPtmeB/A4Eb2FiB4E8BCAL40g3w76AdXhZ3D+KtNR5SUiAvBrAF5m5k8Z/5pU+0YxtmcSwGPYeF3/HMDHx5bHId8PYuPx/WMAL/UyAvh+AC8AeLX7fe9I8n0GGxP6b7GZiT4Ykg3Ax7u2fgXATzYi738A8CcAvoLNQLrUgrwA/gE25vxXALzY/TzWcvvm/MwZgzNmXHCMvRyYMWPGyJhJYMaMC46ZBGbMuOCYSWDGjAuOmQRmzLjgmElgxowLjpkEZsy44JhJYMaMC47/D9E1WlVGbytsAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "visualization = mesh * 2 + cle.detect_label_edges(cells)\n", "\n", "cle.imshow(visualization, color_map='jet')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Analyze and visualize number of touching neighbors\n", "We can also count the touching neighbors and visualize the result as parametric image in colours." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAD8CAYAAADAKumpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAbYUlEQVR4nO3db4wd13nf8e+PlES7idRaJq2wolzJDq2ELmLJJmQXLiy2bmRSLcIEqBMKSaAELpgXUmAjNSqqaRG2hQCqcNI6qBNkU6tWW1d/kNiwEEQUVcVro0AtiTRpWRTDipYVay1WNOs6cpNIMqmnL+6MNFzePzP3nplzzszzAS52dzh779nh3t8+58yZMzIznHMuZWtiN8A552bxoHLOJc+DyjmXPA8q51zyPKicc8nzoHLOJc+DyjnXOUlXSzpSebwo6WMT929rHpWk7cAngbXAfzSzfa28kHMua5LWAt8G3mtmfzZun1YqquKFPwXsALYAN0na0sZrOeey90HgG5NCCuCCll74OuCEmT0DIOleYCfw1LiddfF6Y8OVQV74PRccCvI8oRy/ZHPsJsx09VNPT/y341vSb39qrn5x8vE89Mx7OmzJKmcPnTazDYs8xY9K9pc19jsJR4GXKpuWzGxpwu67gHumPV9bQXU58Fzl6xXgvdUdJO0GdgPw5rfCvzkY5IUPblCQ5wll2w2/HbsJUy1fuwOunrLD2afZdvjBztrTB8sHdkz510Po5yJdtvY9TaxY6vor4JYa+/0LeMnMts7aT9JFwE8Bt0/br63B9HFpcc7/jpktmdlWM9vKJQuFvGvZ8rXT3niuanpIjdh9af0xbULAhTUeDewAvmpmL0zbqa2KagW4ovL1JuD5ll4radVf3G03pFWZNAmg5Wt3eGU1Q52QKpVhFa26mpMIHho3MaPbB+1VVI8DmyVdVZR2u4AHWnot1xGvrMLLrboKWVFJ+mvATwKfm7VvK0FlZmeAW4GHgGPA/WZ2tI3XykmTv7htmzd0PKzCyyms1gBvrPGow8z+0szebGZ/Xud1W2Fmf2xm7zCzt5vZHW29Tm5SCqt5eVidb9H/11zCquz6zXqE5jPTI+hLWHlghWX3KfnAamEwvRYPqkiWD+yIFlghA8bDKvwfnjKwUgwtr6gC2X/D9bGb0EjMwArFw6o9qYVVrIqqrekJnXnw57edt23/Ddez/cCXum/MAsaFVRvTGdoKlfJ5fQpDeKvDKuaUhjKoupZlUI0Lpz5aHV6LBlcXlY/Pt2pfNbi6Di1R/6xeSFl1/R78+W21Qyq3LmAdi3QTu+ye+UB7d7ruGsYao0q2ogpRNeXYBayjDKvUZrqv5t3B/onV9UuuompSNdXRx8qqVLfCil3dhHr92D9HqrqsqgZ/1i90QA1JDmcOFw2Z8vs9rMbrKqwGfdZv86XHgfbW6elrF3C15QM7ku4OzjvQvjqcJoXV0LuYdp9aH1wvL6HpWhJB1YWyC9j3wFo9fpVaBdJk3Kpp28ftP7TwantVhhZWT6j3um2tmd7EO7ZebL99sLuVD/seVq+5LXYDppsWIl0G7CJhlnKX+7yw+p4O1VnMbpqfkKzOMihXwcKvVZXMGJVzsaRWdYbSxriVX+vXoT6fCcxJSgFRzv1qtJhgwtVUqY1rBn0eVYeGMsCeunED7LEDbNbr5zjuZfcJfWjx5xFwYZ3UOLP4a1UNsqIqeWWVhmowxA6pOoY8837NGnjjutmP0AZbUbm0DPWNnxsJLoiQGoMPqtXTFvbOUR7vfShkiwJJ/Iyfy1Ptrl9ggw+q0jwBVf3eJMPKudAErO3+ZQc9RhXSIkEXnFdTri2RLvbzoAK+ojBn/5IIKw8pN0awS6s8qNzCPKTcGEGv/xSwrsajzlNJf0PSH0j6U0nHJP2dSft6UAUWrarykHJjBL9IPWxF9Ulgv5n9GPAuRvcAHWvwg+mhun1VZVh1MsDuAeVWaXUFjUBXJUu6BPgA8EsAZvYK8Mqk/QcfVG1qJbA8mNJxG3Bn7EZEWOm13lm/9ZIOVr5eMrOlytdvA74D/CdJ7wIOAR81s78Y92QeVB1YePqCh5MbI8raY/UrqtMzVk+4AHg38Ktm9qikTwJ7gH85aWfXgbnCygMqfR1XVdEXRiwH0xe3AqyY2aPF13/AKKjGGnRQtTE+FYQHlCOBUBon0BiVmf1vSc9JutrMjgMfBJ6atP9ggypGSE2tqjycXCHJgCqFXeLzV4HPSroIeAb45Uk7DjKoYlZSY8PKQypvgbp/SQdUKWBQmdkRoNYqoEkE1Uk2sq/onu5hX+TWtO+1s4HXRG2GC6lhWGURSpNEuNYviaCqajuwkh2Xctnae6T4pMYJk6wDCqLd3SG5oCrtm3wC4DxDqMJcel4LqOq2KWGVfUjB6FqWFhbGq/Oy2asbal5NuS6Mu4yqFyEF0S5KXugpJT0LfB84C5wxs62SLgXuA64EngV+1sz+72LNnK0aVl5hubaMq6Jm6U1IlTLt+v09Mztd+XoP8IiZ7ZO0p/i60/NakyqsbXhF5ebTNKDKLmDvQirSwnltZONOYFvx+d3AMn4C3mVqngqq1LuQgmwH0w04IMmA3ysuOrzMzE4CmNlJSW8Z942SdgO7Ada9dewuznVqkVAaZ5t2sGw9C6tMg+r9ZvZ8EUYPS/rTut9YhNoSwMVb3xH/vvJuEEKH0Sy9C6tw1/o1slBQmdnzxcdTkj4PXAe8IGljUU1tBE4FaKdzC+k6oHort4pK0g8Ba8zs+8XnNwD/GngAuBnYV3z8QoiGOjcPD6jAcgsq4DLg85LK5/lvZrZf0uPA/ZI+AnwL+PDizXSuvtTCqXfdv5zO+pnZM4zWOV69/f8wWrLBuU6lFlC9FKmi6sXMdOdSD6lt6skt69cAb6jxCCzZa/2cqyv1kOqdnLp+zrkBynAw3S3A16JyWfIxqmHx7orL1toaj8AGFVS9OkXsXAw5LvPi5uddP5clXzivG6lUVd71c1nyimpYvKJy2fKzft1YtgejTcDzgHJZC3jWb9wKwZP2HWRQxeAB5Xoh/Aqfq1cIHmuwQdV2VeXB5HrJJ3z2gweU67X6C+etl3Sw8vVSsVhm1bgVgsfyoArEA8rN0ovlXupXVKenjTkVzlsh2My+PG7HQQfVot0/DyfXVJNbuTW5CW9nAnb9JqwQ7EFVtXxtEVDXvL5t1twmDya3iPfZ9Y32D3F/ylbCLsBg+pQVgscaXFC9FlBjeBC5vqmGXZC7WoarqMauEDxp50EE1bRwcq4LTaupZJUL5y1o0grBk/Q+qDyknAvMF84Lw8NpQO4EPhS7EQPi61GF4SE1EHcWD2DvQ1FbMlNvun0Q7aLkXgWVh9QAVAKqKvWw6pUIC+f1puvnIdVTY0JpknFhtde7hWH5JTTz8YDqsQYhNUkZXrECq1fdPgh21q+pbIPKA6rnAoSUa4f5Wb/pPJwGooWQ2vtQnKrqK/pSr6oqE5z1s36TeUj13J1MHCgPJdaA+1cUZE54GoqgmvUILemKysNpADru4sWqrPrCBGfW1qlvXg36ukkGlQdUjyUw9rS6suoiuPrSBTSJsxfUiY1Xgr5uckHlIeW6Vg0ur7amM8Qray+qsWfYoMpmjMq5trUdUn0YqzLEGdbOfISWVEXl1ZRz6TsbITaSCSoPKefSZ4izEZZPSCKorn7qabg6divckPnYVD2DDqpDP/Ie4FDsZjjnZjDEy9QZTA9r5mC6pLsknZL0ZGXbpZIelvR08fFNlX+7XdIJSccl1f47pY8b+rg1/wliuJNOJii6bnRZTeU+oD6qqC6Y+Qitzlm/zwDbV23bAzxiZpuBR4qvkbQF2AW8s/ie35EU4cqgjnlguQZyD6uzrJ35CG1mUBX32fruqs07gbuLz+8Gfrqy/V4ze9nMvgmcYHQLnNqSr6qmBZIHluu5cowquaCa4DIzOwlQfHxLsf1y4LnKfivFtvNI2i3poKSDvPidc/8t9bCaxQNrPD8m2TOIMo8q9IRPjdk2NnXMbMnMtprZVi7ZcP4T5R5W4ONYGYh1ti/f7l+6Y1TjvCBpI0Dx8VSxfQW4orLfJuD5eRuXXFgtEjgeVm6VHMPKEK9w0cxHXZLWSjos6Y+m7TdvUD0A3Fx8fjPwhcr2XZLWSboK2Aw8NudrAAmG1SKGHFZD/tkj28eeYHdMbmGM6qPAsVk7zazRJN0DbAPWS1oBfgPYB9wv6SPAt4APA5jZUUn3A08BZ4BbzOxsk1aPbUMRVvaJcT3LjoR6o417ntsCPXfKyp/RA+scbayq0Mpt3AvltX4hSNoE/EPgDuDXpu07M6jM7KYJ//TBCfvfUbywq+tOhhFW4IEVUJuBNE3AMah/D/wz4OJZO/rqCXV08aYa2ht3KMFcQ9OxqpBduaYadP3Wl2f1i8fu6vNI+kfAKTOrdUlKEpfQ1KWPW/fdvy4DZEiVlWskVjCt1uASmtNmtnXKv78f+ClJNzK6r80lkv6rmf3CuJ29okrNkCorD+VaUgkpCHcJjZndbmabzOxKRlez/MmkkILMKirouKqKFRpDqqxuY1jhPMHqQfWUwqlq0KsnuDGGFlbggUW6AVUVOqjMbBlYnrZPlkHVyXSFFN40ZRuGFlillv8PUluDKoeQCjk9oYksg6oUZXDddacaXIFDK7WQykU5RtW1rIMKBhJWQ+oGTtJxteXGKy+h6Vr2QQUDCSt3rharLTeZD6anxH/x89JwMN67fYuJMUbVm3lUwS5e9pDK1214F7llKS9FnI2FwyrVkPI3XzNTAsurqcXEWuHTu36QbkC5xfg4VitijFH1qqKCOaoq/wV2rrZXWcPLrJv5CM0rKjcMt8HeI7Eb0Q9+1i+Q2tMVcqimfHwqmL3XVD4/EqsV4y3bg7GbUEus6Qm96/qVZnYBcwgp15pqaLlmcrpdVhYmhlUuIeXVVKtSCatcqil4/Vq/rm+X1cuuX9U53UAPKLfK3mvS6wambHQJTfjB8ll6XVGdxwPAjRGzssqpmoL87pScr9tIewZzqu3quVS6gamL1fXrfVBNPfvnoeAqug6r3Kqpkl9CE1h2Kyp4cEbXVVjtfQiWD+zo5sUC8ktoAqsdUr5mt1ulzflWex869+vlAzvYdkM+lZUv8xJQ40rKw8pNEOqs4OqAqiorqxwCq8HtsoLqVVAt1NWLfYMB7/Yl65zuYIu/H9WuYKqhFWsp4t6MUQUbj/LAcJN0+Ecs5fErn56QCg8rl4AUwyrUYLqkN0h6TNLXJB2V9K+m7Z99UNknlN/ZPZefSEMCqYVVwHlULwN/38zeBVwDbJf0vkk7ZxlUZTi1GlBdVlVewbkpUgurQLd0NzP7f8WXFxaPiSsJZDWY3nnlVHeA3YPGDcSrrAl2uyxJa4FDwI8CnzKzRyftm2VF1blJl9ykfCmOCyeBqSspVVU1u37rJR2sPHavfh4zO2tm1wCbgOsk/e1Jr5lFRZXMGJSH0vAkEFKlFOZbNZiecNrMttZ6TrPvSVoGtgNPjtsn+aBKJqScS0TM+VahZqZL2gD8oAipNwL/gCl/FmZ2/STdJemUpCcr2/ZK+rakI8Xjxsq/3S7phKTjkua6OVEng+XO9UCMLmGgeVQbgS9KegJ4HHjYzP5o0s51KqrPAP8B+M+rtv87M/tEdYOkLcAu4J3A3wT+u6R3mNnZOi0Hr6Cca6rL6wVfDXQJjZk9AVxbd/+ZFZWZfRn4bs3n2wnca2Yvm9k3gRPAdXW+0Sso5+bXXWWV352Sb5X0RNE1fFOx7XLguco+K8W280jaXZ4ReOvKoQWa4ZyDbsIqtxU+fxd4O6MZpSeB3yy2jyuJxk7iMrMlM9tqZls3JD+k71weugirbNajMrMXys8l/T5QDoKtAFdUdt0EPD9365xzjbU5jaG8hKZrc1VUkjZWvvwZXp/78ACwS9I6SVcBm4HHFmuic24ebVRXFmmMauYzSroH2MZopukK8BvANknXMOrWPQv8CoCZHZV0P/AUcAa4pckZP+dc2ka3y0pw4Twzu2nM5k9P2f8O4I5FGuWcCyP01AVDnH01k66fcy4fQbuABmfOrJ35CM3Ptzk3AMsHdow9Jd+UmTh7pvvY8KByztU2Ciq/C41z6Uho5YRkGB5UziWjZyE1GlBffKzKbA2vvLRu8QY15EHlXM8FnfhpgFdUzrmQgs9ON3lQOZeEnnX7gjLgTPernHhQOddTra1Rdaadp53Gg8q5qh5UU60uovcq8FJ7Tz+JB5VzPdL6Sp8G/KDdlxjHg8q5UubVVCfLERsQYZmBJK71O75lc+wmOJe1Tu9Gc6bGI7BkKqpth0cHevnadG606FwOOg0pI8pgehIVVdW2ww++FlrOdcpvMDtbGVRDrahWq4aVV1muM9WwSnjMKtrdkv2s32TeLXRRrK6wIgZXzNu4nydAxSTpCkb3Cv0RRvG3ZGafnLR/FkFV8sByUZXBlXCl1bpw0xPOAP/UzL4q6WLgkKSHzeypcTtnFVQlDywXVceBlVQ1FWh6gpmdZHSrPczs+5KOMboHaH+CquTjWC6qDgIrqZCCJmf91ks6WPl6ycyWxu0o6UpGt3d/dNKTZR1UVdsOP+hh5eKYdbZwziBLLqSgyWD6aTPbOmsnST8M/CHwMTN7cdJ+vQkq8C6hS9QcZxKTDKlSoOkHki5kFFKfNbPPTds3uXlUzvVajblaSYdUoHlUksTotnvHzOy3Zu3fq4rKuSyMC6s7Ew+oUriZ6e8HfhH4uqQjxbZ/bmZ/PG5nDyrnEpBFSEGw6Qlm9j+g/h28PKicc/VFWj3Bg8o5V5/hl9A45xIXafUEDyrnXH2+wqdzLnk+RuVcGPsPX8/2a78Uuxm1ZLn2mi+c59xi9h++/pyPKcsypCItnOdB5XpjdTilHFZZhhS8fq3frEdgHlQRpPwG6psUj3W2IQVeUQ1FTl2TnEw7nikd66xDqpRiUEm6QtIXJR2TdFTSR4vtl0p6WNLTxcc3Vb7ndkknJB2X9KHwzc5TTl2TnPhx7FA5PWHWI7A6FVW5ZOiPA+8DbpG0BdgDPGJmm4FHiq8p/m0X8E5gO/A7ktaGb3o/+JtsMXWPXwrHuRfVVDk9YdYjsJnTE6YsGboT2FbsdjewzOi68J3AvWb2MvBNSSeA64D/GbrxqWvyJsrldHoq5gme8nu6Pta9CKhSpLvQNBqjWrVk6GVFiJVh9pZit8uB5yrftlJsW/1cuyUdlHTwB9/58zmaPlkKvxhN30gp/MXPxaLHqstjncLvYlAJd/2A+kuGMn7pBjtvg9mSmW01s60XbvjrdZuRhXnfCPsPX++B1RE/zgtIsesHE5cMfUHSRjM7KWkjcKrYvgJcUfn2TcDzoRpcV6xliUO8AarP4V1Cl5RUb+k+ZcnQB4Cbi89vBr5Q2b5L0jpJVwGbgcfCNbmZ3G8R73/5zxXyeLR9bHP+vZso4Vu6j10yFNgH3C/pI8C3gA8DmNlRSfczuj/XGeAWM4twGeO5uqiwPFTaldPx7WVIQbqrJ8xYMvSDE77nDuCOBdrVmrbuBdjmm8jPCrZ3fP3YNmTAy92/7KBXT1jkXoBd/3WPdWo9tpyqqFJvqylId4yq7+YZw4r55hnSmcGufs6Qr9PrkIL0pyf0XW6D7m12hVIIw1gVq5sh0sx0D6pVZgVWSr/Qoc+ADf1axEV/3pz+0M0t4Fk/SXdJOiXpyVn7elBNMC6wUnzjhpilncvKA12Y9+cdREiVwk1P+Ayj64Fn8qCaoQyslN+w81yu06R7l/LP3oamP++gQirgwnlm9mXgu3X2HfRZvyb2sYc97IvdjInaDpOuTuPnFIqDCqhS/bN+6yUdrHy9ZGZL876sV1Suti7CMBUptSUp9ceoTpfX8haPuUMKPKhcQ22ebczJIKsp8OkJLh+hQyXVkJrUrsGGVMmnJ7hchAqXVENqksGHFIyqqlmPGiTdw2hBzaslrRTXDY/lQeXmltMCdvOqttFDKiwzu8nMNprZhWa2ycw+PWlfDyq3kBzCZlH7D1/vIRWZT09wC2t6wfQQwq2/4qzz4hWVC6ZOAHlI5S7OynkeVC4ovxyn7+LMT/CgcsGNC6TcQyrlqxK69SrwVzUeYXlQuVZUgyn3kHJVXlG5nklhXauQvKoqdT9G5Wf9nHMN+Fk/55LnVVWcs35eUTnnGvCKKnn72BO7CS4Bw66q4pz184rKOddAnPtleUXVkFdVbtjiTE/wiso510CcisqDyjnXQJzBdA8q51wDRhuD5bP4GNUcfJzKDZePUTnnkudjVM655PkYlXMueV5ROeeS5xWVcy555SU03fKzfs65BsKtniBpu6Tjkk5Imnoq3Ssq51wDYbp+ktYCnwJ+ElgBHpf0gJk9NW5/r6iccw0Eq6iuA06Y2TNm9gpwL7Bz0s4yq3n/5RZJ+g7wF8Dp2G2paT35tBXyam9ObYW82vu3zGzDIk8gaT+jn3mWNwAvVb5eMrOlyvP8Y2C7mf2T4utfBN5rZreOe7Ikun5mtkHSQTPbGrstdeTUVsirvTm1FfJr76LMbHugp9K4p5+0s3f9nHMxrABXVL7eBDw/aWcPKudcDI8DmyVdJekiYBfwwKSdk+j6FZZm75KMnNoKebU3p7ZCfu1NgpmdkXQr8BCwFrjLzI5O2j+JwXTnnJvGu37OueR5UDnnkhc9qJpMo49F0rOSvi7piKSDxbZLJT0s6eni45site0uSackPVnZNrFtkm4vjvVxSR9KpL17JX27OL5HJN2YQnslXSHpi5KOSToq6aPF9mSPb2+ZWbQHo0G0bwBvAy4CvgZsidmmCe18Fli/atu/BfYUn+8B7ozUtg8A7waenNU2YEtxjNcBVxXHfm0C7d0LfHzMvlHbC2wE3l18fjHwv4o2JXt8+/qIXVE1mkafmJ3A3cXndwM/HaMRZvZl4LurNk9q207gXjN72cy+CZxg9H/QmQntnSRqe83spJl9tfj8+8Ax4HISPr59FTuoLgeeq3y9UmxLjQEHJB2StLvYdpmZnYTRLzTwlmitO9+ktqV8vG+V9ETRNSy7Usm0V9KVwLXAo+R5fLMWO6gaTaOP6P1m9m5gB3CLpA/EbtCcUj3evwu8HbgGOAn8ZrE9ifZK+mHgD4GPmdmL03Ydsy2F45u92EHVaBp9LGb2fPHxFPB5RuX8C5I2AhQfT8Vr4XkmtS3J421mL5jZWTN7Ffh9Xu8uRW+vpAsZhdRnzexzxeasjm8fxA6qRtPoY5D0Q5IuLj8HbgCeZNTOm4vdbga+EKeFY01q2wPALknrJF0FbAYei9C+c5Rv+sLPMDq+ELm9kgR8GjhmZr9V+aesjm8vxB7NB25kdDblG8Cvx27PmPa9jdGZnK8BR8s2Am8GHgGeLj5eGql99zDqLv2A0V/0j0xrG/DrxbE+DuxIpL3/Bfg68ASjN/vGFNoL/F1GXbcngCPF48aUj29fH34JjXMuebG7fs45N5MHlXMueR5UzrnkeVA555LnQeWcS54HlXMueR5Uzrnk/X8PuH2bLFm5CgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "neighbor_count_image = cle.touching_neighbor_count_map(cells)\n", "\n", "cle.imshow(neighbor_count_image, color_map='jet', colorbar=True, min_display_intensity=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note, the numbers along the image border may not be accurate. Hence, we should exclude the corresponding cells from the further analysis." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQEAAAD8CAYAAAB3lxGOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAXBklEQVR4nO2df8xddX3HXx9LYQniKDCwlE4qqW6gEZUxiGKmBqGY+GA28VEczSRpUqlR59zKHMZsIIqx/qGBWEdn0WrLgtjGTTdEE9E4fmgQKNgfQrW/1oqFaUyktHz2x3NOOb3Pufeee8+P7/fc834lN/fe73PuPZ/ePp/X+Xy+33PPY+6OEKK7vCB0AEKIsEgCQnQcSUCIjiMJCNFxJAEhOo4kIETHqU0CZnapmW0xs+1mtrKu/QghymF1nCdgZnOArcDFwC7gfuBd7v5o5TsTQpSirkrgfGC7uz/u7geB9cBUTfsSQpTgmJredwGwM/N8F/Dn/TY2M522KET9POnuf9Q7WJcELGfsqEQ3s2XAspr2L4SYzS/yBuuSwC5gYeb5GcCe7AbuvhpYDaoEhAhJXRK4H1hsZouA3cA08O6a9hU1i1edfuTxtr/dM2BLIcJQiwTc/ZCZrQD+C5gDrHH3zXXsSwhRjlqWCEcOYkLbgWwVkKJqQATkx+5+Xu+gzhhsmDwxCBESSaAmBiW7RCBiQhIIxOJVp0sGIgokgcBIBiI0mhisgaqSWpOIomJyJwYlgYqp86guKYiSaHWg7ah1EHUgCbQQyUBUiSRQIU0nZtdksOFVV4YOYSKRBCaALsggFYBEUD2aGKyQGBJxEicPiyb+O3+6ruZIWo9WB+okBgFkmQQZVHHUlxiOQhKoi9gEkNJmETRZ9ndIFFoiFCKPrs8z1HVRESFaRVYEHaoMAEmgNLG2Am0m9JF52P4nTRJqB0RUhBZAETa86spWxFkUSaAEqgKqZZISq02oHejh5D2nFN72AAcBOGn9sXWFMza/WrCTEzfA0++cEzoUETmqBBJO3nPKSALIcmD6YMXRVMeJGw6HDkFEjiRQETGJ4FcLdg7fSIgESYDRWoBBxCCCPAGoGhCDkAQmCFUAYhwkgYoJVQ1IAGJcOi+BqlqBLAemDzYmg18t2FlIAGoJRD86/wWiOiTQS5VLiFUc8WNdNmzjeQItO3sw9wtEOk+gAQ5MHywlgi6U+m0UwKTQ+XagKcZpD4qW+qOi1qA6JkFena4EmmgFxqGJI/+JGw5H0Ra8dsUvAbgJ+Pt7rg0bTEfpbCUQQgCDqoH0qN9k6R+6IkgFkHLTRTcGiqTbdHJiMHQF0Ds/ELrnb7Ii6E38PNpYEbRkglCXF0sJLYGUw/fsCx0CUK8EiiR9HhJBLUgCEI8AIB4JQPUiGDf5s7RFBJe87bYjj0/8l/DzLAOQBEAS6EcVEqgi8XuJWQTZ5M8SsQh0odGYBBAbZSYJX7vil7UIANo5Wfj0de1agi21RGhmO4DfAoeBQ+5+npmdBGwAzgR2AFe4+1PlwhRNMOqyYV2J38tNF90YTUXQ7+jfZkq1A4kEznP3JzNjNwEH3P2TZrYSmOfu/zDkfRppB2KrBGJqB1KGSaCpxM8jpAjGSf4I24LG2oEpYG3yeC1weQ37EDUxqC0IKQAI0xpc8rbbJvLon6VsJfAE8BTgwBfcfbWZPe3uJ2a2ecrd5+W8dhmwLHn62rGDGAFVAsU5657doUPoS10VQR3JHlk1UP3qgJmd7u57zOxU4C7g/cCmIhLoeR+1A5ERswRgfBGEOKpHJILq2wF335Pc7wfuBM4H9pnZfIDkfn+ZfQiRx6itQRfK+nEZWwJmdryZnZA+Bt4CPAJsApYmmy0FNpYNUog8iohAyT+csdsBM3spM0d/mFlq/Kq732BmJwO3A38M/BJ4h7sfGPJeagciI/Z2IEtvaxBj0kfSElR7URF3fxx4Vc74r4E3j/u+QoxKeh5BjMnfBjp1xqCYXGIXQMxnEUoCovXcPW9B6BBajSQgZrHiFT9kyfIdocMQDSEJBCLWScEVr/jhkccSQTeQBAIx56LTQodQCIlg8umUBH59+pPDN+ooK17xw6OqgCwSwWTT6asNhySGdqBf0uexZPkOvnXLmfUFI4LROQn8+vQnozhpaM5FpwUTwSjJLyafTrUDMfGiez/NvGObX9suI4Aly3eoNZhAOimBkHMDL7r307zo3k8feT7v2Nsak0FVFYBEMFl0UgIh6E3+XuqWQdUtgEQwOXTuasNZ6pwbGJTwo/LUwatKvb7OOYAYJgvbcsZgBF8i0l8lboIqkz8lrRDGkYEmAcUwJIGKqCP5exlVBhJAXDx93eEYqoFZdLodgHItQROJP4h+Mmg6+UO3BG1pByB4S6B2II/sSsEwIYRO+l6yE4lXvuysgJGIIsRYBYAkcBT9lg4X3flvDUfSLnQ2YbvREuGEsG7rz4PuX0uGg4m1CgBJoBBPvP1vQofQCiSCdiIJiEppWgRvPOHxRvc3iUgCBVE1UJwmRPDGEx4/IoDrD91T+/7KEHMrAJKAaBnZ5M8SuwhipvPnCYxKzCsFsS0TVrFiULbc/6djLiodQ1kiqgQa+6vEQgDl24Iq+v3rD90TtEqISAB9kQRGRHMDzaAJv+aQBEStjFMN1CGApquBC3av4oLdq/iT98Z1lmkekoConSIiSCf86qwAmhLBBbtXHfU8dhHotGHRCP1OLW667L/+0D21TBb2Jn6b0OrAmMS4ShDb6kAev//Kc6FDmMU4Uhgn6X+25iMjv6Zi9C1CIfLItgnf3Hd/wEjCoDkBIRLqFkCscwOSwJhoqVBMCpKAEB1HEhCCbs4FpAyVgJmtMbP9ZvZIZuwkM7vLzLYl9/MyP7vWzLab2RYzu6SuwIUQ1VCkEvgScGnP2ErgbndfDNydPMfMzgamgXOS19xsZvGfPD0iFy7dz4VL9/Pib7yVF3/jraHDESVpsgqIcXJwqATc/fvAgZ7hKWBt8ngtcHlmfL27P+PuTwDbgfOrCTVeJAMxCrGJYNw5gdPcfS9Acn9qMr4A2JnZblcyNjFcuHR/35+FlkHo6wyKdlL1yUKWM5Z7NqCZLQOWVbz/KEhF8L+X/0fgSOIixrMFxfiVwD4zmw+Q3KeHx13Awsx2ZwB78t7A3Ve7+3l5pzFOCmll0GR1oGpgNEKtCsTUEowrgU3A0uTxUmBjZnzazI4zs0XAYuC+ciHGw6BWYBhNiaAN3x8QM8QigiJLhF8DfgS83Mx2mdnVwCeBi81sG3Bx8hx33wzcDjwKfBu4xt0P1xV822hCBLFWAmoF4kXfIixImSpgGFXOHcRYCXzsBf965PE/3vbegJHMJoaThBr8dqGuMRgrVVYIsVUCWQEAfOKqNXziqjWBohF5SAIFqLMKSOnaeQYSwfOEnhuQBIbQhABSqhBBTO1AbxUg4kQSiIyyIoilHSgiAFUDcSAJDKDJKiBLGRHc+p2tvOnmb1UYzeiMUgFIBDOEbAm0OjCAUBJIKbpqcOt3tvb92Xfft6SqcApRpgUIsXIQw+pASgOrBLmrA5JAH0ILIEuvDAYlfVWMIo+6ev+6pRCTACCcBHSh0ZbRhACAIy3FMBnUOfmXbRWqFkJsAgiJKoE+xFQJwEw10JQA+pEVQuiZ/7JSiFECqgTEQEILAGaqg+++b0lwAUC9VULX0OpADrFVATEReuUhj1HPQoyxCgiJJNCDBNBedEryeEgCGWIVwA0v/GzoEI7iB5+P+2JRg2SgKmA2mhMg3uQX5dC8QTE6XwlIAOMRezUgitN5CYhu8Imr1qgV6EOnJdCGKiC2+YAsbasGpu5//ZFbbIT8s+WdlUAbBCCq47orlh/1PEYRhKKTEmiLAGKuAiaBWEQQsgqAjq0OKPmrZ3r5QtbfsnP4hoHorQB6mbr/9Wz8sx80FE2cdKoS+NHaU4dvJArz+c0zvz7TyxcO2TIMwwSQErIiCF0FQMckADMiSG8x0pYqIBVASmwiKCqAlFhagxB0TgJZYhVB7PQKICUWEYwqgJSmRRBDFQAdl0BstKEK6CeAWBhXAClNiSAWAYCuJwDEM2EYuwSKCiDURGFZAeRR9aRh4OTXHx/ph9oC0Y+qKoOfrflIaAH0RZVAD6GqgpirgHFagCaqgTqO/MP41MMXNr7PClElUARVBUcT6xyABFAdcf4PB0YimKGMAOpcKQghgElGEhC5xFoBhGJSqwCQBPrSZDVw9nXfa2xfRahKAHVUA6oCqqdT3x0YlVQEwyYLY0viMsRaASj560OrAyPQK4Oqk//Kzz5Y6fuNQx0SKLtSEFoAE9QK6M+QVcHV2zfUvo+QMqizEhhVBqGTHyZKAKAlwvaw7kPnBtlv3a3AKHMEMQigKwz9XzezNWa238weyYx93Mx2m9mDye2yzM+uNbPtZrbFzC6pK/BJp2kRxDoXIOqnyP/8l4BLc8Y/6+7nJrf/BDCzs4Fp4JzkNTeb2Zyqgg1NE61AlqZE0KQAilQDqgKaZej/vrt/HzhQ8P2mgPXu/oy7PwFsB84vEZ+omdgqAAmgecr8Bqwws4eSdmFeMrYAyM7+7ErGZmFmy8zsATN7oEQMjdF0FdBFJIAwjCuBW4CzgHOBvcBnknHL2TZ35t/dV7v7eXmzleJ5Qk0S1kleSyABhGMsCbj7Pnc/7O7PAV/k+ZJ/F5D9Hz4D2FMuxPCErgLqEkHIViAVwXVXLJcAAjPWb4GZzc88fTuQrhxsAqbN7DgzWwQsBu4rF6KA6kUQw1yAkj8OiiwRfg34EfByM9tlZlcDN5nZw2b2EPBG4EMA7r4ZuB14FPg2cI27H64t+gYIXQVkqUoEMQhAxMPQ7w64+7tyhm8dsP0NwA1lghL9Wfehc6M4vVhMDjoktJAyFYGqANGLfiNayjgikABEHvqt6AgSgOiHfjNaTJFq4PObXyABjMmEfYOwL/rtmGCU/OPTFQGAJNB6+lUDbRDAtpO/GToEgSQwkCufncvvX/Ke0GEMZRJPLRbNIQn04cpn5x553DYRtKEKSImxGuhSKwCSQGHaIoI2CUDEga4x2EO2AujHH/ziKw1EMjpf+PB1Rx5f8X/vCBhJcb78uedmjf3zu98WIJIZJrwK0IVG+1Ek8bPEKIGsAFJiFkFe8vcSQgaSQCBCSmBUAaTEJII8AfQSixCKJH+WJkUw4QIASSCfcSUAcYigiACyhJDBqInfS1Mi6KoEOj2LVEYAEH6ycFQBANz+h/9eQyT9KSsAgI99dVMFkQzmnL/IuyhWN+hsJVBWAFlCVQTjSADqrwaqSPx+1FEVZAVw1ecuqPz9I0LtQEqVAoAwEhhXACl1iKDO5M9StQj6VQETKARJoOrkz9KUCMomfy9lZdBU4g9jVDGMU/5PgBS6K4E6kz9L3SKoWgAp44ggluTvZZAMqur7WyyDbkqgKQGk1CWCugTQyyAhxJr4efTKoI6JvxbKoHsSaFoAUI8EmhJASp4I2iSAlDtWTzWynxbJIFcCQy802lZCCABmlg1jOH+gDOky4jPX/2XgSNrBbe//n4E/j10SEymBUAJIqVIETVcBonqykohRCJ0+WSh2JIDxaaoVGJVhVUMIJrISiIEy1YCSvxyxCiAlTwQhKwRJoEZGFYGSX4RA7UDNFP1+gQQgQiEJNEDoLxoJMQhJoCEGiUBVgAiJJBAYCUCERhJokN5qQAIQMaDVgYb5/Uvew9q/WhQ6jInmqs9dEOV6fB4xnDykSkBMFMe/8kYgjuQaRiwxSgJiYkgFkBJLkuURU2ySQACO2bg1dAidIaZkS4ktpqESMLOFZvY9M3vMzDab2QeS8ZPM7C4z25bcz8u85loz225mW8zskjr/AW1j3dxnAYmganqrgCwxJV1MsaQUqQQOAR929z8FLgCuMbOzgZXA3e6+GLg7eU7ys2ngHOBS4GYzm1NH8G0jFUCKRFANgwQghjNUAu6+191/kjz+LfAYsACYAtYmm60FLk8eTwHr3f0Zd38C2A6cX3HcE4NEUI6iAojhCBxDDHmMtERoZmcCrwbuBU5z970wIwozOzXZbAGQXZ/ZlYx1jt4jfz+O2biVQ1MvqzmayWKco3+ahE0vH8aa/CmFJwbN7IXAHcAH3f03gzbNGZt1+TAzW2ZmD5jZA0VjKErR5KuTUWNQRVCcsuV/k0kZuwCgoATMbC4zAljn7l9PhveZ2fzk5/OB/cn4LmBh5uVnAHt639PdV7v7eXnXPGs740romI1bJYOGaENyNkWR1QEDbgUec/dVmR9tApYmj5cCGzPj02Z2nJktAhYD91UXcjHWzX02SEVQxT5TGUgIogmKVAKvA/4aeJOZPZjcLgM+CVxsZtuAi5PnuPtm4HbgUeDbwDXufriW6AsQSgZVIREcTZUrAXVXA22pNoZODLr7D8jv8wHe3Oc1NwA3lIirclIR1HkR0jbLpg20aSmwLQKADp4xmFYGVSdsnQIIUQ3EdrnxugTQpmSti05/i3Dd3GfHrgyaPuqnIujaUmKbjv4pbRNL5yqBXsapCkKW/V2aMGxKAFUmbdsEAJLAEdo2gViXCFLJhD4CN73/NiZvVUgCPQyTQUyiqFIEeRVGaBE0TVkRtFUkkkAf8mQQkwBSyopgWHshEdT7uhiY6L9K3CVGmTAcRxy/e/jakV9ThtDyGeX7BS0SQLf+KnHXqHuy8PhX3tiICEIn/yi0KPkHonZAFKbuBI1JAJOS4EWQBMRI1JWoMQmgCJMkCUlAjEzVCRurAPol+iQJACQBMSZVJW6sAujHpAkAJAFRgrIJ3AYBZJN+EgUAWiIUFTDuqkEbJDBh5C4RSgKiMorKQMkfjFwJqB0QlVEkuSWA+JAERKUMSnIJIE4kAVE5eckuAcSLJCBqIZv0EkDcaGJQiO6giUEhxGwkASE6jiQgRMeRBIToOJKAEB1HEhCi40gCQnQcSUCIjiMJCNFxJAEhOo4kIETHkQSE6DiSgBAdRxIQouNIAkJ0HElAiI4jCQjRcSQBITpOLH+a/Engd8l9GziF9sQK7Yq3TbFCu+J9Sd5gFNcYBDCzB/KufxYjbYoV2hVvm2KF9sWbh9oBITqOJCBEx4lJAqtDBzACbYoV2hVvm2KF9sU7i2jmBIQQYYipEhBCBCC4BMzsUjPbYmbbzWxl6HjyMLMdZvawmT1oZg8kYyeZ2V1mti25nxcotjVmtt/MHsmM9Y3NzK5NPustZnZJJPF+3Mx2J5/vg2Z2WQzxmtlCM/uemT1mZpvN7APJeLSf71i4e7AbMAf4OfBS4Fjgp8DZIWPqE+cO4JSesZuAlcnjlcCnAsX2BuA1wCPDYgPOTj7j44BFyWc/J4J4Pw78Xc62QeMF5gOvSR6fAGxNYor28x3nFroSOB/Y7u6Pu/tBYD0wFTimokwBa5PHa4HLQwTh7t8HDvQM94ttCljv7s+4+xPAdmb+DxqjT7z9CBqvu+91958kj38LPAYsIOLPdxxCS2ABsDPzfFcyFhsO/LeZ/djMliVjp7n7Xpj5ZQFODRbdbPrFFvPnvcLMHkrahbS8jiZeMzsTeDVwL+38fPsSWgKWMxbjcsXr3P01wBLgGjN7Q+iAxiTWz/sW4CzgXGAv8JlkPIp4zeyFwB3AB939N4M2zRmL4fMdSGgJ7AIWZp6fAewJFEtf3H1Pcr8fuJOZEm+fmc0HSO73h4twFv1ii/Lzdvd97n7Y3Z8DvsjzJXTweM1sLjMCWOfuX0+GW/X5DiO0BO4HFpvZIjM7FpgGNgWO6SjM7HgzOyF9DLwFeISZOJcmmy0FNoaJMJd+sW0Cps3sODNbBCwG7gsQ31GkCZXwdmY+Xwgcr5kZcCvwmLuvyvyoVZ/vUELPTAKXMTPr+nPgo6HjyYnvpczM+P4U2JzGCJwM3A1sS+5PChTf15gpoZ9l5kh09aDYgI8mn/UWYEkk8X4ZeBh4iJlEmh9DvMDrmSnnHwIeTG6Xxfz5jnPTGYNCdJzQ7YAQIjCSgBAdRxIQouNIAkJ0HElAiI4jCQjRcSQBITqOJCBEx/l/Pqu0tNkgfeYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "cells_ex_border = cle.exclude_labels_on_edges(cells)\n", "\n", "cle.imshow(cells_ex_border, labels=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "After correcting the label image, we can also correct the parametric image." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAD8CAYAAADAKumpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAYG0lEQVR4nO3df+xddX3H8eerBaooKlhlhZKBim64zEoa3MJi2JwWmBH8w6UkW/iDpPsDDGZbRplb2i5hgWW6kUzNvk4m2xyMqARinJUxDTGZQHEVKJVRhUmlo+ucE0XAlvf+uOfC5dv749zvPT8+n3Nej+Tk+72nn+89n+/pva/v5/M5n/O5igjMzFK2qu0KmJnN4qAys+Q5qMwseQ4qM0ueg8rMkuegMrPkOajMrHGS3iJp98j2Q0kfmli+rnlUks4HrgdWA38TEdfWciAzy5qk1cD3gHdExH+OK1NLi6o48MeAC4CzgEsknVXHscwse+8Cvj0ppACOqenA5wD7IuI7AJJuBi4CHhpXWDo+4DU1VcXMBg4ciojXLfIMb5Li6TJHgj3AMyO7liJiaULxzcBN056vrqA6FXh85PF+4B2jBSRtAbYMHr36xW/NrCY7JrZYyvoJcHmJcn8Ez0TExlnlJB0HvA+4elq5uoJKY/a9ZDCsSNclAOkU33BolgEBx1b7lBcA34iIJ6cVqiuo9gOnjTxeDzxR07GSFjt3vPC9Nm1rsSZmixOVh8YlzOj2QX3TE+4FzpR0RtG02wzcXtOxzKwhwxbVrK3Uc0nHA+8GPj+rbC1BFRGHgSuAncBe4JaI2FPHsVI22poa99gsN6uAl5fYyoiIpyPitRHxf2WOW4uI+GJEvDki3hgR19R1nNw4rCxnw67frK1qnplek2mB5LCyXFXZ9ZtHXYPpNsMwrDzAbjmpYTC9FAdVyxxYlpMapieU4qCqwUq6duN+xuFlqXFQdUSV40/Ln8vBZW0T5a/qVclBlRF3E61tHqOy0hxY1pa2un6enlChpqcdxM4dvZrqEBv687umyvOobMX6EFjDkHJYtcvzqGxhsXNHJ7uDy8NpUlhpd/d+99QMb6FpmoOqIqm0aLo0fjVv62lceYdXtTyYnrFUQmpU7q2rqrp4ZZ7HYVaeB9PNWuJxr/I8RmXWotGwcgtrOnf9MpRity93bbdwZh2/z0Em4NgyqXG42uM6qCwpbYdUGcM69jGwVq2Cl68pUdBBlQ63pqqVQ0j1nQTHtJAaDqpltlH+zbJ9U/F1Z02VWcRVgzd+H//qW31Kd/0q5qt+hW3smCukRg0DK0VupVilBKwusVXMQVWRpMLqqrYrYJ3V0s1+Dirm6+5Nk0RYjQkpt6qsMg4qW5hbUlY3AWtKbGWeSnqNpM9K+pakvZJ+eVJZB1XFWmtVOaSsCdW2qK4HvhQRPwe8jcFngI7V+6Cqqts3avumBgPrKkqFlLt/VomKgkrSq4B3Ap8CiIjnIuIHk8p7ekKNapm+sEDLqc8TFevQ2+kf5a7qrZW0a+TxUkQsjTx+A/DfwN9KehtwH3BlRPx43JM5qBqwfdOCYdWDbp1bfJkov87LoYjYOOXfjwHOBj4YEXdLuh7YCvzxuMK97/o1ZUVdwZLdunk5FKrTu3NZ3WD6fmB/RNxdPP4sg+Aaq9ctqjrGpyrRQAsqlW5L797ouato5byI+C9Jj0t6S0Q8DLwLeGhS+d62qNoIqamtqquorQU1Sdsh0fbxbQWqver3QeAzku4HNgB/OqlgL4OqzZbU2LBqcQyq6bCIDTte2LqiS7/LTBUGVUTsjoiNEfGLEXFxRPzvpLK97vq15YWrgRtarUYj+vImTqUr3Yga7uWbpXctqmTHpVpUR5h0rdU0zfbdg60Xr62WbqFxi8oq05dgGtq+++h929jBDjrcslpF6Vtkqj5sb/TiL94KLRIyfWo9ldHp11mOLSpJjwFPAUeAwxGxUdJJwD8BpwOPAb85bZDM0jHvOEsfw2lcK6p3Ml3h81cj4tDI463AnRFxraStxeMezK3uhz6GE8wfUJ3tAg4XzmtYHV2/i4Abi+9vBC6u4RhWk2lB1MeQGg6UWyHHrh8QwJclBfDXxU2HJ0fEAYCIOCDp9eN+UNIWYMvg0asXrIZVqY+BBNUHUidbVS19pvuihzw3Ip4owugOSd8q+4NFqC0BSKfEgvUwK6Xp1lHnwmp4r1/DFgqqiHii+HpQ0q3AOcCTktYVral1wMEK6mm2EHffKpJbi0rSK4BVEfFU8f17gD8BbgcuBa4tvt5WRUXNVsIBVbHcggo4GbhV0vB5/jEiviTpXuAWSZcB3wU+sHg1zcpLLZw61/1r4arfioMqIr7DYJ3j5fv/h8GSDWaNSi2gOqmlFlWvZqZbd6UeUp2Zrb4KeFmJrWK+18+yl3pIdU5OXT8z66EMB9NtAX1Yi8o6yGNU/eLuimVrdYmtYr0Kqk5dIjZrQ6b3+tkKuetnWfLCec1IpVXlrp9lyS2qfnGLyrLlq37N2MG21ibgOaAsaxVe9Ru3QvCksr0MqjY4oKwTql/hc/kKwWP1NqjqblU5mKyTPOGzGxxQ1mnlF85bK2nXyOOlYrHMUeNWCB7LQVURB5TN0onlXsq3qA5NG3MqHLVCcETcNa5gr4Nq0e6fw8l6p8Ku34QVgh1U44z+hZsVWg4mW0T2ramhCgbTp6wQPFbvg2rUpBdSXz+Vxewo1bWoxq4QPKmwg8qsAZ1pTQ0XzlvQpBWCpx3WZpjnY87NOs+rJ5jN6bq2K9AzLd3r56Aqya2qxFzHCyG1fWerNZmpM90+cFCZlTISUKNSD6tOaaHr58H0OWj3Nl8BbNocXbtxYbV9U3VVMXwLjdlRKhh/GoZXW4HVqW4fVHbVbyWHtTl4rKohHiRPVqyevVXNQWXpqSGk2hrD6swHjxZCcOSY2VvVHFSWhuuYOFBeFYdVBRxU1ks1h9Nyvjq4mBAcXr1q5lY1D6avgK/+LSCBsaflYdXEQHsnlngBQuLIMWVi47lKj+ugst4bDS5PZ5guEM+tPq5EyWqDyl0/s0LdIdWFsapAHGb1zK1qDqoV8jQF66sjHDNzq5q7fmZWWiCO1HGPzAwOKjM8NlWWg8rMkheIZykzmF6tmWNUkm6QdFDSgyP7TpJ0h6RHiq8njvzb1ZL2SXpYUif/TsXOHcTOHY3PAbJ6NNmayn1AfdCian6Mqsxg+qeB85ft2wrcGRFnAncWj5F0FrAZeGvxMx+X1Hw7sWkOLJtD7mF1hNUzt6rNDKric7a+v2z3RcCNxfc3AheP7L85Ip6NiEeBfQw+AqczYueUF5kDyzpuOEbVdFCttI12ckQcAIiIA8UHCAKcCnx9pNz+Yt9RJG0BtgwevXqF1UjUMKyuarUW6XGIZy+glnlSs1Q9j0pj9sW4ghGxFBEbB5+menzF1UhEAzfa2mLautqXb/cv3TGqcZ6UtA6g+Hqw2L8fOG2k3HrgiZVXLy1Tu32zOKxsmRzDKhDPcdzMrSxJqyX9u6QvTCu30qC6Hbi0+P5S4LaR/ZslrZF0BnAmcM8Kj9E9fQ6rPv/uHVLDGNWVwN5ZhWa20STdBJwHrJW0H9gGXAvcIuky4LvABwAiYo+kW4CHgMPA5RFxZJ5ap2qh1tSocW/YPoxlDX9HB9ZL5LaqwvBevypIWg/8BnAN8LvTys4Mqoi4ZMI/vWtC+WuKA1tZ19GPsAIHVgdUOAb1l8AfACfMKuibkkuorDU1Td/euH0J5hJyGquao+u3VtKukW3L6PNIei9wMCLuK3Nc30IzQyMhNdSnlpVlaY5baA4NruhPdC7wPkkXMvhcm1dJ+oeI+K1xhd2iSk2fWlYO5exUdQtNRFwdEesj4nQGd7P866SQAgfVVI22pkY5rHonl+5fWzPTHVSp6ltYObCyUXVQRcRXI+K908p4jGqC1lpTo/p2K87y37PmsPYaVPOrcnrCPBxUlq7R4Ko4tBxSKzMco2qagyoHvhrYeGvLxhveQtM0B5XlqcbWlk3mpYgTksT4lJU352x3d/sW4zGqBDikMubbc2rnMaoEJBtSfR+fmteUwHJrajHu+rUo2YCyxXgcqxYOqhY4pMzKe55VPMuaxo/b+6CynrgKtu9uuxLd4BZVw7JoTXl8qjLbN4x8v7utWoyXy+J5bY1R9fZevyxCymozGlo2H9+U3JBsQsqtqVqlEla5tKbgxXv9Zm1V61XXzwFly23fkF43MGWDW2iaH0zvVYtKm/L5y2XNabNllVNrCvL7pORsjYZVki0st6Za4ZZVOW0t89KrFtVybmHZqKZbVrm1poZy+qRkq4NbU61rKqxyDSl3/VqiTdvS7AJaa+qcb5VrQA15HlWL3AW0SapqYe1gW/YhBS9+XNasrWq9b1ENDcOqtdaVu33JGg0r7c4/bBbR1jIvblEt49aVTdL3kBryzPREOKzMxqtqMF3SyyTdI+mbkvZImtqVcdfPrAS3pgYqnEf1LPBrEfEjSccCX5P0zxHx9XGF3aKaoNFWlcenLCMVfaR7RMSPiofHFltMKu8W1RSlB9gdNNYTz7Oqso/LkrQauA94E/CxiLh7UlkHVQkTA8sB1Qvu9r1Uya7fWkm7Rh4vRcTSaIGIOAJskPQa4FZJvxARD457MgfVHLRpG7HBk0P7xCH1UnNMTzgUERtLPWfEDyR9FTgfGBtUHqMys9IqvOr3uqIlhaSXA78OfGtS+ZlBJekGSQclPTiyb7uk70naXWwXjvzb1ZL2SXpYkj+cyKxjKppHtQ74iqT7gXuBOyLiC5MKl2nDfRr4K+Dvlu3/i4j489Edks4CNgNvBU4B/kXSm4u+aPbc7bO+e764hWZREXE/8Pay5We2qCLiLuD7JZ/vIuDmiHg2Ih4F9gHnlK2MmaVO2S3zcoWk+4uu4YnFvlOBx0fK7C/2HUXSFkm7BlcGnl6gGs1wa8qsvWVeVhpUnwDeCGwADgAfKfZrTNmxk7giYikiNg6uDBy/wmqYWdOyWY8qIp4cfi/pk8BwEGw/cNpI0fXAEyuuXSLcmjIbyGopYknrRh6+nxfnPtwObJa0RtIZwJnAPYtV0cxSES2NUc18Rkk3AecxmGm6H9gGnCdpA4Nu3WPA7wBExB5JtwAPAYeBy3O/4ufWlNmLBh+XVf3CeLPMDKqIuGTM7k9NKX8NcM0ilTKzNAXiyPPNd/18C42ZlRdw+LCDyswSFiGOHG4+NhxUZlbaIKjcojJLhldOGCNwUJmlwiE1XsQqnntmTePHdVCZWXkBuEWVFs+hMlsm5KBKiUOqv9ztmyKAw+Nu6a2Xg8rM5nO4+UM6qJZxS6rf3Jqa4XngmeYP66DC4WRWWgA/bf6wvQ8qh5QNuTVVQgAtLDPQ+6Ayszl5jKpZbk2ZzSloJah6+7l+Dilbzq+JEoZBNWurWC9bVH5B2iSjrw2PWY3hq371c0DZPJa/XhxchQpaTJJOY/BZoT/DIP6WIuL6SeV7EVQOKKvC8HXU68CqbnrCYeD3IuIbkk4A7pN0R0Q8NK5w54PKIWVV63VgVTQ9ISIOMPioPSLiKUl7GXwGaP+CyiFldeplYJW/6rd28OHCL1iKiKVxBSWdzuDj3e+e9GSdDSqHlDVl1mutU0FWfjD90ODDhaeT9Ergc8CHIuKHk8p1MqgcUpaSzl1JrGj6gaRjGYTUZyLi89PK9nYelVkbsv8jWtE8Kkli8LF7eyPio7PKd7JFZZaycWGVTUurupnp5wK/DTwgaXex7w8j4ovjCjuozKy8iqYnRMTXgNIr8DmozKw8r55gZskLfAuNmSWupdUTHFRmVp5X+DSz5HmMyqwa2r0tm/lK2UxLGOWF88wWM3zj5xAAOdTxKC0tnOegss5Y/sZPOQhSrttUw3v9Zm0Vc1C1INsXaYZSPNcp1qk0t6j6IaeuSU6mnc+UznVKdVmxFINK0mmSviJpr6Q9kq4s9p8k6Q5JjxRfTxz5masl7ZP0sKRN1Vc7Tzl1TXLi89ig4fSEWVvFyrSohkuG/jzwS8Dlks4CtgJ3RsSZwJ3FY4p/2wy8FTgf+Lik1dVXvRv8JltM2fOXwnlOoQ4LG05PmLVVbOb0hClLhl4EnFcUuxH4KnBVsf/miHgWeFTSPuAc4N+qrnzq5nkT5XI5PRUredMPf6bpc92JgBpq6VNo5hqjWrZk6MlFiA3D7PVFsVOBx0d+bH+xb/lzbZG0a7Bc6dMrqPqUeibwwpi3DinUOReLnqsmz3Xn/l8T7voB5ZcMZfzSDXHUjoiliNg4WK70+LLVyMJKX5zava17L+xE+TwvIMWuH0xcMvRJSesi4oCkdcDBYv9+4LSRH18PPFFVhcvKuZk/+hzuElpSUv1I9ylLht4OXFp8fylw28j+zZLWSDoDOBO4p7oqzyf3VkrOda9Dleej7nPbyf+7hD/SfeySocC1wC2SLgO+C3wAICL2SLqFwedzHQYuj4gWbmN8qSZaWJ18YSYkp/ObU13n0tLqCYo4avio+UrolIAtjR+3ytCq+4XZ9y5gnee3jnObZlDtuK/MR1hNo+M3Bm/aNbvgA1r4WKN6vXrCItMCmn4htjXm1rY03/DT5Vjn0lIdo+q6lYxhtflCzH3MbR5N/Z45jXu1LvXpCV2XWwDUVdfheWj7XLTVYrUZWpqZ7qBaZtabNKUXdNUtgb7fi5jTRNLWVHjVT9INkg5KenBWWQfVBLm8cat4c+USzE1YZLJub1Q3PeHTDO4HnqnXg+ll5PACnPeiwErG5Po0iF/3+cxahff6RcRdxW15MzmoOqKJyYtNhFVOb/qc6lqZ8lf91g7u433BUkQsrfSw7vpZaX2ayZ1SXZJSfozq0PBe3mJbcUiBg8rmVOfVxpzkVt/KeHqC5aLqN2mqb/pJ9Uq1vo3x9ATLRVVv1tze9LnVtxZRYitB0k0MFtR8i6T9xX3DYzmobMX6MO9otI451DcnEXFJRKyLiGMjYn1EfGpSWQeVLaQPb94UZur3nacn2MLmvWHab/qctbPOi1tUVpkyAeSQyl07K+c5qKxSvh2n69qZn+CgssqNCySHVFc8D/ykxFYtB5XVwlfLuqqdFpUH0602Dqiuan6JTweVmc2hnat+Diozm0M7i6Y7qMxsDm5RmVnyhlf9muWgMrM5uOtnZslz18/MkucWlZklzy0qM0te4MF0M0ucW1RmljyPUZlZ8tyiMrPkuUVlZslzi8rMktfOLTReOM/M5lDdmumSzpf0sKR9krZOK+sWlZnNoZqun6TVwMeAdwP7gXsl3R4RD40r7xaVmc2hshbVOcC+iPhORDwH3AxcNKlwIi2qA4dgx4+BQ23XpKS15FNXyKu+OdUV8qrvzy7+FAd2wva1JQq+TNKukcdLEbE08vhU4PGRx/uBd0x6siSCKiJeJ2lXRGxsuy5l5FRXyKu+OdUV8qvvoiLi/IqeSuOeflJhd/3MrA37gdNGHq8HnphU2EFlZm24FzhT0hmSjgM2A7dPKpxE16+wNLtIMnKqK+RV35zqCvnVNwkRcVjSFcBOYDVwQ0TsmVReERO7hWZmSXDXz8yS56Ays+S1HlTzTKNvi6THJD0gafdwboikkyTdIemR4uuJLdXtBkkHJT04sm9i3SRdXZzrhyVtSqS+2yV9rzi/uyVdmEJ9JZ0m6SuS9kraI+nKYn+y57ezIqK1jcEg2reBNwDHAd8EzmqzThPq+Riwdtm+PwO2Ft9vBa5rqW7vBM4GHpxVN+Cs4hyvAc4ozv3qBOq7Hfj9MWVbrS+wDji7+P4E4D+KOiV7fru6td2immsafWIuAm4svr8RuLiNSkTEXcD3l+2eVLeLgJsj4tmIeBTYx+D/oDET6jtJq/WNiAMR8Y3i+6eAvQxmVCd7fruq7aAaN43+1JbqMk0AX5Z0n6Qtxb6TI+IADF7QwOtbq93RJtUt5fN9haT7i67hsCuVTH0lnQ68HbibPM9v1toOqrmm0bfo3Ig4G7gAuFzSO9uu0Aqler4/AbwR2AAcAD5S7E+ivpJeCXwO+FBE/HBa0TH7Uji/2Ws7qOaaRt+WiHii+HoQuJVBc/5JSesAiq8H26vhUSbVLcnzHRFPRsSRiHge+CQvdpdar6+kYxmE1Gci4vPF7qzObxe0HVRzTaNvg6RXSDph+D3wHuBBBvW8tCh2KXBbOzUca1Ldbgc2S1oj6QzgTOCeFur3EsM3feH9DM4vtFxfSQI+BeyNiI+O/FNW57cT2h7NBy5kcDXl28CH267PmPq9gcGVnG8Ce4Z1BF4L3Ak8Unw9qaX63cSgu/RTBn/RL5tWN+DDxbl+GLggkfr+PfAAcD+DN/u6FOoL/AqDrtv9wO5iuzDl89vVzbfQmFny2u76mZnN5KAys+Q5qMwseQ4qM0ueg8rMkuegMrPkOajMLHn/DzuU+VrXGWTbAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "neighbor_count_image_ex_border = neighbor_count_image * (cells_ex_border != 0)\n", "\n", "cle.imshow(neighbor_count_image_ex_border, color_map='jet', colorbar=True, min_display_intensity=0)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now, we can measure the number of neighbors. We can either just read those numbers and put them in a list ..." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "cl.Array([[0., 5., 6., 7., 6., 6., 7., 6., 6., 5., 5., 6., 6., 7., 6., 6.,\n", " 6., 6., 6., 6., 6.]], dtype=float32)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "cle.read_intensities_from_map(cells_ex_border, neighbor_count_image_ex_border)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "... we can also read these values together with all other statistics and put them in a pandas DataFrame." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
labelnumber_of_neighborscentroid_xcentroid_y
015.074.59269723.731028
126.0173.74137935.582230
237.049.46091551.928082
346.091.78376852.666199
456.0142.72821062.962471
567.0189.97320668.526794
676.073.69918188.797356
786.0115.32108390.665649
895.0163.09773399.824242
9105.058.206852122.341309
10116.098.847565134.280228
11126.0147.731705127.326218
12137.0200.522903112.817368
13146.067.880905162.991592
14156.0122.943306157.515717
15166.0171.010025153.333328
16176.043.426266189.492981
17186.096.272728192.326477
18196.0140.815964194.446030
19206.0193.817596193.019318
\n", "
" ], "text/plain": [ " label number_of_neighbors centroid_x centroid_y\n", "0 1 5.0 74.592697 23.731028\n", "1 2 6.0 173.741379 35.582230\n", "2 3 7.0 49.460915 51.928082\n", "3 4 6.0 91.783768 52.666199\n", "4 5 6.0 142.728210 62.962471\n", "5 6 7.0 189.973206 68.526794\n", "6 7 6.0 73.699181 88.797356\n", "7 8 6.0 115.321083 90.665649\n", "8 9 5.0 163.097733 99.824242\n", "9 10 5.0 58.206852 122.341309\n", "10 11 6.0 98.847565 134.280228\n", "11 12 6.0 147.731705 127.326218\n", "12 13 7.0 200.522903 112.817368\n", "13 14 6.0 67.880905 162.991592\n", "14 15 6.0 122.943306 157.515717\n", "15 16 6.0 171.010025 153.333328\n", "16 17 6.0 43.426266 189.492981\n", "17 18 6.0 96.272728 192.326477\n", "18 19 6.0 140.815964 194.446030\n", "19 20 6.0 193.817596 193.019318" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "statistics = cle.statistics_of_labelled_pixels(neighbor_count_image_ex_border, cells_ex_border)\n", "\n", "table = pd.DataFrame(statistics)\n", "\n", "# rename a column\n", "table = table.rename(columns={\"mean_intensity\": \"number_of_neighbors\"})\n", "\n", "# only filter out a subset of all columns; only what we care\n", "table = table[[\"label\", \"number_of_neighbors\", \"centroid_x\", \"centroid_y\"]]\n", "\n", "table" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Exercise\n", "Analyze a larger field of view with more cells and vary the parameters `random_sigma_x` and `random_sigma_y` of the `artificial_tissue_2d` function. Use a touching-neighbor-count map to count the number of touching neighbors before and after applying a median filter to the map." ] }, { "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 }