Skip to content
GitLab
Menu
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
Main Repositories
otb
Commits
25656a2f
Commit
25656a2f
authored
Jan 27, 2021
by
Julien Osman
Browse files
Merge branch '2040-RPC_model' of gitlab.orfeo-toolbox.org:orfeotoolbox/otb into 2040-RPC_model
parents
2c879376
3eb0a12b
Pipeline
#6638
failed with stages
in 89 minutes and 14 seconds
Changes
2
Pipelines
2
Hide whitespace changes
Inline
Side-by-side
Data/Baseline/OTB/Files/otbGCPsToRPCSensorModelImageFilterWithoutDEMOutput.txt
View file @
25656a2f
{"LineOffset": "2224.5", "SampleOffset": "2224.5", "LatOffset": "43.6089", "LonOffset": "1.42524", "HeightOffset": "0", "LineScale": "225.5", "SampleScale": "225.5", "LatScale": "1", "LonScale": "1", "HeightScale": "1", "LineNum": [ "0.0017
4891
", "0.4698
79
", "-182.989", "-
4.38776
e-1
5
", "0.0001
83514", "3.261
03e-23", "
1.880
11e-36", "-0.00013
3312
", "-0.0029
5
522", "0", "0", "4.137
24
e-05", "-2.6611
7
e-05", "0", "-0.0037833", "-0.0039722
9
", "0", "0", "0", "0", ], "LineDen": [ "1", "0.11
9555
", "-0.7
88492
", "0", "-0.001762
4
5", "0", "0", "-0.69207", "-0.7251
7
5", "0", "0", "4.
36024
e-06", "-1.2
1576
e-05", "0", "-1.6
0709
e-05", "-1.
09766
e-05", "0", "0", "0", "0", ], "SampleNum": [ "-4.0
2927
e-05", "140.123", "0.7205
52
", "-
6.37125
e-1
5
", "-0.00030
6972
", "
3.63979
e-25", "
7.44761
e-39", "7.06
041
e-05", "0.000393
491
", "0", "0", "0.0051405", "0.0016987
2
", "0", "3.047
58
e-06", "3.0597
6
e-05", "0", "0", "0", "0", ], "SampleDen": [ "1", "-0.01406
19
", "0.0561
058
", "0", "-0.0016542
6
", "0", "0", "-0.720
302
", "-0.23804
8
", "0", "0", "-2.11
371
e-07", "-2.52
599
e-06", "0", "1.287
11
e-06", "7.082
05
e-06", "0", "0", "0", "0", ], }
Residual ground error:
0.
12432
13845
{"LineOffset": "2224.5", "SampleOffset": "2224.5", "LatOffset": "43.6089", "LonOffset": "1.42524", "HeightOffset": "0", "LineScale": "225.5", "SampleScale": "225.5", "LatScale": "1", "LonScale": "1", "HeightScale": "1", "LineNum": [ "0.0017
6507
", "0.4698
26
", "-182.989", "-
3.97557
e-1
4
", "0.0001
67599", "-4.114
03e-23", "
4.642
11e-36", "-0.00013
2256
", "-0.0029
9
522", "0", "0", "4.137
19
e-05", "-2.6611
6
e-05", "0", "-0.003783
2
3", "-0.0039722", "0", "0", "0", "0", ], "LineDen": [ "1", "0.11
6648
", "-0.7
95797
", "0", "-0.0017625", "0", "0", "-0.6920
5
7", "-0.72515
9
", "0", "0", "4.
24205
e-06", "-1.2
217
e-05", "0", "-1.6
2292
e-05", "-1.
11751
e-05", "0", "0", "0", "0", ], "SampleNum": [ "-4.0
3023
e-05", "140.123", "0.7205
49
", "-
3.96237
e-1
4
", "-0.00030
7018
", "
-6.10744
e-25", "
3.08505
e-39", "7.06
292
e-05", "0.000393
857
", "0", "0", "0.005140
4
5", "0.0016987
1
", "0", "3.047
72
e-06", "3.0597
5
e-05", "0", "0", "0", "0", ], "SampleDen": [ "1", "-0.01406
52
", "0.0561
119
", "0", "-0.0016542
8
", "0", "0", "-0.720
295
", "-0.23804
6
", "0", "0", "-2.11
567
e-07", "-2.52
7
e-06", "0", "1.287
29
e-06", "7.082
13
e-06", "0", "0", "0", "0", ], }
Residual ground error: 1
.
2432
282311
Modules/Filtering/Projection/src/otbRPCSolver.cxx
View file @
25656a2f
...
...
@@ -184,7 +184,7 @@ void computeCoefficients(const std::vector<PrecisionType> & f,
// compute the residual
auto
residual
=
M
.
transpose
()
*
w2
*
(
M
*
coeffs
-
r
);
auto
residualValue
=
inner_product
(
residual
,
residual
);
res
=
inner_product
(
residual
,
residual
);
}
outCoeffs
.
assign
(
coeffs
.
begin
(),
coeffs
.
end
());
...
...
@@ -266,7 +266,7 @@ void RPCSolver::Solve(const GCPsContainerType& gcpContainer, PrecisionType& rmsE
{
const
auto
&
imagePoint
=
gcp
.
first
;
minc
=
std
::
min
(
imagePoint
[
0
],
min
l
);
minc
=
std
::
min
(
imagePoint
[
0
],
min
c
);
maxc
=
std
::
max
(
imagePoint
[
0
],
maxc
);
minl
=
std
::
min
(
imagePoint
[
1
],
minl
);
...
...
@@ -385,7 +385,7 @@ void RPCSolver::Solve(const GCPsContainerType& gcpContainer, PrecisionType& rmsE
+
(
gcp
.
first
[
1
]
-
outPoint
[
1
])
*
(
gcp
.
first
[
1
]
-
outPoint
[
1
]);
}
rmsError
=
std
::
sqrt
(
rmseAcc
)
/
numberOfPoints
;
rmsError
=
std
::
sqrt
(
rmseAcc
/
numberOfPoints
)
;
}
}
\ No newline at end of file
Write
Preview
Supports
Markdown
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