@pyimport, nifty.
Thanks!
On Monday, September 21, 2015 at 11:03:35 AM UTC-7, Tom Breloff wrote:
>
> This seemed to match:
>
>
>> using PyPlot, PyCall
>> @pyimport matplotlib.colors as COL
>> @pyimport numpy as np
>> xaxis = linspace(1, 20, 101)
>> yaxis = linspace(2, 5, 101)
>> data = Float64[x+y for x in xaxis, y in yaxis];
>> X, Y = np.meshgrid(xaxis, yaxis)
>> plt = pcolormesh(X, Y, data', norm = COL.LogNorm(vmin=minimum(data),
>> vmax=maximum(data)))
>> colorbar()
>> savefig("test.png")
>
>
>
>
> On Mon, Sep 21, 2015 at 1:34 PM, Jan Strube <[email protected]
> <javascript:>> wrote:
>
>> A small python snippet that does roughly what I want is below.
>> I'm using asymmetric data to remind myself whether to transpose the
>> matrix or not.
>> Thanks for your help.
>>
>>
>> from matplotlib import pyplot as plt
>> from matplotlib.colors import LogNorm
>> import numpy as np
>> import sys
>> xaxis = np.linspace(1, 20, 101)
>> yaxis = np.linspace(2, 5, 101)
>> data = np.zeros((100, 100))
>>
>> for x in range(100):
>> for y in range(100):
>> data[x,y] = xaxis[x] + yaxis[y]
>> X, Y = np.meshgrid(xaxis, yaxis)
>> # I can figure out how to make this figure in Julia
>> #p = plt.pcolormesh(X, Y, data.T)
>>
>> # This is what I would like instead
>> p = plt.pcolormesh(X, Y, data.T, norm=LogNorm(vmin=np.min(data), vmax=np.
>> max(data)))
>> plt.colorbar(p)
>> plt.savefig("test.png")
>>
>>
>>
>>
>>
>> On Sunday, September 20, 2015 at 8:54:26 PM UTC-7, Tom Breloff wrote:
>>>
>>> Yes you should probably use `collect`.
>>>
>>> With regards to plotting... can you post the pyplot code that generates
>>> the graph that you want? We may be able to either show you how to do it in
>>> julia, or it will help in future development by pinpointing a deficiency.
>>> Thanks.
>>>
>>> On Sun, Sep 20, 2015 at 11:11 PM, Jan Strube <[email protected]> wrote:
>>>
>>>> I'm trying to write some data and axis definitions to HDF5 for later
>>>> plotting in pyplot. (Because I haven't figured out how to do lognorm on
>>>> pcolormesh in PyPlot.jl)
>>>> Writing the data - a 2D array - is no problem.
>>>> Writing the axes - linspace(min, max, 100) - doesn't work, because I
>>>> just found out that linspace creates a LinSpace object, not an array, and
>>>> HDF5 doesn't know how to write that.
>>>> My question is: What is an idiomatic way to turn LinSpace into an
>>>> Array? Is collect the recommended way to do this?
>>>>
>>>>
>>>
>
