「第16回燐字書き取り」分割

第16回ですかー

分割した画像はこれ。

1. leastsquare.py

私(skarsna.meloviliju)にとって二度目の作業、13回のログをたよりに行う。

原稿の四隅(1枚目)に位置する4字の（マスの左上の）座標をGIMPで取ると、

      y1,   y2       s1 s2
左上:   58,   58 --->  0, 0
右上:  921,   58 --->  9, 0
左下:   58, 1314 --->  0,19
右下:  920, 1314 --->  9,19

だったので、これをleastsquare.pyのsとyに突っ込んでやると下の理論が動いて必要な行列がでるらしい。科学の力ってすごいね。

うまいパラメータt₁₀, t₁₁, t₁₂, t₂₀, t₂₁, t₂₂を持ってきて、

となるようにしてやれば良い。ただしs₀たちは全部1。

Y ≒ T₁Sで最小二乗法をやるのでYS^T=T₁SS^Tを解けばよく、T₁=YS^T(SS^T)^-1なのでT₁は

[[ 2.81250000e+03  1.43000000e+02 -5.26315789e-02]
[ 3.17000000e+02  0.00000000e+00  9.92631579e+01]]

ほーん。

2. split.py

linziListに燐字の表を打ち込む。これが一番疲れた。どういう対応かは画像とファイル見て察してほしい。

split.pyを走らせる。

SystemError: tile cannot extend outside image

は？

存在しない画像範囲にアクセスして死んでるらしいが、前回と違って出力画像が見えない。

j.v.の力を借りて30分ほど粘ったが、プログラムには問題がなさそうということでファイル形式を確認したところ、jpgファイルを使用していることが「判明」。

jpgは座標でなんか悪さをするらしいのでGIMPでpng化してsplit.pyを走らせるとうまくいった。よかったね。