小编mod*_*rog的帖子

我正在尝试使用R和raster包更改光栅投影。输入栅格投影为兰伯特方位角；参数在这里：

Coordinate System:
Lambert_Azimuthal_Equal_Area
False_Easting: 4321000,000000
False_Northing: 3210000,000000
Central_Meridian: 10,000000
Latitude_Of_Origin: 52,000000
GCS_ETRS_1989
Datum: D_ETRS_1989
Prime Meridian: 0


PROJCS
 ["ETRS_1989_LAEA",
   GEOGCS ["GCS_ETRS_1989",
           DATUM ["D_ETRS_1989",
                  SPHEROID ["GRS_1980",6378137.0,298.257222101]],
           PRIMEM["Greenwich",0.0],
          UNIT["Degree",0.0174532925199433]],
   PROJECTION["Lambert_Azimuthal_Equal_Area"],
   PARAMETER["False_Easting",4321000.0],
   PARAMETER["False_Northing",3210000.0],
   PARAMETER["Central_Meridian",10.0],
   PARAMETER["Latitude_Of_Origin",52.0],
   UNIT["Meter",1.0]]

我需要将它们转换为 ESRI ASCII 格式的简单栅格，使用经度和纬度坐标，使用墨卡托式投影，像元大小为 0.1 度（我希望我能很好地解释自己，因为我没有足够的 GIS技能，抱歉）。我需要的是格式.ASC中的栅格，其中栅格的每个值对应于 size 的单个单元格N x N，其中N以度为单位（例如 0.1 度），栅格坐标以经度/纬度为单位。

我尝试在 中使用raster库R，并遵循为该projectRaster函数找到的示例。但是经过多次尝试使用多个参数后，我无法正确制作。我想我没有使用正确的参数进行投影、数据或类似的东西。

这是我尝试过的。我在 R 中加载栅格，然后使用以下方法设置其投影：

>crs(r)<-"+proj=laea +lat_1=52 +lon_0=-10 +ellps=GRS80"

然后我定义输出投影并尝试转换并保存：

>newproj <- "+proj=lonlat +lat_1=52 +lon_0=-10 +ellps=WGS84"
>pr2 …