scale = 3e8/(interpolated_x*interpolated_x)
scale_matrix = scale*np.ones(spec.shape)
spec_scaled = scale_matrix*spec
interpolated_spec = scipy.interpolate.interp1d(well, spec_scaled, 'linear', axis=1)(interpolated_x)[/code]
ohne Probleme.
Wie sieht die komplette Fehlermeldung inklusive Traceback aus?
`scale_matrix` ist überflüssig, es geht einfach `spec_scaled = scale * spec`
Habe das in scale * spec geändert. Dachte anfangs dass die Skalierung vor der Interpolierung evt. effizienter sei.
Du hast Recht, es funktioniert so. Problem jetzt ist der Skalierungsfaktor.
spec ist das Ergebnis einer F.T. von Energiewerten über der Zeit in Energiewerten über Frequenz ... und die sollen ja in Energiewerte über Wellenlänge umgerechntet werden.
D.h.
(i) äquidistantes Wellenlängeninterval erzeugen - ist erledigt
(ii) Energiewerte skallieren - Skalierung stimmt nicht
(iii) Energiewerte auf Intervall (i) interpolieren - Code ist da
Aber soweit funktioniert es.
Vielen Dank für die Hilfe