Tweaked versioning, added stats, parallel image stuff

Author: flip phillips

.vscode/tasks.json

@@ -17,6 +17,12 @@
             "type": "shell",
             "command": "/usr/bin/env wolframscript -f install.wls",
             "group": "build"
+        },
+        {
+            "label": "release",
+            "type": "shell",
+            "command": "./release.sh",
+            "group": "build"
         }
     ]
 }

CHANGELOG.md

@@ -16,7 +16,14 @@
 -->
 
 ## [Unreleased changes]
+### Added
+- `Statistics.wl` - routines from the pre v5 stats packages
+- `Overlap` - where do two lists of events overlap. There's probably a time-series function that does this but I just have data and don't want to have to cram it into a TS. (e.g., my data is spatial so TS is the wrong metaphor)
+- `ParallelImageApplyIndexed`
+- `GetCompilerFunctions`
 
+### Fixed
+- missing trailing ; in `Sugar.wl`
 
 * * *
 

FPTools/Data.wl

@@ -30,3 +30,5 @@ ImportAsDataset[f_,
 
    True,
    Dataset[raw]]]
+
+Overlap[d1_,d2_]:= {Max[{Min[d1], Min[d2]}],Min[{Max[d1], Max[d2]}]}

FPTools/FPTools.wl

@@ -14,7 +14,7 @@ Module[{files},
 
   (* load in any session *)
   files={"Introspection.wl","Entity.wl","Image.wl","Internet.wl","Sugar.wl",
-         "Versioning.wl","Globals.wl","Graphics.wl","Data.wl","Signals.wl","GA.wl"};
+         "Versioning.wl","Globals.wl","Graphics.wl","Data.wl","Signals.wl","GA.wl","Statistics.wl"};
 
   Map[Get[FileNameJoin[{DirectoryName[$InputFileName], #}]] &, files];
 

FPTools/Image.wl

@@ -164,3 +164,22 @@ LensDistortionCorrection[i_,{tx_,ty_},\[Theta]_,{p1_,p2_},kVec_,opt:OptionsPatte
 	rt=ut@*RotationTransform[\[Theta]];
 
 	ImageTransformation[i,rt[LensDistortPoint[tt[#],{p1,p2},kVec]]&,PlotRange->{{0.1,0.9},{0.1,0.5}},opt]]
+
+(* this is still in development-ville, right now it is for a single image plane *)
+
+ParallelImageApplyIndexed[f_, img_] := 
+ Module[{rows, cols, ix, id, decode, pd, res},
+  {rows, cols} = ImageDimensions[img];
+  
+  decode[i_] := {Quotient[i - 1, rows] + 1, Mod[i - 1, rows] + 1};
+  
+  (*ix=Table[{r,c},{r,1,cols},{c,1,rows}];*)
+  
+  ix = Table[i, {i, 1, rows*cols}];
+  id = Flatten[ImageData[img]];
+  
+  pd = Transpose[{id, ix}];
+  
+  res = ParallelMap[f[#[[1]], decode[#[[2]]]] &, pd];
+  Image[Partition[res, rows]]
+  ]

FPTools/Introspection.wl

@@ -124,3 +124,5 @@ GetSymbolsDefinedInMX[mxfile_] := Module[
         ];
         If[res=={{"System`Private`ConvertersPrivateDumpSymbol"}},None,res]
     ];
+
+GetCompilerFunctions[] := Sort[Compile`CompilerFunctions[]]

FPTools/PacletInfo.m

@@ -5,8 +5,8 @@
 	Publisher -> "Skidmore Vision Lab",
 	Copyright -> "©2018- Flip Phillips & Misc Contributors",
 	License -> "MIT",
-	Version -> "1.0.4",
-	BuildNumber -> "160",
+	Version -> "1.0.6",
+	BuildNumber -> "185",
 	MathematicaVersion -> "11.2+",
 	URL -> "https://github.com/flipphillips/FPTools",
 	Thumbnail -> "Documentation/icon.png",
@@ -21,7 +21,7 @@
 			"FPTools`$FPToolsVersion","FPTools`$FPToolsDebug",
 
 			(* data *)
-			"FPTools`ImportAsDataset", 
+			"FPTools`ImportAsDataset", "FPTools`Overlap",
 
 			(* Entity *)
 			"FPTools`EntityPropertiesDataset",
@@ -38,6 +38,7 @@
 			"FPTools`BlurFaces","FPTools`FindFaceImages",
 			"FPTools`AddAlphaChannel",
 			"FPTools`LensDistortPoint", "FPTools`LensDistortionCorrection",
+			"FPTools`ParallelImageApplyIndexed",
 
 			(* ffmpeg *)
 			"FPTools`$FFMpegPath", "FPTools`ImportMP4Frame", "FPTools`ImportMP4Frames", "FPTools`ImportWebFrame",
@@ -48,13 +49,25 @@
 			(* introspection *)
 			"FPTools`Spelunk","FPTools`PropertiesAvailability","FPTools`PropertiesDataset",
 			"FPTools`NewKernelEvaluate","FPTools`GetSymbolsDefinedInMX",
+			"FPTools`GetCompilerFunctions",
 
 			(* signals *)
 			"FPTools`FindZeroCrossings","FPTools`UnitizedGaussianPDF",
 			"FPTools`FFTShift",
 			"FPTools`InverseFFTShift",
 			"FPTools`RealSphericalHarmonicY",
 
+			(* statistics *)
+			"FPTools`LocationReport",
+			"FPTools`ShapeReport",
+			"FPTools`DispersionReport",
+			"FPTools`StandardDeviationMLE",
+			"FPTools`VarianceMLE",
+			"FPTools`StandardErrorOfSampleMean",
+			"FPTools`VarianceOfSampleMean",
+			"FPTools`CoefficientOfVariation",
+			"FPTools`ZeroMean",
+			
 			(* sugar *)
 			"FPTools`ApplyIf",
 			"FPTools`ToCartesianRules",

FPTools/Statistics.wl

@@ -0,0 +1,34 @@
+
+(* some routines I used to use from the old statistics packages *)
+
+LocationReport[data_] := {Mean -> Mean[data], 
+  HarmonicMean -> HarmonicMean[data], Median -> Median[data]} /; VectorQ[data]
+
+ShapeReport[list_] := {
+	Skewness -> Skewness[list] ,
+	QuartileSkewness -> QuartileSkewness[list] ,
+	Kurtosis -> Kurtosis[list]} /; VectorQ[list]
+
+DispersionReport[list_] :=
+	{
+	Variance -> Variance[list] ,
+	StandardDeviation -> StandardDeviation[list] ,
+	MinMax -> MinMax[list] ,
+	MeanDeviation -> MeanDeviation[list] ,
+	MedianDeviation -> MedianDeviation[list] ,
+	QuartileDeviation -> QuartileDeviation[list]
+	}   /; VectorQ[list]
+
+StandardDeviationMLE[data_] := 
+ StandardDeviation[data] Sqrt[(Length[data] - 1)/Length[data]]
+
+VarianceMLE[data_] := Variance[data] (Length[data] - 1)/Length[data]
+
+StandardErrorOfSampleMean[data_] := 
+ StandardDeviation[data]/Sqrt[Length[data]]
+
+VarianceOfSampleMean[data_] := Variance[data]/Length[data]
+
+CoefficientOfVariation[data_] := StandardDeviation[data]/Mean[data]
+
+ZeroMean[data_] := data - Mean[data]

FPTools/Sugar.wl

@@ -7,9 +7,9 @@ ApplyIf[f_, g_][x_] := ApplyIf[f, g, x];
 
 ToCartesianRules[] := 
   MapThread[
-   Rule, {{r, \[Theta], \[Phi]}, ToSphericalCoordinates[{x, y, z}]}]
+   Rule, {{r, \[Theta], \[Phi]}, ToSphericalCoordinates[{x, y, z}]}];
 
 ToSphericalRules[] := 
   MapThread[
    Rule, {{x, y, z}, 
-    FromSphericalCoordinates[{r, \[Theta], \[Phi]}]}];
+    FromSphericalCoordinates[{r, \[Theta], \[Phi]}]}];

FPTools/Usage.wl

@@ -20,6 +20,8 @@ PropertiesDataset::usage = "PropertiesDataset[object] returns a TableForm of the
 NewKernelEvaluate::usage = "NewKernelEvaluate[expr] evaluates expr in a fresh, clean kernel and returns the result.";
 GetSymbolsDefinedInMX::usage = "GetSymbolsDefinedInMX[mxfile] returns a list of symbols defined in mxfile.";
 
+GetCompilerFunctions::usage = "GetCompilerFunctions[] returns the list of functions that are compilable.";
+
 (* ::Section:: *)
 (*Entity*)
 
@@ -48,6 +50,7 @@ LensDistortPoint::usage = "LensDistortPoint[{x,y},{p1,p2},{k1,...}] projects poi
 
 LensDistortionCorrection::usage = "LensDistortionCorrection[image,{tx,ty},theta,{p1,p1},{k1,k2,k3...}] performs a lens distortion correction along with a generalized translation (t) and rotation (theta). p1 and p2 specify 1st order astygmatism and k1-6 the different orders of distortion beyond that. Note you only need to supply k's up to the desired order of correction.";
 
+ParallelImageApplyIndexed::usage = "ParallelImageApplyIndexed[f,im] maps f over im using parallel kernels. At this point, its only a single image plane";
 
 (* ::Section:: *)
 (*Internet*)
@@ -87,6 +90,8 @@ ImportAsDataset::usage = "ImportAsDataset[f] imports f and uses the first row as
 
 UnitizedGaussianPDF::usage = "UnitizedGaussian[mu,sigma] gives a pure function of a PDF whose value at the mean = 1.0.";
 
+Overlap::usage = "Overlap[l1,l2] gives the intersection of the range of two lists.";
+
 (* ::Section:: *)
 (*Signals*)
 
@@ -102,4 +107,22 @@ RealSphericalHarmonicY::usage = "RealSphericalHarmonicY[l,m,th,phi] returns real
 
 GARepopulate::usage = "GARepopulate[pop,fitfunc] creates a new population using fitfunc.";
 GAMutate::usage = "GAMutate[chromo] mutates a chromosome.";
-GACrossover::usage = "GACrossover[c1,c2] performs a random crossover of c1 & c2. Returns two offspring.";
+GACrossover::usage = "GACrossover[c1,c2] performs a random crossover of c1 & c2. Returns two offspring.";
+
+(* ::Section:: *)
+(*Statistics*)
+LocationReport::usage = "LocationReport[list] gives the Mean, HarmonicMean, and Median location statistics for list.";
+DispersionReport::usage = "DispersionReport[list] gives the Variance, StandardDeviation, SampleRange, MeanDeviation, MedianDeviation, and QuartileDeviation dispersion statistics for list.";
+ShapeReport::usage = "ShapeReport[list] gives the Skewness, QuartileSkewness, and KurtosisExcess shape statistics for list.";
+
+StandardDeviationMLE::usage = "StandardDeviationMLE[list] gives the standard deviation of the entries in list.  Division by n is used, giving a maximum likelihood estimate of the population standard deviation.";
+
+VarianceMLE::usage = "VarianceMLE[list] gives the variance of the entries in list. Division by n (rather than n-1) is used, giving a maximum likelihood estimate of the population variance (use Variance[list] for an unbiased estimate).";
+
+StandardErrorOfSampleMean::usage = "StandardErrorOfSampleMean[list] gives an unbiased estimate of the standard error (standard deviation) of the sample mean, using the entries in list as a sample from the population.";
+
+VarianceOfSampleMean::usage = "VarianceOfSampleMean[list] gives an unbiased estimate of the variance of the sample mean, using the entries in list as a sample from the population.";
+
+CoefficientOfVariation::usage = "CoefficientOfVariation[list] gives the coefficient of variation, defined as the ratio of the standard deviation to the mean of the entries in list. If the entries are not all positive, then $Failed is returned."
+
+ZeroMean::usage = "Zeroize the data about the mean.";