pcolor + contour

[Puckja]

With 3.0, the following code seems to overlap contour on pcolor correctly:
octave:> pcolor(z);shading flat; hold on
octave:> contour(z);

However, I have a real data set, but just cannot show the contour line with the following code:
h1=pcolor(d1); shading flat; hold on;
[c2,h2]=pcolor(eps,1); % with colormap(jet)

where d1 is a 60x120 intensity map with magnitude ranging [0 5] and eps is same size map with bottom half being 1 and upper half being 2.3 (just like a half plane so I only need one contour to see the separation).

I can not see the contour even by using
set(h2,"edgecolor","white");

BUT, if I do [c3,h3]=contour(d1), then immediately both contour for d1 and eps show up. Then, I can set the edgecolor of h2 to whatever I want.

This may be a workaround if I can delete the contour of d1 (h3).  I tried to delete(h3), but it didn't work either.  So I am stuck.

Thanks for your help!!

Puck

[David Bateman]

Pukeja wrote:

> With 3.0, the following code seems to overlap contour on pcolor correctly:
> octave:> pcolor(z);shading flat; hold on
> octave:> contour(z);
>
> However, I have a real data set, but just cannot show the contour line with
> the following code:
> h1=pcolor(d1); shading flat; hold on;
> [c2,h2]=pcolor(eps,1); % with colormap(jet)
>
> where d1 is a 60x120 intensity map with magnitude ranging [0 5] and eps is
> same size map with bottom half being 1 and upper half being 2.3 (just like a
> half plane so I only need one contour to see the separation).
>
> I can not see the contour even by using
> set(h2,"edgecolor","white");
>
> BUT, if I do [c3,h3]=contour(d1), then immediately both contour for d1 and
> eps show up. Then, I can set the edgecolor of h2 to whatever I want.
>
> This may be a workaround if I can delete the contour of d1 (h3).  I tried to
> delete(h3), but it didn't work either.  So I am stuck.
>
> Thanks for your help!!
>
> Puck
>  

gnuplot doesn't do proper hidden line removal in this case, though there
appears to be a patch applied to the CVS version to address this. So
this is really a gnuplot issue and we can't do much about it in octave..
Try the latest CVS version of gnuplot and see if that helps..

D.


--
David Bateman                                David.Bateman@...
Motorola Labs - Paris                        [phone] (Ph)
Parc Les Algorithmes, Commune de St Aubin    [phone] (Mob)
91193 Gif-Sur-Yvette FRANCE                  [phone] (Fax)

The information contained in this communication has been classified as:

