From b9b51dad5177df6f8f5d23dd98dbf6a3824d0735 Mon Sep 17 00:00:00 2001
From: "Knollmueller, Jakob (kjako)" <jakob@knollmueller.de>
Date: Mon, 17 Dec 2018 14:48:21 +0100
Subject: [PATCH] now with correct response of the human eye

---
 nifty5/frequency_sensitivity.npy | Bin 0 -> 8320 bytes
 nifty5/plot.py                   |  21 +++++++++++----------
 2 files changed, 11 insertions(+), 10 deletions(-)
 create mode 100644 nifty5/frequency_sensitivity.npy

diff --git a/nifty5/frequency_sensitivity.npy b/nifty5/frequency_sensitivity.npy
new file mode 100644
index 0000000000000000000000000000000000000000..21cd6145ac2d8f83e3d600a4cb5973d1a3bc6b4d
GIT binary patch
literal 8320
zcmbR27wQ`j$;eQ~P_3SlTAW;@Zl$1ZlV+i=qoAIaUsO_*m=~X4l#&V(cT3DEPSsIR
zFV09TNL9B|&@eJJ(@`+d)KREaAQNzfYAWAM*c)uWXY&fba@`htlNHO}&zE$uXMlq3
z@7W$7KM`!dWw%G##h_OEzul*=&UJBx%S&#1^+xYXu>I=8FKV?b+U$?#f7m@c-wiHb
zU32c+)ce8qi%)k<ZJpb0KQUNz(ac%yaQVAVn;8DS3bvnpS!3<}Lml=}=2jikk9xr6
z?K!7$rF;vvuhDsUBR;pwUf*d4+tzw7xctKNl_IAZL+rEjru+TU?XiEZQ?zkPpbuQ0
zd0MQT5pRfn)CS)pQ{MI3&llH;DN*)=%csUTYRr@fvG@AKY^A@x-#(n_8C%F#f4KZ{
zB~>E^l@NOwsk_&oY@cA?@$LMLnU4eD@|rI>Y|`~Z>_x112pP?oX#eE3&brPsfpGcm
z^>19yS%ui&$j~eFpEb!|cGHcSmUn~T@^4Eo_?x<g*iVzXes<l}$@Yt9z7MlB41vr0
z+iZ!P6BuG2a;9F3PjITeMkAYuWosy0ek0S2Waiiqd%-%5_FSK7_UkfaFIwCTgUbt_
ztj^8K2(dq}fAR5|_UZN}0VaJGx)E^s(*F4Ji^U=KJr64#bPmt3-(~LYZP6GBm%kqB
z+-OlBVlS<H+(+yGO#7dOulg9TMZx7QWz@UpcZJxOh3#x|-8$P|(xN2t{?Qn?{Om_u
zQ`x44*k9;ey(1xeuD!*-zsh^J#=+(Pu6Q>ucVURV_TIVgsukzk>pXHUKO2+)mygfC
zwBpK|5PRnz?vqYuFR-^c+7YwID-kY#&~)1->m4EX9|RqBZ@MkC_gQh)akWztTweLl
z%sq<^huF{4wmR`laglvYw~+cQt7N!*>!I3XoaaOAXRKbyC2h6XzBu{Fd6CW(xcu{u
z#0v$tLhN}9GnmBEme{A|R4nTGoC=rs^mn~={b`7OvzhuYfk{j4qe^v3S>n>+@~cHP
zpV)m2vH$V==|}E;%k2H?9=O+>$bif9-r;$-^iPO=-oYg=*<LNTcWR#^^~ER)E}y^X
z!*?EzQ2WPi$qyODR@$3Q@OYa&HybX0G2;qjiEya>8#dO?g638BOq)3rGKzBH^2Yi*
zxNpma+OvhdWnz(BV}B~t<?W7+e7OA7ud_v+v_kC_mR_xyzHY6(*B|@4E$0j1^51vY
z%da#IwYU1R_nVf>di!-}?>eu&RS1`lY)aM?a16Cyr5J2+KWBrz@cimkSDzNa<#&6y
z8<+cp+M9d2e_dg|(Y~TUQt;o0Vz|64zqakY@KF0_tdl$qwKv(%SoY@l-JnvqeEn5E
zch}@l`^5#yNsU&U?I-mAPQLrT3@-n0_NRa~d7<_Rx3^SX^xR_KQNhJ<x2FOw@0@fk
zO1LW2UdGLBvSQ>``-VixId?Rx;POi~cBNFchT7lUaD8=Q#y0y3FAa-3$ExA-?C<C1
zJ)96~KTSIG$iedM_65e~S8gZQ!sRo!Hdc7d3ANwpqT00V;12sF>xesZJL}={r>j$&
z)-MmWxBJTcre1xQy{2pPfxsUPaCu!P&mOVOq4r0XJ_?QO-EF@x=hmbv+D&lzeh%Gf
zwfjQtJ;M&}H~PHCUUG_IYIuAzT>j%Zfdx-ahT5NFTcF7kxzB#$?hf-?lUm^N!BfAi
z^u8KuAF{n}^3#3$?KvLb=Zrnp3YXs+ePi>c2ch=$3MR2^2M*Zx&fKqYBfT9iFRrv_
zujHFh`&07@-=A21&|Y}HvB^TIPPlyKi}}YJzK7b2$*yng?muKdbIEf1tg~Hk`P=K7
zFFt1qvrn4+M%b_Xu)X{$X>a3(9=N=1Y5HwHzA*cBl3L3I;*Z!bUEd$}$EXi3Ki}H>
z=~l@w`#CeExs&{k+8^B2St|FhA1=?xr2k$<HOyY`WJ7La`Z4?8Qyi9jQ=bTzPd*{^
zv)Lfbe*cfB=NB{_x1Z*-hxcB<B)I(1zHdygt;6gCG`x)t%sgRlz+zE;v~n_BUOn^{
zPmp_<{loALQ|@m%X@6?><_BvQO@YgIO79cf85Cw;)Gqq{-|18Ksm1zUQ%+BX%fEiO
zP(dLs%zma4+v#&(Puu@<{bp3?JsmFZyP{RAJu}Qc^oNMn0o}9q-{vMdy%(7YmtUWk
zY4WZl%wF)Ea_7d>bM|k)o(asoHw!K=XzFVh+7M=ccCpcmMRU*FKlQape7$@QT)z0X
zfybWiF#Dr>7AT*}x?q22>EeRS!g+A{s|Q7bl&6Q;Z*E?BV;1{Gd+t|!^+)9wz~#+5
zenfXI3bSXvpuW-W@J0K1zOO3!qZh*EXZqbq{kS&FK1gO+=cBqy_J$je7DY~91egCU
zav(oqXPEt>43(rdtIPIBL>A^KUtSECkG;LPa{rMq`;S`|S!sQ{Y#&}xllGQxDO`U4
z!uDph3t{#SjEchNHeRuRbvry^yVo+fykbUnZ}06ed#^8Vlq;IA+CR8m>~iz$a=3i6
zp8xbO&%*32ez))xCttJgS<SX3s(vM0{^=K^h0&kF?A4iMzven$x3{RCseH#^HC*0(
zx7g}Kf5Ys>IJe%-P`+XR*h^<}?4vbs`IU{owrFyO+p`Mz9!qAvY2PokpXvVOb#QrZ
zk9+$jh=kj}{hD+7`^%g5flo@lFk5VZ%jfbPI`LgT+}?mI=kcD0x9oME-{D*pyAdvb
z{>suzaoXYbC34xn8*kpWXM5!$AJDuBE^j!y^UhJTaQhubS;FoY@7O<hv&Q7xip_BO
z$%(nobezKN|F~sp^PIeEzw5n-*PL@(;PPM913ygm4Yz;y>SRsAfqV8p0>187{IL}-
zAO6nd*RP0hd*$Na(sS3`w@;S05z%Vi4wv7#MS>+MCER{Z?X<f&J`e2AzirAe;@b(A
zm#+T9don-V{<x#QZug-F_EvjN?y;J=3oc*l_(0sCI^6yt>xr6`77y*0bqNSMtL%Zx
z-)BFfIK3_0{?EyTGso9Gv=@pDPxW5E7cTF3Zn^fqiQ)El*Ux?aT=S9r47r6S=8N{*
zYqT&;Pd>NaZt;|E)6}`)_Q{VpBi1Z^Wbe;5H^S=T0sGADoGRLH4%o3q=i8rI5pMrd
zr0hh2;$!=u+LvOH7Z2K>I2)9-_sL1Sbj4s#<1OL#dG3!aOy@qf5BXwxwfy2Cdy!AR
z7nbK;v^)9SEO^%baQiQnek*@UJh2bYZkjyp;$eH|+M@R#UR|}*UMCgHcsksEM+46{
z=P6I@SswfqX5~0)|K%}5tWd`dyWWz2X&KkT?R}dnQ;PYY+OsWsu9oI?%s%(yZ?-eL
zZ`r-KdQ@=!VYvOH=4Y1{bv?D`NV#Zvyyv+6vwz<<w7$P%7sz<5%Is~peMM`q<0aN-
z_S|Z_eKgLWuutIm6lJe`-){5ql`V6BgxmABA6U=W_{^U7#gh0=!Bh5cddiw=RzI*4
z?d$DhVTrK!%CTcn`S#3y-Kk4;d6Q4u|C*uTdBy9IU3p01jBNf0`-pRIW{PD!x0hJB
z<j!%$v-S>qb?qZqp4i=#3SD$bD#AY1V75Nzz32AL<z-fi8_wCcJ+v~tv-hc;^+Ss_
zmTD39MUy;_{tkR$|6NVR$5#G={b^PY?S%5@c5|1@Y+GO$VPF5{d(^w5FYI-Mcz<U;
zyI?P<6{`5e;HBOFT!sVewh{K%IjhY7YQ407HhslU(+d~vy**YjuDbQwF45%4$vlq;
z`}y-qOaJt}wC8Ga>#y8%$$qky%C5et@9Yl$I)3?TaD;uDapC6QpI+MQ<vdZ@yW+Ba
z9-IG8%M~B&R1dDcYZD(~uX`fr<FBw+_Q7FIyBMZkv7hw4TXYTUC%g9c{^yIcBJ2f^
zEfv1G^Oe2Z-&_9*TdvyQbMad%SM}NMrEl>^uCfUGy~%7blZ9W~n@cZPHK+HQeWk#j
zO*>Y6weuDU|6SM^VV}vEzQi~8wY`RAYS{5<*X>2_XzmEQ{@rfvEi2ZWJrVXI^N+t_
zKmOWYGF*%4)BGFu>z7a8`-b(G9sfc(euo(m_6KY<Ob;o%vF9p%zgA}DO?%H`*28Tk
zf9wj=nIxAlj<El@%{Q^6<c+;xRN<U?J8#)5l((H(b^o8;W!+~=eCs3Z6E|_k%-H+J
zK48vX`=8ft+fSW3MK8gK(cbjSNuAPN5%w3>o($W=_SU}d9_QoWFL&%`y?^@c)fXoF
z>ATjN-8mXz@368x@PYqZ`!fn*ZR-T@*)Q-k-(R|&#r|i*L<i@K5%wFFxcacnd~460
zyhxK%|Gxdw#U1SoRc!Xr?xkL<?nKyErs{QkzxdYvT%yKADYXaoZ_4}5e*Mp3zn3>M
zMDRs~eV+nT=WgbA_Ve~$H(Ivkf&G2PZy}5BaNEmYv5Bkr9AUrs_SMdM(|7iHlI>Rv
zogdm?>W;nRx0BDlah5{H{eKbm+o!GUa!!9|Z(C|sw)f&g`^R06HvT;&VE;IQrO1sd
z(%!yx-|{un-r3t$?J%5|@W_7JW<J*m{zCSy>MyF-ibmQiMY8g4S^v&{?$n|!%RfD`
zKg8B-Sh`%qe)-$ets;t%_N=--MY~VDv;Vt6$9!YUWBa#}-({ozi`jE-S>Io+6KVgP
zVe<7u58v6xo_V%&m)H|~$+g*B_TiHD*;SKgJ~EHAKXv7%&8c7S?DsxjU~^=}6MGki
z2j7&}OWU7yC|~U99BHq^$v?48_PxE^7J-)soSxcSe4hM*<%yiVK6~`q4StdKnuRwy
zXSlz&ulJKl*#GFMy`4+4>W6qm`-x}mwu?ta+N<4M*tDwfz5R|C>eu%+J+pt%S99&I
zud@B;$w~+7QX}n^9kZ+V%zkhGvB}tPkJ5AdrkhjSPg|?nheolTep(P|Z?_@LY{{Yb
z_HE@gyM+!sxA*GXT-EQVZoggO)fJzbNc#?5bHTqC-rI9;(4Tn5<AwdMbIa}2>$U79
zpP#+AxjoYU(kk&cq4(e0Pk8Ytt?l9qdma8cO#9dC+E=gL_(E!Oq`kDpk7MiJzPA@u
zoMz(?_0s;Fx<|5hfr0(q;;Ek+=SAA5F1@jk{m*;*S&^PXzaG7`XZLV&E$1+@x3{YN
z^I~PBJ;VC<43R(H+w<JwDxDGa%06Wh>xGlaCiV*%V%Y+=M%tgznLY9O+xPbTTQ~Nc
zIsM9BD^TQ?9f!F+^Kl1(?FS<5XRLBlDtr9iUU)*p5(eAX_8<9`pKM-ZX`kAwA|-bw
z(mqw=(<YXi@9o91|L-<g_u78@Im34vfj0KXgE>@MuSeRq9NfHV=85<AvpL_Z+laie
z4_W_U>eN;C_L@>}bl*IRw3j-QoL0O3z5T^$F54Ni-q>#wbCh1n>SW*j;JkV8yGZ->
z*My}Kr@go5P?_a+eE%E!gI6bJ|2K89f4gk6<F22P_AU<(pLMByZ*RBp{nZ}PxAteI
zJlY+Z=4S7oJKbB6HOk(;DnY?F`Mtf?RlDgIGT+*poi}8gd(Ojtqj7a;hd`8lyN&CQ
zsV?vB^V*Y-@7ns-euCTA4@+}=?1g{D#lM%1vOoV%^TIvV_x77CJH9Vrd}seP?oMtv
zYk+;}0jJC`^(cFZgM2I0+27lLd%Z@z-}jyUA}PDa(N_ZPueYfc?=_0DPiX#JU-06c
zz2&@%;gvJr+52CcACa&)*xu5Yt5(G>%HC<mx2w(D-r1Wiy!Cs-;dk~Cg3+_=W`x<#
z7Jk>(?HOgCEc;yKdh<K`PkmX-roVY-fA5@N-1Hfd_J40)obV|m%D$!PdazdDJNxxD
zzXIw+-`gkmoKXBQBgQ^{!PZ%k2~qaTo*$o7EBDSmJFDw=irstr7oPP`J!Zz+A55FM
z<Un?my+A;^^3+dn?FByF?mm(B-ah|X__XIDN%qRRwd>T&qwHg!u8L$i@z%c9!dFVD
z|GmBP4)MRYTT<*>KPT+yYl^a8S@o{Jc<x*KS9SYVL~nm@f4%d1@`W#H_Rn{^9QxWD
zW&f4m=-8v8xAx&j_05;we{bItJ^j>?luUci28}Z@Go$SH3f}+bqWjkV#}3c<=TG0;
zyI4G2s{1p?ezhCVwZluI?6)16o>PARjs3)ab^W^!-`nedm*L<kF0kk2{cvAvLzMmc
z@_<JxC%&=QPGkLl>F#@b+aOQd#>^u7{L7bLPTU=3ze0;M;*Rqh`y)1w&z-pW-rj3%
zY0<a%68noYw}1X|EXw}big3U0H(%R(>$pxeU;p0TZSG&E81r&_<Am9N<1a<oJL@la
zWu5-oUO81`a(l&ldk(p4g{`4g_EXjB*^l3ivR|<EMO(v{SN6Z=a!<VO_}>2HymbrK
z71!E-f0HDr_cF@<pS7X#(Z*Nym#)3*Q|Eqff2{xT^Vc2?_K}<2q^Ev~vKLoaeu|6l
zm3^ocd!yEgclJuB{_NsyY_i{7sipeof0X@c$Jlz&SugD+CEH^{o8H-%N#0xP`k=*L
z)`3qynLFCPK3UjALG7je-S2D`4L0xWPiO3(P;1|AU(fc*;*?mly-DlI7h0QM*e|@+
zCcOIFTYH5a3#zxy>$HD(=BksSQnY>d1TLFM^%wSbPQq7P``+4XI~>^BrrK-oJb9<j
z4Bcq^KimsUC(e3qFZF%;v^?dv_6naF58O_gV81kSZrFc|X#2c&Q~i^?&+WgjmWc@8
z@y1^3n6l~@qe=Gc3XKVAF46X9!*zV#);+VIrWTjF!1IlLgZIPgDWX&CGoPhoo%N5l
z_Ye5o^z`Oa`;zOX=C)g2+h=g=aeUu9&HnTnj}nvUX!{EuyLR3+d1@cj!OM9|?zMfu
z>8dG_X*2D0i*@Q|r$yTbIW~T{+W*8}-}33Lvd&lb7Tw1c_VCTI@3$1_U@D5XUm+*7
zIOG0fd%o@2)m{%?+81ZD?3>Fy-~QwO&l59iqwQy?ScGtjKDJ+1J7Jp;*Gv1Tn^A6g
zDGTj`k6oXAp(EP9Lo1GTmDeMCUze>PZyCR^|5dYjEwA%p`>j2@mzqzBwy!W~+2+;o
z(EcsgEuF=&&+TU^@QAKeTxu^KJb(SX`O)@@yBkh<Zhl~Y@V7M2`r2pqD}s6gS>u=6
zS4uYTWL+I?ui{o+wfn?<d)2bq{R^i(wU4ZF^UHm>(*E}Sw8J^uqV3tlzQ#Y+x@X_{
z;a<^{^-t`%7cqEWO<rSfyUgqC<%7}oFP`SPa`)b`f1Bl?)_&r#{oa$i+-;w&v!9=%
zcirl2w0%_FKXo3y+xE+j<THspdt~n+Ut7!@yV0J}Sm@!x8`1XJ%(uP-Pr7N(&~DBy
z$MVR2*WH*;FCsSEC;$BVisMPNeRapygR{kM*hhK(=hjeqXs<oTrtnC}R{Nv-Z+^*t
zA8pT6C^j!+)iwL4?om^v+#lF)>RMv*pkTYbdfVQA*M3FYKXQ6}?RDx^d-kc1AD8^T
zXa8yrpKsIso%WqR3pwoAV(brzY~asVykhS%OK9$vjJx)V>nl=K{Px(t7HSn*Di~wG
z=vj@&r#qMI8|OQQeYtwuKH}h`+GBh6+56tikl~SyvFDst#=Gz6MSBIghO$}Zx9oXC
zC$B2mc)))Bd>^%9jTn22Nuied2QSzcpP#OEMenA4N#5G5&-{n%1=9=+ZW+hepOwix
z<xqLv{%FJcd%VxC+uv$Gu2qzG*uGdt#LCe=#@<<H!-Kevv-Z8rZ~xY<zGnYxW<u=p
zr$_9se){gb!Yjs}?b@b#ccwG;_jHa}O<RA}KH-D+rNbwW*_-dY?aLn)V}GK4tDHg2
zDSQ2}+2w0DU$H;#X72uA?Fsvt_4~ui5@YPUHFxZl_;$j+p}uI#k?oi5_2w$=`7`m9
z{cqRBiFb2i>_Z-TotSa*xc!&e4%H0bF4@m%&$)ET=&XG#PkXjYMU4H)|9^FN96V-!
z{r1tt3_h3a_jAdLrVCuK-+w8qbaiu#y`AjF$=7!twcjP}eTHGtMf;5xnGdG4U$j@8
z=~pk*A7kIZF60xm`-uHVjex>4zc1MLta|9~DSO$zIo_zVa#oDJfcL&d^RFDX&)@o_
z%Q60f{eS=8PPZAa+CNnloBUvDjQ#w+te<b*9<smH{K+G?^}K!YR}s71XV>iA-~628
zzA?sL?~+)675gFkYjdP)!hO!!Kbxgz`Tx`n`<0vSE?c)J#{S^RNw)RI2klEXH(uSS
zaMu24>58Be8n^AaD-UiEJrQF+`SjD>cM=cSukz;$>-cusK74s()3?cY?Q`vy?y9*Q
zV}E>BVM&R?e)|WnXI+asbIShZ4X#Zp|L@zMXYD-l_+E_tR$=qhr4#nrtF3j^*}mkY
zz0J%^nO|)l+8dt9Ip_5%#y;+E|J~dCd+axM9)0((?u0$tmn`+zhDY|3Ck5Qt_%+6!
z?fHvPq5fU=)@hy_-6D_MU+`modgSn9`>zouk0cml?XCW7ue`~>(|*Ro9P2dQWA<s)
z!fA%@p4f-WOT4M)iM9WID@X9`l<oE@k^L=FACK6JYbSsBfAy*T&Zob>J`<0%-?LM3
z<pIfU_68R>o;|nau>Hwfo2t(3eP%Dc=Kg<Q<yd<whUBCF=WnrJV&pnYzw(g%)i|5d
zPm7-0*A^Y-+^QFAFMG=KL1yG;d(T4=Ck-MF+V?jdeY;2Qh5dbt<-#&nvGxp841Ql0
z+GPJOFz1J%`2l;cWqxVxvtHOc{_m1)c8#@v6e})fcV>gV`r-e6X*&DtSx;TR7by19
ze(}*f^;ZG0_EVf$7uY^oZ$ImmlabuSJ@#SOxF&5o@zS2PC&(}`Cf2?(%Wd0*_;vPS
zr2$u_O6|6v;pIC^rTvwCdXSm*j`UdjxQXp&l5elE7h;Gj<6W}T{?@+!MUz}#+n<z>
za+NQRwSQ;)B>2_7)%MQPPKVER?y%ptu43m+tvB}C_y78})y3Lh+&wAn>y1_RQI@AR
zeVf1C-g;r?ga(zj_Pt9VM!fBewcqOOSjq5trF~&|&?32Q+w6Btj!IXTduRVX`)E?g
z)L8pp%hr7SF22%!3&(z2-&0%d7qY8AxN!8HeW20GoZSmz?Z4eL)2?<~VJ~!Ht=*wJ
zTkIRImfyRq`rdx?kDfB6HL>=ecn+rgD_L$|IDL-2`TNcGaf?sfxz_sLUUYv!L+AEb
z`&oW|5$wyC*?)4GoBNSrv;97iN0T4#d~aXg8q)RQP^|sG8_W*jcbD4h$nYOHA-Ku@
zgSSOe-;Ve8H@z*UgrAGGPfXsazFc^zed^~1p#@4C?UkDQ3@o?5x3?COnYZs|ti9#2
zf2C{mme|{F(E7PVaf5x2>lds0A@A+y-uS;<^=Yhq$$MdDH>Jh)zP&fHuk)_6|5;h^
zF){6({r`E7H}-spwO=oxIpMzOB722h(N9xuthP^G_h$K%^0)ShsmFJJ{vB&yzM!|p
zk!PX(&66+ouUfUzUWY+D=SRgG`@`C+k4CY_+5eh-y?+PW0{bIZZ&<r@EVqy1uKNEz
z@wL6`r@r$Ch2rdg={dW_GtIMK|M0~2yQxd<JH<{;dhYegzI{j0O%1s?`!~ghgEIN%
z*z@nXwO+GuvAxia&}Hm`FYRB}g+1=qjI&?*Q$Oyd`7HbNBWWAAZ(C^pc!6`ro>kB7
zy<M%|elv-)kFB4RnVC4lzJBJ@)b1bi?YDLsW)-?TvtP?C_dV7j&Yo>?x#;rM)9jzM
z)!d$QbB_H57iHZG*Pqz)Ut(lD>K$h<RMxZCM0BdXP;uGm$2_y_tC!E1BA@!$zHr7<
zE}if=ds)wgsfX$&+uOz$?ACIaZm-J4_b2buL;K5dCq*VF#o2%Q@=H+V$V7Y9X*&0&
z&zxc}_2+1hu<HYRQ`NO{zjEX3k8H15KKWU{eb~cTkqWL!_I4#&TLPZnv7i2Wf<|Ix
zoc-cfX6^efJ@x{Ta#J^a>$CUxaKf+3<);15jU`4WTH@>%7wo<pc)7#ginV$!&*U!q
z#UT<hsn@RBM^{AH=ue2VzxRLh(+vf!_F_)nvu!im?2kF~1+T2QWWU$W#%<c{ID5(d
zHE;R(n(RN;Nzcr)Z?eCV${}*w^qjpst3trvWpVa1oRhDlhSu4yKD^B$BBIV-U{dy>
zK+co)ji;HTQZ~ieYww=@L4IC^z0SdQ<z>MY_DgT{d}H`}*#7av=gFt{#@YA8Ul)tF
zF1Ba1_+I3*zR3PWU5->j@;-alh|{@7C*$llU3vQSj9tF{T2}?8ol1H3y(^yF%W>IZ
zzg%v8`OGVE_CBv<=G|t?u#b=qlDqjV%|7T-%;HMbP4=8mCp9uWh_hFATq@x$oMgZM
zaJn_m{{;KQD;_478&=t8uP*COe;sGfRPVR<?ayfYqdj*ntBFV1H@TeY_6=HMf3`4c
z>bY-m_Is;?JYCNQ+w)9%X<%v)Wbc|e|5(5K9D98WyZNR}@%GB+g&wX{_O<`#VX|JQ
z(#w9wp___3S5C2?_)l@=9Nu_)`M6$|>Nr>XwYK{O`&^yvZ4Sx2t|{oV|9pgXGqYs8
ez3DCffEh;C_B&Fhytt=qVbAfuby9n8vpoPN6E_zC

literal 0
HcmV?d00001

diff --git a/nifty5/plot.py b/nifty5/plot.py
index 6fd18c940..3c3d485b0 100644
--- a/nifty5/plot.py
+++ b/nifty5/plot.py
@@ -61,16 +61,17 @@ def _mollweide_helper(xsize):
 
 
 def _rgb_data(spectral_cube):
-    def _eye_sensitivity(energy_bins):
-        a = np.arange(0, 1, 1 / energy_bins)
-        rgb = np.empty((3, energy_bins))
-        rgb[0] = np.exp(-(a - 5 / 12) ** 2 / (2 * (2.5 / 12) ** 2))
-        rgb[1] = np.exp(-(a - 6.5 / 12) ** 2 / (2 * (2 / 12) ** 2))
-        rgb[2] = np.exp(-(a - 10 / 12) ** 2 / (2 * (1 / 12) ** 2))
-        rgb[0] /= rgb[0].max()
-        rgb[1] /= rgb[1].max()
-        rgb[2] /= rgb[2].max()
-        return rgb
+    def _eye_sensitivity(energy_bins, spacing=None):
+        from scipy.ndimage import zoom
+        rgb_high = np.load('frequency_sensitivity.npy')
+
+        # if spacing != None:
+        #     spacing = np.arange(0, 1, 1 / energy_bins)
+
+        rgb = zoom(rgb_high.T[1:],(1,energy_bins/len(rgb_high.T[0])))
+
+
+        return np.clip(rgb,1e-15, rgb.max())
     rgb = _eye_sensitivity(spectral_cube.shape[-1])
     rgb_data = np.tensordot(spectral_cube, rgb, axes=[-1, -1])
     rgb_data = np.log(rgb_data)
-- 
GitLab