Skip to content
Projects
Groups
Snippets
Help
Loading...
Help
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
I
imagej-elphel
Project
Project
Details
Activity
Releases
Cycle Analytics
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Charts
Issues
3
Issues
3
List
Board
Labels
Milestones
Wiki
Wiki
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Charts
Create a new issue
Commits
Issue Boards
Open sidebar
Elphel
imagej-elphel
Commits
a9f8afde
Commit
a9f8afde
authored
Jun 26, 2019
by
Andrey Filippov
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
committing w/o breaking - 2
parent
61820deb
Changes
3
Show whitespace changes
Inline
Side-by-side
Showing
3 changed files
with
71 additions
and
30 deletions
+71
-30
Aberration_Calibration.java
...com/elphel/imagej/calibration/Aberration_Calibration.java
+2
-0
MatchSimulatedPattern.java
.../com/elphel/imagej/calibration/MatchSimulatedPattern.java
+53
-30
EProperties.java
src/main/java/com/elphel/imagej/common/EProperties.java
+16
-0
No files found.
src/main/java/com/elphel/imagej/calibration/Aberration_Calibration.java
View file @
a9f8afde
...
@@ -375,6 +375,8 @@ public static MatchSimulatedPattern.DistortionParameters DISTORTION =new MatchSi
...
@@ -375,6 +375,8 @@ public static MatchSimulatedPattern.DistortionParameters DISTORTION =new MatchSi
// additionally multiply by Hamming
// additionally multiply by Hamming
256, // FFTSize (was 128)
256, // FFTSize (was 128)
32, // FFTSize_lwir
32, // FFTSize_lwir
32, // FFTOverlap (was 32)
4, // FFTOverlap_lwir
0.5, //fftGaussWidth
0.5, //fftGaussWidth
0.0, //phaseCorrelationFraction
0.0, //phaseCorrelationFraction
1.5, // 2.5, //6.0, // 2.0, // 0.0, // correlationHighPassSigma, - pixels in frequency domain
1.5, // 2.5, //6.0, // 2.0, // 0.0, // correlationHighPassSigma, - pixels in frequency domain
src/main/java/com/elphel/imagej/calibration/MatchSimulatedPattern.java
View file @
a9f8afde
...
@@ -725,19 +725,19 @@ public class MatchSimulatedPattern {
...
@@ -725,19 +725,19 @@ public class MatchSimulatedPattern {
/* returns array of 3 arrays: first two are 3-element wave vectors (x,y,phase), last - 3-rd order correction coefficients */
/* returns array of 3 arrays: first two are 3-element wave vectors (x,y,phase), last - 3-rd order correction coefficients */
public
double
[][]
findPatternDistorted
(
public
double
[][]
findPatternDistorted
(
double
[][]
bayer_pixels
,
// pixel array to process (no windowing!), two greens will be used
double
[][]
bayer_
mono_
pixels
,
// pixel array to process (no windowing!), two greens will be used
PatternDetectParameters
patternDetectParameters
,
PatternDetectParameters
patternDetectParameters
,
double
min_half_period
,
double
min_half_period
,
double
max_half_period
,
double
max_half_period
,
boolean
greens
,
// this is a pattern for combined greens (diagonal), adjust results accordingly
boolean
greens
,
// this is a pattern for combined greens (diagonal), adjust results accordingly
String
title
){
// title prefix to use for debug images
String
title
){
// title prefix to use for debug images
if
(
bayer_pixels
==
null
)
return
null
;
if
(
bayer_
mono_
pixels
==
null
)
return
null
;
if
(
bayer
_pixels
.
length
<
4
)
return
null
;
if
(
bayer_mono
_pixels
.
length
<
4
)
return
null
;
if
(
bayer_pixels
[
0
]==
null
)
return
null
;
if
(
bayer_
mono_
pixels
[
0
]==
null
)
return
null
;
if
(
bayer
_pixels
[
3
]==
null
)
return
null
;
if
(
bayer_mono
_pixels
[
3
]==
null
)
return
null
;
int
size2
=
bayer
_pixels
[
0
].
length
;
int
tile_size2
=
bayer_mono
_pixels
[
0
].
length
;
int
size
=(
int
)
Math
.
sqrt
(
size2
);
int
tile_size
=(
int
)
Math
.
sqrt
(
tile_
size2
);
int
hsize
=
size
/
2
;
int
hsize
=
tile_
size
/
2
;
int
hsize2
=
hsize
*
hsize
;
int
hsize2
=
hsize
*
hsize
;
double
[]
quarterHamming
=
initWindowFunction
(
hsize
,
patternDetectParameters
.
gaussWidth
);
double
[]
quarterHamming
=
initWindowFunction
(
hsize
,
patternDetectParameters
.
gaussWidth
);
double
[]
patternCorr
=
new
double
[
6
];
// second order non-linear pattern correction (perspective+distortion)
double
[]
patternCorr
=
new
double
[
6
];
// second order non-linear pattern correction (perspective+distortion)
...
@@ -747,25 +747,25 @@ public class MatchSimulatedPattern {
...
@@ -747,25 +747,25 @@ public class MatchSimulatedPattern {
double
[][]
quarter_pixels
=
new
double
[
9
][];
double
[][]
quarter_pixels
=
new
double
[
9
][];
double
[][][]
quarter_patterns
=
new
double
[
9
][][];
double
[][][]
quarter_patterns
=
new
double
[
9
][][];
int
[]
quarterIndex
={
0
,
// top left
int
[]
quarterIndex
={
0
,
// top left
size
/
2
,
// top right
tile_
size
/
2
,
// top right
size
*
size
/
2
,
// bottom left
tile_size
*
tile_
size
/
2
,
// bottom left
(
size
+
1
)*
size
/
2
,
// bottom right
(
tile_size
+
1
)*
tile_
size
/
2
,
// bottom right
(
size
+
1
)*
size
/
4
,
// center
(
tile_size
+
1
)*
tile_
size
/
4
,
// center
size
/
4
,
// top
tile_
size
/
4
,
// top
size
*
size
/
4
,
// left
tile_size
*
tile_
size
/
4
,
// left
(
size
+
2
)*
size
/
4
,
// right
(
tile_size
+
2
)*
tile_
size
/
4
,
// right
(
2
*
size
+
1
)*
size
/
4
};
// bottom
(
2
*
tile_size
+
1
)*
tile_
size
/
4
};
// bottom
int
i
,
j
,
iq
;
int
i
,
j
,
iq
;
int
index
,
qindex
;
int
index
,
qindex
;
if
(
this
.
debugLevel
>
2
)
SDFA_INSTANCE
.
showArrays
(
bayer_
pixels
,
size
,
size
,
title
+
"-bayer"
);
if
(
this
.
debugLevel
>
2
)
SDFA_INSTANCE
.
showArrays
(
bayer_
mono_pixels
,
tile_size
,
tile_
size
,
title
+
"-bayer"
);
for
(
iq
=
0
;
iq
<
9
;
iq
++)
{
for
(
iq
=
0
;
iq
<
9
;
iq
++)
{
index
=
quarterIndex
[
iq
];
index
=
quarterIndex
[
iq
];
qindex
=
0
;
qindex
=
0
;
for
(
i
=
0
;
i
<
hsize
;
i
++)
{
for
(
i
=
0
;
i
<
hsize
;
i
++)
{
for
(
j
=
0
;
j
<
hsize
;
j
++)
{
//quarter_pixels[iq][qindex++]=input_pixels[index++];
for
(
j
=
0
;
j
<
hsize
;
j
++)
{
//quarter_pixels[iq][qindex++]=input_pixels[index++];
green0
[
qindex
]=
bayer
_pixels
[
0
][
index
];
green0
[
qindex
]=
bayer_mono
_pixels
[
0
][
index
];
green3
[
qindex
++]=
bayer
_pixels
[
3
][
index
++];
green3
[
qindex
++]=
bayer_mono
_pixels
[
3
][
index
++];
}
}
quarter_pixels
[
iq
]=
combineDiagonalGreens
(
green0
,
green3
,
hsize
,
hsize
);
quarter_pixels
[
iq
]=
combineDiagonalGreens
(
green0
,
green3
,
hsize
,
hsize
);
index
+=
hsize
;
// jump to the next line
index
+=
hsize
;
// jump to the next line
...
@@ -11554,6 +11554,9 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11554,6 +11554,9 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
// additionally multiply by Hamming
// additionally multiply by Hamming
public
int
FFTSize
;
public
int
FFTSize
;
public
int
FFTSize_lwir
;
public
int
FFTSize_lwir
;
public
int
FFTOverlap
;
// 32 used for aberration kernels, former FFT_OVERLAP
public
int
FFTOverlap_lwir
;
// 4
public
double
fftGaussWidth
;
public
double
fftGaussWidth
;
public
double
phaseCorrelationFraction
=
1.0
;
// 1.0 - phase correlation, 0.0 - just cross-correlation
public
double
phaseCorrelationFraction
=
1.0
;
// 1.0 - phase correlation, 0.0 - just cross-correlation
public
double
correlationHighPassSigma
;
public
double
correlationHighPassSigma
;
...
@@ -11628,6 +11631,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11628,6 +11631,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
int
zeros
,
int
zeros
,
int
FFTSize
,
int
FFTSize
,
int
FFTSize_lwir
,
int
FFTSize_lwir
,
int
FFTOverlap
,
int
FFTOverlap_lwir
,
double
fftGaussWidth
,
double
fftGaussWidth
,
double
phaseCorrelationFraction
,
double
phaseCorrelationFraction
,
double
correlationHighPassSigma
,
double
correlationHighPassSigma
,
...
@@ -11695,6 +11700,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11695,6 +11700,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
this
.
zeros
=
zeros
;
this
.
zeros
=
zeros
;
this
.
FFTSize
=
FFTSize
;
this
.
FFTSize
=
FFTSize
;
this
.
FFTSize_lwir
=
FFTSize_lwir
;
this
.
FFTSize_lwir
=
FFTSize_lwir
;
this
.
FFTOverlap
=
FFTOverlap
;
this
.
FFTOverlap_lwir
=
FFTOverlap_lwir
;
this
.
fftGaussWidth
=
fftGaussWidth
;
this
.
fftGaussWidth
=
fftGaussWidth
;
this
.
phaseCorrelationFraction
=
phaseCorrelationFraction
;
this
.
phaseCorrelationFraction
=
phaseCorrelationFraction
;
this
.
correlationHighPassSigma
=
correlationHighPassSigma
;
this
.
correlationHighPassSigma
=
correlationHighPassSigma
;
...
@@ -11764,6 +11771,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11764,6 +11771,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
this
.
zeros
,
this
.
zeros
,
this
.
FFTSize
,
this
.
FFTSize
,
this
.
FFTSize_lwir
,
this
.
FFTSize_lwir
,
this
.
FFTOverlap
,
this
.
FFTOverlap_lwir
,
this
.
fftGaussWidth
,
this
.
fftGaussWidth
,
this
.
phaseCorrelationFraction
,
this
.
phaseCorrelationFraction
,
this
.
correlationHighPassSigma
,
this
.
correlationHighPassSigma
,
...
@@ -11832,6 +11841,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11832,6 +11841,8 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
properties
.
setProperty
(
prefix
+
"zeros"
,
this
.
zeros
+
""
);
properties
.
setProperty
(
prefix
+
"zeros"
,
this
.
zeros
+
""
);
properties
.
setProperty
(
prefix
+
"FFTSize"
,
this
.
FFTSize
+
""
);
properties
.
setProperty
(
prefix
+
"FFTSize"
,
this
.
FFTSize
+
""
);
properties
.
setProperty
(
prefix
+
"FFTSize_lwir"
,
this
.
FFTSize_lwir
+
""
);
properties
.
setProperty
(
prefix
+
"FFTSize_lwir"
,
this
.
FFTSize_lwir
+
""
);
properties
.
setProperty
(
prefix
+
"FFTOverlap"
,
this
.
FFTOverlap
+
""
);
properties
.
setProperty
(
prefix
+
"FFTOverlap_lwir"
,
this
.
FFTOverlap_lwir
+
""
);
properties
.
setProperty
(
prefix
+
"fftGaussWidth"
,
this
.
fftGaussWidth
+
""
);
properties
.
setProperty
(
prefix
+
"fftGaussWidth"
,
this
.
fftGaussWidth
+
""
);
properties
.
setProperty
(
prefix
+
"phaseCorrelationFraction"
,
this
.
phaseCorrelationFraction
+
""
);
properties
.
setProperty
(
prefix
+
"phaseCorrelationFraction"
,
this
.
phaseCorrelationFraction
+
""
);
properties
.
setProperty
(
prefix
+
"correlationHighPassSigma"
,
this
.
correlationHighPassSigma
+
""
);
properties
.
setProperty
(
prefix
+
"correlationHighPassSigma"
,
this
.
correlationHighPassSigma
+
""
);
...
@@ -11891,6 +11902,7 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11891,6 +11902,7 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
properties
.
setProperty
(
prefix
+
"legacyMode"
,
this
.
minUVSpan
+
""
);
properties
.
setProperty
(
prefix
+
"legacyMode"
,
this
.
minUVSpan
+
""
);
}
}
public
void
getProperties
(
String
prefix
,
Properties
properties
){
public
void
getProperties
(
String
prefix
,
Properties
properties
){
// EProperties properties = (EProperties) pproperties;
if
(
properties
.
getProperty
(
prefix
+
"correlationSize"
)!=
null
)
if
(
properties
.
getProperty
(
prefix
+
"correlationSize"
)!=
null
)
this
.
correlationSize
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"correlationSize"
));
this
.
correlationSize
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"correlationSize"
));
if
(
properties
.
getProperty
(
prefix
+
"correlationSizeLwir"
)!=
null
)
if
(
properties
.
getProperty
(
prefix
+
"correlationSizeLwir"
)!=
null
)
...
@@ -11905,6 +11917,17 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
...
@@ -11905,6 +11917,17 @@ error=Sum(W(x,y)*(F^2 + 2*F*(A*x^2+B*y^2+C*x*y+D*x+E*y-Z(x,y)) +(...) )
this
.
FFTSize
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTSize"
));
this
.
FFTSize
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTSize"
));
if
(
properties
.
getProperty
(
prefix
+
"FFTSize_lwir"
)!=
null
)
if
(
properties
.
getProperty
(
prefix
+
"FFTSize_lwir"
)!=
null
)
this
.
FFTSize_lwir
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTSize_lwir"
));
this
.
FFTSize_lwir
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTSize_lwir"
));
if
(
properties
.
getProperty
(
prefix
+
"FFTOverlap"
)!=
null
)
this
.
FFTOverlap
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTOverlap"
));
if
(
properties
.
getProperty
(
prefix
+
"FFTOverlap_lwir"
)!=
null
)
this
.
FFTOverlap_lwir
=
Integer
.
parseInt
(
properties
.
getProperty
(
prefix
+
"FFTOverlap_lwir"
));
// finally shortened :
// this.FFTOverlap= properties.getProperty(prefix+"FFTOverlap", this.FFTOverlap);
// this.FFTOverlap_lwir= properties.getProperty(prefix+"FFTOverlap_lwir",this.FFTOverlap_lwir);
if
(
properties
.
getProperty
(
prefix
+
"absoluteCorrelationGaussWidth"
)!=
null
)
if
(
properties
.
getProperty
(
prefix
+
"absoluteCorrelationGaussWidth"
)!=
null
)
this
.
absoluteCorrelationGaussWidth
=
Boolean
.
parseBoolean
(
properties
.
getProperty
(
prefix
+
"absoluteCorrelationGaussWidth"
));
this
.
absoluteCorrelationGaussWidth
=
Boolean
.
parseBoolean
(
properties
.
getProperty
(
prefix
+
"absoluteCorrelationGaussWidth"
));
if
(
properties
.
getProperty
(
prefix
+
"zeros"
)!=
null
)
if
(
properties
.
getProperty
(
prefix
+
"zeros"
)!=
null
)
...
...
src/main/java/com/elphel/imagej/common/EProperties.java
0 → 100644
View file @
a9f8afde
package
com
.
elphel
.
imagej
.
common
;
import
java.util.Properties
;
public
class
EProperties
extends
Properties
{
private
static
final
long
serialVersionUID
=
-
425120416815883045L
;
public
int
getProperty
(
String
key
,
int
value
){
return
Integer
.
parseInt
(
getProperty
(
key
,
""
+
value
));
}
public
double
getProperty
(
String
key
,
double
value
){
return
Double
.
parseDouble
(
getProperty
(
key
,
""
+
value
));
}
public
boolean
getProperty
(
String
key
,
boolean
value
){
return
Boolean
.
parseBoolean
(
getProperty
(
key
,
""
+
value
));
}
}
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment