I am using leaflet and the R library sp. The CRS from my leaflet map is EPSG3857 - WGS84 Web Mercator. I am retrieving the coordinates of my current map view with getBounds(). For further calculations I want to create a polygon from these coordinates. That works fine until I want to assign the CRS.
When I use this
+proj=longlat +datum=WGS84 +nodefs +ellps=WGS84 +towgs84=0.0.0 +init=epsg:3857 my Polygon is a square, which fits my current map view. I check that visually through plotting the Polygon. The plot is a square. That fits since my map view is a square too. When I use the proj4 I obtain from http://spatialreference.org/ref/sr-org/7483/ my polygon is shaped rectangular.
+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs How comes the two are different? Shouldnt when I use the epsg:3857 it all be the same?! Am I missing something crucial here?
The values of the corners of my Polygon (which come back from the getBounds() function) are
> coordinates [,1] [,2] [1,] 2.570801 46.89023 [2,] 17.424316 55.15377 [3,] 17.424316 46.89023 [4,] 2.570801 55.15377
أكثر...
When I use this
+proj=longlat +datum=WGS84 +nodefs +ellps=WGS84 +towgs84=0.0.0 +init=epsg:3857 my Polygon is a square, which fits my current map view. I check that visually through plotting the Polygon. The plot is a square. That fits since my map view is a square too. When I use the proj4 I obtain from http://spatialreference.org/ref/sr-org/7483/ my polygon is shaped rectangular.
+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs How comes the two are different? Shouldnt when I use the epsg:3857 it all be the same?! Am I missing something crucial here?
The values of the corners of my Polygon (which come back from the getBounds() function) are
> coordinates [,1] [,2] [1,] 2.570801 46.89023 [2,] 17.424316 55.15377 [3,] 17.424316 46.89023 [4,] 2.570801 55.15377
أكثر...