[PYTHON] Finden Sie das Verhältnis der Fläche des Biwa-Sees nach der Monte-Carlo-Methode

Discraimer Dieser Beitrag enthält Inhalte, die so gut oder besser sind als die technischen Inhalte, die Shiga-Bürger der Präfektur Shiga lieben. Wenn Sie reine technische Artikel lesen möchten, wenden Sie sich bitte an andere. Es ist technisch nichts Erstaunliches passiert.

Shiga Bürger

Hallo. Ich bin Bürger der Präfektur Shiga. Woran denkst du, wenn du die Nummer 1/6 hörst?

Tatsächlich ist dies das Verhältnis der Fläche des Biwa-Sees zur Gesamtfläche der Präfektur Shiga. Diese Zahl ist ein Wert, den jeder in der Präfektur, der in der fünften Klasse der Grundschule oder höher ist, kennen sollte, und ein guter Teil der Bürger der Präfektur weiß es.

Weil in der Präfektur Shiga ein Geschäft namens ** Floating School ** entwickelt wird und Fünftklässler in der Präfektur auf einem Lernboot übernachten können, das auf dem Biwa-See schwimmt (es ist als Präfekturschule positioniert). Wir werden Schulungsaktivitäten durchführen, einschließlich. Dieses Schiff heißt ** Uminoko ** und im Jahr 2018 hat sich das seit über 30 Jahren beliebte Schiff zurückgezogen und ist jetzt ein neues Schiff, das erneuert wurde.

Es ist eine ziemlich gute Erfahrung, also ist es eines der guten Dinge an Shiga-Bürgern.

Bitte schauen Sie einmal. : https://uminoko.jp

Übrigens ist es interessant, dass viele Grundschüler, die an Bord dieses Seeigels gehen, sagen, dass sie noch das Lied "Ship of Hope" singen können, an das sie sich vor dem Training erinnern. Natürlich bin ich einer von ihnen.

Auf diesem Schiff beobachten Grundschüler in der Präfektur Plankton, erleben Seilarbeiten und führen Aktivitäten mit einem Cutter (kleines Schiff) durch. Aus diesem Grund lerne ich vor und nach dem Tod, aber einer der Werte, die häufig auftreten, ist dieser 1/6.

Vor kurzem hat Kokuyo einen ** Rokubunnoichi Nocho ** als Antwort auf diese 1/6 Namenserkennung veröffentlicht. Es wird Papier verwendet, das aus Yoshi an der Küste des Biwa-Sees hergestellt wurde, und es ist ziemlich gut. Kaufen Sie es also bitte.

https://www.kokuyo-shop.jp/sp/ProductDetail.aspx?sku=4560107742501

Kommen wir nun zum Hauptthema.

Wie ich bisher erklärt habe, ist es 1/6, aber leider sagen das alle Leute von Kyoto und Osaka. ** "Ist Shiga oder der halbe Biwa-See?" **

Es ist natürlich, so zu denken. Sie sehen wahrscheinlich nicht die reichen und riesigen Wälder, die den Biwa-See umgeben.

Wir können es jedoch nicht wie erwartet akzeptieren.

Die Leute von Osaka und Kyoto können das also. Wir schlagen vielmehr ein Flächenverhältnis-Experiment vor, das überall auf der Welt durchgeführt werden kann.

Monte-Carlo-Methode

Wenn Sie diesen Artikel lesen, wissen Sie. Es ist auch eine sehr beliebte Methode zum Ermitteln des Umfangsverhältnisses. Bei der Berechnung des Umfangsverhältnisses ** das Verhältnis der Fläche des Kreises, die im Quadrat enthalten sein soll **.

Ja, das Verhältnis der Fläche in einer bestimmten Fläche kann probabilistisch bis zu einer bestimmten Genauigkeit berechnet werden. Daher wird die Fläche des Biwa-Sees in der Präfektur Shiga nach der Monte-Carlo-Methode berechnet.

Holen Sie sich die Karte

Bereiten Sie zunächst eine Karte vor, die als Grundlage für die Bestimmung des Gebiets des Biwa-Sees und des Gebiets der Präfektur Shiga dient. Dies ist vorerst in Ordnung.

https://www.freemap.jp/itemFreeDlPage.php?b=shiga&s=shiga 1.png

Kartenverarbeitung

