#DATA

Download from [this website][1] the following ZIP file: **NUTS_2006_60M_SH.zip** and save all the **NUTS_RG_60M_2006 files**. In what follows I refer to this as "NUTS3 regions".

Download from [this website][2] the following nc-file: **cru_ts3.23.1991.2000.pet.dat.nc.gz** This is a raster with data all over the world.

#Desired output
The nc-file is a raster of the world with values in different colors. 
[![raster data][3]][3]

The goal is now to translate these data to a NUTS3 polygon regions and save it back as a **nuts3-shapefile with the data added to the polyons**.
[![Add data to polygons][4]][4]



#Software
I can work with R, but since I couldn't make it work in R, I switched to ArcGis. If the data are added to the NUTS3 shapefile, I can continue working with it in R.

Any solution in R or ArcGis is welcome, as long as the output is 1 polygon shapefile with all the values added.


#Problem
I can open the nc-file in ArcGis (toolbox, multidimension tools, make netCDF Raster Layer). I can also overlap the nuts3-file with the raster. But I can't merge the data together in one shapefile (by taking the average of the underlying rasters that overlap with one polygon) and drop all the data that I don't need.

I found some examples on how to do this on Google, but only one of them explained how to start from an nc-file. As far as I understood, I first have to right click on the data-layer to export the data as a raster (as it is now only an nc-file and not yet a shapefile). However, this is not working since I get the following error: "Output file format with specified pixel type or number of bands or colormap is not supported."

  [1]: http://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts
  [2]: https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_3.23/cruts.1506241137.v3.23/pet/
  [3]: https://i.sstatic.net/UrfOD.jpg
  [4]: https://i.sstatic.net/ieVO2.jpg