Pablo Romero wrote:
> Eric,
> Its Still not working, still getting white/empty areas where Z=0.0
> Here's my ipython output:
>
> [6] ga-> Znp=np.array(Z.filled())
> [7] ga-> Znp
> Out[7]:
> array([[ 0. , 0. , 0. , ..., 0. ,
> 0. , 0. ],
> [ 0. , 0. , 0. , ..., 0. ,
> 0. , 0. ],
> [ 0. , 0. , 0. , ..., 0. ,
> 0. , 0. ],
> ...,
> [ 1.023 , 1.551 , 2.5079999 , ..., 0.72600001,
> 0.85799998, 1.023 ],
> [ 1.58399999, 1.74899995, 1.9799999 , ..., 0.92400002,
> 1.41899991, 1.58399999],
> [ 0. , 0. , 0. , ..., 0. ,
> 0. , 0. ]])
> [8] ga-> Lv=(0,1,3,5,6,7,8,9,10,12,14,16,18,20,25,30,35,40,50,75)
> [9] ga-> norm = mpl.colors.BoundaryNorm(Lv,256)
> [10] plt.contourf(Znp, Lv, norm=norm)
> Out[10]:
> [11] Znp[0,0]
> Out[11]: 0.0
> [12] Znp[-1,0]
> Out[12]: 0.0
> [13] Znp[1,0]
> Out[13]: 0.0
> [14] Znp[0,1]
> Out[14]: 0.0
> [15] Znp[0,-1]
> Out[15]: 0.0
You missed the [-1,-1] corner. I thought at least one corner would be
non-zero, but maybe not. At this point, it doesn't matter.
>
> This is very frustrating. :(
>
> Is there anything else Im missing here?
> anything else I can try?
Now I am wondering whether there is any difference between your
numpy/matplotlib installation and mine that is making this
difference--although that really does not make any sense, either. But
let's check it. Please use the same method as above to make Xnp and
Ynp, and then use np.savez to write out a file with Xnp, Ynp, and Znp:
np.savez('XYZ.npz', Xnp=Xnp, Ynp=Ynp, Znp=Znp)
and ftp XYZ.npz to ftp://currents.soest.hawaii.edu/pub/incoming. Then I
can easily try to reproduce exactly what you are doing, and we will see
if I get the same result, or a different one.
>
> Im beginning to think that negative/"less than zero" values are actually
> being plotted in those land areas, even though my array is only showing
> zeros, '0.0'.
I don't see how this could be.
>
> Is there any way I can plot the output in a grid format?
> i.e. print the text value of each point in the array at each point on the
> plot?
> This way, I can see if a '0.0' value is really being plotted in those
> white/blank areas.
>
Not trivially; among other things, I think you have too many points. And
you have already verified that the values print as 0.0. I am sure they
really are 0.0.
> Or, is there another/better way to keep diagnosing the problem?
>
I can't think of anything now other than letting me try it on my system,
as suggested above.
You might also see what happens with a workaround:
Znp0 = Znp.copy()
Znp0[Znp==0] = 0.1
Then try plotting Znp0.
Eric
------------------------------------------------------------------------------
Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are
powering Web 2.0 with engaging, cross-platform capabilities. Quickly and
easily build your RIAs with Flex Builder, the Eclipse(TM)based development
software that enables intelligent coding and step-through debugging.
Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
