From 3ac7685e4b3b16e31eefd04381718a3d39eb4258 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Micha=C5=82=20Drza=C5=82?= <michal.drzal@gmail.com>
Date: Sun, 26 Jun 2022 22:39:17 +0200
Subject: [PATCH] Added scaling using zscale

---
 rfplot.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/rfplot.py b/rfplot.py
index 7070b4b..8373bbc 100755
--- a/rfplot.py
+++ b/rfplot.py
@@ -18,6 +18,7 @@ import h5py
 import json
 import requests
 from astropy.time import Time
+from  astropy.visualization import ZScaleInterval
 C = 299792.458 # km/s
 
 from modest import imshow
@@ -112,7 +113,7 @@ def main():
         site_location = wgs84.latlon(site["lat"], site["lon"], site["height"])
 
     # Create plot
-    vmin, vmax = np.percentile(s.z, (5, 99.95))
+    vmin, vmax = ZScaleInterval().get_limits(s.z)
 
     # Time limits
     tmin, tmax = mdates.date2num(s.t[0]), mdates.date2num(s.t[-1])