I am new to all things GIS so apologies for a basic questions. I am trying to specify the CRS of a climatic raster layer in R. It is an equal area projection based on Europe Equidistant Conic and the datum is D European 1950. Can anyone help me specifiy the argument in R?

projection(myraster) = CRS("+proj=eqdc +ellps=???+ datum=???")

Many thanks if you can help!

  • Is the map area in the south Permian Basin? – Mintx Dec 12 '13 at 17:36
  • In your other answer you refer to ECA&D E-OBS - where did you find that these data are provided in European 1950 datum? Please indicate the URL or reference... – markusN Dec 12 '13 at 20:22
  • I assume tsutsume wants to make a map in that projection. The data looks like degrees. – AndreJ Dec 12 '13 at 20:44

"Europe Equidistant Conic" is also coded as ESRI:102031 and its Proj.4 definition is:

+proj=eqdc +lat_0=0 +lon_0=0 +lat_1=43 +lat_2=62 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs
| improve this answer | |

QGIS reports European Datum ED50 as

+proj=longlat +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +no_defs

But there are als definitions

IGNF:ED50G +proj=longlat +towgs84=-84.0000,-97.0000,-117.0000 +a=6378388.0000 +rf=297.0000000000000 +units=m +no_defs
ED50(ED77) EPSG:4154 +proj=longlat +ellps=intl +towgs84=-117,-132,-164,0,0,0,0 +no_defs

It mainly depends on what part of Europe your data is from. But I doubt if it has a great impact on the analysis of your data.

| improve this answer | |
  • Wow, thanks for the quick response. It is the whole extent of Europe. The original climatic data set is from the E-OBS gridded climatic data set. I will try your suggestions. – tsutsume Dec 12 '13 at 17:41

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.