今度は、実写のステレオ写真（右目、左目用の写真がペアになったもの）から、奥行きデータを作成、凸凹写真（仮）にしてみました

拡大表示

画像はEmguのサンプルより。
コードは下に。
—-
この凸凹写真（仮）のアイディアは、 富士フイルム　FinePixREAL　3D　W1 の発売当初から考えていた、、、のですが、やっと（実験的にですが）　形にすることができました。

decoboco.py

import clr
clr.AddReference("Emgu.CV")
import Emgu.CV as cv
import System
import Emgu.CV.Structure as S
clr.AddReference("System.Drawing")
import System.Drawing as D
import Emgu.CV.CvEnum as E

fnleft="left.jpg"
fnright="right.jpg"
rev=0.
rev=1.

b=4
iml=cv.Image[S.Gray,System.Byte](fnleft)
imr=cv.Image[S.Gray,System.Byte](fnright)
im2=cv.Image[S.Gray,System.Int16](iml.Size)
im2r=cv.Image[S.Gray,System.Int16](iml.Size)
#bm=cv.StereoBM(E.STEREO_BM_TYPE.CV_STEREO_BM_BASIC,0)
#bm.FindStereoCorrespondence(iml,imr,im2)
gc = cv.StereoGC(8, 2)
gc.FindStereoCorrespondence(iml,imr,im2,im2r)
#im2 = im2 * (-16)
minvalue,maxvalue,minpoint,maxpoint=im2.MinMax()
dz=(maxvalue[0]-minvalue[0])
dz2=dz/2.0

if rev :
    #scale=-128*rev+128./(maxvalue[0]-minvalue[0])bm
    #im3=im2.Add(S.Gray(-dz2+minvalue[0]))
    im2-= minvalue[0]
    #im3._Mul(-1)
    #minvalue,maxvalue,minpoint,maxpoint=im3.MinMax()
    scale=130./(dz)
    im2*=scale
    im2=255-im2
    #im2+=128
    #im4=im2.ConvertScale[System.Byte](scale,64.)
#else :
#    scale=64./dz2
#    im3=im2.ConvertScale[System.Byte](scale,0.)

#imb=im3.SmoothBlur(16,16)

im2=im2.SmoothBlur(b,b).SmoothBlur(b,b)
im2.Save("000_z.jpg")

#for i in [iml,imr,im2,im2r,im3,imb]:
for i in [iml,imr,im2,im2r]:
    i.Dispose()

import shutil
if fnleft.endswith(".jpg"):
    shutil.copy(fnleft,"000_tex.jpg")
else:
    im=cv.Image[S.Bgr,System.Byte](fnleft)
    im.Save("000_tex.jpg")
    im.Dispose()

Related posts

• Pythonで名前つきリスト その2 (0)
• 【Py】webページの先頭数バイトだけ取得 (0)
• wget.py webをまとめて取得　などなど (0)
• CUIでプログレスバーもどき (0)
• ドロネー三角形 matplotlibで (0)

プリンタのインクリボン待ち中。

• 無料のMathematica Player + IronRuby でグラフィック出力（hirax.net::inside out）
  当blogの記事を取り上げて頂きました。
  Ruby(Rails)も憶えたいなあと思っているのですが、なかなか手が出せないでおります（＾＾；

その後、IronPythonから、MathematicaPlayerのフロントエンドを呼び出す方法が判ったのでメモ。

　Plot3Dなどを使って描画した３Ｄグラフ等を、マウスで自由にグリグリと向きを変えて見ることが出来るようになります。

まずは　いろいろと初期化
カーネルMathKernel.exeのパスを指定して起動

import clr
clr.AddReferenceToFile("Wolfram.NETLink")
from Wolfram.NETLink import *

kernelLink=MathLinkFactory.CreateKernelLink(
"-linkmode autolaunch -linkname 'C:/Program Files/Wolfram Research/Mathematica Player/7.0/mathkernel.exe'")

kernelLink.WaitAndDiscardAnswer()

ここまでは従前どおり。

次に、フロントエンドMathematicaPlayer.exeの絶対パスを指定して起動します。

#主にkernellinkオブジェクトのEvaluatetoInputFormメソッドを使っていきます
＞＞＞　kernelLink.EvaluateToInputForm("1+1",0)
2  # 1+1なので、２が返ってきました

#ConnectToFrontEnd などの、フロントエンド関連命令を使えるようにします。
kernelLink.EvaluateToInputForm('Needs["JLink`"]',0)

#フロントエンドを起動します。
kernelLink.EvaluateToInputForm('ConnectToFrontEnd["C:/Program Files/Wolfram Research/Mathematica Player/7.0/MathematicaPlayer.exe"]',0)

　

これで、フロントエンドを使う準備が完了しました。

Plot3D関数を使い、3次元グラフを作成

UseFrontEnd関数、CreatekernelLink.EvaluateToInputFormcumentを使いフロントエンドに表示します。

p=kernelLink.EvaluateToInputForm("p=Plot3D[Cos[x*y],{x,-3,3},{y,-3,3}]",0)

kernelLink.EvaluateToInputForm("UseFrontEnd[CreateDocument[ p ]]",0)

グラフが、MathematicaPlayerのウィンドウに表示されたら成功です

マウスでドラッグすると、グリグリと向きが変化します。

使い終わったら、カーネルとの接続を切断します。

#フロントエンド MathematicaPlayer.exeを終了します
kernelLink.EvaluateToInputForm('CloseFrontEnd[]',0)

#MathKernel.exeを終了します。
kernelLink.Close()

残念ながら、インタラクティブにグラフを変化させたりとか、そういうのは判りませんでした。

Related posts

• 【iPy】凸凹写真（仮）　ステレオ・ペア写真から奥行きデータ作成編 (0)
• 【もうすぐW杯】サッカーボールな地球儀のペーパークラフト (0)
• 地震雲？ (0)
• 怪獣算 (0)
• Mathematicaでメタボール (0)

こちらの続きです

IronPythonから無料のMathematica Playerのカーネルを呼び出して、グラフなどの画像を出力してみました。

>>> from netlink import *
>>> m=Mathematica()
>>> m.open()
>>> im=m.im("Plot[Sin[X],{x,0,2Pi}]",100,"RGB")
>>> im.Save("sin.png")
>>> im.Dispose()

（前略） 上記の続きで
m.do("""solution=
  NDSolve[
    {Derivative[1][x][t]== -3*(x[t]-y[t]),
     Derivative[1][y][t]==
       26.5*x[t]-y[t]-x[t]*z[t],
     Derivative[1][z][t]==x[t]*y[t]-z[t],
     x[0]==z[0]==0,
     y[0]==2},{x,y,z},{t,0,17}];""")

im=m.im("""ParametricPlot3D[
  Evaluate[{x[t],y[t],z[t]}/.solution],
   {t,0,17},PlotPoints ->1000,ViewPoint->{1,0,1}]""",500,"RGBA")

im.Save("r.jpg")

im=m.im("""ListLinePlot[{Re[#], Im[#]} & /@
   Accumulate[Exp[I 10.^4 EulerGamma Sqrt[Range[5000]]]],Axes->False ]""",500,"RGBA")

m.im("Plot[Sin[x],{x,0,2Pi}]",100,"RGBA")

netlink.pyとして保存

#-*- coding:utf-8 -*-
"netlink.py"
import clr , System
from System import Array,Byte
clr.AddReferenceToFile("Wolfram.NETLink")
from Wolfram.NETLink import *
import sys
import image

exe1="mathkernel.exe"
exe2="MathematicaPlayer.exe"
mlpath="-linkmode autolaunch -linkname "
#Full Path To MathKernel.exe
mlpath+=" 'C:\\Program Files\\Wolfram Research\\Mathematica Player\\7.0\\%s' "

def fe():
  NL=MathLinkFactory.CreateMathLink(mlpath%exe2)
  return NL

def ev(txt):
  return eval(txt.replace("{","(").replace("}",")"))

class Mathematica(object):
  def open(self):
    "init"
    self.KL=MathLinkFactory.CreateKernelLink(mlpath%exe1)
    self.KL.WaitAndDiscardAnswer()

  def do(self,command):
    "evaluate to input form"
    return self.KL.EvaluateToInputForm(command, 0)

  def get(self,name):
    "get value in symbol"
    #変数の値を、.NETのオブジェクトとして読み出し
    #変数のみ。　式の評価結果などを直接読み出すことは出来ません。
    self.KL.PutFunction("EvaluatePacket",1)
    self.KL.PutSymbol(name)
    self.KL.EndPacket()
    self.KL.WaitForAnswer()
    try:
      ret=self.KL.GetObject()
    except Exception:
      self.KL.ClearError()
      self.do("")
      self.KL.Flush()
      while not self.KL.Ready:
        pass
      while self.KL.Ready:
        self.KL.WaitAndDiscardAnswer()
      ret=self.do(name)
    return ret

  def imagedata(self,txt,w,mode="RGB"):
    "evaluate to image data array"
    if mode=="RGBA":
      bg=",Background->None"
      self.do("$ipycolortrans={{0,0,1,0},{0,1,0,0},{1,0,0,0},{0,0,0,1}};")
    else :
      bg=""
      self.do("$ipycolortrans={{0,0,1},{0,1,0},{1,0,0}};")

    tmpl='$ipyimagedata=Rasterize[%(txt)s,"Image",RasterSize->%(w)s%(bg)s];'
    self.do(tmpl%locals())

    if self.do('ImageQ[$ipyimagedata]')!="True":
      return
    self.do('$ipyimagedata=ImageData[$ipyimagedata,"Byte"];')
    self.do('$ipyimagedata2=$ipycolortrans.#&[#]&/@#&/@$ipyimagedata;')
    return self.get("$ipyimagedata2")    

  def im(self,txt,w,mode="RGB"):
    "evaluate to image"
    dat=self.imagedata(txt,w,mode)
    if dat :
      return image.fromArray(dat)

  def close(self):
    self.KL.EvaluateToInputForm('Quit[]', 0)
    self.KL.Close()

Related posts

• Pythonで名前つきリスト その2 (0)
• 【Py】webページの先頭数バイトだけ取得 (0)
• wget.py webをまとめて取得　などなど (0)
• CUIでプログレスバーもどき (0)
• ドロネー三角形 matplotlibで (0)

IronPythonで画像ファイルを扱うための簡単なモジュール。
憶えやすいように、関数名をPILっぽい名前にしてあります。

あと、IronPython2.0から、PILやNumpyなどの、Cで書かれたモジュールを呼び出す方法もメモ。

動作確認 IronPython RC2

image.pyとして保存

# -*- coding: utf8 -*-
"image.py　　　small pil"
import clr
import System
from System import Array,Byte
clr.AddReferenceByPartialName("System.Drawing")
from System.Drawing import Bitmap,Rectangle,RotateFlipType,Color
from System.Drawing.Imaging import PixelFormat,ImageLockMode
from System.Runtime.InteropServices import Marshal

def new (mode,size):
  w,h=size
  mode2,n=_mode(mode)
  bmp = Bitmap(w, h, mode2)
  return bmp

def _mode(mode):
  mode=mode.lower()
  if mode=="rgb":
    mode2=PixelFormat.Format24bppRgb
    nchannel=3
  elif  mode in "rgba" :
    mode2=PixelFormat.Format32bppArgb
    nchannel=4
  else :
    raise 'mode = "RGB" or "RGBA"'
  return mode2,nchannel

def fromArray(*args):
  if len(args)==3:
    mode,size,dat=args
    return fromArray1D(mode,size,dat)
  elif len(args)==1 :
    dat=args[0]
    if isinstance(dat,Array):
      h,w,n=[dat.GetLength(i) for i in xrange(3)]
      a=Array[Byte](dat)
      mode=[None,None,None,"RGB","RGBA",][n]
      return fromArray1D(mode,(w,h),a)
  return None

def fromArray1D(mode,size,dat):
  #http://d.hatena.ne.jp/okazuki/20060828/1156726184
  w,h=size
  mode2,n=_mode(mode)
  bmp=new(mode,(w,h))
  rct=Rectangle(0, 0, w, h)
  bd=bmp.LockBits( rct, ImageLockMode.WriteOnly, mode2);
  Marshal.Copy(dat, 0, bd.Scan0, w*h*n);
  bmp.UnlockBits(bd);
  bmp.RotateFlip(RotateFlipType.Rotate180FlipXY)
  return bmp

def frombuffer(mode,size,buf):
  dat=Array[Byte](buf)
  return fromArray(mode,size,buf)

def iterInt(txt):
    pat=re.compile(r"(\d+)")

    return (int(x.group(1)) for x in pat.finditer(txt) )

def toArray(bmp):
  u"Bitmapの内部データにアクセスする"
  #"http://d.hatena.ne.jp/Guernsey/20090226"
  # http://blogs.yahoo.co.jp/spike_spike690/4735887.html
  # Pixelデータにアクセスする領域を設定する
  w,h=bmp.Width,bmp.Height
  rect = Rectangle(0, 0, w,h)
  # BitmapDataクラスのインスタンスを生成し、LockBits関数でBitmapをLockする。
  bmpData = bmp.LockBits(rect,ImageLockMode.ReadOnly, PixelFormat.Format24bppRgb)
  # BitmapData内の各種項目を抜き出す。
  ptr = bmpData.Scan0
  imageStride =bmpData.Stride
  # コピーするデータ量を設定．PixelFormat に関わらず stride * 高さ
  bytes = imageStride * h
  imageData = Array[Byte](Length(bytes))
  # 確保したデータ配列にBitmap内のデータをコピー
  Marshal.Copy(ptr, imageData, 0, bytes)
  # Bitmapのロックを解除する。

  bmp.UnlockBits(bmpData)
  return imageData

def fromMstring(mode,size,txt):
  "from string to Bitmap"
  w,h=size
  mode2,n=_mode(mode)
  dat=Array[Byte](Length(w*h*n,iterInt(txt)))
  im=fromArray(mode,size,dat)
  return im
def toMstring(x):
  if isinstance(x,Bitmap) :
    bmp=x
    w=bmp.Width
    n=bmp.GetPixelFormatSize(bmp.PixelFormat)//8
    a=toArray(bmp)
    return 'Partition[Partition[{%s},%s],%s]'%(
                        ",".join((str(i) for i in a)),n,w)

#  #http://sonic64.com/2005-08-04.html
def fromstring(mode,size,txt):
  "fromstring PIL style"
  w,h=size
  mode2,n=_mode(mode)
  dat=Encoding.Unicode.GetBytes(txt)[::2]
  im=fromArray(mode,size,dat)
  return im

def tostring(x):
  "tostring PIL style"
  bmp=x
  a=toArray(bmp)
  return "".join((chr(i) for i in a))

def load(fname):
  return Bitmap(fname)
import re

class Length(object):
  def __init__(self,length,itr=None):
    self.length=length
    self.itr=itr
  def __len__(self):
    return self.length
  def __iter__(self):
    if self.itr :
      length=self.length
      def counter(var=[0]):
        var[0]+=1
        if var[0]>length :
          raise StopIteration
        return var[0]
      return (counter() and i for i in self.itr )
    else :
      return ( 0 for i in xrange(self.length))

Related posts

• Pythonで名前つきリスト その2 (0)
• 【Py】webページの先頭数バイトだけ取得 (0)
• wget.py webをまとめて取得　などなど (0)
• CUIでプログレスバーもどき (0)
• ドロネー三角形 matplotlibで (0)