https://www.daiphys.com/portal/ 2026-05-01T09:45:12+00:00 https://www.daiphys.com/portal/ https://www.daiphys.com/portal/_media/wiki/dokuwiki.svg text/html 2022-12-14T15:50:09+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/fits.html?rev=1671033009&do=diff FITS Sample import pyfits hdulist = pyfits.open('sample.fits') hdulist.info() hdulist[0].header['name'] scidata = hdulist[1].data scidata = hdulist['TABLE'].data hdulist.close() text/html 2024-08-08T03:22:33+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/gdown.html?rev=1723087353&do=diff Gdown Google Drive Public File Downloader * An alternative way to curl/wget for working around the security notice on Google Drive when downloading a large file References Gdown <https://github.com/wkentaro/gdown> text/html 2023-03-10T05:55:57+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/index.html?rev=1678427757&do=diff Python * FITS * Gdown * Jupyter * List/Array * Logging * Matplotlib * Minuit * Misc * Mmap * Numpy * OpenCV * OpenPyXL * Pandas * Pickle * PIL/Pillow * Pip * Polars * Pytest * Quandl * ROOT * Scapy * Scipy * Seaborn * Socket * Streamlit * Threading * Tifffile * Uv * Venv * Yt text/html 2025-07-14T05:36:15+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/jupyter.html?rev=1752471375&do=diff Jupyter Tips * Jupyter Notebook上で通常の右クリックによるメニュー(名前を付けて画像を保存など)を表示したい場合はShift+右クリック Install <https://jupyter.org/install> ## Jupyter Lab pip install jupyterlab jupyter-lab jupyter-lab --port=8888 --ip=0.0.0.0 jupyter-lab --port=8888 --ip='*' --NotebookApp.token='' --NotebookApp.password='' --no-browser text/html 2023-11-02T17:21:54+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/list-array.html?rev=1698945714&do=diff List/Array Sample Map def f(value) : return value * 2a = [1, 2, 3]list(map(f, a)) list(map(lambda x: (x * 2), a)) [f(x) for x in a] [(x * 2) for x in a]## [2, 4, 6] Map for Lists a = [1, 2, 3] b = [4, 5, 6]list(map(lambda x, y: (x * y), a, b))## [4, 10, 18] text/html 2023-06-29T14:46:50+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/logging.html?rev=1688050010&do=diff Logging Sample import logging logging.basicConfig(level=logging.DEBUG)import logging import traceback logging.basicConfig(filename=logfile, level=logging.DEBUG, format="[%(asctime)s] pid=%(process)d:%(levelname)s:%(name)s:%(message)s") logging.getLogger().addHandler(logging.StreamHandler()) logging.error(traceback.format_exc()) text/html 2024-02-22T02:32:14+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/matplotlib.html?rev=1708569134&do=diff Matplotlib rcParams plt.rcParams['figure.figsize'] = [12, 8] #plt.rcParams.update(plt.rcParamsDefault) Plot plt.figure(figsize = (12, 8)) plt.plot(x, y, 'o') plt.errorbar(x, y, xerr=(xl, xh), fmt='o', color='black', ecolor='red', capsize=5) plt.xlim(0, 1) plt.xticks([]) plt.yticks([]) plt.xticks(np.arange(0, 1, 0.1)) plt.xlabel('X') plt.grid() plt.legend() plt.savefig('figure.png', bbox_inches='tight') plt.show() text/html 2023-10-13T03:55:36+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/minuit.html?rev=1697169336&do=diff Minuit Tutorials Google Colaboratory <https://colab.research.google.com/drive/1jufurZXLxzkC1FzyOXpq6I7_-oDQZvCt> Sample import numpy as np from iminuit import Minuitdx = np.arange(10) dy = ((50.0 * dx) + 100.0) * (1.0 + (0.05 * np.random.randn(len(dx)))) dy_error = np.sqrt(dy) model = lambda x, a, b : ((a * x) + b) chisq = lambda a, b : sum(((dy - model(dx, a, b)) ** 2.00) / (dy_error ** 2.00)) text/html 2025-10-22T17:10:46+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/misc.html?rev=1761153046&do=diff Misc Built-in Functions <https://docs.python.org/3/library/functions.html> Dir <https://docs.python.org/3/library/functions.html#dir> import sysdir()## ['__annotations__', '__builtins__', '__doc__', '__loader__', '__name__', '__package__', '__spec__', 'sys']dir(sys)## ['__breakpointhook__', '__displayhook__', '__doc__', '__excepthook__', '__interactivehook__', '__loader__', '__name__', '__package__', '__spec__', '__stderr__', '__stdin__', '__stdout__', '__unraisablehook__', '_base_executab… text/html 2022-12-02T17:14:41+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/mmap.html?rev=1670001281&do=diff Mmap Server import queue import mmap mapfile = '.malloc.mmap' with open(mapfile, 'wb') as f : f.write(b'\x00') q = queue.Queue() n = 0 with open(mapfile, 'r+b') as f : mm = mmap.mmap(f.fileno(), 0) while True : mm.seek(0) flag = mm.readline() if flag == b'\x01' : n += 1 q.put(n) mm[:] = b'\x00' else : if not q.empty() : print(q.get()) text/html 2023-05-29T03:21:47+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/numpy.html?rev=1685330507&do=diff Numpy Runtime Error np.seterr(divide='ignore') np.seterr(divide='raise') np.seterr(invalid='ignore') np.seterr(invalid='raise') Print Option np.set_printoptions(suppress=True) References text/html 2025-08-25T03:56:16+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/opencv.html?rev=1756094176&do=diff OpenCV Install pip install opencv-pythonpython import cv2 Troubleshooting - Missing libGL ## Traceback (most recent call last): ## ImportError: libGL.so.1: cannot open shared object file: No such file or directorysudo apt install libgl1 Sample Reading 8-bit/16-bit/32-bit TIFF text/html 2023-05-11T05:18:57+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/openpyxl.html?rev=1683782337&do=diff OpenPyXL Sample import openpyxl## Workbook wb = openpyxl.load_workbook('sample.xlsx') wb = openpyxl.Workbook() wb.save('sample.xlsx') wb.sheetnames## Worksheet ws = wb['Sheet'] ws = wb.worksheets[0] ws.title## Cell A1 c1 = ws['A1'] c1 = ws.cell(1, 1) c1.coordinate c1.row c1.column c1.value = 100 ## Cell A2 c2 = ws['A2'] c2 = ws.cell(2, 1) c2.value = "=SUM(A1:A1)" text/html 2025-08-23T10:30:01+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/pandas.html?rev=1755945001&do=diff Pandas Cheat Sheet <https://github.com/pandas-dev/pandas/tree/main/doc/cheatsheet> Display pd.get_option('display.max_rows') pd.get_option('display.max_columns') pd.set_option('display.max_rows', None) pd.set_option('display.max_columns', None) I/O df = pd.read_csv('sample.csv') Misc df.replace({'n':0, 'y':1}) df[df['a'] == 'y'].sort_values(by='b') df[(df['a'] == 'y') & (df['b'].notnull())].sort_values(by='c') df.loc[df['a'] == 0, ['b']] = 1 df['a'].notnull() * 1 text/html 2023-02-20T07:02:20+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/pickle.html?rev=1676876540&do=diff Pickle I/O import pickle data = 1 pickle.dump(data, open('sample.pkl', 'wb')) data = pickle.load(open('sample.pkl', 'rb'))with open('sample.pkl', 'wb') as f : pickle.dump(data, f)with open('sample.pkl', 'rb') as f : data = pickle.load(f) text/html 2025-08-25T09:21:30+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/pil.html?rev=1756113690&do=diff PIL/Pillow Install pip install pillowpython from PIL import Image Sample Loading Image img = Image.open("input.tif") img_array = np.array(img) Troubleshooting - Decompression Bomb Error <https://pillow.readthedocs.io/en/stable/reference/Image.html> <https://stackoverflow.com/questions/56174099/how-to-load-images-larger-than-max-image-pixels-with-pil> ## PIL.Image.DecompressionBombError: Image size (400000000 pixels) exceeds limit of 178956970 pixels, could be decompression bomb DOS attack. … text/html 2024-07-23T13:54:01+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/pip.html?rev=1721742841&do=diff Pip Requirements pip freeze > requirements.txt pip install -r requirements.txt GitHub pip install git+https://github.com/numpy/numpy.git Upgrade pip install --upgrade pip setuptools Troubleshooting - No Matching Distribution Upgrade setuptools to fix the following error: text/html 2024-02-28T02:22:18+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/polars.html?rev=1709086938&do=diff Polars I/O df = pl.read_csv('sample.csv') References Polars <https://pola.rs/> <https://qiita.com/_jinta/items/fac13f09e8e8a5769b79> <https://zenn.dev/bee2/articles/e8623a603752ff> text/html 2023-07-20T04:53:59+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/pytest.html?rev=1689828839&do=diff Pytest Test Results * PASSED * FAILED * SKIPPED * XFAIL: Expected Failure * XPASS: Unexpected Passing * ERROR References <https://zenn.dev/imaki/scraps/8e70706f595cdb> text/html 2022-12-15T04:01:08+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/quandl.html?rev=1671076868&do=diff Quandl Sample import quandl quandl.ApiConfig.api_key = "abcdefg..." data = quandl.get("NIKKEI/INDEX") data = quandl.get("BOE/XUDLJYD") data = quandl.get("CUR/JPY") text/html 2022-12-16T02:15:53+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/root.html?rev=1671156953&do=diff ROOT Preparation import ROOT st = tgraph.FindObject('stats') ROOT.gROOT.SetBatch(1) Read Data data = np.genfromtxt("input.dat", unpack=True) xval = np.array(data[0]) yval = np.array(data[1]) xerr = np.array(data[2]) yerr = np.array(data[3]) xerl = np.array(xerr) xerh = np.array(xerr) yerl = np.array(yerr) yerh = np.array(yerr) xmin = np.min(xval) xmax = np.max(xval) ymin = np.min(yval) ymax = np.max(yval) text/html 2025-09-22T03:15:16+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/scapy.html?rev=1758510916&do=diff Scapy Sample import pandas as pd from scapy.all import rdpcap from scapy.all import IP packets = rdpcap("input.pcap") records = [] for i, pkt in enumerate(packets): if IP in pkt: records.append({"idx": i, "time": pkt[IP].time, "src": pkt[IP].src, "dst": pkt[IP].dst, "id": pkt[IP].id, "len": pkt[IP].len}) else: records.append({"idx": i, "time": pkt[IP].time, "src": None, "dst": None, "id": None, "len": None}) df = pd.DataFrame(records) text/html 2022-12-14T15:45:43+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/scipy.html?rev=1671032743&do=diff Scipy Newton-Raphson ## derive x at f(x) = 5 def g(x) : value = (f(x) - 5.00) return(value) scipy.optimize.newton(g, 0.00) text/html 2023-10-22T15:30:02+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/seaborn.html?rev=1697988602&do=diff Seaborn References Seaborn <https://seaborn.pydata.org/> <https://qiita.com/4m1t0/items/76b0033edb545a78cef5> text/html 2022-12-02T18:12:54+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/socket.html?rev=1670004774&do=diff Socket TCP Server import socket import time def tcpserver(host, port, nlis=1) : sock = socket.socket(socket.AF_INET, socket.SOCK_STREAM) sock.setsockopt(socket.SOL_SOCKET, socket.SO_REUSEADDR, 1) sock.bind((host, port)) sock.listen(nlis) print('## waiting for connection at %s:%s') % (str(host), str(port)) return(sock) sock = tcpserver('localhost', 54321) conn, addr = sock.accept() print('## connected by %s:%s') % (str(addr[0]), str(addr[1])) n = 0 while True : n +… text/html 2025-07-02T03:17:37+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/streamlit.html?rev=1751426257&do=diff Streamlit References Streamlit <https://streamlit.io/> text/html 2024-07-10T01:28:14+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/threading.html?rev=1720574894&do=diff Threading References <https://qiita.com/castaneai/items/9cc33817419896667f34> text/html 2025-08-25T09:07:02+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/tifffile.html?rev=1756112822&do=diff Tifffile Install pip install tifffilepython import tifffile Sample Loading Image img_array = tifffile.imread("input.tif") Extracting Alpha Channel img_array = tifffile.imread("input.tif") img_array_alpha = img_array[..., 1] References <https://github.com/cgohlke/tifffile> <https://pypi.org/project/tifffile/> text/html 2026-04-21T09:08:49+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/uv.html?rev=1776762529&do=diff Uv Sample References <https://docs.astral.sh/uv/> text/html 2023-11-28T06:38:00+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/venv.html?rev=1701153480&do=diff Venv Sample sudo apt install python3-pip sudo apt install python3-venv python3 -m venv venv source venv/bin/activate pip install -r requirements.txt deactivate Venv in Jupyter See Jupyter for more details on how to use virtual environments in Jupyter Notebook text/html 2022-12-16T01:10:52+00:00 Anonymous (anonymous@undisclosed.example.com) https://www.daiphys.com/portal/python/yt.html?rev=1671153052&do=diff Yt Sample import yt.mods filename = "sample.hdf5" pf = yt.mods.load(filename) pc = yt.mods.PlotCollection(pf, [0.0,0.0,0.0]) pc.add_slice("Density", "z", use_colorbar=True, periodic=False)pc.set_width(width=5.00E+12, unit='cm') pc.set_xlim(-1.0E12, 1.0E12) pc.set_ylim(0.0, 1.0E12) pc.set_zlim(1.0E+00, 1.0E+30, nticks=5) pc.autoscale() pc.set_cmap(hot)