diff --git a/Code/Radiometry/otbSarRadiometricCalibrationToImageFilter.txx b/Code/Radiometry/otbSarRadiometricCalibrationToImageFilter.txx index 211b39e05e524dac280ba32a518da5d9458f2360..95a8f4ebeb08a914bfcdc5d45a288858297f13a9 100644 --- a/Code/Radiometry/otbSarRadiometricCalibrationToImageFilter.txx +++ b/Code/Radiometry/otbSarRadiometricCalibrationToImageFilter.txx @@ -76,11 +76,23 @@ SarRadiometricCalibrationToImageFilter<TInputImage, TOutputImage> std::cout<<"SarRadiometricCalibrationToImageFilter starts trouble"<<std::endl; incidenceAngle = function->GetIncidenceAngle(); - std::cout<<"function->GetIncidenceAngle():"<<std::endl; - std::cout<<function->GetIncidenceAngle()<<std::endl; incidenceAngle->SetPointSet(imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()); std::cout<<"imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle():"<<std::endl; - std::cout<<imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()<<std::endl; + typename ParametricFunctionType::PointType point; + point.Fill(0); + typename ParametricFunctionType::PointSetType::PixelType pointValue; + pointValue = itk::NumericTraits<ParametricFunctionType::PointSetType::PixelType>::Zero; + unsigned int nbRecords = imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()->GetNumberOfPoints(); + + // Fill the linear system + for (unsigned int i = 0; i < nbRecords; ++i) + { + imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()->GetPoint(i, &point); + imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()->GetPointData(i, &pointValue); + //std::cout << "point = " << point << std::endl; + std::cout << pointValue << " / "<<point<<std::endl; + } + //std::cout<<imageMetadataInterface->GetRadiometricCalibrationIncidenceAngle()<<std::endl; incidenceAngle->SetPolynomalSize(imageMetadataInterface->GetRadiometricCalibrationIncidenceAnglePolynomialDegree()); std::cout<<"imageMetadataInterface->GetRadiometricCalibrationIncidenceAnglePolynomialDegree():"<<std::endl; std::cout<<imageMetadataInterface->GetRadiometricCalibrationIncidenceAnglePolynomialDegree()<<std::endl;