Fixing the design: In various part of the FireVoxel there is an inconsistent treatment of time values.  Certain structure hold time in "Absolute format" (i.e. seconds since Jan 1, 1900), other structures keeps them in seconds since the start of acquisitions.  This causes several glitches, f.e. during transfer of the Input function, from Automatic IDIF finder into Parametric Model calculation model.
Presently converting all the internal structure into the "Absolute format".  Relative times would only be constructed temporarily for specific purposes such as Display and during the Parametric Model calculations.

Defect Corrected:  When saving the TAC such as Input Function from ROIStats4D dialog box, the TAC name was improperly formatted and later incorrectly displayed inside the Parametric Models dialog.

Updated the layout of the Dynamic Experiment\Parametric Map calculation dialog.  Added new relevant buttons.

DCDI Liver mode: Performing yet another rearrangement of the variables.  Principal part is to have the S= ETS (extracellular transit speed) = 1/EMTT= (Fa+Fv+Ki)/Ve

Goal is the 100x times speed improvement during the Grid-search part of the optimization.  This can be achieved by using the ETS as the top variable of the nested loop. In that case all the exponential tables can
be precalculated once for each value of ETS.  Exponent calculations take about 99% of calculation time.

This is very relevant to plans to calculate the maps for the whole liver.

Performing rearrangement of Sourbron's DCDI equation, moving to new
optimization variables:


1. ETR = Ve/(Fa+Fv),
2. F = (Fa+Fv) ,
3.  fa = Fa/(Fa+Fv),
4. H = Ki/(Fa+Fv)

This might have definite advantages, due to simplified parameter ranges.
ETR is closely related to EMTT (extracellular mean transit time) and is in [0,100] sec interval.
F  is the total Inflow.

fa and H are just coefficients and should be in [0,1] range.

This seems to be a definite improvement over the {Ve,Fa,Fv,Ki }  variable set where all 4 nontrivial
intervals are required.

Calculated the Hepatic Uptake rate (Ki) for 2 patients. This is a very early and preliminary result.  Both maps use the same color scale. The units of measurement are "1/100/min" as per Sourbron's paper.

Implemented derived parameters as per Sourbron's paper.
Primary optimization parameters are Fa,Fv,Ve,Ki.
So far there was no need for the "Arterial Delay" as Henry had foreseen.

Performed DCDI calculation for two patients LP and KA.  Below are the "Total Inflow" maps (=Fa+Fv).  Color mapping is using the same color scale (see color bar on the right)

Liver DCDI modelling function is fast enough to process large ROIs on the voxel-wise basis
Below is the parametric map Arterial Flow (Fa) in  "ml/min/100ml" units as in Sourbron's paper.
Total calculation time is about 3 sec.

Liver DCDI model: 
a) Fully parallelized the "Grid search" preceding the iterative Amoeba  (Simplex) optimization
b) Integrated the Parallel Amoeba variant

Obtained first fitting result for the Liver DCDI model (Sourbron).
In this example 2 input functions are shown in Signal intensity units.
The target ROI is in red.
Fitted Concentration Curve (black) is overlaid on top of the Data Concentration Curve.
Optimal parameter values are shown at the bottom of the diagram. They are normalized
to Sourbron's unit from the paper and on a quick comparison are of the same order of magnitude.
Note:  Arterial delay of "0" was used.

Running tests of the Liver DCDI model, debugging and targeting the visualization of the model fit as the next waypoint.

DCDI (Dual Compartment Dual Input) model:  Implementing a fitting model, within the FireVoxel Dynamic Framework...

Implemented a support function, that in Play4D -  ALL the documents that have a compatible dynamic information will be played simultaneously.  This is very helpful for the displaying results of 4D registration,

In preparation to release Build 130:  Restoring several parts of the Automatic Build Test that use the .FVX file format. Format was changed recently due to the use of the new 3D Volumentric primitives.

Dialog "Split ROI by Threshold":
Implemented a feature, so user is able to enter the desired threshold manually.

Portal Vein Input function (obtained automatically) was compared between:
a)  Unregistered DCDI dataset
b) DCDI dataset registered for the whole liver (URAL-measure\Affine Transform)

difference is substantial, in line with the substantial motion across the axial plain which is eliminated as the result of the registration.

Compared the Portal Vein Input Function  derived from
a) Whole registered liver     VS.
b)  (starting from wholly registered liver) Local registration with the small ROI enveloping the Portal Vein.

There is small gross misregistration (during b.) on the initial timepoints due to the absence of features (contrast edges on early timepoint).  But besides that there seems to be little difference between a) and b).

So in the initial DCDI modelling we will restrict to the Whole liver registration only.

Developed a function, when given a 4D volume and an ROI (3D or 4D),   FVX advances to timepoint where amount of information (or entropy) over the given ROI is maximum. 

This is an extremely valuable function for all 4D registrations.  It allows user to automatically select an "Anchor" point for 4D registration.

Prior to this development, every individual workflow required an individually crafted recommendation which timepoint to use as an "anchor" for the registration.

Re-implementing the frequent {VolumeCropSlice, VolumeCropTimepoint} operations which are frequently used in processing.  Idea is to perform it through copying compressed blocks, instead of VolumeGetBox-VolumeSetBox sequence that require   Decompress\Compress.

Pair of this function is used extensively in Registrations and in over 70 locations inside the FireVoxel.