sbosse/simndt2b
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fri 16 May 23:06:54 CEST 2025

Browse Source
This commit is contained in:
sbosse 2025-05-16 23:08:29 +02:00
parent fc2237aed0
commit cea974c335
1 changed files with 18 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/SimNDT/core/simulation.py
View File

@ -3,6 +3,9 @@
""" """
simulation.py simulation.py
"Numercial" Material grid: Simulation.Im derived from scenario.Iabs in create_numericalModel,
just a resizing of Iabs (labeled material matrix)
Created by Miguel Molero on 2013-09-30. Created by Miguel Molero on 2013-09-30.
Copyright (c) 2013 MMolero. All rights reserved. Copyright (c) 2013 MMolero. All rights reserved.
""" """
@ -81,7 +84,7 @@ class Simulation:
self.t = self.dt*np.arange(0,self.TimeSteps) self.t = self.dt*np.arange(0,self.TimeSteps)
def create_numericalModel(self, scenario): def create_numericalModel(self, scenario):
print("Creating numerical simulation model (material labeled)")
#Spatial Scale #Spatial Scale
Pixel_mm = float(scenario.Pixel_mm) Pixel_mm = float(scenario.Pixel_mm)
@ -91,6 +94,20 @@ class Simulation:
self.Im = imresize(scenario.Iabs, self.Rgrid, interp='nearest') self.Im = imresize(scenario.Iabs, self.Rgrid, interp='nearest')
self.MRI, self.NRI = np.shape(self.Im) self.MRI, self.NRI = np.shape(self.Im)
def update_numericalModel(self, scenario, Im, Iabs):
print("Updating numerical simulation model (material labeled)")
#Spatial Scale
Pixel_mm = float(scenario.Pixel_mm)
if Im is not 0:
self.Im = Im
self.MRI, self.NRI = np.shape(self.Im)
if Iabs is not 0:
scenario.Iabs = Iabs
Mp = np.shape(scenario.Iabs)[0]*1e-3/Pixel_mm/self.dx
self.Rgrid = float(Mp/np.shape(scenario.Iabs)[0])
self.TapGrid = np.around(scenario.Tap * self.Rgrid)
self.Im = imresize(scenario.Iabs, self.Rgrid, interp='nearest')
self.MRI, self.NRI = np.shape(self.Im)
def rotate_model(self, theta, scenario): def rotate_model(self, theta, scenario):
Theta = theta * (180.0/np.pi) - 270. Theta = theta * (180.0/np.pi) - 270.