0

我有一个我下载的光栅，无论出于何种原因，在读入它时都没有“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
4

1 回答 1

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文件吗？

于 2015-03-01T16:43:41.000 回答