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