[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary

_______________________________________________
Help-octave mailing list
Help-octave@...
https://www.cae.wisc.edu/mailman/listinfo/help-octave

[Puckja]

Puckja wrote:

With 3.0, the following code seems to overlap contour on pcolor correctly:
octave:> pcolor(z);shading flat; hold on
octave:> contour(z);

However, I have a real data set, but just cannot show the contour line with the following code:
h1=pcolor(d1); shading flat; hold on;
[c2,h2]=pcolor(eps,1); % with colormap(jet)

where d1 is a 60x120 intensity map with magnitude ranging [0 5] and eps is same size map with bottom half being 1 and upper half being 2.3 (just like a half plane so I only need one contour to see the separation).

I can not see the contour even by using
set(h2,"edgecolor","white");

BUT, if I do [c3,h3]=contour(d1), then immediately both contour for d1 and eps show up. Then, I can set the edgecolor of h2 to whatever I want.

This may be a workaround if I can delete the contour of d1 (h3).  I tried to delete(h3), but it didn't work either.  So I am stuck.

Thanks for your help!!

Puck

OK, I find the way to delete the h3 contour accidentally.  h3 is actually an array so by doing delete(h3(1)) the unwanted contour of d1 is deleted and leave the one I want.

It may be solve more gracefully by upgrade to new gnuplot, but if I don't have time to try various gnuplot to find out yet.

Anyone has found out which version of gnuplot solves this issue?  

Thanks!

Puck

[David Bateman]

Puckja wrote:

>
> Puckja wrote:
>  
>> With 3.0, the following code seems to overlap contour on pcolor correctly:
>> octave:> pcolor(z);shading flat; hold on
>> octave:> contour(z);
>>
>> However, I have a real data set, but just cannot show the contour line
>> with the following code:
>> h1=pcolor(d1); shading flat; hold on;
>> [c2,h2]=pcolor(eps,1); % with colormap(jet)
>>
>> where d1 is a 60x120 intensity map with magnitude ranging [0 5] and eps is
>> same size map with bottom half being 1 and upper half being 2.3 (just like
>> a half plane so I only need one contour to see the separation).
>>
>> I can not see the contour even by using
>> set(h2,"edgecolor","white");
>>
>> BUT, if I do [c3,h3]=contour(d1), then immediately both contour for d1 and
>> eps show up. Then, I can set the edgecolor of h2 to whatever I want.
>>
>> This may be a workaround if I can delete the contour of d1 (h3).  I tried
>> to delete(h3), but it didn't work either.  So I am stuck.
>>
>> Thanks for your help!!
>>
>> Puck
>>
>>    
>
> OK, I find the way to delete the h3 contour accidentally.  h3 is actually an
> array so by doing delete(h3(1)) the unwanted contour of d1 is deleted and
> leave the one I want.
>
> It may be solve more gracefully by upgrade to new gnuplot, but if I don't
> have time to try various gnuplot to find out yet.
>
> Anyone has found out which version of gnuplot solves this issue?  
>
> Thanks!
>
> Puck
>
>  

I think no version, except maybe the CVS of gnuplot solve this issue as
its to do with the hidden surface removal capabilities of gnuplot. There
was a gnuplot patch discussed on the gnuplot lists in December that
addressed this and perhaps its been applied, though I don't follow the
gnuplot lists and can't say.

D.


--
David Bateman                                David.Bateman@...
Motorola Labs - Paris                        [phone] (Ph)
Parc Les Algorithmes, Commune de St Aubin    [phone] (Mob)
91193 Gif-Sur-Yvette FRANCE                  [phone] (Fax)

The information contained in this communication has been classified as:

[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary

_______________________________________________
Help-octave mailing list
Help-octave@...
https://www.cae.wisc.edu/mailman/listinfo/help-octave

[Puckja]

David:

Could you please elaborate how to get Gnuplot CVS version a little bit more?  I tried to get a version by the following commands (from http://gnuplot.sourceforge.net/development/index.html#DownloadCVS)
  cvs -d:pserver:[e-mail]:/cvsroot/gnuplot login
  cvs -z3 -d:pserver:[e-mail]:/cvsroot/gnuplot co -P gnuplot

But when I tried to compiled it (first of all there is no ./configure in the directory), I got some error and cannot proceed:
Linux:/src/gnuplot# ./prepare
./prepare: 47: aclocal: not found

Some part of the preparation process failed.
Please refer to INSTALL for details.

How did I fix it?

Also, I have heard "mouse zooming" is possible with new gnuplot.  Is it true?  I definitely would like to have this capability (not just scale the image of the plot, but the axis scaled in as well.)

Thanks!

Puck

David Bateman wrote:

Pukeja wrote:
> With 3.0, the following code seems to overlap contour on pcolor correctly:
> octave:> pcolor(z);shading flat; hold on
> octave:> contour(z);
>
> However, I have a real data set, but just cannot show the contour line with
> the following code:
> h1=pcolor(d1); shading flat; hold on;
> [c2,h2]=pcolor(eps,1); % with colormap(jet)
>
> where d1 is a 60x120 intensity map with magnitude ranging [0 5] and eps is
> same size map with bottom half being 1 and upper half being 2.3 (just like a
> half plane so I only need one contour to see the separation).
>
> I can not see the contour even by using
> set(h2,"edgecolor","white");
>
> BUT, if I do [c3,h3]=contour(d1), then immediately both contour for d1 and
> eps show up. Then, I can set the edgecolor of h2 to whatever I want.
>
> This may be a workaround if I can delete the contour of d1 (h3).  I tried to
> delete(h3), but it didn't work either.  So I am stuck.
>
> Thanks for your help!!
>
> Puck
>  
gnuplot doesn't do proper hidden line removal in this case, though there
appears to be a patch applied to the CVS version to address this. So
this is really a gnuplot issue and we can't do much about it in octave..
Try the latest CVS version of gnuplot and see if that helps..

D.


--
David Bateman                                [e-mail]
Motorola Labs - Paris                        [phone] (Ph)
Parc Les Algorithmes, Commune de St Aubin    [phone] (Mob)
91193 Gif-Sur-Yvette FRANCE                  [phone] (Fax)

The information contained in this communication has been classified as:

[x] General Business Information
[ ] Motorola Internal Use Only
[ ] Motorola Confidential Proprietary

_______________________________________________
Help-octave mailing list
[e-mail]
https://www.cae.wisc.edu/mailman/listinfo/help-octave

[John W. Eaton]

On 31-Jan-2008, Puckja wrote:

| Could you please elaborate how to get Gnuplot CVS version a little bit more?

Isn't that really a question for the gnuplot developers?

jwe
_______________________________________________
Help-octave mailing list
Help-octave@...
https://www.cae.wisc.edu/mailman/listinfo/help-octave

[zaphodD2]

Puckja wrote:

I got some error and cannot proceed:
Linux:/src/gnuplot# ./prepare
./prepare: 47: aclocal: not found

Some part of the preparation process failed.
Please refer to INSTALL for details.

How did I fix it?

Theres a couple of reasons this error might be appearing. A first guess would be, you don't have automake installed. Give installing that a go and see how you get on.

Maybe post back if your problem persists?