自从更新sp包以来，我通过调用sp:：CRS定义得到了一个警告

自从更新sp软件包后，我收到一条以前没有收到的警告：

 1: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded ellps unknown in CRS definition: +proj=stere +lat_0=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +R=6370040 +units=km +no_defs
2: In showSRID(uprojargs, format = "PROJ", multiline = "NO") :
  Discarded datum unknown in CRS definition

这是我用来定义CRS的代码

 sp::CRS(
      "+proj=stere +lat_0=90 +lat_ts=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +a=6370040 +b=6370040 +to_meter=1000 +no_defs"
    )

我尝试创建一个投影，但它以前工作过：） 但自从更新后，我得到了这个警告。 我做错了什么？ 提前谢谢大家,

致以最良好的祝愿，

---

Andreas

在最近的版本中，sp和其他空间软件包发生了一些重大变化。例如，见：

及

---

与您的问题相关的是，crs和crs之间存在差异。在你的定义中，你有参数a和b，它们指的是地球的半径。因为它们是相同的，所以您使用的是球体。PROJ库的新版本要求，如果使用proj4符号，则必须使用WGS84（或密切相关的）基准。否则，您需要提供EPSG代码或WKT表示法。比如说：

p4 <- "+proj=stere +lat_0=90 +lat_ts=90 +lon_0=10 +k=0.93301270189 +units=km +datum=WGS84"
sp::CRS(p4)
#CRS arguments:
# +proj=stere +lat_0=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +units=km  +datum=WGS84 +units=m +no_defs 

---

p4谢谢您的回复。我还发现，有一个重大的更新，我只是不能掌握我现在的问题。我想也许有人会马上看到这个问题。嗨，罗伯特，谢谢你的回复。正在读取Radolan数据：。to_meter=1000是一个单元的大小。我认为a和b定义了总面积。如果我使用你的建议，我会得到以下警告：在rgdal：：rawTransform（projto_int，projfrom，nrow（xy），xy[，1]，：3376620投影点由于我省略了
units=km
，所以没有完成提示。我在“是”中添加了这一点，这似乎删除了警告。奇怪的是，我现在在绘图方面得到了不同的结果。要么是我以前做错了，要么是现在出了问题。
wkt <- 'PROJCS["unknown",
GEOGCS["unknown",
    DATUM["unknown",
        SPHEROID["unknown",6370040,0]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]]],
PROJECTION["Polar_Stereographic"],
PARAMETER["latitude_of_origin",90],
PARAMETER["central_meridian",10],
PARAMETER["scale_factor",0.93301270189],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["kilometre",1000,
    AUTHORITY["EPSG","9036"]],
AXIS["Easting",SOUTH],
AXIS["Northing",SOUTH]]'

But that still does not work because of the datum being unknown 

CRS(SRS_string=wkt)

#CRS arguments:
# +proj=stere +lat_0=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +R=6370040 +units=km +no_defs 
#Warning messages:
#1: In showSRID(SRS_string, format = "PROJ", multiline = "NO") :
#  Discarded ellps unknown in CRS definition: +proj=stere +lat_0=90 +lon_0=10 +k=0.93301270189 +x_0=0 +y_0=0 +R=6370040 +units=km +no_defs
#2: In showSRID(SRS_string, format = "PROJ", multiline = "NO") :
#  Discarded datum unknown in CRS definition