| @@ -687,7 +687,7 class SpectraProc(ProcessingUnit): | |||||
| 687 | for ip in range(new_cinterfid): |
|
687 | for ip in range(new_cinterfid): | |
| 688 | ind = junkspc_interf[:, new_interfid[ip]].ravel().argsort() |
|
688 | ind = junkspc_interf[:, new_interfid[ip]].ravel().argsort() | |
| 689 | jspc_interf[new_interfid[ip] |
|
689 | jspc_interf[new_interfid[ip] | |
| 690 | ] = junkspc_interf[ind[nhei_interf / 2], new_interfid[ip]] |
|
690 | ] = junkspc_interf[ind[nhei_interf // 2], new_interfid[ip]] | |
| 691 |
|
691 | |||
| 692 | jspectra[ich, :, ind_hei] = jspectra[ich, :, |
|
692 | jspectra[ich, :, ind_hei] = jspectra[ich, :, | |
| 693 | ind_hei] - jspc_interf # Corregir indices |
|
693 | ind_hei] - jspc_interf # Corregir indices | |
| @@ -747,17 +747,17 class SpectraProc(ProcessingUnit): | |||||
| 747 |
|
747 | |||
| 748 | ind = numpy.abs(jcspc_interf[mask_prof]).ravel().argsort() |
|
748 | ind = numpy.abs(jcspc_interf[mask_prof]).ravel().argsort() | |
| 749 |
|
749 | |||
| 750 | median_real = numpy.median(numpy.real( |
|
750 | median_real = int(numpy.median(numpy.real( | |
| 751 | junkcspc_interf[mask_prof[ind[list(range(3 * num_prof / 4))]], :])) |
|
751 | junkcspc_interf[mask_prof[ind[list(range(3 * num_prof // 4))]], :]))) | |
| 752 | median_imag = numpy.median(numpy.imag( |
|
752 | median_imag = int(numpy.median(numpy.imag( | |
| 753 | junkcspc_interf[mask_prof[ind[list(range(3 * num_prof / 4))]], :])) |
|
753 | junkcspc_interf[mask_prof[ind[list(range(3 * num_prof // 4))]], :]))) | |
|
|
754 | comp_mask_prof = [int(e) for e in comp_mask_prof] | |||
| 754 | junkcspc_interf[comp_mask_prof, :] = numpy.complex( |
|
755 | junkcspc_interf[comp_mask_prof, :] = numpy.complex( | |
| 755 | median_real, median_imag) |
|
756 | median_real, median_imag) | |
| 756 |
|
757 | |||
| 757 | for iprof in range(num_prof): |
|
758 | for iprof in range(num_prof): | |
| 758 | ind = numpy.abs(junkcspc_interf[iprof, :]).ravel().argsort() |
|
759 | ind = numpy.abs(junkcspc_interf[iprof, :]).ravel().argsort() | |
| 759 | jcspc_interf[iprof] = junkcspc_interf[iprof, |
|
760 | jcspc_interf[iprof] = junkcspc_interf[iprof, ind[nhei_interf // 2]] | |
| 760 | ind[nhei_interf / 2]] |
|
|||
| 761 |
|
761 | |||
| 762 | # Removiendo la Interferencia |
|
762 | # Removiendo la Interferencia | |
| 763 | jcspectra[ip, :, ind_hei] = jcspectra[ip, |
|
763 | jcspectra[ip, :, ind_hei] = jcspectra[ip, | |
General Comments 0
You need to be logged in to leave comments.
Login now