|
setup diary | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||
これまで，pdfを分割したり，結合したりする必要があるときには，imagemagickを使ったり，img2pdfやpdfuniteを使ったりしてきた．一旦画像に変換するかしないかなどの違いはあるが，これらを使う場合には，以下のようにする．
convert *.jpg out.pdf img2pdf -o out.pdf *.jpg pdfunite *.pdf output.pdf
gsを使って，似たようなことができるということなので，やってみた．まず，pdfから一部を取り出すには，次のようにする．
gs -dNOPAUSE -dQUIET -dBATCH -sOutputFile=out.pdf -dFirstPage=1 -dLastPage=200 -sDEVICE=pdfwrite in.pdf
次に，pdfを結合するには，以下のようにする．
gs -q -dNOPAUSE -dBATCH -sDEVICE=pdfwrite -sOutputFile=out.pdf in1.pdf in2.pdf in3.pdfpdfのまま処理するので，ファイルサイズはあまり変わらないかとも思ったら，いろいろとやっていたら二倍ぐらいのサイズになってしまった．サイズが変わらない変換方法があれば良いのに．
\documentclass[a4j,dvipdfmx]{jarticle}
\usepackage{color}
\usepackage{epic,eepic}
\usepackage{tikz}
まず，picture環境では，以下のように単位を指定してから，putコマンドで，座標と配置するものを定義すると，図が描ける．
\unitlength 5mm
\begin{picture}(10,10)
\put(5,5){\circle{10}}
\put(5,5){\makebox(0,0){$O$}}
\put(0,0){\color{red} \blacken\path(0,0)(2,2)(2,0)(0,0)}
\put(7,5){\color{cyan} \circle*{2}}
\put(7,6){\color{magenta} \shade\circle{2}}
\end{picture}
一方，tikzでは，tikzpicture環境を使い，デフォルトでは10mm単位なので，必要に応じてscaleを変更して，コマンド，座標，配置するものと指定して，最後にセミコロンをつけて，図を描く．
\begin{tikzpicture}[scale=0.5]
\draw (5,5) circle [radius=5];
\draw (5,5) node {$O$};
\filldraw[red] (0,0)--(2,2)--(2,0)--cycle;
\shade[left color=yellow, right color=cyan] (7,5) circle [radius=1];
\filldraw[draw=magenta,fill=magenta,fill opacity=0.5] (7,6) circle [radius=1];
\end{tikzpicture}
ここで示した例は，そのほんの一部であるが，グラデーションで色を塗ったり，半透明にしたりと，picture環境ではできないことが，tikzではできる．
それぞれのファイルに書かれている数値の最大値をファイル毎に表示したくて，やり方を考えてみた．awkを使ってやろうとしたけど，絶対値や最大値を求めるのが面倒で，ファイル毎の処理も大変だったので，結局rubyを使ってやってしまった．例えば，こんな感じで一行にした．
ruby -e "puts Dir.entries('.').grep(/harm/).sort.map{|fn|open(fn){|f|f.readlines}[1..-1].map{|l|l.split[3].to_f.abs}.max}"
ファイルの処理を書くのは，rubyでは長くなってしまうし，列を取り出す処理も少し面倒なので，別の言語の方が良いかも知れないけど．
texで表を作っていて，左右とつなげたいときには，
\multicolumn{2}{c}{contents}
などとすれば良いが，上下をつなげたいときは，multirowというパッケージを使わないといけないようだ．一方，htmlだと簡単に上下左右をつなげることができる．左右と上下をつなげるには，それぞれ以下のようにする．
<td colspan=2> <td rowspan=2>
これを使って，光速(299792458m/s)の語呂合わせを考える表を作ってみた．
| 2 | 肉球,憎く,磁極,布告,復刻 | 肉,服,不休,普及 | ||
| 9 | ||||
| 9 | 食うな,困難,救難,急な | 苦無く | ||
| 7 | 無く,泣く,急く | 泣くに,石窟 | ||
| 9 | 工夫,国,苦痛,苦に,給餌,給付,クジ | |||
| 2 | 強いツヤ,によいつえー | 不死,不信,ニシン,自信 | ||
| 4 | 死後は | よ,し | ||
| 5 | 御飯,小屋,碁盤,合板,ゴーヤ,ファイヤー | |||
| 8 | ||||
これらの語句をつなぎ合わせると，いろいろなものができるが，その中でましなのは，これかな．
肉食うな ，工夫し御飯(脳梗塞)
他の定数もこのような表を作って考えてみようかな．
Taylor展開の係数を計算するのが面倒だったので，maximaを使ってやってみた． どの変数について何次まで展開するかを指定すると，計算してくれる． 具体的な例はこんな感じ．
f(x,r):=1/sqrt(1-2*x*r+r^2); taylor(f(x,r),r,0,3);
公式をつかって，微分と比較するなら，このようにして確かめられる．
diff(f(x,r),r,3)/3!;
PCと装置とのインターフェースは，時代とともに徐々に変わって来ているが，長く使われているものとして，RS232Cがある．その通信方法にはいろいろな流儀があるので，特定の装置と通信するのが意外に難しい場合がある．
昔のwindowsからRS232Cを介して制御していた古い装置があるのだが，linux上のwineから動かして使っていた．しかし，wineのversionによっては動かなかったり，クセのある装置である．以前のPCは2017年6月にセットアップしたようだが，そのときも苦労していたようだ．そのシステムもずいぶん古くなってしまったので，別のPCに新しくlinuxを入れて，それから制御しようとしたが，うまくいかない．
今回はオシロを使って，どのような信号が流れているのかを調べてみた．制御できるPCでは，通信の最初にDTR線の電圧がパルス的に変化していたが，うまくいかないPCでは，DTRは最初だけステップ関数的に変化して，二度目以降は動いていないことが分かった．調べてみると，linux kernelは，DTR/DSRフロー制御をサポートしていないという記述を見つけた．この装置を最新のlinux上のwineから動かすのは難しいのかも知れない．
maximaで%pi/5を含む三角関数をいじっていて，少し苦労したので，メモをしておく．表示を一行でできるようにdisplay2d:false;とした後で， exp(%i*%pi/3);とすると，(sqrt(3)*%i)/2+1/2という答えが出てくる． ところが，exp(%i*%pi/5);では%e^((%i*%pi)/5)と，ほぼそのままの結果しか得られない．そこで，
trigrat(exp(%i*%pi/5)+exp(-%i*%pi/5));
とすると，
%i*(sin((2*%pi)/5)-sin((3*%pi)/5))-cos((3*%pi)/5)+cos((2*%pi)/5)+1
と，三角関数で表してはくれるが，複雑な式になる．少なくとも素直に2*cos(%pi/5)として欲しいのに，なかなかそうなってくれない． これを改善するためには，load(ntrig);と，パッケージを入れる必要がある． すると，(sqrt(5)+1)/2と，その値を計算してくれて，黄金比が出てくる． maximaで三角関数などを思ったように変形するのは，意外に難しい．
装置についたカメラをwebから監視したいと思っていたのだが，linuxからビデオキャプチャできる良いものが無かったので，なかなか実現できなかったが，ようやくできるようにすることができた．使ったのはPCA-DAV2というUSBビデオキャプチャで，入手が困難だったが売っているところを見つけてなんとか買うことができた．それをdebian bullseyeからmotionを使ってwebに流すまでの設定についてメモをしておきたい．
まず，そのままだと認識されないので，/etc/modprobe.d/em28xx.confの中にoptions em28xx card=29と記述して，再起動すると認識するようになる．カード番号については，13,40,68という情報もあったが，自分の環境では29が正しいようだ．
次に，そのままだとPALが標準設定になっているので，これをNTSCに変更する必要がある．v4l-utilsをインストールして，v4l2-ctl --set-standard=0とすると，NTSCになる．すると，cheeseなどで画像が見れるようになる．
そして，web配信はmotionを使った．これをインストールしたら，/etc/motion/motion.confを変更して必要な設定を行う．daemonをoffにして，NTSCにするためにnormを1にするのが重要で，あとは自分の環境に合わせて設定する．今回はpicture_outputもmovie_outputもoffにして，保存は行わないようにした．root権限でmotionを実行すると，8081ポートから動画が見れるようになる．
最後に，motionが自動で起動するように，まずはchown motion:motion /var/log/motion/motion.logとしてログファイルの権限を変更する．そして，systemctl start motionとすると，自動起動できるようになる．
まだ，装置が実際に動いているときの動作を確認していないので，なんとも言えないが，便利になると良いな．
あと，cookieで日本語を使っている場合に，CGI.unescape(ci[0])などとしなければならなかった．
qda=ARGV[0]
dat=qda.sub(/\.qda$/i,"_.dat")
d=open(qda,"rb"){|f|f.read}
n=d[0,16].unpack("n*")
c=n[8] # byte size of the last comment
n=n[1] # the number of data
sz=d[0x200,2*n].unpack("n*") # size of each data
tp=d[0x200+2*n,2*n].unpack("n*") # type
names=d[0x200+4*n,40*n].unpack("a40"*n).map(&:strip)
r=d[(0x200+4*n+40*n)..-1]
sz=sz.zip(tp).map{|s,t|
a,b="g","" # big endian float
a,b="N","Z21" if t==1 # string
a,b="G","" if t==2 # big endian 64bit float
*data,r=r.unpack(a*s+"n"*s+b*s+"a*")
msk=data[s,s] # mask
str=data[2*s,s] # string
data=data[0,s] # numeric
if t==2 # date or time
td=r[0,8].unpack("C*")[3]
data=data.map{|l|(Time.gm(l>3e9?1904:2000)+l).getutc.strftime("%d-%m-%Y")} if td==0 #date
data=data.map{|l|Time.at(l).getutc.strftime("%H:%M:%S")} if td==1 #12:34:56
data=data.map{|l|["%02d"]*2*":"%l.divmod(60)} if td==2 #12:34
end
r=r[136..-1]
data=str if t==1
data
}
mx=sz.map(&:size).max
sz=sz.map{|l|l.fill(0,l.size...mx)}.transpose
open(dat,"w"){|f|
f.puts "# "+r.gsub(/\x0d\x0a?/,"\n# ") if r.size>0
f.puts names*", "
sz.each{|l|f.puts l*", "}
}
まず，引数として与えたファイルを読み込むが，2バイト目からdataの数を，14バイト目からがコメントのサイズを表すようなので，それを取り込む．200hバイトから，それぞれのdataの長さがあり，その次に型だと思われるdataがある．その後に，40バイトずつのdataの名前があり，そのあとにdata本体とマスクだと思われるdataと文字列の場合には21バイドずつの文字列があり，その後にはdataの名前と良くわからないdataが合計136バイトあるようだ．ちなみに，それぞれはbig endianである．さらに，日付の表記は日月年の順番で，内部では秒単位の倍精度浮動小数点で表されているようである．時刻も同様である．ただ，文字列のdataの場合に，data本体には何が保存されているかはよく調べなかった．まあ，特定のversionのKaleidaGraphでセーブしたdataしか変換できないかも知れないが，今回の目的は達成したので，良しとしよう．