🤖

QRコードをQR分解する

2021/04/10に公開

Python

math

tech

概要

QRコードをQR分解します。

QRコードをQR分解する

QRコードを良く目にすると思います。例えばこんなのです。

これを見ると、まるで疎行列のように見えてきますね。なので、これを行列だと思ってQR分解したくなりますね。

QR分解とは、正方行列 $A$ を、直交行列 $Q$ と上三角行列 $R$ の積、

A = QR

と分解することです。なぜQR分解が必要かはその辺にいるガチ勢に聞いてください。ではさっそくQRコードをQR分解してみましょう。以下、Google Colabで実行することを想定していますが、何を使ってもかまいません。

まずはQRコードを作るのに必要なqrcodeをインストールしておきましょう。

!pip install qrcode

必要なライブラリをimportしましょう。

import qrcode
import numpy as np
import scipy.linalg as linalg
import matplotlib.pyplot as plt
from PIL import Image, ImageDraw

まずはQRコードを作ります。

BOX_SIZE = 10
qr = qrcode.QRCode(box_size=BOX_SIZE,border=0)
qr.add_data('Hello QR Code')
qr.make()
img = qr.make_image()
img

いきなりqrcode.make("Hello QR Code")としてQRコードを作ることもできますが、そうするとまわりに「枠」ができてしまうので、それを消すためにQRCodeクラスのコンストラクタでboder=0を指定しています。

ここから行列 $A$ を作りましょう。要素を見て、BOX_SIZE倍縮小するだけです。

LX, LY = img.size
data = np.array(img.getdata()).reshape(LX, LY)
X, Y = LX//BOX_SIZE, LY//BOX_SIZE
A = np.zeros((X, Y))
for ix in range(X):
  for iy in range(Y):
    A[ix][iy] = 255-data[ix*BOX_SIZE][iy*BOX_SIZE]

あとで使うので、行列を可視化する関数を作っておきましょう。

def get_image(M):
  M = np.abs(M)
  M = M/np.max(M)
  im = Image.new("L",(LX,LY), "white")
  draw = ImageDraw.Draw(im)
  for ix in range(X):
    for iy in range(Y):
      c = 255-int(M[iy][ix]*255)
      sx = ix*BOX_SIZE
      sy = iy*BOX_SIZE
      draw.rectangle((sx, sy, sx+BOX_SIZE, sy+BOX_SIZE), fill=c)
  return im