クロフトンの公式 CTなどのボリュームデータにおいて、何らかの方法でセグメンテーションを行ったあと、そのラベルデータの表面積を計算したいときがあります。このとき、クロフトンの公式というものを使うと正確に計算できます。これ…

方針 簡単のために単位球を考えます。乱数tで緯度を決め、乱数uで経度を決めます。ただし、緯度ごとに、緯線の長さが異なるため、球面上に均一に分布させるために、この緯線の長さに比例した頻度でtを生成する必要があります。 なお…

Pythonでの生成方法を書いていきますが、根底にある考え方は他の言語でも共通です。 手始めに、一様乱数から 確立密度関数で表現すると、[0, 1)において、P(x) = 1。 # [0, 1)の乱数を生成 import…

Pythonでは3Dポリゴンを扱うことができます。STLというフォーマットはポリゴンのXYZ座標を列挙するだけなので、ライブラリを使わずとも扱えますが、numpy-stlというライブラリを使うともっと簡単です。インストー…

プログラミングを始めようとMacを手にしたあなたのために、プログラミングを始めるために必要な環境構築の方法をここにまとめておきます。 Pythonのインストール こちらのページにあるとおりですが、まずはPythonをイン…

概要 今回は、すでにchainerがインストールされているDockerイメージに少しだけ手を加えて、OpenCVとJupyter Notebookが使えるようなイメージを作成します。GPUをDocker上で使うためにはd…

PythonでOpenCVを使うには import cv2 import numpy as np PythonとOpenCVのインストールの方法はこちらです。インストールが終わったなら早速始めてみましょう。OpenCVを…

マルチページTIFFをPythonで読み込むためにはPIL(Python Image Library)を用います。次のコードではマルチページTIFFの30枚目から100枚目までを読み込み、512*512ピクセルにリサイズ…

Pythonで連番画像を読み込む方法の紹介です。画像の読み込み自体はOpenCVのimread関数で行いますが、この関数にはファイル名の文字列を渡す必要があります。今回はPythonで連番ファイルの文字列を生成する方法の…

今回はPythonで厳密な重心を求める手法を紹介します。まずは重心の定義から見ていきましょう。C++版はこちら。 幾何学的には、ある図形の、そのまわりでの一次モーメントが 0 であるような点のこと。図形 D （およびその…