2024-01-01から1年間の記事一覧

QGISでGoogle map satelliteの画像をシェープファイルで切り抜く

普通にラスタツールでやるとエラーが出る ラスタ計算ではなく、マスクして見えなくするアプローチを取る。 やり方は以下 gis.stackexchange.com マスクするレイヤーのシンボロジーで、「Inverted Polygon」を選択し、塗りつぶしを白色、枠線なしで指定する。…

QGISで目的の座標に対応するシェープファイルの情報を抽出する

ケニアのWardのシェープファイル https://data.humdata.org/dataset/administrative-wards-in-kenya-1450Google map上でLag, Lngを取得(この場合、-1.993590, 36.867553) QGISプラグインでLag Lon toolsをインストール Plugins>Lat Lon Tools>Zoom to Coor…

温度と相転移時の圧力pkの関係(ファンデルワールスの状態方程式より)

周囲温度をTを上げていった際の、気体→液体への相転移時の圧力Pk(臨界温度Tc以上では臨界圧力pcに一致)をみたい。 ファンデルワールスの状態方程式で、周囲温度Tを上げていく。 ファンデルワールスの状態方程式でP-V曲線を計算し、変曲点をMatlabで計算する…

ファンデルワールスの状態方程式と理想気体の状態方程式の違い

clear; close all; a = 3.59; b= 0.0427; R = 0.08206; T = 300; v0 = 0.05; vf = 2; v = linspace(v0,vf,1000);%m3 p = R*T./(v-b) - a./v.^2; p_ideal = R*T./v; semilogx(v,p,'b'); hold on; semilogx(v,p_ideal,'r--'); ylim([60 80]); legend('Van der …

ギブスの自由エネルギーと圧力pのグラフ

ファンデルワールスの状態方程式 上記をP-V線図にプロットすると以下になる。 ギブスの自由エネルギーgA-gMは体積Vを用いて以下のように表される。 いま、ギブスの自由エネルギーの最下点(極小点)gMのエネルギーを仮にゼロとすると(任意の定数なのでゼロ…

ファンデルワールスの状態方程式でp-v関係を求める

clear; close all; a = 3.59; b= 0.0427; R = 0.08206; T = 300; v0 = 0.05; vf = 2; v = linspace(v0,vf,1000);%m3 p = R*T./(v-b) - a./v.^2; semilogx(v,p); ylim([60 80]); もういっちょ clear; close all; %Parameters for Van der Waals equation (exa…

Pressure transmitterとPressure transducerの違い

Pressure transmitter:圧力を電流に変換する Pressure transducer(圧力変換器、圧力センサ):圧力を電圧(and 電流)に変換する

MATLABの色設定

QGISの起動時のプラグイン読み込みのプロファイリング

gis.stackexchange.comGo to View --> Panels --> Debugging/Development Tools Select Project Load from the category list

GDALのインストール(Intel Mac, M1 Mac)

環境はIntel Mac, macOS 15.1 Sequoia. brew install gdalでGDAL-3.10がインストールされた。 pip install gdalでgdal-3.10のpipをインストールしようとするとビルドに失敗する。 ソースからのインストールでなんとかできたが、いくつか不可解な点は残った。…

Ansys用の計算サーバ選定

www.ansys.com

urdf-vizのインストール

github.comデフォルトのUbuntu20のrustcだと動かなかったので、一度アンインストールしたあと、curlで再度インストールした。 phoenixnap.com curl https://sh.rustup.rs -sSf | sh 便利である。View post on imgur.comimgur.com

KUKAで与えられた座標にロボットを動かす

まずはシミュレータ上で動かす 1)RViz上で指定した位置に動かす 2)座標を指定して動かす(1点) 3)座標を指定して動かす(スクリプトを用いて連続した値を送り続ける:G-codeを模擬)まずはRvizから。 逆運動学のエンジンにはtrack_ikを用いる。 まず…

STKでSentinel-2の通過時刻を計算する

sorabatake.jp

ソファーの構造検討

ベースライン 設計:35mm角 30cm感覚で配置 荷重:100kgf 応力 変位 12.7mm 設計目標 100kgfの荷重で最大変位が1mm以下 50mm角に変更 応力 最大変位 3.6mm 75mm角に変更 応力 最大変位 0.92mm 厚みは50mmに決定 配置間隔を20mmと25mmにしたときの比較 結局は…

