设置自定义 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 不会定义网格单元格的大小。您可能必须通过设置网格的范围来做到这一点。实际上，我在网格单元所在的坐标系中看不到有关 where 的任何信息。 extent(ndvi) 告诉你什么？您确定没有您没有的配套 .hdr 文件吗？

设置自定义 proj4string

Setting custom proj4string

r

gis

arcgis

raster

rgdal