質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

88.05%

ラスター画像からベクター画像への変換について

解決済

回答 2

投稿 編集

  • 評価
  • クリップ 1
  • VIEW 3,329

score 8

現在、図面の画像データをCADデータにしようとしています。

スクリプトでは、potraceというソフトとpythonを連携してラスター画像をベクター画像に変換しています。
ソフト:http://iphone.moo.jp/app/?p=386
アルゴリズム:http://sesamecake.blog84.fc2.com/blog-entry-155.html

おおよそ上手くできているのですが、以下のことをルールベースで改善したいと思っています.
・アウトラインの修正
potraceの特性上、ラスター画像での1本線はベクター画像では外枠の2本線で表現されるため、隣接する2本線を1本に統合したいです。
おすすめのアルゴリズムがあれば,ご教授いただきたいです.
イメージ説明

・修正方法
細かく見ると、雑な箇所もありこれらを特定の閾値を使って修正したいと思っていますが,手法が思いつきません.
また、上記のような箇所を特定するために、ラスベク変換において定量的に精度を示すものが有れば教えてください。
イメージ説明

よろしくお願い致します。

(0302 追記)
もう一点伺いたいことがあります.
ご回答いただいスクリプトで別画像のラスベク変換をした際,直角に曲がる線分が大きな曲線となる箇所がありましたが,原因は推測できるでしょうか.

よろしくお願いいたします.
イメージ説明

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 2

checkベストアンサー

+1

アウトラインの修正

Skeletonizeという手法が使えそうです。
これにて閉領域の中心を通る線を得られます。
さらにSkeleton Networkにてこれをグラフ化できます。
グラフのEdgeから線、あるいはBezier曲線を得ることができます。
以上を組み合わせると以下のようなSVGを出力するコードになります。

import numpy as np
from numpy import array, linalg, matrix
from scipy.special import comb
import cv2
from skimage.morphology import skeletonize
import sknw
import svgwrite

#
# Bézier curve fitting with SciPy
# https://stackoverflow.com/questions/12643079/b%C3%A9zier-curve-fitting-with-scipy
#
Mtk = lambda n, t, k: t**(k)*(1-t)**(n-k)*comb(n,k)
bezierM = lambda ts: matrix([[Mtk(3,t,k) for k in range(4)] for t in ts])

def lsqfit(points, M):
    M_ = linalg.pinv(M)
    return M_ * points

def get_bezier(pts):
    points = array(pts)
    ts = array(range(points.shape[0]), dtype='float')/(points.shape[0]-1)
    M = bezierM(ts)
    control_points = lsqfit(points, M)
    return control_points.tolist()


if __name__ == '__main__':

    fname = 'test'
    src = cv2.imread(fname + '.png', cv2.IMREAD_GRAYSCALE)
    _, src = cv2.threshold(src, 192, 255, cv2.THRESH_BINARY)

    # Skeleton化
    # https://scikit-image.org/docs/dev/auto_examples/edges/plot_skeleton.html
    ske = skeletonize(~(src != 0))
    ske_gray = (ske * 255).astype(np.uint8)
    ske_rgb = cv2.cvtColor(ske_gray, cv2.COLOR_GRAY2RGB)
    cv2.imwrite(fname + '_ske.png', ske_rgb)

    # Skeleton Networkを作成
    # https://github.com/Image-Py/sknw
    graph = sknw.build_sknw(ske.astype(np.uint16), multi=True)

    dwg = svgwrite.Drawing(fname + '.svg', profile='tiny')

    # Edge
    for (s,e) in graph.edges():

        pt_s = graph.node[s]['o'].tolist()
        pt_e = graph.node[e]['o'].tolist()

        for g in graph[s][e].values():

            # 開始 + 中間点 + 終点
            pts = g['pts'].tolist()
            pts = [pt_s] + pts + [pt_e]

            # 点群にフィットするBezierのパラメータを取得
            params = get_bezier(pts)

            # Bezierとして描画
            d = 'M{},{} C{},{} {},{}, {},{}'.format( params[0][0], params[0][1], params[1][0], params[1][1], params[2][0], params[2][1], params[3][0], params[3][1])
            p = dwg.path( d=d, stroke='#000', fill='none', stroke_width=5)
            dwg.add(p)

            # 線をそのまま描画
            #for i in range(len(pts)-1):
            #    dwg.add(dwg.line(pts[i], pts[i+1], stroke='#000', stroke_width=5))

    dwg.save()


結果は、まあそれなりですね。

元画像
イメージ説明
Skelton化した画像
イメージ説明
SVGをpng化した画像
イメージ説明

修正前

変換結果の画像を見ると、アウトラインのひとつひとつはもともとは2次元ポリゴンで構成されているかと思われます。
であればMedial axis、限定(直線)的にはStraight skeletonが求めたいラインに近いものになるかと思います。
pythonでのそのものずばりの実装は見つかりませんでしたが、以下が参考になるかと思います。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/03/03 00:58

    わざわざ,コードを書いてくださりありがとうございます!!

    おかげで,ラスベク変換がうまくいきました.
    skimageというモジュールがあることもわかり勉強になりました.

    よろしくお願い申し上げます。

    キャンセル

+1

隣接する2本線を1本に統合したいです

アップロードされたラスター画像は線が楕円形の丸になっているような箇所がありますので、ここの連続した楕円の連続から接続された線に復元するのは、恐らくラスター→ベクター変換より大変だろうと察します。

というのも、これをやろうとすると、楕円から長手方向の向きの検出、楕円ではない部分の長手方向の向きの検出、これらの接続が必要になるためです。この処理は、ラスター画像からベクター画像を起こす処理と同じだと思います。

そう考えると、タイトルにある(=本当にやりたいことの)

ラスター画像からベクター画像への変換

をやろうとするには、ラスター画像からそのまま「線画のベクター」を吐き出すソフトを探した方が筋が良さそうです。CADにしたい、とのことですが、使っているソフトがpotraceだそうですので、SVGで出力できれば十分と判断しました。

オフラインであればこれが使えそうです。inkscapeのモジュールのようですが、pythonで書かれていますので、必要な機能を抜き出して使えそうな気がします。

イメージ説明
※https://github.com/fablabnbg/inkscape-centerline-traceより借用


オンラインでもOKであればこんなのもありますね。

[ともにSVGに変換後PNGにして1/4に縮小 左:センターライン 右:アウトライン]
イメージ説明イメージ説明

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 88.05%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る