python – Czit notes

Konwersja formatu macierzowego deal.ii na zgodny z Matlab

przez kszyh | 10 sierpnia 2015 - 10:47 |10 sierpnia 2015 C++, Linux, Python

Format wyjściowy biblioteki deal.ii

(0,0) 6.73636e-09
(0,1) 2.76187e-09
(0,2) -3.27721e-09
(0,3) -3.6236e-11
(0,4) 1.13877e-08
(0,5) -5.36419e-10

(0,0) 6.73636e-09

(0,1) 2.76187e-09

(0,2) -3.27721e-09

(0,3) -3.6236e-11

(0,4) 1.13877e-08

(0,5) -5.36419e-10

Matlab wymaga numeracji od 1 i linii ze spacjami zaraz przecinków:

1 1  6.73636e-09
1 2  2.76187e-09
1 3  -3.27721e-09
1 4  -3.6236e-11
1 5  1.13877e-08
1 6  -5.36419e-10

1 1 6.73636e-09

1 2 2.76187e-09

1 3 -3.27721e-09

1 4 -3.6236e-11

1 5 1.13877e-08

1 6 -5.36419e-10

Skrypt dokonujący konwersji:

#!/usr/bin/python
import sys
import re

with open(sys.argv[1]+".new", 'a+') as w:
    with open(sys.argv[1]) as file:
        for line in file:
            m = re.match(r'\((.*?),(.*?)\)(.*?)$', line)
            w.write('{0} {1} {2}\n'.format(int(m.group(1))+1, int(m.group(2))+1, m.group(3)))

#!/usr/bin/python

import sys

import re

with open(sys.argv[1]+".new", 'a+') as w:

with open(sys.argv[1]) as file:

for line in file:

m = re.match(r'$(.*?),(.*?)$(.*?)$', line)

w.write('{0} {1} {2}\n'.format(int(m.group(1))+1, int(m.group(2))+1, m.group(3)))

Uruchomienie:

./matrix_convert.py A.dat

1	./matrix_convert.py A.dat

Rendering powierzchni NURBS w ParaView

przez kszyh | 7 marca 2014 - 16:16 |7 marca 2014 Lifehacking, Skrypty

Dodaj komentarz

try: paraview.simple
except: from paraview.simple import *
paraview.simple._DisableFirstRenderCameraReset()

points = {
1:[0,0,2],
10:[2.103215,2.103215,2],
11:[-7.55E-08,1.728084,2]
[...]
}

cm = [
[1,2,9],
[2,3,10],
[3,4,11],
[...]
]



for key, value in points.iteritems():

    Sphere4 = Sphere()
    Sphere4.Radius = 0.15
    Sphere4.Center = value
    Sphere4.PhiResolution = 20
    Sphere4.ThetaResolution = 20
    
    RenderView1 = GetRenderView()

    DataRepresentation4 = Show()
    DataRepresentation4.EdgeColor = [0.0, 0.0, 0.5000076295109483]
    DataRepresentation4.SelectionPointFieldDataArrayName = 'Normals'
    DataRepresentation4.AmbientColor = [0.0, 0.0, 0.0]
    DataRepresentation4.CubeAxesColor = [0.0, 0.0, 0.0]
    DataRepresentation4.ScaleFactor = 0.1

    RenderView1.CameraPosition = [0.12249341845958465, 1.6783114555009098, -2.82751287919202]
    RenderView1.CameraClippingRange = [1.8715334193481947, 5.084339235287565]

    DataRepresentation4.DiffuseColor = [1.0, 0.0, 0.0]

    Render()

    
    

for value in cm:
    for i in range(1, len(value)):
        Line2 = Line()

        # toggle the 3D widget visibility.
        active_objects.source.SMProxy.InvokeEvent('UserEvent', 'ShowWidget')
        RenderView1 = GetRenderView()
        RenderView1.CameraClippingRange = [0.04313226219643915, 43.13226219643914]

        fromi = value[0]
        toi = value[i]
        frompt = points[fromi]
        topt = points[toi]
        Line2.Point1 = frompt
        Line2.Point2 = topt
        
        DataRepresentation2 = Show()
        DataRepresentation2.EdgeColor = [0.0, 0.0, 0.5000076295109483]
        DataRepresentation2.SelectionPointFieldDataArrayName = 'Texture Coordinates'
        DataRepresentation2.AmbientColor = [0.0, 0.0, 0.0]
        DataRepresentation2.CubeAxesColor = [0.0, 0.0, 0.0]
        DataRepresentation2.ScaleFactor = 0.10752979516983033


        DataRepresentation2.Representation = 'Wireframe'
        DataRepresentation2.LineWidth = 2.0

        Render()

try: paraview.simple

except: from paraview.simple import *

paraview.simple._DisableFirstRenderCameraReset()

points = {

1:[0,0,2],

10:[2.103215,2.103215,2],

11:[-7.55E-08,1.728084,2]

[...]

}

cm = [

[1,2,9],

[2,3,10],

[3,4,11],

[...]

]

for key, value in points.iteritems():

Sphere4 = Sphere()

Sphere4.Radius = 0.15

Sphere4.Center = value

Sphere4.PhiResolution = 20

Sphere4.ThetaResolution = 20

RenderView1 = GetRenderView()

DataRepresentation4 = Show()

DataRepresentation4.EdgeColor = [0.0, 0.0, 0.5000076295109483]

DataRepresentation4.SelectionPointFieldDataArrayName = 'Normals'

DataRepresentation4.AmbientColor = [0.0, 0.0, 0.0]

DataRepresentation4.CubeAxesColor = [0.0, 0.0, 0.0]

DataRepresentation4.ScaleFactor = 0.1

RenderView1.CameraPosition = [0.12249341845958465, 1.6783114555009098, -2.82751287919202]

RenderView1.CameraClippingRange = [1.8715334193481947, 5.084339235287565]

DataRepresentation4.DiffuseColor = [1.0, 0.0, 0.0]

Render()

for value in cm:

for i in range(1, len(value)):

Line2 = Line()

# toggle the 3D widget visibility.

active_objects.source.SMProxy.InvokeEvent('UserEvent', 'ShowWidget')

RenderView1 = GetRenderView()

RenderView1.CameraClippingRange = [0.04313226219643915, 43.13226219643914]

fromi = value[0]

toi = value[i]

frompt = points[fromi]

topt = points[toi]

Line2.Point1 = frompt

Line2.Point2 = topt

DataRepresentation2 = Show()

DataRepresentation2.EdgeColor = [0.0, 0.0, 0.5000076295109483]

DataRepresentation2.SelectionPointFieldDataArrayName = 'Texture Coordinates'

DataRepresentation2.AmbientColor = [0.0, 0.0, 0.0]

DataRepresentation2.CubeAxesColor = [0.0, 0.0, 0.0]

DataRepresentation2.ScaleFactor = 0.10752979516983033

DataRepresentation2.Representation = 'Wireframe'

DataRepresentation2.LineWidth = 2.0

Render()

Archiwa tagu: python

Konwersja formatu macierzowego deal.ii na zgodny z Matlab

Rendering powierzchni NURBS w ParaView

Ostatnie wpisy

Najnowsze komentarze

Kategorie

Meta