[R-sig-Geo] Scaling coarse raster to finer raster using overlapping area weights

[Dijk, Michiel van]

Dear all,

Just to inform you that the installation problem I had before is fixed (thanks Jeroen Ooms!). I have a setwd() line in my Rprofile.site which seemed to break the installation (probably should change that). Jeroen somehow created a fix in the package and now the installation works on my machine too (even with the setwd() in  Rprofile.site).

I will try the package after Christmas.
Good holidays!
Michiel 

-----Original Message-----
From: Dijk, Michiel van 
Sent: 21 December 2016 17:19
To: Edzer Pebesma; r-sig-geo at r-project.org
Subject: RE: [R-sig-Geo] Scaling coarse raster to finer raster using overlapping area weights

Hi Edzer,

Thanks for the swift reply. Unfortunately, already with installing the package something goes wrong. I get the following message:

*** arch - i386
rm -f RcppExports.o bbox.o gdal.o gdal_geom.o gdal_read.o gdal_write.o geos.o proj.o wkb.o sf.dll mkdir -p ../inst "C:/PROGRA~1/R/R-33~1.2/bin/i386/Rscript.exe" "../tools/winlibs.R"
cp -r ../windows/gdal2-2.1.1/share/gdal ../inst/
cp: cannot stat '../windows/gdal2-2.1.1/share/gdal': No such file or directory
make: *** [winlibs] Error 1
Warning: running command 'make -f "Makevars.win" -f "C:/PROGRA~1/R/R-33~1.2/etc/i386/Makeconf" -f "C:/PROGRA~1/R/R-33~1.2/share/make/winshlib.mk" CXX='$(CXX1X) $(CXX1XSTD)' CXXFLAGS='$(CXX1XFLAGS)' CXXPICFLAGS='$(CXX1XPICFLAGS)' SHLIB_LDFLAGS='$(SHLIB_CXX1XLDFLAGS)' SHLIB_LD='$(SHLIB_CXX1XLD)' SHLIB="sf.dll" OBJECTS="RcppExports.o bbox.o gdal.o gdal_geom.o gdal_read.o gdal_write.o geos.o proj.o wkb.o"' had status 2
ERROR: compilation failed for package 'sf'
* removing 'C:/Users/vandijkm/R/library/3.3.2/sf

It seems it is unable to find the gdal files. I have gdal installed but not in the .../Windows/ folder. Perhaps that is the problem?
Best,
Michiel

-----Original Message-----
From: R-sig-Geo [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Edzer Pebesma
Sent: 21 December 2016 16:23
To: r-sig-geo at r-project.org
Subject: Re: [R-sig-Geo] Scaling coarse raster to finer raster using overlapping area weights

Hi Michiel,

in my first reply to this list yesterday I mentioned the function sf::st_interpolate_aw (area-weighted) which has an argument "extensive"
that has to be set to TRUE or FALSE; when TRUE, amounts are redistributed proportionality to the target area. It seems you have an extensive variable (amount), so you could use that.

For this you need to install package sf from source, by

devtools::install_github("edzer/sfr")

Both datasets need to be polygon. Let me know if it works, or when you have trouble getting it to work (more likely!)

Best regards,

On 21/12/16 16:06, Dijk, Michiel van wrote:
> Hi,
> 
> I hope you can help me with the following. I have a coarse raster (0.5 x 0.5 degree) with information on global water demand (WD), which I want to combine with a higher resolution (5 arc min) global raster (which I have both in raster and polygon format). The higher resolution raster contains the locations of simulation units (SIMUs) that are related to a model. I would like to redistribute the  global water demand raster over the SIMUs. I first started with the 'resample' function in the raster package but this function allocates the total value for water demand to the cells with smaller resolution, resulting in much too large total water demand. What I am looking for is some sort of weighted disaggregation of values in the water demand raster to the SIMU raster using the share in overlapping area between cells in the two rasters. 
> 
> For example:
> Suppose we have cell A in the water demand raster and this cell contains SIMUa, SIMUb and SIMUc (which can have different sizes). I would like to allocate the following value to SIMUa: Value in cell A x (Area SIMUa)/(Area cell A). Off course the first problem is that it is very well possible that only a part of SIMUa overlaps with Cell A so before I can start I need to (1) find the overlap between the Water demand raster and the SIMU raster and (2) calculate the area of SIMUa that overlaps with cell A. Here I already run into problems. 
> 
> I looked on the web for solutions and came across similar issues but none of it really worked. This post seems to come close:
> http://gis.stackexchange.com/questions/155701/overlaying-grids-with-a-
> polygon-to-determine-the-average-value-of-the-overlappin
> 
> however when I run the gInteresects command to find the overlap  I get the following error: 
> 
> Error in RGEOSBinPredFunc(spgeom1, spgeom2, byid, func) : 
>   rgeos_binpredfunc_prepared: maximum returned dense matrix size 
> exceeded
> 
> I assume that this is a memory issue but I could not find any information on the internet to solve this. Any recommendations to help me solve this particular issue or suggestions on other/better solutions to down scale the water demand raster are welcome. 
> 
> Many thanks,
> Michiel
> 
> 
> M. (Michiel) van Dijk, PhD
> Research scholar | Ecosystems Services and Management (ESM) | 
> International Institute for Applied Systems Analysis (IIASA) Senior 
> researcher (out of office) | International Policy Division (IB) | 
> Wageningen Economic Research
> 
> |Schlossplatz 1 - A-2361 Laxenburg, Austria, room s-138
> S +43 2236 807 537
> Skype: michiel.van.dijk
> S https://www.researchgate.net/profile/Michiel_Van_Dijk/
> / http://nl.linkedin.com/pub/michiel-van-dijk/1a/805/346
> 
> www.iiasa.ac.at
> www.wur.eu/economic-research
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

--
Edzer Pebesma
Institute for Geoinformatics  (ifgi),  University of Münster Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:   http://www.jstatsoft.org/
Computers & Geosciences:   http://elsevier.com/locate/cageo/