Als nächstes färben Sie die weiße Karte. Dies ist ein Prozess, mit dem die Software das ** Nicht-Shiga-Präfekturgebiet ** ** Shiga-Präfektur-Landgebiet ** ** Shiga-Präfektur-Biwako-Gebiet ** auf der Karte identifizieren kann. Bitte füllen Sie es entsprechend aus. Dieses Mal benutzte ich Seashore, um ** rgb (255,38,0) im Landgebiet der Präfektur Shiga ** und rgb (0,0,255)` im Gebiet des Biwa-Sees in der Präfektur Shiga **. Ich habe es gefärbt.

map.png

Puyoguramingu

Im Grunde ist es dasselbe wie das Umfangsverhältnis zu finden, einen zufälligen Punkt in dem Gebiet festzulegen (Quadrat oder Präfektur Shiga) und zu bestimmen, ob es sich in dem angegebenen Gebiet befindet (Kreis oder Biwa-See). Wiederholen. Ich glaube nicht, dass es dir wirklich wichtig ist, weil es dasselbe ist wie das, was du mit einem Kreis machst, aber was du dir ein wenig Sorgen machen musst, ist, dass der Nenner ** nicht die Epochenzahl ist. ** **. Teilen Sie die Anzahl der Personen, die den Biwa-See ** betreten haben, durch die Anzahl der Personen, die die Präfektur Shiga betreten haben **. Es ist farbig, um dies zu identifizieren. Wie Sie alle wissen, können Sie die Farbdaten einer bestimmten Koordinate im Bild erhalten, indem Sie mit OpenCV auf die Adresse "[y, x]" des Arrays zugreifen, das von "imread" angezeigt wird. Ich glaube nicht, dass es noch etwas gibt, worüber man sich Sorgen machen muss. Derzeit ist die Anzahl der Versuche auf 50.000 festgelegt.


import cv2
import numpy as np

print("Biwako Area inspector")
print("©︎ Yusei ito 2020")
print("=========================")


print("Importing image..")
img = cv2.imread('./map.png') #Processed map image file
height, width, channels = img.shape[:3]
print("Image Imported.")
print("\tWidth:"+str(width))
print("\tHeight:"+str(height))

count_max=50000
count_shiga=0
count_biwako=0

for i in range(0,count_max):
    px=img[np.random.randint(height),np.random.randint(width)]#Be attention. Array is formed as [y,x] and return is formed [blue,green,red]
    if (px[0]==0 and px[1]==38 and px[2]==255):
        count_shiga+=1
    if (px[0]==255 and px[1]==0 and px[2]==0):
        count_shiga+=1
        count_biwako+=1
    if (count_max%10==0):
        devisionBy=1 if count_shiga==0 else count_shiga
        print(str(i)+"\t"+str(count_biwako/devisionBy))


print("All Process succeed.\r\n\r\n")    
print("Attempt:"+str(count_max))
print("Shiga:"+str(count_shiga))
print("Biwako:"+str(count_biwako))

print("Overall Score:"+str(count_biwako/count_shiga))


Da das Ausführungsergebnis redundant ist, fügen Sie es erst ganz am Ende ein. Da es sich um Rokubunnoichi handelt, beträgt der erwartete Wert 1/6 = 0,1666666 ....


Attempt:50000
Shiga:21143
Biwako:3493
Overall Score:0.16520834318687036

Es ist gut. Es ist sehr gut.

Deshalb ist die Präfektur Shiga ein guter Ort. Bitte komm einmal.

Ein vollständiger Satz von Artikeln zu diesem Beitrag ist auf GitHub zusammengefasst. Passende siehe bitte.

https://github.com/YuseiIto/Biwako_Rate

Recommended Posts

Finden Sie das Verhältnis der Fläche des Biwa-Sees nach der Monte-Carlo-Methode
Geschwindigkeitsvergleich jeder Sprache nach der Monte-Carlo-Methode
Erhöhen Sie die Geschwindigkeit der Monte-Carlo-Methode zur Implementierung von Cython-Ausschnitten.
Mit Reiskörnern bestreuen, um das Umfangsverhältnis zu ermitteln (Monte-Carlo-Methode)
Schätzung von π nach der Monte-Carlo-Methode
Suchen Sie den Bereich des Summensatzes überlappender Rechtecke
Ich konnte pypy3.6-7.3.1 nicht mit macOS + pyenv installieren, aber ich konnte pypy3.6-7.3.0 installieren. Ich fühlte den Wind von Pypy nach der Monte-Carlo-Methode.
So finden Sie den Bereich des Boronoi-Diagramms
Dominion-Kompressionsspiel nach Monte-Carlo-Methode analysiert
#Monte Carlo-Methode zum Ermitteln des Umfangsverhältnisses mit Python
Die Kopiermethode von pandas.DataFrame ist standardmäßig Deep Copy
Verstehen Sie die Metropolitan Hasting-Methode (eine der Methoden in der Markov-Ketten-Monte-Carlo-Methode) mit der Implementierung
Ermitteln Sie den Durchmesser des Diagramms anhand der Suche nach Breitenpriorität (Python-Speicher).
Einführung in die Monte-Carlo-Methode
[Wissenschaftlich-technische Berechnung nach Python] Monte-Carlo-Simulation nach der Metropolenmethode der Thermodynamik des 2D-Anstiegsspinsystems
Python: Ermitteln Sie die erforderliche Plattendicke für den Grenzknickdruck des freiliegenden Rohrs nach der Brent-Methode
cv2.Canny (): Macht die Anpassung der Kantenerkennung mit der Canny-Methode angenehm
Finden Sie die mysteriöse Veränderung der Pokemon-Bilderbuchbeschreibung nach Levenstein-Entfernung heraus
Finden Sie die Definition des Wertes von errno
Simulieren Sie die Monte-Carlo-Methode in Python
Stellen Sie das Verhältnis von Topcoder, Codeforces und TOEIC nach Bewertung grafisch dar (Pandas + Seaborn).
Ich habe versucht, die optimale Route des Traumlandes durch (Quanten-) Tempern zu finden
Ermitteln Sie den Mindestwert der Funktion mithilfe der Partikelgruppenoptimierungsmethode (PSO).