Idl Moving Average Filter

Der Moving Average als Filter Der gleitende Durchschnitt wird oft für die Glättung von Daten in Anwesenheit von Rauschen verwendet. Der einfache gleitende Durchschnitt wird nicht immer als der Finite Impulse Response (FIR) - Filter erkannt, der es ist, während er tatsächlich einer der gebräuchlichsten Filter in der Signalverarbeitung ist. Wenn man sie als Filter betrachtet, kann man sie beispielsweise mit gefensterten Filtern vergleichen (siehe Artikel zu Tiefpaß-, Hochpass - und Bandpass - und Bandsperrfiltern für Beispiele). Der Hauptunterschied zu diesen Filtern besteht darin, daß der gleitende Durchschnitt für Signale geeignet ist, für die die Nutzinformation im Zeitbereich enthalten ist. Von denen Glättungsmessungen durch Mittelung ein Paradebeispiel sind. Window-sinc-Filter, auf der anderen Seite, sind starke Künstler im Frequenzbereich. Mit Ausgleich in der Audioverarbeitung als typisches Beispiel. Es gibt einen detaillierteren Vergleich beider Arten von Filtern in Time Domain vs. Frequency Domain Performance von Filtern. Wenn Sie Daten haben, für die sowohl die Zeit als auch die Frequenzdomäne wichtig sind, dann möchten Sie vielleicht einen Blick auf Variationen auf den Moving Average werfen. Die eine Anzahl gewichteter Versionen des gleitenden Durchschnitts zeigt, die besser sind. Der gleitende Durchschnitt der Länge (N) kann so definiert werden, wie er üblicherweise implementiert ist, wobei der aktuelle Ausgabeabtastwert der Durchschnitt der vorhergehenden (N) Abtastwerte ist. Als Filter betrachtet, führt der gleitende Durchschnitt eine Faltung der Eingangsfolge (xn) mit einem rechteckigen Puls der Länge (N) und der Höhe (1N) durch (um den Bereich des Pulses und damit die Verstärkung des Filters zu erzeugen , eins ). In der Praxis ist es am besten, (N) ungerade zu nehmen. Obwohl ein gleitender Durchschnitt auch unter Verwendung einer geraden Anzahl von Abtastwerten berechnet werden kann, hat die Verwendung eines ungeradzahligen Wertes für (N) den Vorteil, daß die Verzögerung des Filters eine ganzzahlige Anzahl von Abtastwerten ist, da die Verzögerung eines Filters mit (N) Proben genau ((N-1) 2). Der gleitende Durchschnitt kann dann exakt mit den ursprünglichen Daten ausgerichtet werden, indem er um eine ganze Zahl von Abtastwerten verschoben wird. Zeitdomäne Da der gleitende Durchschnitt eine Faltung mit einem rechteckigen Puls ist, ist sein Frequenzgang eine sinc-Funktion. Dies macht es ähnlich dem Dual des Fenstersynchronfilters, da es sich hierbei um eine Faltung mit einem Sinc-Puls handelt, der zu einem rechteckigen Frequenzgang führt. Es ist diese sinc Frequenzantwort, die den gleitenden Durchschnitt ein schlechter Darsteller im Frequenzbereich macht. Allerdings führt es sehr gut im Zeitbereich. Daher ist es perfekt, um Daten zu löschen, um Rauschen zu entfernen, während gleichzeitig eine schnelle Sprungantwort beibehalten wird (1). Für das typische Additiv-Weiß-Gauß-Rauschen (AWGN), das oft angenommen wird, bewirkt die Mittelung (N) - Proben, dass das SNR um einen Faktor (sqrt N) erhöht wird. Da das Rauschen für die einzelnen Proben unkorreliert ist, gibt es keinen Grund, jede Probe unterschiedlich zu behandeln. Daher wird der gleitende Durchschnitt, der jeder Probe das gleiche Gewicht gibt, die maximale Menge an Rauschen für eine gegebene Sprungantwortschärfe beseitigen. Implementierung Da es sich um ein FIR-Filter handelt, kann der gleitende Durchschnitt durch Faltung implementiert werden. Es hat dann die gleiche Effizienz (oder das Fehlen davon) wie jedes andere FIR-Filter. Sie kann aber auch rekursiv und effizient umgesetzt werden. Es folgt unmittelbar aus der Definition, daß diese Formel das Ergebnis der Ausdrücke für (yn) und (yn1) ist, dh, daß die Veränderung zwischen (yn1) und (yn) ein zusätzlicher Term (xn1N) ist Das Ende, während der Term (xn-N1N) von Anfang entfernt wird. In praktischen Anwendungen ist es oft möglich, die Division durch (N) für jeden Term auszulassen, indem die resultierende Verstärkung von (N) an einer anderen Stelle kompensiert wird. Diese rekursive Umsetzung wird viel schneller als Faltung. Jeder neue Wert von (y) kann mit nur zwei Additionen anstelle der (N) Additionen berechnet werden, die für eine einfache Implementierung der Definition erforderlich wären. Eine Sache, mit der Sie nach einer rekursiven Implementierung Ausschau halten, ist, dass Rundungsfehler akkumulieren. Dies kann ein Problem für Ihre Anwendung sein oder auch nicht, aber es bedeutet auch, dass diese rekursive Implementierung tatsächlich mit einer Integer-Implementierung besser funktionieren wird als mit Gleitkommazahlen. Dies ist sehr ungewöhnlich, da eine Gleitkomma-Implementierung gewöhnlich einfacher ist. Der Schluss davon muss sein, dass Sie die Nützlichkeit des einfachen gleitenden Durchschnittsfilters in Signalverarbeitungsanwendungen nie unterschätzen sollten. Filter Design Tool Dieser Artikel wird mit einem Filter Design Tool ergänzt. Experimentiere mit verschiedenen Werten für (N) und visualisiere die resultierenden Filter. Versuchen Sie es jetztMarshall Perrins IDL-Routinen Liste Letzte Aktualisierung: Mo Sep 17 17:04:02 2012. Routinen nach Kategorie Astronomische Utility FINDINGCHART - eine Suche nach einem Objekt auf dem Bildschirm oder PDF. LASERTIMES - Berechnen Sie die Start - und Stoppzeiten für den Lick Laser in LST. Diese entsprechen 11 Uhr und 5 Uhr Ortszeit (entweder PST oder PDT) LSTNOW - Rücklaufstrom lst, für ein gegebenes Observatorium MAGERRPROP - Wandelt Flüsse in Größenordnungen mit Fehlerfortpflanzung um. MULTIFINDINGCHART - rufen Sie findingchart für eine Reihe von Quellen in einer Datei aufgeführt. SUNTIMES - Berechnen Sie APPROXIMATE Sonnenauf - und - untergang mal ZD2AIRMASS - Berechnen Sie die Luftmasse gegen den Zenitabstand. Data Cubes CUBEDISPGRIDSPEC - Overplot ein Gitter von Spektren auf ein Graustufenbild eines Würfels Dies wurde von einer Figur in (ich glaube) einer von Tracey Becks Papiere inspiriert. CUBECOLLAPSE - Einen Würfel entlang der Wellenlängenachse kollabieren. Der Standardwert ist, den Cube zu durchschnittlich, aber es gibt Optionen für Median oder Summe kombinieren stattdessen. Es gibt einige Unterstützung für die Verwendung von Qualitätsflaggen, um schlechte Pixel anzuzeigen, die beim Zusammenbrechen überspringen. CUBEGETWAVEAXIS - Rückgabe der Achse eines Würfels als Wellenlänge CUBEREAD - Einlesen eines Datenwürfels in eine standardisierte Struktur CUBEREBIN - Rebin einen Würfel in den X - und Y-Dimensionen, wobei die Z-Dimension unverändert bleibt. CUBESUBTRACTCONTINUUM - Grobe Kontinuums-Subtraktion für IFS-Datacubes Fast jeder Algorithmus ist besser als dieser EXTAST3 - Extrahieren Sie ASTrometry-Parameter aus einem FITS-Image-Header. (3D-Version für Datacubes) PUTAST3 - Fügt WCS astrometry-Parameter in einen gegebenen FITS-Header. (3D Version für Datacubes) XYZ2ADL - Berechnen R. A. Dec und Wellenlänge von X, Y und Z und einer FITS-Astrometriestruktur ISMYSQLPRESENT - Gibt es hier eine mysql-Datenbank? Diese Routine ermöglicht es Ihnen, den Code entweder mit oder ohne mysql, abhängig von der lokalen Maschinenfunktionalität, auszuführen. MYSQLCHECK - überprüfen Sie, ob eine Verbindung zu mysql geöffnet offen ist, wenn es nicht MYSQLQUERY2 - Submit MySQL-Abfrage und erhalten Antwort als ein Array von Strukturen. QUERY2MASS - Wenn Sie einen Namen angeben, fragen Sie den 2MASS-Katalog auf Vizier nach seinen J-, H - und Ks-Größen (und deren Fehlern) ab. QUERYIRAS - Wenn Sie einen Namen haben, fragen Sie den IRAS-Katalog auf Vizier nach seinen J-, H - und Ks-Größen (und ihren Fehlern) ab. QUERYUSNOB - APRINT - drucken Sie etwas aus, das als IDL-Array formatiert ist (dh in der IDL-Quellcode-Syntax). FITSLOADER - Lädt einen bestimmten Satz von FITS-Dateien, sucht nach geeigneten Bias - und Flat-Field-Dateien für sie, wendet diese an und wählt optional die Skalierung Ganz auf eine feste Belichtungszeit eingestellt. HDRCONCAT - Verkettet FITS-Header in ein 2D-String-Array, wobei intelligente Dinge mit Achsengrößen bei Bedarf durchgeführt werden. HDRCOPY - kopiert Schlüssel von einem Header zu einem anderen PICKFITS - wähle eine oder mehrere passt Dateien mit der GUI und dann laden. READTEXT - eine Textdatei vollständig in ein Array einlesen SELECTFITS - Der Zweck dieses Programms besteht darin, dem Benutzer die Auswahl einer FITS-Bilddatei zum Lesen zu ermöglichen. Die Bilddaten werden als Ergebnis der Funktion zurückgegeben. Das beste Merkmal dieses Programms ist die Möglichkeit, das Bild zu durchsuchen, bevor Sie es lesen. SIGDISP - Wie die IRCAL sigdisp, macht eine grobe n-Sigma lineare Strecke eines Bildes. SIXTYSTRING - Wie die Goddard-Bibliotheken sechzig, aber es gibt Strings des Formats DD: MM: SS zurück. Dies ist die Umkehrung des Zehnfadens. STRC - Formatiert etwas als String, entfernt die zusätzlichen Leerzeichen und unnötige Nullen. Funktioniert auf Skalaren oder Arrays. SXPARARR - Wie sxparr, aber für eine Reihe von passt Kopfzeilen. TENSTRING - Tut das gleiche wie die goddard idl astro bibliothek zehn prozedur (konvertiert sexagesimal coords in dezimal), sondern arbeitet auf string argumente der form DD: MM: SS oder DD MM SS. TEXPRINT - drucken Sie ein IDL-Array in der LaTeX-Tabellensyntax aus, das zum Einfügen in Ihr Papier geeignet ist. Bildanzeige ALOGSCALE - intelligent logarithmisch skalieren ein Bild für die Anzeige. ASINHSCLCOLOR - RGB-Asinh-Skalierung nach der Art von Lupton et al. ASINKSCLCONTOURS - schaffen Konturfarben für Overplotting. ASINHDEMO - zeigt die Verwendung von Asinh-Skalierungsroutinen Bildverarbeitung AORADNORM - Wenn Sie ein Bild (vermutlich ein AO-Bild) erhalten, erstellen Sie ein unscharfes maskiertes Bild. Dann ein Radialgeräuschprofil für sie berechnen und durch das Profil teilen. APPHOT - Programm zum Durchführen einer Kreisblende - Photometrie auf einer bildgebenden Blende, Größe und Innen - und Außenradius für den Himmelsubtraktionsring. AUTOREGISTER - Automatische Registrierung einer Reihe von Bildern. AVGMED - Median durchschnittlich ein Stapel von FITS-Bildern, mit optionaler medianer Gewichtung. BIAS - schnelles Skript zum Erstellen eines Bias-Frames und Anzeige von nützlichen Informationen CALCPHOTNOISE - Geben Sie ein Bild und eine Abbildung der Belichtungszeiten für dieses Bild an, berechnet das Photonenrauschen als Funktion der Position. CENPSFSUB - zentrierende psf Subtraktionsroutine. CROPBADMASK - Gegeben ein Bild und eine Maske, zuschneiden Sie das Bild auf die minimale Größe, die die Maske enthält. FFTREBIN - Rebins ein Bild, indem nötigenfalls Nullen zwischen den FFT-Komponenten hinzugefügt werden. FFTSHIFT - verschiebt ein Bild durch (dx, dy) Pixel unter Verwendung von Fourier-Transformationen. FFTSHIFTCUBE - FFTFILT - Hoch - oder Tiefpassfilter ein Bild im Fourierraum. FINDMAXSTAR - bei einem Bild findet man die x, y-Koordinaten des hellsten Gegenstands. FITPLANE - passt auf ein Flugzeug zu einem Bild FIXNANFITS - mit dem Namen einer fits-Datei öffnen, öffnen Sie die Datei, fixieren Sie alle NaN-Pixel in ihr durch Interpolation von den nächsten Nachbarn, dann speichern Sie die Datei FIXNANS - fix NANs in einem Bild, z. B Vor Kreuzkorrelation oder so. FIXPIX - gegeben ein Bild oder Stapel von Bildern und eine schlechte Pixelmaske, wird in schlechte Pixel füllen, indem die NPIX am nächsten guten Pixel, werfen die höchsten und niedrigsten der Haufen, und dann arithmetisch durchschnittlich. FWCENTROID - Implementiert den robusten Schwebefenster-Schwerpunkt-Algorithmus, der für JWST-Zielerfassungen verwendet wird. HIMCUT - Wie IMCUT, schneiden Sie eine Subregion aus einem Bild, und aktualisieren FITS-Header. Sehr ähnlich zu HEXTRACT, aber mit unterschiedlichen Argumenten IMAGESHIFTMASKEDGES - Beim Verschieben eines Bildes manipulieren Sie die fehlerhafte Pixelmaske so, dass die korrekten Kanten des verschobenen Bildes als schlecht erkannt werden. IMCENTERF - Programm, um den Massenmittelpunkt eines Bildes um den Punkt (x, y) zu berechnen, geben Sie die Antwort in (xcen, yc) zurück. IMCUT - Funktion zum Ausschneiden eines quadratischen Unterabschnitts aus einem Bild, wobei der Benutzer die Mitte des Unterabschnitts spezifiziert. INDICES - wie Pythons-Indices-Befehl, gibt Koordinatenindizes für ein Array zurück MATCHFILTER - Faltet ein Bild mit einer normalisierten Version von sich selbst zusammen. Schnell und dreckig. MATRIXDFT - Matrix Diskrete Fourier Transformation, nach Soummer et al. Dies ist nicht so schnell wie eine FFT, sondern erlaubt Ihnen, willkürlich wählen Sie die Stichprobe und Bereich in der Fourier-Domäne abgedeckt. MEANS - gibt die Mittel für einen Stapel von Bildern zurück MEDARR2 - Wie medarr, außer du kannst einen Satz von Belichtungszeiten angeben, um die Bilder zu skalieren, bevor sie ihre Mediane nehmen. MKDOMEFLAT - macht eine Kuppel flach durch Laden einer Reihe von Bildern, summiert sie, und dann Flattening. Gegebenenfalls geht es um Polarisierungsprobleme. MKSKYFLAT - Median-Kombinieren Sie eine Reihe von himmelblauen Bildern in einem Master Himmel flach. MKTWIFLAT - bildet ein flaches Feld aus einer Reihe von Bildern des Dämmerungshimmels unter Verwendung einer iterativen Anpassungstechnik für jedes Pixel, um das relative Flachfeld zu messen und jede konstante (d. h. thermische) Komponente zu berücksichtigen. Die für die Herstellung von Infrarot-Twiflaten verwendet werden sollen. MOSF - Routine zu verschieben und Stapel einer Reihe von Bildern, um ein endgültiges Mosaik zu schaffen, mit Masken, um schlechte Pixel oder kosmische Strahlen auszuschließen NEWSKYFLAT - Median-kombinieren eine Reihe von himmelblauen Bildern in einem Master-Himmel flach. OPTPSFSUB - Optimierung der PSF-Subtraktionsroutine. OPTSUB - PSF Subtraktionsoptimierer. Wenn zwei Bilder bereits ausgerichtet, was Skalierungsfaktor optimiert ihre Subtraktion OPTSUBWITHOFFSET - PSF Subtraktion Optimierer. Wenn zwei Bilder bereits ausgerichtet sind, was Skalenfaktor optimiert ihre Subtraktion Diese Version löst auch für eine konstante Verschiebung zwischen den beiden Bildern. RADGEN - Programm zur Erzeugung eines 2-d-Bildes mit 1-d-Profil (radial oder elliptisch) mit einer bestimmten Mitte und radialen Ausdehnung. RADNOISE - Radiallüfterkennungsprofil RADPLOTF - Programm zur Berechnung des radialelliptischen Profils eines Bildes bei gegebener Blende, Größenbereich und Innen - und Außenradius für Himmelsubtraktionsring. RECENTER - eine Reihe von Bildern und anfängliche Vermutungen zur Registrierung berechnen endgültige Offsets mit IMCENTERF REGISTERDEMO - Demonstrieren Sie, wie subreg. pro und mosf. pro zu registrieren und Mosaikbilder zu verwenden SETSKY - Addieren Sie Offsets zu Bildern, um den Hintergrundhimmelspegel gleich zu machen Der Hintergrund des ersten Bildes. SKYSUB - Subtrahieren eines 2D-Himmelsrahmens (oder eines Biasrahmens) aus einem 3D-Array von Bildern standardmäßig, es werden die Eingabebilder gelöscht, außer wenn STAT gespeichert wird - gibt Median-, Mittel-, Min-, Max - und Std-Abweichung für ein gegebenes Bild STDDEVARR - berechnen Die Standardabweichung für jedes Pixel in einem Würfel STDDEVS - Rückgabe der Standardabweichungen für einen Stapel von Bildern SUBMEDIAN - Subtrahieren des Medians von jedem Slice eines Datenwürfels SUBREG - Subpixelregistrierung von Bildern SUBREGSHIFTSTOPEAKS - Angesichts eines durch subreg erzeugten Schichtarrays, konvertiere dies Auf ein Array von Spitzenpixelpositionen in einem Bild, wie es für die Eingabe in die Mosf-Funktion benötigt wird. TUNEREGISTER - Bildregistrierung testen. Nimmt einen bereits registrierten Stapel von Bildern (eventuell von Hand gemacht), dann schneidet eine Subregion in es und läuft subreg auf dieser Subregion. IRCAL Pipeline GETIRCALFILTER - Angenommen, ein FITS-Header gibt einen String zurück, der den aktuellen IRCAL-Filter beschreibt. IRCALADDWCS - Hinzufügen von WCS-Koordinaten zu einem FITS-Header für IRCAL IRCALSTREHL - Berechnen Sie Strehl für ein IRCAL-Bild. IRCALBADPIXELS - Markieren Sie anhaltend schlechte ircal Pixel als NANs Verwendung: ircalbadpixels, imgs Dies ist nicht all die schlechten Pixel, nur diejenigen, die ich sagen kann, sind schlecht, aber die automatisierten Sachen nicht abholen aus irgendeinem Grund. IRCALBADPIXFROMLIST - Lesen Sie schlechte Pixel-Liste aus Datei und gelten für einen Bildstapel IRCALDEGHOST - entfernen Sie lästige negative Ghosts, die durch Kanal-Crosstalk verursacht werden. IRCALDEWARP - Verzerrung aus einem IRCAL-Bild entfernen IRCALDEWARPSHIFTS - Konvertieren einer Liste von Bildverschiebungen von rohen IRCAL-Koordinaten zu gedämpften Koordinaten IRCALFIXPIX - Dies ist ein Frontend zu nsfixpix. pro, das die IRCAL-Fehler-Ecken kennt und diese in der schlechten Pixelmaske bewahrt . IRCALGETJ2000 - Die IRCAL-Koordinaten der aktuellen Epoche in J2000-Koordinaten umwandeln IRCALREGRID - IRCAL-Bilder wiederholen, um die anamorphotische Vergrößerung zu kompensieren IRCALSATMASK - Bestimmen der gesättigten Regionen in IRCAL-Daten und der Maske an NaN Behandelt IRCAM - und NIRC2-Bilder. IRCALZEROPT - Frontend-Routine zur photometrischen Kalibrierung von IRCAL-Daten REDIRCALSKY - Version des IRCAL-Codes zum Erstellen von Himmelsdateien zur Eingabe von redircalsub. REDIRCALSUB - Reduzierung der IRCAL-Daten mit Subpixel-Registrierung und Mosaik. Mathematik GAUSSIAN2D - Berechnen Sie die 2-d-Gaußsche Funktion und optional das Derivat auf der Grundlage von Goddard IDL Astros Gaussian. pro RESTRICTION: Im Moment nur kreisförmig symmetrische Gaussianer GRIDEVAL - Wertet eine Funktion aus, die als String für alle Punkte eines gelieferten Rasters in x geliefert wird Und y. MAXES - die Maximalwerte jedes Bildes für einen Stapel von Bildern zurückgeben MODE - Berechnen Sie den Modus (häufigste Element in einem Array) optional mit Binning MONTECARLOMEAN - verwendet einen Satz von Zahlen, verwendet monte-carlo Methoden, um sowohl den Mittelwert als auch den Standard zu bestimmen Abweichung des Mittelwertes. MONTECARLOMEDIAN - verwendet einen Satz von Zahlen, verwendet Monte-Carlo-Methoden, um sowohl den Median und die Standardabweichung des Median zu bestimmen. POISS - Berechnet die Poisson-Verteilung als Funktion von X und M. PRODUCTERRPROP - Implementiert die Fehlerausbreitungsformel für das Produkt aus zwei Zahlen, von denen angenommen wird, dass sie keine Korrelation haben ATV - Interaktive Anzeige von 2D - oder 3D-Bildern. ATVMAKEMOVIE - erstellen Sie einen Film aus einem Bildstapel im ATV. FINDLOCALMAX - Finden Sie lokales Maximum in der Nähe einer gegebenen Position im Bild FINDCLOSEST - Gegeben ein Array, finden Sie den Index, dessen Wert in diesem Array am nächsten zu einer bestimmten Zahl ist. FINDFWXM - Finden Sie die volle Breite bei Xmax. D. h. X 0,5 für FWHM. MRECENTER - Finden Sie genaues Zentrum für einen Stern, wie RECENTER, aber mit MPFITPEAK, um die Spitze zu tun. WHEREIS - Wenn der 1-d-Index eines Pixels in einem Array angegeben ist, werden die x - und y-Koordinaten, die diesem Pixel entsprechen, zurückgegeben. WHEREISMAX - Angeben eines Arrays, gibt Position und Wert des maximalen Pixels zurück ALLSKYPLOT - Plot einige Punkte auf dem ganzen Himmel ARROWS2 - überblendet einen Pfeil auf einem Bild an einer gegebenen Position und einem Winkel. DRAWSCALEBAR - Zeichnen einer Skalierungsleiste auf einem Bild. DRDPIX - Zeigt interaktiv die X-Position, die Y-Position und den Pixelwert des Cursors an. GETWHITE - Farbtabellenindex für Weiß, entweder im X - oder im Postscript-Modus. HOR - Zeigt eine horizontale Linie auf einem Diagramm mit dem angegebenen y-Wert an. IMCONTOUR - Erstellen Sie ein Konturdiagramm mit astronomischen Koordinaten beschriftet. IMDISPGETAXES - Erstellen von Arrays mit Achsenwerten für die Verwendung mit IMDISP IMDISPWITHCONTOURS - Imdisp ein Bild, und Overplot Konturen mit subtilen Farben a la Tufte LABELOPLOT - Wie oplot, sondern beschriftet die Handlung mit einigen Text-String sowie. LOGHIST2D - Wrapper für hist2d, der logarithmische Bins und einige andere Optionen erlaubt MULTIPLOT - Erstellen Sie mehrere Plots mit gemeinsamen Achsen. PLOTDRAWZOOMBOX - Zoomlinien für einen Plotausschnitt zeichnen. PLOTRECTANGLE - Zeichnet ein Rechteck in Datenkoordinaten, wahlweise gedreht PLOTRESTORECOORDS - Gegenstück zu plotsavecoords PLOTSAVECOORDS - Speichern des aktuellen Achsensatzes. PPLOT - produzieren ein hübsches Plot für Postscript. RADECGRID - zeichnen Sie ein Gitter von Linien mit konstantem RA und Dec. RAINBOW - Geben Sie einen Vektor von Farbindizes von fsccolor zurück. SAVEPLOT - Speichern eines Plots auf einem Datenträger im PNG-Format VER - Zeichnen Sie eine vertikale Linie auf einem Diagramm mit dem angegebenen x-Wert. WIN - erstellen IDL-Plot-Fenster in Zeilen in einer schönen gekachelt dargestellt. Programmierung CHECKDIR - Prüft, ob ein gegebenes Verzeichnis existiert GETMYNAME - Gibt den Namen und den Dateipfad zum Quellcode der aufrufenden Prozedur zurück GETNUM - druckt einen String und erlaubt es dem Benutzer, eine Zahl wie READ einzugeben wenn der Benutzer wieder prompt tut eine Reihe GETYN eingeben - erhalten eine Antwort auf eine Frage yesno, so dass der Benutzer zu schlagen den Standard respone HASNANS zu wählen - ist ein Array haben NaNs darin STACKPOP - erscheint eine Variable aus einem Stapel. ANMERKUNG Wenn der Stapel nur ein Element enthält, wird es nach dem Abspielen des Elements undefiniert. Siehe auch stackpush. pro STACKPUSH - schiebt eine Variable auf einen Stack. Wenn der Stapel undefiniert ist, wird er erzeugt. STRUCTMERGE - Gegeben zwei Strukturen, A und B, die jeweils ein Array sein können und die unterschiedliche Tag-Namen haben können, erstellen Sie eine neue Struktur C, die die Vereinigung von A und B ist. TIMEIT - Testen Sie, wie lange ein gegebener IDL-Befehl ist Zum Ausführen. WHICH - Druckt vollständige Dateinamen in IDL-Pfadsuchreihenfolge für eine bestimmte Routine. Routine Beschreibungen Kategorie: Astronomische Utility Liste der Routinen ZD2AIRMASS - Berechnen Luftmasse vs Zenith Entfernung. LSTNOW - Rücklaufstrom lst, für eine gegebene Sternwarte SUNTIMES - Berechnen Sie APPROXIMATE Sonnenauf - und - untergangszeiten LASERTIMES - Berechnen Sie die Start - und Stoppzeiten für den Lick Laser in LST. Diese entsprechen 23 Uhr und 5 Uhr Ortszeit (entweder PST oder PDT) FINDINGCHART - machen Sie eine Findliste für ein Objekt auf dem Bildschirm oder PDF. MULTIFINDINGCHART - rufen Sie findingchart für eine Reihe von Quellen in einer Datei aufgeführt. MAGERRPROP - Wandelt Flüsse in Größenordnungen um, mit Fehlerfortpflanzung. ZD2AIRMASS FWCENTROID Kategorie: IRCAL Pipeline Liste der Routinen IRCALBADPIXELS - Mark anhaltend schlechte ircal Pixel als Nans Nutzung: ircalbadpixels, imgs Dieses ist nicht alle schlechten Pixel, nur diejenigen, die ich sagen kann, sind schlecht, aber die automatisierte Sachen aus irgendeinem Grund abholen tut. IRCALBADPIXFROMLIST - Lesen Sie schlechte Pixel-Liste aus Datei und gelten für einen Bildstapel IRCALDEGHOST - entfernen Sie lästige negative Ghosts, die durch Kanal-Crosstalk verursacht werden. IRCALDEWARP - Entfernen Verzerrung von einem IRCAL Bild IRCALDEWARPSHIFTS - konvertiert eine Liste von Bildverschiebungen von rohen IRCAL coords coords zu dewarped IRCALGETJ2000 - Konvertieren Sie die aktuellen Epoche IRCAL coords in J2000 coords IRCALREGRID - resamplen IRCAL Bilder für anamorphe Vergrößerung IRCALZEROPT zu kompensieren - Front-End Routine zur photometrischen Kalibrierung von IRCAL-Daten IRCALSTREHL - Berechnung von Strehl für ein IRCAL-Bild. GETIRCALFILTER - Gibt einen FITS-Header an, gibt einen String zurück, der den aktuellen IRCAL-Filter beschreibt. IRCALFIXPIX - Dies ist ein Front-End zu nsfixpix. pro, die über die IRCAL schlechten Ecken weiß und bewahrt sie in der schlechten Pixel-Maske. IRCALSATMASK - Bestimmen Sie gesättigte Bereiche in IRCAL-Daten und Maske zu NaN Behandelt auch IRCAM - und NIRC2-Bilder. IRCALADDWCS - Hinzufügen von WCS-Koordinaten zu einem FITS-Header für IRCAL REDIRCALSKY - Version des IRCAL-Codes zum Erstellen von Sky-Dateien für die Eingabe von redircalsub. REDIRCALSUB - Reduzierung der IRCAL-Daten mit Subpixel-Registrierung und Mosaik. IRCALBADPIXELS IRCALBADPIXFROMLIST IRCALDEGHOST IRCALDEWARP IRCALDEWARPSHIFTS IRCALGETJ2000 IRCALREGRID IRCALZEROPT IRCALSTREHL GETIRCALFILTER IRCALFIXPIX IRCALSATMASK IRCALADDWCS REDIRCALSKY REDIRCALSUB Kategorie: Mathematik Liste der Routinen GAUSSIAN2D - Berechnen Sie die 2-D-Gauß-Funktion und das Derivat gegebenenfalls auf Basis von Goddard IDL Astros Gaussian. pro Restriktion: Gerade jetzt nur kreissymmetrische Gaussians GRIDEVAL - Wertet eine Funktion aus, die als String geliefert wird, für alle Punkte eines gelieferten Rasters in x und y. MAXES - die Maximalwerte jedes Bildes für einen Stapel von Bildern zurückgeben MODE - Berechnen Sie den Modus (häufigste Element in einem Array) optional mit Binning POISS - Berechnet die Poisson-Verteilung als Funktion von X und M. MONTECARLOMEAN - bei gegebenem Zahlensatz , Verwendet monte-carlo-Methoden, um sowohl den Mittelwert als auch die Standardabweichung des Mittelwerts zu bestimmen. MONTECARLOMEDIAN - verwendet einen Satz von Zahlen, verwendet Monte-Carlo-Methoden, um sowohl den Median und die Standardabweichung des Medians zu bestimmen. PRODUCTERRPROP - Implementiert die Fehlerfortpflanzungsformel für das Produkt von zwei Zahlen, nahm keine Korrelation GAUSSIAN2D MONTECARLOMEAN MONTECARLOMEDIAN PRODUCTERRPROP Kategorie haben: Misc Liste der Routinen FINDCLOSEST - ein Array gegeben, den Index, dessen Wert in diesem Array zu finden, ist am nächsten zu einer bestimmten Zahl. FINDFWXM - Finden Sie die volle Breite bei Xmax. D. h. X 0,5 für FWHM. FINDLOCALMAX - Finden Sie lokales Maximum in der Nähe von gegebenen Position im Bild MRECENTER - Finden Sie präzise Zentrum für einen Stern, wie RECENTER, aber mit MPFITPEAK, um die Spitze zu tun. WHEREISMAX - Angeben eines Arrays gibt Position und Wert von maximalem Pixel zurück WHEREIS - Wenn der 1-d-Index eines Pixels in einem Array angegeben wird, werden die x - und y-Koordinaten, die diesem Pixel entsprechen, zurückgegeben. ATV - Interaktive Darstellung von 2D - oder 3D-Bildern. ATVMAKEMOVIE - erstellen Sie einen Film aus einem Bildstapel im ATV. FINDCLOSESTFrequency Reaktion des Lauf Average Filter Der Frequenzgang eines LTI-System ist die DTFT der Impulsantwort, die Impulsantwort eines L - sample gleitenden Durchschnitt ist, da das Filter mit gleitendem Mittelwert ist FIR, die Frequenzantwort auf die endliche Summe reduziert Wir Kann die sehr nützliche Identität verwenden, um den Frequenzgang zu schreiben, wo wir ae minus jomega haben lassen. N 0 und M L minus 1. Wir können an der Größe dieser Funktion interessiert sein, um zu bestimmen, welche Frequenzen durch den Filter ungedämpft werden und welche gedämpft werden. Unten ist ein Diagramm der Größe dieser Funktion für L 4 (rot), 8 (grün) und 16 (blau). Die horizontale Achse reicht von Null bis pi Radiant pro Probe. Man beachte, daß der Frequenzgang in allen drei Fällen eine Tiefpaßcharakteristik aufweist. Eine konstante Komponente (Nullfrequenz) im Eingang durchläuft das Filter ungedämpft. Bestimmte höhere Frequenzen, wie z. B. pi 2, werden durch das Filter vollständig eliminiert. Wenn es aber die Absicht war, ein Tiefpassfilter zu entwerfen, dann haben wir das nicht sehr gut gemacht. Einige der höheren Frequenzen werden nur um einen Faktor von etwa 110 (für den 16-Punkte-gleitenden Durchschnitt) oder 13 (für den vier-Punkte-gleitenden Durchschnitt) gedämpft. Wir können viel besser als das. (1-exp (-iomega)) H8 (18) (1-exp (- & omega; & sub4; (1-exp (-iomega)) (1-exp (-iomega)) (1-exp (& ndash; H16)) Achse (0, pi, 0, 1) Copyright - 2000 - Universität von Kalifornien, Berkeley