Skip to content

Commit

Permalink
fixed issue #48
Browse files Browse the repository at this point in the history
  • Loading branch information
zitmen committed Sep 4, 2016
1 parent 727f448 commit 79ca75f
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 7 deletions.
2 changes: 1 addition & 1 deletion pom.xml
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,7 @@

<groupId>cz.cuni.lf1.lge</groupId>
<artifactId>ThunderSTORM</artifactId>
<version>dev-2016-07-27-b1</version>
<version>dev-2016-09-04-b1</version>
<packaging>jar</packaging>

<name>ThunderSTORM</name>
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -53,14 +53,13 @@ protected void drawPoint(double x, double y, double z, double dx, double dz) {
//3D gaussian blob integrated in over voxel
//mathematica function for definite integral:
//Integrate[1/(Sqrt[(2*Pi)^3*s1^2*s2^2*s3^2]*E^((1/2)*((x - x0)^2/s1^2 + (y - y0)^2/s2^2 + (z - z0)^2/s3^2))), {z, a, b}, {x, ax, ax + 1}, {y, ay , ay + 1}]
final int MAXRADIUS = (int) (defaultDX / resolution * 5);
if(isInBounds(x, y)) {
x = (x - xmin) / resolution;
y = (y - ymin) / resolution;
dx = dx / resolution;
int u = (int) x;
int v = (int) y;
int actualRadius = (this.radius < 0) ? (int) Math.min(ceil(dx * 3), MAXRADIUS) : this.radius;
int actualRadius = (this.radius < 0) ? (int) ceil(dx * 3) : this.radius;
double sqrt2dx = Math.sqrt(2) * dx;

int w = threeDimensions ? ((int) ((z - zFrom) / zStep)) : 0;
Expand Down Expand Up @@ -101,8 +100,4 @@ protected void drawPoint(double x, double y, double z, double dx, double dz) {

}
}

private double squareDist(double x1, double y1, double x2, double y2) {
return (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2);
}
}

0 comments on commit 79ca75f

Please sign in to comment.