Question

我有一个我下载的光栅，无论出于何种原因，没有＆＃39; proj4string＆＃39;当它被读入时。

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

文档说它是地理参考的，我想知道我是否可以创建一个自定义＆＃39; proj4string＆＃39;定义数据。我也可以访问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

1 个答案: