[PYTHON] Draw a "breast curved surface" in a 3D graph (1)

"Draw a graph in the programming language Julia" I succeeded in moving the sample for the time being, but the data created by my own program (in Julia language) I decided to find out how to draw a graph using.

This time, referring to "Practical technique of Julia language" learned from "Breast curved surface", "Breasts" in 3D space I made a sample to draw a "curved surface". As usual, the calculation result obtained by Julia is displayed as a graph in PyQtGraph. (Please create an environment and move it)

• Oppai.jl (Julia source code) passes an array to a function called graph_run () in the Python source. Let Python draw the graph for you. ..
• Oppai.py (Python source code) draws a 3D graph based on the array passed by Julia (51 line graphs).
• When calling Python code using PyCall, C: \ julia-0.2.1-win32 \ bin is included in the load path, so both Oppai.jl and Oppai.py are C: \ julia- Save it under 0.2.1-win32 \ bin.
• When saving the source code, the character encoding should be utf-8.
• To launch the sample, double-click Oppai.jl (Julia source code) from Explorer.

Note: For the environment construction, refer to "Drawing a graph in the programming language Julia".

↑ Drag with the mouse or press the cursor key to rotate the 3D graph. You can zoom in / out by turning the mouse wheel.

jl:Oppai.jl(C:\julia-0.2.1-win32\Save to bin)

# -*- coding: utf-8 -*-

function bust(x,y)
1/8* (6*exp(-((2/3*abs(x) - 1)^2 + (2/3 *y)^2) - 1/3*(2/3*y + 1/2)^3)+ 2/3 *exp(-2.818^11*((abs(2/3*x) - 1)^2+ (2/3 *y)^2)^2) + 2/3*y - (2/3*x)^4)
end

nx = 100
ny = 51
x = linspace(-3, 3, nx)
y = linspace(-3, 3, ny)
z = zeros(Float64, length(y), length(x))

for iy = 1:length(y)
    for ix = 1:length(x)
        z[iy, ix] = bust(x[ix], y[iy])
    end
end

using PyCall
@pyimport Oppai
Oppai.graph_run(x, y, z, 10)

python:Oppai.py(C:\julia-0.2.1-win32\Save to bin)

# -*- coding: utf-8 -*-

from pyqtgraph.Qt import QtCore, QtGui
import pyqtgraph.opengl as gl
import pyqtgraph as pg
import numpy as np

def graph_run(plot_data_x, plot_data_y, plot_data_z, distance):

    app = QtGui.QApplication([])

    w = gl.GLViewWidget()
    w.resize(600,400)
    w.opts['distance'] = distance
    w.show()
    w.setWindowTitle(u'Julia+PyQtGraph sample')

    x = plot_data_x
    y = plot_data_y
    for i in range(len(y)):
        yi = [y[i]]*len(x)
        z = plot_data_z[i,:]
        pts = np.vstack([x,yi,z]).transpose()
        plt = gl.GLLinePlotItem(pos=pts, color=pg.glColor((i,len(y)*1.3)), width=(i+1)/10., antialias=True)
        w.addItem(plt)

    app.exec_()