<div dir="ltr"><div><div><div>So, the problem here is that what you are asking for is in conflict with the settings. You are giving bounding box that starts at llcrnrlon=lon_inf and goes to urcrnrlon=lon_sup, and that is the extent of the data you are providing, too. But you want to view it in a different bounds. It is possible to get what you want, but it starts to get tricky at this point.<br><br></div>Before venturing further, I should point out to you that Basemap is effectively deprecated. Its support will be terminated in 2020. The CartoPy project, also built on top of matplotlib, is a far superior library, and I think would be much easier for you to use. It doesn't do everything that basemap does yet, but the things that it does do, it does better. I suggest checking it out before we go any further here.<br><br></div>Cheers!<br></div>Ben Root<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 10, 2017 at 2:57 PM, Mauro Cavalcanti <span dir="ltr"><<a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear Ben,<br><br></div>Thanks a lot again for your help and patience. <br><br>After modifying the call to Basemap according to your suggestion:<span class=""><br><br>m = Basemap(llcrnrlat=lat_inf, llcrnrlon=lon_inf, urcrnrlat=lat_sup, urcrnrlon=lon_sup)<br><br></span></div>the grid is correctly displayed where it belongs (thr west coast of Africa). But now I got a map centered at the Pacific basin, with the grid appearing on the lower left corner of it! In fact, See the attached image. I would like a conventional map centered around the Atlantic basin (ie., with center coordinates at lat_0=0 and long_0=0).</div><div><br></div><div>After using the superb Matplotlib for almost a decade (in the beginning, I even got some help here from the legendary John Hunter!), I nonetheless feel somewhat ashamed of having found a potential bug in the library...<br></div><div><br></div>All the best,<br></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-10-10 12:20 GMT-03:00 Benjamin Root <span dir="ltr"><<a href="mailto:ben.v.root@gmail.com" target="_blank">ben.v.root@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Ah, the problem was that by default, the limits for a Basemap goes from -180 to 180, and -90 to 90. If you do:<br><br>m = Basemap(llcrnrlat=lat_inf, llcrnrlon=lon_inf,<br>            urcrnrlat=lat_sup, urcrnrlon=lon_sup)<br><br></div>Then things line up correctly, and you don't need the origin keyword argument, the transpose, or even the extent argument.<br><br></div>In fact, there might even be a bug, as I would have expected specifying the extent should have worked regardless of the original bounds. Might need to look into that.<br><br></div>Cheers!<br></div>Ben Root<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 9, 2017 at 7:54 PM, Mauro Cavalcanti <span dir="ltr"><<a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi,<br><br></div>Thanks for you reply and suggestions.<br><br></div>I changed the imshow call to:<br><br>im = m.imshow(grilla_salida, cmap='summer', extent=(lon_inf, lat_inf, lon_sup, lat_sup), aspect='auto', interpolation='none')<br><br></div>However, the figure is stil wrong (see attachment).<br><br></div>Maybe if instead of imshow, should I use meshgrid/pcolormesh?<br><br></div>Best regards,<br></div><div class="m_2806886750797440327m_3298287048520949079HOEnZb"><div class="m_2806886750797440327m_3298287048520949079h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-10-09 11:29 GMT-03:00 Benjamin Root <span dir="ltr"><<a href="mailto:ben.v.root@gmail.com" target="_blank">ben.v.root@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>First, you shouldn't need to transpose your image... that'll effectively rotate the data by 90 degrees. Second, you didn't specify the extents of your image, so Basemap is putting everything starting at coordinate 0,0 in the default projection.<br><br></div>If you specify the extent keyword argument to imshow as the (lon1, lat1, lon2, lat2) tuple for the lower-left and upper right corners, you won't even need the origin='lower', and you definitely won't need the transpose.<br><br></div>Cheers!<br></div>Ben Root<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Oct 8, 2017 at 2:33 PM, Mauro Cavalcanti <span dir="ltr"><<a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Dear ALL,</div><div><br></div><div>I have a simple dataset of longitudes/latitudes (see the attached csv file). <br><br>From such data, I want to generate a grid like this:<br><br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 1 0 0 0 0 0 0 <br>0 0 0 0 0 0 1 1 2 0 0 0 0 <br>0 0 0 0 0 1 1 1 1 0 0 0 0 <br>0 0 0 1 0 1 0 0 0 0 0 0 0 <br>0 0 0 2 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 1 3 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 1 0 0 0 <br>0 0 0 0 0 0 0 0 0 1 0 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 0 1 0 0 <br>0 0 0 0 0 0 0 0 0 0 1 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0 <br>0 0 0 0 0 0 0 0 0 0 0 0 0<br><br>which gives the number of data records in each cell of the grid, using one of the variables in the dataset ("spp") as a categorical (grouping) factor.<br><br>From this grid, I then want to create a heat map, superimposed on a Matplotlib/Basemap.<br><br>I wrote some code which does what I want (see the attachments).<br><br>It (mostly) works, but te problem is that the grid image is not being displayed correctly: as shown in the attached figure, it appears too small, and in the lower left corner of the map, instead of where it should be (the West coast of Africa, along the Gulf of Guinea). <br><br>Thanks in advance for any assistance you can provide.<br><br></div>Best regards,<span class="m_2806886750797440327m_3298287048520949079m_457455388731623633HOEnZb"><font color="#888888"><span class="m_2806886750797440327m_3298287048520949079m_457455388731623633m_-5485353596679859117HOEnZb"><font color="#888888"><br clear="all"><div><br>-- <br><div class="m_2806886750797440327m_3298287048520949079m_457455388731623633m_-5485353596679859117m_1000018557826917125gmail_signature">Dr. Mauro J. Cavalcanti<br>E-mail: <a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a><br>Web: <a href="http://sites.google.com/site/maurobio" target="_blank">http://sites.google.com/site/m<wbr>aurobio</a><br>"Life is complex. It consists of real and imaginary parts."</div>
</div></font></span></font></span></div><span class="m_2806886750797440327m_3298287048520949079m_457455388731623633HOEnZb"><font color="#888888">
<br>______________________________<wbr>_________________<br>
Matplotlib-users mailing list<br>
<a href="mailto:Matplotlib-users@python.org" target="_blank">Matplotlib-users@python.org</a><br>
<a href="https://mail.python.org/mailman/listinfo/matplotlib-users" rel="noreferrer" target="_blank">https://mail.python.org/mailma<wbr>n/listinfo/matplotlib-users</a><br>
<br></font></span></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><span class="m_2806886750797440327HOEnZb"><font color="#888888"><br>-- <br><div class="m_2806886750797440327m_3298287048520949079m_457455388731623633gmail_signature" data-smartmail="gmail_signature">Dr. Mauro J. Cavalcanti<br>E-mail: <a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a><br>Web: <a href="http://sites.google.com/site/maurobio" target="_blank">http://sites.google.com/site/m<wbr>aurobio</a><br>"Life is complex. It consists of real and imaginary parts."</div>
</font></span></div>
</div></div></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><br>-- <br><div class="m_2806886750797440327gmail_signature" data-smartmail="gmail_signature">Dr. Mauro J. Cavalcanti<br>E-mail: <a href="mailto:maurobio@gmail.com" target="_blank">maurobio@gmail.com</a><br>Web: <a href="http://sites.google.com/site/maurobio" target="_blank">http://sites.google.com/site/<wbr>maurobio</a><br>"Life is complex. It consists of real and imaginary parts."</div>
</div>
</div></div></blockquote></div><br></div>