I have a raster that I downloaded and for whatever reason there is no 'proj4string' when it is read in.
ndvi= raster("av_TOTND1996v4.bsq")
> proj4string(ndvi)
[1] NA
The documentation says that it is georeferenced, and I am wondering if I could create a custom 'proj4string' to define the data. I also have access to ARCGIS, however I am not really familiar with it, and can't figure out how to reset the center of the projection.
The georefernce information provided is:
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
Your best bet would be to try and find the code for the data from whoever supplied you the data. Otherwise you will have to construct the proj4 string, assuming you have all the details...
You can search for spatial reference codes on http://spatialreference.org/ so if your data is for a know geographic area you could perhaps search for that location (eg if its a project on Madagascar, maybe it uses a Madagascan grid system).
That site is also a good source of example proj4 strings that might teach you how to construct one from your info. For example, here's a European Lambert Equal Area conical projection that might be a bit like yours:
http://spatialreference.org/ref/epsg/3035/proj4/
Although your centre longitude of -100 tells me its not this (nor is it Madagascar...)
The HTML version of the projection is a bit more enlightening as to what the parameters mean:
http://spatialreference.org/ref/epsg/3035/html/
Note that a proj4string won't define the size of your grid cells. You'll probably have to do that by setting the extent of your grid. Actually I don't see any information about where in your coordinate system the grid cells are. What does extent(ndvi) tell you? Are you sure there's not a companion .hdr
file you've not got?