Skip to content

Commit 575ef93

Browse files
authored
Small updates to applyRtc() (#82)
* disable polarimetric symmetrization by default * revert changes to `symmetrize_cross_pol_channels` * Update GCOV and GSLC specification XMLs * Revert changes to the GCOV and GSLC specification XMLs * consider valid the pixels with zero values * pass `rtc_min_value` instead of `rtc_min_value_db` to `_applyRtc()`
1 parent bf6f02c commit 575ef93

File tree

1 file changed

+4
-10
lines changed

1 file changed

+4
-10
lines changed

cxx/isce3/geometry/RTC.cpp

Lines changed: 4 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -73,12 +73,6 @@ void _clip_min_max(std::complex<T>& radar_value, float clip_min, float clip_max)
7373
*/
7474
using T_real = typename isce3::real<T>::type;
7575

76-
// no data (complex)
77-
if (std::abs(radar_value) == 0) {
78-
radar_value *= std::numeric_limits<T_real>::quiet_NaN();
79-
return;
80-
}
81-
8276
// clip min (complex)
8377
if (!std::isnan(clip_min) && std::abs(radar_value) < clip_min)
8478
// update magnitude without changing the phase
@@ -375,20 +369,20 @@ void applyRtc(const isce3::product::RadarGridParameters& radar_grid,
375369
if (input_raster.dtype() == GDT_Float32 ||
376370
(input_raster.dtype() == GDT_CFloat32 && flag_complex_to_real))
377371
_applyRtc<float>(input_raster, *rtc_raster, output_raster,
378-
rtc_min_value_db, abs_cal_factor, clip_min, clip_max, info,
372+
rtc_min_value, abs_cal_factor, clip_min, clip_max, info,
379373
flag_complex_to_real);
380374
else if (input_raster.dtype() == GDT_Float64 ||
381375
(input_raster.dtype() == GDT_CFloat64 && flag_complex_to_real))
382376
_applyRtc<double>(input_raster, *rtc_raster, output_raster,
383-
rtc_min_value_db, abs_cal_factor, clip_min, clip_max, info,
377+
rtc_min_value, abs_cal_factor, clip_min, clip_max, info,
384378
flag_complex_to_real);
385379
else if (input_raster.dtype() == GDT_CFloat32)
386380
_applyRtc<std::complex<float>>(input_raster, *rtc_raster, output_raster,
387-
rtc_min_value_db, abs_cal_factor, clip_min, clip_max, info,
381+
rtc_min_value, abs_cal_factor, clip_min, clip_max, info,
388382
flag_complex_to_real);
389383
else if (input_raster.dtype() == GDT_CFloat64)
390384
_applyRtc<std::complex<double>>(input_raster, *rtc_raster,
391-
output_raster, rtc_min_value_db, abs_cal_factor, clip_min,
385+
output_raster, rtc_min_value, abs_cal_factor, clip_min,
392386
clip_max, info, flag_complex_to_real);
393387
else {
394388
std::string error_message =

0 commit comments

Comments
 (0)