| map_set without any grid lines [message #88783] |
Wed, 18 June 2014 16:36  |
audrey.schaufelberger
Messages: 10
Registered: June 2014
|
Junior Member |
|
|
I have a very simple question: I want to use map_set but I don't want it to display any grid lines. I tried setting glinethick to 0 and or spacing the glines very far apart, but nothing seems to work. Does anyone have any suggestions?
best,
Audrey
|
|
|
|
|
|
| Re: map_set without any grid lines [message #88785 is a reply to message #88784] |
Wed, 18 June 2014 17:55  |
robert.vorburger
Messages: 2
Registered: June 2014
|
Junior Member |
|
|
I didn't... here's the corresponding excerpt from my code:
map_set, /orthographic, -53, 191, 0, $ azimuthal, lambert
/grid, latdel=30, londel=30, /isotropic
clr_0= bytscl(map, max=max(map), min=min(map), top=254)
for i=0,size_map[1]-1 do begin
for j=0,size_map[2]-1 do begin
polyfill, i + [-0.5,0.5,0.5,-0.5], $
j-90 + [-0.5,-0.5,0.5,0.5], $
color=clr_0[i, j]
endfor
endfor
map_grid, latdel=10, londel=10, glinethick=1, glinestyle=2
map_grid, latdel=30, londel=30, glinethick=4, glinestyle=0
map_set, /orthographic, 0, 0, 0, $
glinethick=0, /isotropic, /advance
the second map_set plots a grid even without the /grid keyword...
p.s: the idea behind this is to overplot an ellipse for which I need to reset the map, because I want to plot the ellipse on the southern hemisphere of the globe but don't want it to be 'distorted'...
p.p.s: I am using the ellipse function you give to one of the questions on your website. thank you for that!
|
|
|
|
| Re: map_set without any grid lines [message #88786 is a reply to message #88785] |
Wed, 18 June 2014 18:19 |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
robert.vorburger@gmail.com writes:
>
> I didn't... here's the corresponding excerpt from my code:
>
> map_set, /orthographic, -53, 191, 0, $ azimuthal, lambert
> /grid, latdel=30, londel=30, /isotropic
>
> clr_0= bytscl(map, max=max(map), min=min(map), top=254)
>
>
> for i=0,size_map[1]-1 do begin
> for j=0,size_map[2]-1 do begin
> polyfill, i + [-0.5,0.5,0.5,-0.5], $
> j-90 + [-0.5,-0.5,0.5,0.5], $
> color=clr_0[i, j]
> endfor
> endfor
>
> map_grid, latdel=10, londel=10, glinethick=1, glinestyle=2
> map_grid, latdel=30, londel=30, glinethick=4, glinestyle=0
>
> map_set, /orthographic, 0, 0, 0, $
> glinethick=0, /isotropic, /advance
>
> the second map_set plots a grid even without the /grid keyword...
>
> p.s: the idea behind this is to overplot an ellipse for which I need to reset the map, because I want to plot the ellipse on the southern hemisphere of the globe but don't want it to be 'distorted'...
>
> p.p.s: I am using the ellipse function you give to one of the questions on your website. thank you for that!
Well, then, don't set the GLINETHICK keyword!
This code doesn't draw a grid on the second plot:
window, xs=800, ys=400
!p.Multi=[0,2,1]
map_set, /orthographic, -53, 191, 0, $
/grid, latdel=30, londel=30, /isotropic
map_grid, latdel=10, londel=10, glinethick=1, glinestyle=2
map_grid, latdel=30, londel=30, glinethick=4, glinestyle=0
map_set, /orthographic, 0, 0, 0, /isotropic, /advance
map_continents
!P.Multi=0
END
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thue. ("Perhaps thos speakest truth.")
|
|
|
|
|
|
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