Ansysでのエラー(At least one body has been found to have only 1 element in at least 2 directions along with reduced integration)

At least one body has been found to have only 1 element in at least 2 directions along with reduced integration. This situation can lead to invalid results or solver pivot errors. Consider changing to full integration element control or me…

Kerasで多出力の回帰モデル

machinelearningmastery.com # use mlp for prediction on multi-output regression from numpy import asarray from sklearn.datasets import make_regression from keras.models import Sequential from keras.layers import Dense # get the dataset def …

IMEの文字種切り替えの吹き出しを消す

scrapbox.io 以下を実行 $ defaults write kCFPreferencesAnyApplication TSMLanguageIndicatorEnabled 0

XacroをUrdfに変換するシェルスクリプト

/usr/local/bin/xacro2urdf #!/bin/bash if [ "$#" -eq 1 ]; then filename="$1" extension="${filename##*.}" if [[ ${extension} == "xacro" ]]; then robot=${filename%.*} echo rosrun xacro xacro -o ${robot}.urdf ${robot}.xacro rosrun xacro xacro …

Ubuntu 20.04 LTS (arm64)のダウンロード

M1 Mac用のイメージ(arm64)がほしい Ubuntu Focal (Daily Build)のページに行く ftpmirror.your.org

OctaveでGAを用いた最適化

ローゼンブロック関数の最適化 pkg load ga; function y = rosenbrock(x) % Rosenbrock function, commonly used for testing optimization algorithms y = 100 * (x(2) - x(1)^2)^2 + (1 - x(1))^2; end nvars = 2; lb = [-5, -5]; % Lower bounds ub = [5,…

CadQueryを用いてPythonでパラメトリック3Dモデリング

import cadquery as cq import numpy as np # NACA 4-digit wing profile generator (e.g., NACA 2412) def naca_4_digit(m, p, t, chord_length, num_points): x = np.linspace(0, 1, num_points) * chord_length yt = 5 * t * (0.2969 * np.sqrt(x/chord_l…

Macで動画(ファイル)の日付を変更する

Created dateを変更するとき あるファイルから日付(Created date)を抽出 以下のフォーマットを取得(%m/%d/%Y %H:%M:%S) $ stat -f %Sm -t "%m/%d/%Y %H:%M:%S" IMG_1984.JPG(ファイル名) 03/30/2012 00:15:32SetFIleでCreate dateを変更 上記のフォーマ…

Octaveでスタートアップスクリプトを設定する

以下のbig.mスクリプトを常に参照できるようにする。 big.m set(findall(gcf,'-property','FontSize'),'FontSize',20); set(findall(gcf,'-property','MarkerSize'),'MarkerSize',14); set(findall(gcf,'Type','Line'),'LineWidth',2); grid; 1. big.mを任意…

m1 macにOctaveをインストールする

Homebrewだとエラーが出たので、以下で配布されているapp bundleを利用する。 octave-app.org discussions.apple.com ダウンロードしたOctave.appを起動すると壊れていると言われるので、以下を実行する必要がある。 $ xattr -c /Applications/Octave-9.2.app

Ubuntu20にOpenRaveのインストール

robots.uc3m.es sudo apt install rapidjson-dev sudo apt install git # probably already installed sudo apt install libboost-filesystem-dev libboost-system-dev libboost-thread-dev libboost-iostreams-dev sudo apt install libxml2-dev libode-dev…

電子観望

Samさんの動画 NEWTONY Player One Ceres-C SkyWather AZ-GTi www.youtube.com その2 Skywatcher EVOGUIDE 50ED SIGHTRON SV305-SJ SkyWather AZ-GTi www.youtube.comSamさんのブログ hoshizolove.blog.jpevScope equinox 36万は高すぎ www.youtube.com

Zoteroの設定

Linked Attachment Base Directoryに~/Dropbox/Zoteroを指定 Data Directoryに~/Zoteroを指定 ZotMoovのDirectory to Move Files toに~/Dropbox/Zoteroを指定

ガウス過程

qiita.comzenn.devqiita.comwww.hello-statisticians.com

多目的最適化の問題

tech-oh.com 上記の例を解こうと思ったが、RBFによる応答曲面で計算した出力yと、制御因子Xのベクトルを掛け算するところで詰まった。進め方としては以下のようになっている。 import numpy as np from pymoo.core.problem import Problem import math POP_S…