誰よりも速くNumpy数値計算のコードを最適化しろ！

• Python 3.8.2
• numpy 1.18.2
• ipython 7.13.0 （計測用）

このページは単なる文句というか意見をさらしておくだけのものにつき予めご了承ください。タグもOpinionになっているしあくまで意見だということで。

特にドロー系の画像編集ソフトなどでは、指定された2次元領域に関係しているオブジェクトを素早く特定して処理を行うために空間インデックスといわれる類のデータ構造をしばしば用いている。その内の一つである領域四分木というデータ構造とアルゴリズムを見ていき、Pythonで素朴に実装してみる。具体的なチンパンコード付き。

地理院地図はウェブページから普通に閲覧は出来るが、ウェブマップサービスとしてタイル座標からそれぞれの画像にアクセスすることもできる。これを使ってPythonから画像を取得して1枚の画像にするプログラムを置いておく。

国土地理院が提供する標高タイルをダウンロードしてGeoTiff形式にする簡単なPythonプログラムをメモ用に置いておく。