Barry,
On Wed, Apr 23, 2014 at 4:37 PM, Barry Rowlingson
b.rowling...@lancaster.ac.uk wrote:
On Wed, Apr 23, 2014 at 2:42 PM, Agustin Lobo alobolis...@gmail.com wrote:
I think that the output of get_map commands, i.e.
gmap - get_map(location = rprobextLL@bbox, maptype = hybrid, source
=
(follow up of previous message, I pressed enter as if I were in the R
session and sent the message...)
On Wed, Apr 23, 2014 at 4:37 PM, Barry Rowlingson
b.rowling...@lancaster.ac.uk wrote:
- I imagine the web services expect EPSG:4326.
What they expect, yes. But what about what we retrieve?
EPSG:4326, or WGS84, is only one way to refer to places on the earth by
longitudes and latitudes; it seems that google maps come in a Mercator
projection, EPSG:3857, read e.g. this:
http://docs.openlayers.org/library/spherical_mercator.html
or
But the fact that the output of get_map() includes the bounding box in
geographic
coordinates (epsg:4326) does not imply that the matrix is on the same
projection.
This should be specified.
I can provide you a utm raster along with bounding box coordinates in
geographic coordinates (epsg:4326).
So do you mean (and I've not tried your code yet, sorry) that the corners
may well be at the given lat-long points but half way along any edge might
not be at the halfway-point in lat-long?
if you convert the bbox coords to epsg:3857 then maybe you have a correctly
projected raster in those
Dear all,
I'd like to know if someone knows a easy way of putting the legend inside the
sp::bubble plot like we can do with sp::spplot.
Best,
Alessandro Samuel-Rosa
---
Graduate School in Agronomy - Soil Science
Federal Rural University of Rio de Janeiro
Seropédica, Rio de Janeiro, Brazil
---
A second issue with the web mercator projection epsg:3857 is that some
software seem to store an incorrect proj.4 string for the projection,
leading to shifts when passing objects from one software to another. See:
Dear Rgeos,
Id like to know if there is some workaround to import the lines contained
by a shape file, avoiding memory problems when the only thing you need is
the lines it contains.
I mean, reading the lines of the world's coast lines can be done in a
instant with:
wcoastline -
On Thu, Apr 24, 2014 at 4:47 PM, Agus Camacho agus.cama...@gmail.comwrote:
Dear Rgeos,
Id like to know if there is some workaround to import the lines contained
by a shape file, avoiding memory problems when the only thing you need is
the lines it contains.
I mean, reading the lines of
A possibility is, if you are interested in some specific area of the mata
atlantica, to crop the file externally with OGR just for your study area, and
read the cropped file into R.
Javier
---
__
From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on
behalf of Barry
Thanks for the fast response,
@ Barry
The .shp file is 420,241KB, the .shx is 6,994KB and the .dbf is 1,729,941KB.
@ Javier, id like to use the whole region.
Might be any easy way to reduce the number of vertex previously to import?
These shape files are available at:
On Thu, 24 Apr 2014, Agus Camacho wrote:
Thanks for the fast response,
@ Barry
The .shp file is 420,241KB, the .shx is 6,994KB and the .dbf is 1,729,941KB.
@ Javier, id like to use the whole region.
But you'll never be able to render the detail in this coastline anyway.
Why not just use
Thanks Roger,
I didnt get what you mean with the intermediate binary dataset, you mean
the mata_atlantica.shp file (420,241KB)?
I tried maptools::Rgshhs and also getRgshhsMap but got the error:
Error in Rgshhs(fn, xlim = xl.west, ylim = ylim, shift = shift, level =
level, :
Data not same
I think i might have not made myself clear in first place:
I want the lines from mata_atlantica.shp, but those are much more than
shorelines.
They are inland polygons.
thanks
2014-04-24 14:07 GMT-03:00 Agus Camacho agus.cama...@gmail.com:
Thanks Roger,
I didnt get what you mean with the
On Thu, 24 Apr 2014, Agus Camacho wrote:
I think i might have not made myself clear in first place:
I want the lines from mata_atlantica.shp, but those are much more than
shorelines.
They are inland polygons.
Then you'll need a machine with a lot of memory. The URL you gave does not
lead to
On 04/24/2014 03:52 PM, Barry Rowlingson wrote:
So do you mean (and I've not tried your code yet, sorry) that the
corners may well be at the given lat-long points but half way along any
edge might not be at the halfway-point in lat-long?
Indeed:
library(ggmap)
Dear Agus,
Thanks for sending me the shapefile off list. It uncovered some unicode
bugs in my code.
If you want to try some experimental code (banned even!), you can try
https://github.com/thk686/rgdal2
Here's some output:
x = openOGR(/Users/tkeitt/Desktop/mata_atlantica/)
x
INFO: Open of
Hi everyone!
I have been using makegrid and it worked fine. Until tonight. I upgraded to
R-3.1.0 on my Ubuntu, and am using plotKML and gstat from R-Forge.
Here is my new snag:
u2 - SpatialPoints(makegrid(texb,n=800))
Warning message:
In makegrid(texb, n = 800) : NaNs produced
Error in
18 matches
Mail list logo