Радиус сферы в строке proj4 6371116, а после showWKT 6370997.
Делал так:
Код: Выделить всё
library(sp)
library(rgdal)
OurProj <- "+proj=laea +lat_0=90 +lon_0=0 +ellps=sphere +a=6371116 +b=6371116 +units=m +no_defs"
print(showWKT(OurProj))
PROJCS[\"Lambert_Azimuthal_Equal_Area\",GEOGCS[\"GCS_Normal Sphere (r=6370997)\",DATUM[\"D_unknown\",SPHEROID[\"sphere\",6370997,0]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.017453292519943295]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_origin\",90],PARAMETER[\"central_meridian\",0],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1]]"