Tag: numpy

Interpolation und Extrapolation für große Arrays

Ich habe ein großes Array y auf einem ungleichmäßigen, geordneten Raster x . Die Länge des Arrays ist typischerweise N ~ 2 ^ 14 bis N ~ 2 ^ 18. Ich möchte eine Spline-Interpolation (oder quadratisch) des Arrays bekommen. Das Problem, dem ich gegenüberstehe, ist, dass auch für niedrigere Werte von N die Interpolation sehr […]

Auswählen von Indizes für ein 2d-Array in numpy

Das funktioniert ganz gut in 1 Dimension: # This will sort bar by the order of the values in foo (Pdb) bar = np.array([1,2,3]) (Pdb) foo = np.array([5,4,6]) (Pdb) bar[np.argsort(foo)] array([2, 1, 3]) Aber wie mache ich das in zwei Dimensionen? Argsort arbeitet gut, aber die Auswahl funktioniert nicht mehr: (Pdb) foo = np.array([[5,4,6], [9,8,7]]) […]

Anpassung einer Verteilung mit dem Histogramm mit scipy

Ich möchte eine Verteilung mit scipy (in meinem Fall, mit weibull_min) zu meinen Daten passen. Ist es möglich, dies mit dem Histogramm zu tun, und nicht die Datenpunkte? In meinem Fall, weil das Histogramm Integer-Bins von Größe 1 hat, weiß ich, dass ich meine Daten auf folgende Weise extrapolieren kann: import numpy as np orig_hist […]

Numpy.polyfit gibt nützliche Passform, aber unendliche Kovarianzmatrix

Ich versuche, ein Polynom zu einem Satz von Daten zu passen. Manchmal kann es passieren, dass die Kovarianzmatrix, die von numpy.ployfit nur aus inf besteht, obwohl die Passung nützlich zu sein scheint. Es gibt keine numpy.inf oder 'numpy.nan' in den Daten! Beispiel: import numpy as np # sample data, does not contain really x**2-like behaviour, […]

-9999 als fehlender Wert mit numpy.genfromtxt ()

Lets sagen, ich habe eine dumme Textdatei mit dem Inhalt: Year Recon Observed 1505 162.38 23 1506 46.14 -9999 1507 147.49 -9999 -9999 wird verwendet, um einen fehlenden Wert zu bezeichnen (nicht fragen). Also, ich sollte in der Lage sein, dies in eine Numpy Array mit lesen: import numpy as np x = np.genfromtxt("file.txt", dtype […]

Wie kann ich diese Art von ASCII-Dateien effizient mit Python laden?

Ich habe große Fortran-generierte ASCII-Dateien formatiert wie folgt: xyz num_line index 1 float 2 float … num_line float x2 y2 z2 num_line2 index2 1 float 2 float … num_line2 float … Die Anzahl der Blöcke kann bis zu Tausenden und die Anzahl der Zeilen in jedem Block zu Hunderten gehen. Lassen Sie uns ein Beispiel […]

Wie zu entfernen, zufällig, Zeilen aus einem Dataframe aber aus jedem Label?

Dies ist für ein maschinelles Lernprojekt. Ich habe ein Dataframe mit 5 Spalten als Features und 1 Spalte als Label (Abbildung A). Ich möchte zufällig 2 Zeilen entfernen, aber von jedem Label. So, da es 12 Reihen gibt (4 für jedes Etikett); Ich werde am Ende mit 6 Reihen (2 von jedem Label) (Abbildung B). […]

3D-Extrapolation in Python (grundsätzlich scipy.griddata erweitert um extrapolieren)

Ich benutze die griddata Funktion in scipy, um 3 und 4 dimensionale Daten zu interpolieren. Es funktioniert wie ein Champion, außer dass es ein paar NaNs zurückgibt, weil einige der Punkte, die ich brauche, außerhalb des Bereichs der Eingangsdaten sind. Angesichts der Tatsache, dass Nd-Daten nur mit der "linearen" Modus-Interpolation sowieso arbeiten, sollte es ein […]

Wie kann ich einen Vektor auf eine bestimmte Länge mit numpy auflösen und / oder abschneiden?

Ich habe ein paar Listen: a = [1,2,3] b = [1,2,3,4,5,6] Die von variabler Länge sind. Ich möchte einen Vektor der Länge fünf zurückgeben, so dass, wenn die Eingangslistenlänge <5 ist, dann wird es mit Nullen auf der rechten Seite gefüllt, und wenn es> 5 ist, dann wird es am 5. Element abgeschnitten werden. Zum […]

Vektorisierung für Schleife mit wiederholten Indizes in Python

Ich versuche, ein Snippet zu optimieren, das viel gelernt wird (Millionen von Zeiten), so dass jede Art von Geschwindigkeitsverbesserung (hoffentlich das Entfernen der For-Loop) wäre großartig. Ich berechne eine Korrelationsfunktion von einigen j'th Partikel mit allen anderen C_j (| r-r '|) = sqrt (E ((s_j (r') – s_k (r)) ^ 2)) gemittelt über k. Meine […]

Python ist die beste Programmiersprache der Welt.