Python
Pythonで実装してみた 元画像を読み込むところ。BGRに分割し、分割された領域ごとの平均輝度値を求める import cv2 import glob import numpy as np orgfile = "original02.png" orgimg = cv2.imread(orgfile) P=20 #分割する縦横のピクセル数 Nx = int(orgi…
可能な限り全てのツイートを検索して取得する。ただし、サーチバン等により検索で出て来ないツイートは取得できない。 import tweepy tweets = [] query = "物理の駅" for page in tweepy.Cursor(api.search_tweets,q=query, tweet_mode='extended',result_t…
try: raise "ERROR" except: try: import traceback traceback.print_exc() except:pass
https://docs.tweepy.org/en/latest/api.html#tweepy.API.search_tweets API.search_tweets(q, *, geocode, lang, locale, result_type, count, until, since_id, max_id, include_entities) Returns a collection of relevant Tweets matching a specified …
LISE++ファイルは、テキストファイルで作られているのでPythonで容易に読み取ることができる。 [] で囲まれた文字列をカテゴリ、その中で記述された = で区切られたKeyとValueを辞書型で詰めている。; 以降はコメントとして無視している。 def read_lise_lpp…
Jupyter lab等でグラフをインタラクティブに動かしたいとき、例えば特定の値を変えたグラフをセルを実行することなく次々と表示させたい時、ipywidgets.interactを使うと良い。 使用例 %matplotlib inline import matplotlib.pyplot as plt import ipywidget…
文字列→datetime, unixtime, 文字列 import datetime tstr = '2020/9/17 9:00:01' tdt = datetime.datetime.strptime(tstr, '%Y/%m/%d %H:%M:%S') d_week = {'Sun':'日','Mon':'月','Tue':'火','Wed':'水','Thu':'木','Fri':'金','Sat':'土'} print(type(tdt…
session を使い回すところがポイント import requests session = requests.Session() url = "https://test.test/login" data = { "email": "mail@mail.mail", "password": "password", } ret = session.post(url=url, params=data) assert(ret.status_code==…
日本語訳は下に (Japanese translation below) This problem was fixed in version v2.2.2 and later. Please download the latest version of the slackdump. There is an excelent software called slackdump that allows you to back up all your visible…
物理では、例えば、β と Energy [MeV/u] のように、同じ量だが異なる表記をすることがある。 これを1つのグラフに同時に表示したいとき、以下のようにする。 import numpy as np import matplotlib.pyplot as plt Zs = np.arange(1, 100, 1) data1 = 0.5+Zs*…
dumpsで桁数を制御することはできないが、loadsで制御することが可能。速度は遅くなるが、dumps→loads→dumps にすると良い。 import json print(json.dumps([23.67, 23.97, 23.87])) print(json.dumps([23.67000000000002, 23.96999999999999, 23.8700000000…
Jupyter Lab で以下のようなエラーが出た。tornado がエラーだとか。 Exception in callback <TaskWakeupMethWrapper object at 0x000002767EB13B80>(<Future finis...a24"\r\n\r\n'>) handle: <Handle <TaskWakeupMethWrapper object at 0x000002767EB13B80>(<Future finis...a24"\r\n\r\n'>)> Traceback (most recent call last): File "c:\user…</future></handle></future></taskwakeupmethwrapper>
いつからかSlackのチャンネル名に日本語が使えるようになったが、SlackからエクスポートしたZIPファイルのファイル名、ディレクトリ名はcp437 コードページ437 - Wikipedia でエンコードされているため、UTF-8にするには、cp437としてエンコードする必要があ…
旧ワークスペースの一部のチャンネルのみ、新ワークスペース(新WS)にCSVファイルで移行するためのPythonスクリプトを書いたので紹介する。 準備することは、 旧ワークスペースからエクスポートしたデータ (アーカイブしたZIPファイル) ワークスペースのチャ…
C:\Users\Masahiro\Downloads\slack>jupyter lab Traceback (most recent call last): File "c:\users\masahiro\appdata\local\programs\python\python37\lib\runpy.py", line 193, in _run_module_as_main "__main__", mod_spec) File "c:\users\masahiro\a…
この記事の Awkward のバージョンは1.X.Xです 例外処理はちゃんとしてないので注意せよ。 import numpy as np import awkward as ak import matplotlib.pyplot as plt # 配列作成 obj={} rs = np.random.RandomState(0) obj["memx"]=rs.rand(1000) obj["memy…
まず、各定数の確認する。 import numpy as np elementary_charge = 1.602176634*10**-19 # 電気素量 C c_light_m = 299792458 # 光速 m/s electron_mass_kg = 9.1093837015*10**-31 # 電子質量 kg permittivity_vacuum = 8.8541878128*10**-12 # 真空の誘電…
standard_weights = [1.008, 4.003,\ 6.941, 9.012, 10.81, 12.01, 14.01, 16.00, 19.00, 20.18,\ 22.99, 24.31, 26.98, 28.09, 30.97, 32.07, 35.45, 39.95,\ 39.10, 40.08, 44.96, 47.87, 50.94, 52.00, 54.94, 55.85, 58.93, 58.69, 63.55, 65.38, 69.72,…
WindowsでもC++を直接呼ぶ方法があるらしいが、めんどくさいので引数で情報を渡して、標準出力で結果を得ることにする。 import subprocess exe = 'puroguramu.exe' arg = "hikisu" #引数が1つの例 proc = subprocess.Popen([exe, arg], stdout=subprocess.P…
phst.hateblo.jp の続編 mzfit が使っている tensorflow が重すぎなのと、モジュール読み込みに時間がかかるのと、機能的に不足する点があったのと、マルチプラットフォームであるべきという思想に反するので、lmfitを使ってマルチガウシアンでのフィッティ…
名古屋大学の障害情報一覧にあるネットワーク障害の発生年月日、発生月をPythonで解析してみる。2022/05/14 時点でのHTML構造に対応しているので、ホームページの構造が更新されると使えなくなることをご承知おきください。 データをダウンロード import url…
qiita.com NumPyでは、発生する乱数をシードを使ってコントロールすることが可能で、グローバルシードを使う方法と、乱数状態を保持したクラスを使う方法がある。 グローバルシードを変更する np.random.seed(seed) と、乱数状態クラスを取得するための np.r…
この記事の Awkward のバージョンは1.X.Xです 単一原子の基底状態のイオンの質量をPythonで求める。データは nds.iaea.org から引っ張ってきた。APIの説明。 def get_ion_mass(z, a, q): import os gs_filename = "ground_states.csv" if not os.path.exists…
Pythonはインタプリタなので、文字列から動的なimportが可能である。 Pythonのコードの中で、NewModule.py を作って、それを読み込んでみよう。 NewModule.py に書くテキスト print(0) a = 10 + 0 Pythonコード。import importlib が重要。 module_name = "N…
2つめの配列を条件にして、他の配列も削除したり、他の配列の値も変えたりする。 後者では、変える値は1つめから3つめまで全て同じ値になってしまう。異なる値にする方法はあるのだろうか。 import numpy as np a1 = np.array([10, 20, -30, 99, 33]) a2 = n…
qiita.com lmfit.github.io pip3 install lmfit で lmfitをインストールする コード import numpy as np import lmfit from lmfit.lineshapes import gaussian2d, lorentzian import matplotlib.pyplot as plt npoints=10000 np.random.seed(0) x = np.rando…
phst.hateblo.jp の二次関数版 import numpy as np # データ生成 x = np.array([10.09659323, 10.40523609, 10.72069954, 11.03063238, 11.33665493, 11.6391395, 11.9381814, 12.23925842, 12.5349865, 12.82765475, 13.11407145, 13.40273189, 13.69487096…
分離度(Resolution, R)=2.5のヒストグラムがどういうのか想像できなかったので作ってみた。分離度はそれぞれの分布を正規分布としたとき、隣り合うピークの平均値の差 と、各ピークの標準偏差を使い、以下の式で定義している。 半値全幅 FWHM は なので、 と…
periodic_table = ["H", "He",\ "Li", "Be", "B", "C", "N", "O", "F", "Ne",\ "Na", "Mg", "Al", "Si", "P", "S", "Cl", "Ar",\ "K", "Ca", "Sc", "Ti", "V", "Cr", "Mn", "Fe", "Co", "Ni", "Cu", "Zn", "Ga", "Ge", "As", "Se", "Br", "Kr",\ "Rb", "Sr",…
参考にしたサイト qiita.com わざわざ書き直すまでもないが、一応メモとして残しておく。 import requests from requests.auth import HTTPDigestAuth import datetime from time import sleep username = "user" password = "password" url = 'http://192.1…
