Tuesday 18 June 2019

coordinate system - USA Albers Equal Area Conic USGS version (LANDFIRE) projection not read by R?



I am working with the North America Forest Dynamics dataset. As mentioned it uses the USA Albers Equal Area USGS version projection, which is primarily different from the USA Albers Equal Area projection in the measure of the latitude of origin (former is 23 and latter is 37.5).


The datum is clearly read by Arc as D_NorthAmerica 1983. which is the datum. However, this datum information is not read by R (using raster package). Hence, resulting analysis and written rasters have unknown datum when processing is done in R.


Interestingly, the famous spatial reference website does not have a Proj4 string info for this, which I could manually input into the crs() command. But on the same website, the .Prj file can be downloaded and on opening this ile as a text pad file a datum i.e. North America 1983 datum and hence transformation is mentioned.


Why is R not recognizing the datum around which the transformation has to happen for this projection?




No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...