设置自定义 proj4string

Question

我下载了一个光栅，但由于某种原因，在读入时没有“proj4string”。

ndvi= raster("av_TOTND1996v4.bsq")
> proj4string(ndvi)
[1] NA

文档说它是地理参考的，我想知道是否可以创建一个自定义的“proj4string”来定义数据。我也可以访问 ARCGIS，但是我对它不是很熟悉，也不知道如何重置投影的中心。

提供的地理参考信息是：

Spatial_Reference_Information:
Horizontal_Coordinate_System_Definition:
Planar:
Map_Projection:
Map_Projection_Name: Lambert Azimuthal Equal Area
Lambert_Azimuthal_Equal_Area:
Longitude_of_Projection_Center: -100.000000
Latitude_of_Projection_Center: 45.000000
False_Easting: 0.000000
False_Northing: 0.000000
Planar_Coordinate_Information:
Planar_Coordinate_Encoding_Method: row and column
Coordinate_Representation:
Abscissa_Resolution: 1000.000000
Ordinate_Resolution: 1000.000000
Planar_Distance_Units: meters
Geodetic_Model:
Horizontal_Datum_Name: D_Sphere_ARC_INFO
Ellipsoid_Name: Sphere_ARC_INFO
Semi-major_Axis: 6370997.000000
Denominator_of_Flattening_Ratio: infinity

Answer 1

最好的办法是尝试从向您提供数据的人那里找到数据的代码。否则你将不得不构造 proj4 字符串，假设你有所有的细节......

您可以在http://spatialreference.org/ 上搜索空间参考代码，因此如果您的数据是针对已知地理区域的，您或许可以搜索该位置（例如，如果它是马达加斯加的一个项目，它可能使用马达加斯加网格系统）。

该站点也是 proj4 字符串示例的良好来源，可以教您如何根据您的信息构建字符串。例如，这是一个可能与您的有点相似的欧洲兰伯特等面积圆锥投影：

http://spatialreference.org/ref/epsg/3035/proj4/

虽然你的中心经度 -100 告诉我不是这个（也不是马达加斯加......）

投影的 HTML 版本更能说明参数的含义：

http://spatialreference.org/ref/epsg/3035/html/

请注意，proj4string 不会定义网格单元的大小。您可能必须通过设置网格的范围来做到这一点。实际上，我在坐标系中看不到任何有关网格单元格所在位置的信息。范围（ndvi）告诉你什么？你确定没有你没有的伴侣.hdr文件吗？