Skip to content

Instantly share code, notes, and snippets.

Milad Pourrahmani Miladiouss

Block or report user

Report or block Miladiouss

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
Miladiouss /
Last active Oct 29, 2019
Demonstrates that evaluating wcs.pixel_to_world is efficient for large vectors but very inefficient for single calls.
from time import sleep, time
import numpy as np
from astropy import wcs
w = wcs.WCS(naxis=2)
w.wcs.ctype = ["RA---SIN","DEC--SIN"]
w.wcs.cdelt = (6.2/3600, 6.2/3600)
w.wcs.crval = (30, 60)
w.wcs.crota = (0, 15)
# ----------------------------------------
Miladiouss /
Created Oct 11, 2019
Scientific Plotting Template in Python using matplotlib.
Author: Miladious
Latest Modification Date: Oct 11, 2019
Creating a publication quality figure using matplotlib requires a lot of tweaks.
In this gist, I show the main tweaks for creating publication quality plots.
This is an example how one would plot astronomical images in python.
Feel free to modify to your field's standards.
Miladiouss /
Last active Oct 25, 2019
Usual imports and function for an IPython / Jupyter environment, specially for data science (e.g. matplotlib, numpy, time, ...).
# IPython
from IPython.core.interactiveshell import InteractiveShell
from IPython.display import clear_output
InteractiveShell.ast_node_interactivity = "all"
from IPython import display
# Scientific
Miladiouss / HSC-Cat-Selection.sql
Last active Jun 15, 2019
An HSC-PDR2 SQL query to extract RA, Dec, tract, patch of clean galaxies with a magnitude cut.
View HSC-Cat-Selection.sql
-- An HSC-PDR2 SQL query to extract RA, Dec, tract, patch of clean galaxies with a magnitude cut.
-- Useful Links:
-- /,%22d%22%3A-0.019535944017082533,%22fovy%22%3A0.0002503311005562077,%22roll%22%3A0%7D,%22sspParams%22%3A%7B%22type%22%3A%22SDSS_TRUE_COLOR%22,%22filter%22%3A%5B%22HSC-I%22,%22HSC-R%22,%22HSC-G%22%5D,%22simpleRgb%22%3A%7B%22beta%22%3A22026.465794806718,%22a%22%3A1,%22bias%22%3A0.05,%22b0%22%3A0%7D,%22sdssTrueColor%22%3A%7B%22beta%22%3A26003.369421497522,%22a%22%3A1,%22bias%22%3A0.14439999999999995,%22b0%22%3A0%7D%7D%7D
, tract
Miladiouss /
Last active Jul 1, 2019
Calculate and appends ordinality to a dataframe
def appendOrdinality(df, columns, ascending, strictlyIncreasing=False):
Appends absolute and relative ordinality to a dataframe.
Oridinality is the normalized position of a row in a sorted dataframe.
Pandas DataFrame
List of column names for to be used for sorting (e.g. ['prob_1'])
Miladiouss /
Last active Mar 11, 2019
Print a stamp in Python with date, time, device name, and user name (e.g. 'Fri Mar 08 2019 at 00:41:02 on myPC by User')
from datetime import datetime
timeFormat = '%a %b %d %Y at %H:%M:%S'
import socket
deviceName = socket.gethostname()
import getpass
userName = getpass.getuser()
t1 =
print(t1.strftime(timeFormat) + ' on ' + deviceName + ' by ' + userName)
Miladiouss / GPU Detection Test in Python and
Created Mar 8, 2019
Ensure Python and PyTorch can detect all the GPUs in your system.
View GPU Detection Test in Python and
# GPU Detection Tests
from torch.cuda import device_count
import gpustat
If pytorch raises "RuntimeError: cuda runtime error (30)" after suspension,
run the following commands in linux terminal:
sudo rmmod nvidia_uvm
sudo rmmod nvidia
Miladiouss /
Created Feb 22, 2019
Simple code for bifurcation of logistic map
# Library import for plotting
import matplotlib.pyplot as plt
%matplotlib inline
# Define Logistic Map
def f(x, r):
return r * x * (1 - x)
# Define history and parameters
history = [0.5]
Miladiouss /
Last active Jul 12, 2019
Easy astronomy FITS file handling for Python. It includes easy cutout and save function. Additionally, a percentile normalization method is provided which is ideal for scaling FITS files to better visualization (similar to MinMax of DS9).
import numpy as np
from pathlib import Path
# AstroPy
from astropy import units as u
from import fits
from astropy.wcs import WCS
from astropy.coordinates import SkyCoord, Angle
from astropy.nddata import Cutout2D
Miladiouss /
Last active Mar 12, 2019
Image Histogram Transformation.
import numpy as np
from matplotlib import pyplot as plt
import torch
class HistogramTransform(object):
Transforms the distribution of the input tensor to match that
of the list of template histograms corresponding to each channel.
You can’t perform that action at this time.