Pythonのnumpyで四元数してみた

Pythonで四元数を使う方法について調べていたところ、numpy-quaternion · PyPIに辿り着きました。このページでは、numpy-quaternionのインストール方法について紹介していきます。

numpy-quaternionをインストールする

①powershellを管理者権限で起動して「pip install numpy-quaternion」と入力して実行します。

②インストールできたら、以下のコードを保存してPythonで実行してみます。

import numpy as np
import quaternion

q1 = np.quaternion(1,2,3,4)
q2 = np.quaternion(5,6,7,8)
print("q1 * q2 ", q1 * q2)

③使ってみると、こんな警告が表示されました。どうやらnumbaというものをインストールしたほうがいいよってことだそうです。

numbaをインストールする

①管理者権限でpowershellを起動して「pip install numba」と入力して実行します。

②先ほどのソースコードの「import quaternion」の行の下に「import numba」と入力して実行してみます。すると、今度は警告が出なくなりました。

numpy-quaternionの使い方

下記のページのUsageには、numpy-quaternionの使い方が幾つかのサンプルで紹介されていますので、ご覧になってみてください。

github.com

ありがとうございました

QGIS3 GeoPackageの使い方

QGIS3で「GeoPackage」というファイル形式を使ってみたので、このページで紹介します。GeoPackageについて、まだ詳しいことが分かっていません。正確な情報はこちらのページに書いてあるのですが、正直何を言っているのか...。

ファイルを保存する

①今回はシェープファイル、Tabファイル、アスキーグリッドをQGISで読み込んで、それらをGeoPackageで保存してみます。

ブラウザの「GeoPackage」を右クリックして、「Create Database」を選択します。

③保存先を指定して、OKボタンを押します。

④GeoPackageの下に保存先に指定したファイルが登録されました。

⑤GeoPackageに保存したいファイルをQGISに読み込みます。

⑥読み込んだレイヤーを、保存先のGeoPackageにドラッグします。成功すると画像のようなダイアログが表示されます。

 

⑦保存先のGeoPackageの下に、ドラッグしたレイヤーが表示されました。

ファイルを読み込む

①本当にまとめられたか確認をしてみます。ブラウザのGeoPackageの、保存先に指定したファイル名を右クリックして「接続を削除」を選択します。一度QGISを終了します。

②GeoPackageの下に何も表示されていないことを確認し、ブラウザの「GeoPackage」を右クリックして、「New Connection」を選択します。

③ファイル選択ダイアログが表示されるので、先ほど保存したGeoPackageを指定します。するとGeoPackageの下に指定したファイルが表示されます。

④GeoPackageの下の、ファイル名のさらに下に、GeoPackageにまとめたレイヤーが表示されています。これらをダブルクリックすると、下の赤丸のレイヤーウィンドウにレイヤーが表示されます。

複数のファイルだったのが、単一のファイルで済むようになりましたが、もっと良い点や使い方が見つかりましたら、追記していきたいと思います。

ありがとうございました。

DotSpatialをインストールしてみました。

『DotSpatial』というライブラリを最近知りました。GISを作るためのライブラリだそうです。インストールまでは出来るようになったので、その手順を紹介いたします。コンパイルするためには、VisualStudioが必要です。ここではVisualStudio2017を使ってコンパイルをしています。

ソースをダウンロード

①https://github.com/DotSpatial/DotSpatialからソースをダウンロードします。

②ダウンロードしたファイルを解凍します。解凍して出来たフォルダ内の『Source』フォルダに移動して下さい。

コンパイル

③『Compile.cmd』を編集します。私の環境では『%ProgramFiles(x86)%\MSBuild\14.0\Bin\MSBuild.exe』を
『%ProgramFiles(x86)%\Microsoft Visual Studio\2017\Community\MSBuild\15.0\Bin\amd64\MSBuild.exe』に変更しました。

 ④編集したCompile.cmdを実行します。

⑤コンパイルが終わると、binフォルダが出来ているので、フォルダ内に移動します。

デモを動かしてみる

⑥DemoMap.exeを起動してみます。

⑦ウィンドウが表示されました。何か適当なshpファイルを表示させます。赤丸部分のボタンを押すと、ファイル選択ダイアログが表示されるので、shpファイルを選択して下さい。

⑧選択したshpファイルが表示されました。

プログラム作成

⑨自分のプログラムにDotSpatialを使ってみます。VisualStudioを起動して、C#のプロジェクトを新規作成します。

⑩『参照の追加』でDotSpatialのDLLを読み込みます。

⑪先ほどコンパイルしたdotspatialのbinフォルダ内のDLLを選択します。選択したDLLにチェックがされているか確認して、OKボタンを押します。

⑫ツールボックスにDotSpatialのコントロールを追加します。ツールボックスを右クリックして、メニューの『タブの追加』を選択して下さい。

⑬新規のタブを『DotSpatial』としました。

⑭DotSpatialを右クリックして、メニューの『アイテムの選択』をクリックします。

⑮表示されるダイアログの、参照ボタンをクリックします。

⑯コンパイルして出来たbinフォルダ内のDLLを選択します。

⑰このようなエラーが出ますが、気にせずOKボタンをクリックして下さい。

⑱読み込まれたらOKボタンをクリックします。

⑲DotSpatialの下に、コントロールが追加されました。『Map』を選択して、フォームの適当な場所に貼り付けて下さい。他にCommandButtonもMapに重ならないように貼り付けて下さい。貼り付けたボタンをダブルクリックして、次のコードを入力して下さい。

----------------------------------------------------------------------------------------------

map1.ClearLayers();
map1.AddLayer(@"D:\適当なシェープファイル.shp"); //←適当なshpファイルを指定する
IMapFeatureLayer imp = map1.Layers.Last() as IMapFeatureLayer;
imp.SelectionSymbolizer = new PolygonSymbolizer();
PolygonSymbolizer poly = imp.Symbolizer as PolygonSymbolizer;
poly.SetFillColor(System.Drawing.Color.FromArgb(0)); //ポリゴン塗りつぶし色無し

----------------------------------------------------------------------------------------

ソースコードの上の方に、次のコードも入力して下さい。

----------------------------------------------------------------------------------------

using DotSpatial.Controls;
using DotSpatial.Symbology;

---------------------------------------------------------------------------------------

⑳プログラムをコンパイルして実行します。ボタンを押すと、指定したshpファイルが表示されれば、インストールは上手くいっています。

お疲れさまでした。