| @@ -616,8 +616,8 class SpectraProc(ProcessingUnit): | |||
|
|
616 | 616 | # hei_interf |
|
|
617 | 617 | if hei_interf is None: |
|
|
618 | 618 | count_hei = num_hei / 2 # Como es entero no importa |
|
|
619 | hei_interf = numpy.asmatrix(list(range(count_hei))) + num_hei - count_hei | |
|
|
620 | hei_interf = numpy.asarray(hei_interf)[0] | |
|
|
619 | hei_interf = numpy.asmatrix(list(range(int(count_hei)))) + num_hei - count_hei | |
|
|
620 | hei_interf = numpy.asarray(hei_interf)[0].astype('int') | |
|
|
621 | 621 | # nhei_interf |
|
|
622 | 622 | if (nhei_interf == None): |
|
|
623 | 623 | nhei_interf = 5 |
| @@ -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 |
| @@ -728,7 +728,7 class SpectraProc(ProcessingUnit): | |||
|
|
728 | 728 | # Remocion de Interferencia en el Cross Spectra |
|
|
729 | 729 | if jcspectra is None: |
|
|
730 | 730 | return jspectra, jcspectra |
|
|
731 | num_pairs = jcspectra.size / (num_prof * num_hei) | |
|
|
731 | num_pairs = int(jcspectra.size / (num_prof * num_hei)) | |
|
|
732 | 732 | jcspectra = jcspectra.reshape(num_pairs, num_prof, num_hei) |
|
|
733 | 733 | |
|
|
734 | 734 | for ip in range(num_pairs): |
| @@ -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