Trouvez l'aire de l'ensemble somme des rectangles qui se chevauchent

Overview

a1.png

Comme ça, il y a des moments où vous voulez trouver la zone de la partie de l'ensemble de somme de rectangles qui se chevauchent au hasard.

Vous voyez, ce genre de chose qui est populaire ces jours-ci. (Image empruntée ici: https://note.nkmk.me/python-opencv-face-detection-haar-cascade/) img

Je peux gérer cela en dessinant une ligne auxiliaire et en la calculant manuellement, mais il est assez difficile de la calculer correctement https://mathcommunication.hatenablog.com/entry/2016/10/11/003213#%E5%B9%B3%E9%9D%A2%E5%9B%B3%E5%BD%A2%E3%81%AB%E5%AF%BE%E3%81%99%E3%82%8B%E5%8F%AF%E6%B8%AC%E6%80%A7%E3%81%A8%E9%9D%A2%E7%A9%8D

NP difficile? C'est impossible pour les gens ordinaires, alors

politique

Dessinez directement et comptez les pixels.

Ce ne sera pas aussi beau que lorsqu'il est calculé à l'aide d'une formule, et il y aura une erreur dans la résolution, Je pense qu'il n'y a pas beaucoup de cas où cela devient un problème.

En particulier,

Binariser

b1.png

Compter les pixels noirs (couleur de remplissage)

Puisque opencv a une fonction appelée countNonZero, j'utilise ceci. https://docs.opencv.org/2.4/modules/core/doc/operations_on_arrays.html#countnonzero

code

Allons-y.

requirements.txt


cycler==0.10.0
kiwisolver==1.2.0
matplotlib==3.2.1
numpy==1.18.4
opencv-python==4.2.0.34
Pillow==7.1.2
pyparsing==2.4.7
python-dateutil==2.8.1
six==1.14.0

Puisqu'il est aléatoire, la position et la taille du rectangle changeront à chaque exécution.

a.py


from PIL import Image, ImageDraw
import numpy as np
import cv2

import random
import matplotlib


rects = [
    (random.randrange(200),
        random.randrange(200),
        random.randrange(50, 120),
        random.randrange(50, 120)
    ) for x in range(20)
]

colors = random.sample(matplotlib.colors.cnames.keys(), len(rects))

#Pour le dessin au trait
im1 = Image.new('RGB', (200, 200), (255, 255, 255))
d1 = ImageDraw.Draw(im1)
#Pour la confirmation de conversion binaire
im2 = Image.new('RGB', (200, 200), (255, 255, 255))
d2 = ImageDraw.Draw(im2)

#Dessin rectangulaire
for r, c in zip(rects, colors):
    print(c, r)
    d1.rectangle(r, outline=c)
    d2.rectangle(r, fill=(0, 0, 0))

im1.save("a1.png ")
im2.save("a2.png ")

#Conversion binaire
nim = np.array(im2, dtype=np.uint8)
gim = cv2.cvtColor(nim, cv2.COLOR_RGB2GRAY)
ret, bwim = cv2.threshold(gim, 0, 255, cv2.THRESH_OTSU)

#Image après conversion binaire
cv2.imwrite("b1.png ", bwim)

whitePixels = cv2.countNonZero(bwim)
blackPixels = bwim.size - whitePixels

print("-"*50)
print(whitePixels, blackPixels, bwim.size)
print("filled area rate:", blackPixels / bwim.size)

Courir

bash


$ python a.py
gainsboro (188, 158, 61, 63)
palegreen (174, 182, 118, 53)
chartreuse (75, 34, 82, 81)
olive (170, 10, 93, 111)
mistyrose (183, 65, 78, 52)
indianred (25, 55, 51, 109)
mediumaquamarine (102, 155, 59, 89)
navy (169, 127, 85, 96)
darkviolet (103, 152, 76, 51)
yellowgreen (84, 6, 73, 118)
deeppink (42, 61, 95, 90)
bisque (47, 101, 107, 105)
dimgrey (95, 179, 91, 70)
maroon (146, 182, 88, 62)
whitesmoke (189, 118, 83, 112)
lavender (71, 149, 104, 117)
burlywood (36, 82, 105, 82)
darkslategrey (122, 128, 103, 52)
mediumslateblue (163, 121, 84, 92)
silver (144, 126, 61, 58)
--------------------------------------------------
18580 21420 40000
fill area rate: 0.5355

fin.

référence

https://note.nkmk.me/python-pillow-imagedraw/ https://techtech-sorae.com/pythonopencv%E3%81%A7%E4%BA%8C%E5%80%A4%E7%94%BB%E5%83%8F%E3%81%8B%E3%82%89%E7%99%BD%E3%81%A8%E9%BB%92%E3%81%AE%E9%9D%A2%E7%A9%8D%E6%AF%94%E3%82%92%E7%AE%97%E5%87%BA/

Recommended Posts

Trouvez l'aire de l'ensemble somme des rectangles qui se chevauchent
Comment trouver la zone du diagramme de Boronoi
Trouvez la définition de la valeur de errno
Trouvez le ratio de la superficie du lac Biwa par la méthode de Monte Carlo
Version Migemo de la commande: find ,: mfind
Trouvez le coefficient du polypole le moins carré
Trouver la main de "Millijan" par l'optimisation des combinaisons
Trouvez le nombre de jours dans un mois
Découvrez la fraction de la valeur saisie en python
Trouvez la solution de l'équation d'ordre n avec python
Découvrez le jour par date / heure
Projecet Euler 12 Trouvez le nombre de fractions sans diviser.
Calculer la somme des valeurs uniques par tabulation croisée des pandas
Connaissez l'emplacement du fichier de définition de classe Python.
Découvrez la version de la langue que vous utilisez
Découvrez l'emplacement des packages installés avec pip
Le début de cif2cell
Trouvez le maximum de Python
Le sens de soi
Union Find sur networkX
le zen de Python
L'histoire de sys.path.append ()
La vengeance des types: la vengeance des types
J'ai essayé de trouver l'entropie de l'image avec python
Comment trouver le nombre optimal de clusters pour les k-moyennes
Maya | Découvrez le nombre de polygones dans l'objet sélectionné
Découvrez la largeur apparente d'une chaîne en python
J'ai essayé de trouver la moyenne de plusieurs colonnes avec TensorFlow
[Python] Fonctionnalisation de la formule de Heron et calcul de la surface maximale
Hériter de la bibliothèque standard pour trouver la valeur moyenne de Queue
Python --Trouvez le nombre de groupes dans l'expression regex
Trouver l'index de la valeur maximale (valeur minimale) d'un tableau multidimensionnel
Comment trouver le coefficient de mise à l'échelle d'une ondelette bipolaire
Trouver le diamètre du graphique par recherche de priorité de largeur (mémoire Python)
Trouvez l'écart moyen / standard des valeurs de luminosité dans l'image
Trouver les valeurs propres d'une vraie matrice symétrique en Python