[PYTHON] Trouvez le ratio de la superficie du lac Biwa par la méthode de Monte Carlo
Discraimer Ce message contient un contenu aussi bon ou meilleur que le contenu technique que les citoyens de Shiga adorent dans la préfecture de Shiga. Si vous souhaitez lire des articles purement techniques, veuillez contacter les autres. Rien d'étonnant techniquement n'est arrivé.
Citoyens de Shiga
Bonjour. Je suis citoyen de la préfecture de Shiga. À quoi pensez-vous lorsque vous entendez le chiffre 1/6?
En fait, c'est le rapport entre la superficie du lac Biwa et la superficie totale de la préfecture de Shiga. Ce nombre est une valeur que tout le monde dans la préfecture qui est en cinquième année du primaire ou au-dessus devrait connaître, et une bonne proportion des citoyens de la préfecture le savent.
Car, dans la préfecture de Shiga, une entreprise appelée ** Floating School ** est en cours de développement, et les élèves de cinquième de la préfecture peuvent passer la nuit sur un bateau d'apprentissage flottant sur le lac Biwa (il se positionne comme une école préfectorale). Nous réaliserons des activités de formation notamment. Ce navire s'appelle ** Uminoko **, et en 2018, le navire qui était populaire depuis plus de 30 ans a pris sa retraite et est maintenant un nouveau navire qui a été renouvelé.
C'est une assez bonne expérience, c'est donc l'une des bonnes choses à propos des citoyens de Shiga.
Veuillez jeter un œil une fois. : https://uminoko.jp
En passant, il est intéressant que de nombreux élèves du primaire qui montent à bord de cet oursin disent qu'ils peuvent encore chanter la chanson "Ship of Hope" dont ils se souviennent avant la formation. Bien sûr, je suis l'un d'entre eux.
Sur ce bateau, les élèves des écoles élémentaires de la préfecture observeront le plancton, expérimenteront le travail de la corde et réaliseront des activités avec un cutter (petit bateau). Pour cette raison, je fais avant et après l'apprentissage de la mort, mais l'une des valeurs qui apparaît souvent est ce 1/6.
Récemment, Kokuyo a publié un ** Rokubunnoichi Nocho ** en réponse à cette reconnaissance de nom 1/6. Il utilise du papier produit à Yoshi sur la côte du lac Biwa, et c'est assez bon, alors achetez-le.
https://www.kokuyo-shop.jp/sp/ProductDetail.aspx?sku=4560107742501
Maintenant, passons au sujet principal.
Comme je l'ai expliqué jusqu'à présent, il est de 1/6, mais malheureusement, les habitants de Kyoto et d'Osaka le disent tous. ** "Shiga ou la moitié du lac Biwa?" **
Il est naturel de le penser. Ils ne voient probablement pas les forêts riches et vastes qui entourent le lac Biwa.
Cependant, nous ne pouvons pas l'accepter comme prévu.
Ainsi, les gens d'Osaka et de Kyoto peuvent le faire. Nous proposons plutôt une expérience de rapport de surface qui peut être effectuée n'importe où dans le monde.
Méthode de Monte Carlo
Si vous lisez cet article, vous le savez. C'est également une méthode très populaire pour trouver le rapport de circonférence. Lors du calcul du rapport de circonférence, ** le rapport de l'aire du cercle à contenir dans le carré **.
Oui, le rapport de la zone dans une certaine zone peut être calculé de manière probabiliste jusqu'à une certaine précision. Par conséquent, la superficie du lac Biwa dans la préfecture de Shiga est calculée par la méthode de Monte Carlo.
Obtenez la carte
Tout d'abord, préparez une carte qui servira de base pour déterminer la superficie du lac Biwa et la région de la préfecture de Shiga. C'est bien pour le moment.
https://www.freemap.jp/itemFreeDlPage.php?b=shiga&s=shiga
Traitement de la carte
Ensuite, coloriez la carte blanche. Il s'agit d'un processus qui permet au logiciel d'identifier la ** zone non-préfecture de Shiga ** ** la zone terrestre de la préfecture de Shiga ** ** la zone de Biwako de la préfecture de Shiga ** sur la carte. Veuillez le remplir de manière appropriée. Cette fois, j'ai utilisé Seashore pour ** rgb (255,38,0) «dans la zone terrestre de la préfecture de Shiga ** et» rgb (0,0,255) «dans la région du lac Biwa dans la préfecture de Shiga **. Je l'ai colorée.
Puyoguramingu
Fondamentalement, cela revient à trouver le rapport de circonférence, à définir un point aléatoire dans la zone (carré ou préfecture de Shiga) et à déterminer s'il se trouve dans la zone spécifiée (cercle ou lac Biwa) Répéter.
Je ne pense pas que vous vous en souciez vraiment parce que c'est la même chose que ce que vous faites avec un cercle, mais ce dont vous devez vous inquiéter un peu, c'est que le dénominateur ** n'est pas le nombre d'époque. ** **
Divisez le nombre de ** entrés dans le lac Biwa ** par le nombre de ** entrés dans la préfecture de Shiga **.
Il est coloré pour identifier cela.
Comme vous le savez tous, les données de couleur d'une certaine coordonnée dans l'image peuvent être obtenues en accédant à l'adresse [y, x] du tableau qui apparaît par ʻimread` en utilisant OpenCV.
Je ne pense pas qu'il y ait autre chose à craindre.
Pour le moment, le nombre d'essais est fixé à 50 000.
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))
Le résultat de l'exécution étant redondant, collez-le uniquement à la toute fin. Puisqu'il s'agit de Rokubunnoichi, la valeur attendue est 1/6 = 0,1666666 ....
Attempt:50000
Shiga:21143
Biwako:3493
Overall Score:0.16520834318687036
C'est bon. C'est plutôt bon.
C'est pourquoi la préfecture de Shiga est un bon endroit. Veuillez venir une fois.
Un ensemble complet d'articles sur cet article est résumé sur GitHub. Correspondance s'il vous plaît voir.
https://github.com/YuseiIto/Biwako_Rate
Recommended Posts