Unable to open sensor products (with RPC) in Monteverdi
Mantis Issue 1136, reported by mgrizonnet, assigned to salbert, created: 2016-01-27
I am not able to load any Level 1A products with Monteverdi (QB, pleiades...)
There is the "classic" warning message in ossimRpcModel::lineSampleHeightToWorld
ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate. WARNG> "Ce format de fichier permet d'utiliser les pyramides multi-résolution." WARNING ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate. WARNING ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate. WARNING ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate. WARNING ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate. WARNING ossimRpcModel::lineSampleHeightToWorld: Max number of iterations reached in ground point solution. Results are inaccurate.
I've tried on 2 different platforms with otb/ice/monteverdi develop branches (27/02) and also reproduce the issue with monteverdi 3.0.0 package from ubuntugis-unstable repository (both configurations use ossim library 1.8.20 provided in ubuntugis-unstable.
I'll be glad that the issue come from my local config but I can't find my mistake for now. Does anyone can try to reproduce?
Even stranger, I suspect that the issue is not related with ossim as I'm able to open the same images with otbicewiewer without problem. Note that otbiceviewer use the same ossim library from ubuntugis-unstable.
Does anyone else
Note that changing the default "Proj" of the viewer from "Sensor" (default value) to "None" solve the issue!
1453978009 - salbertIssue seems to be related to issue #1122 (closed) (https://bugs.orfeo-toolbox.org/view.php?id=1122). The "./otbcli_ConvertSensorToGeoPoint -in ~/OTB-Data/mv2-test/QB_TLSE_512x512_R_NIR.tif -input.idx -200 -input.idy 128" test could be made to check it (see #1122 (closed)).
OSSIM 1.8.20-3 has been used to solve this issue.
ossimRpcModel::lineSampleHeightToWorld() is called when transforming a point from one projection to another. If Projection mode is set to None, no transform will occur.
You could also try to reproduce with otbIceViewer by loading QUICKIRD/TOULOUSE and move until x and y would be negative. Then, the problem is reproduced if a NaN is displayed in Ice's HD.
1453987580 - julienCan you try with
$ LC_NUMERIC=C monteverdi
It looks like a locale problem to me. I always launch monteverdi that way.
1453991464 - grizonnetmLC_NUMERIC=C monteverdi solves the issue. Thanks.
I confirm that there is no relationship here with bug 0001122 as ConvertSensorToGeoPoint application