[R-sig-Geo] projectRaster across the antimeridian

Bacou, Melanie mel at mbacou.com
Wed Jul 6 03:03:49 CEST 2016


Same behavior on a Linux server with R 3.3.0 and rgdal 1.1-8. Session 
details below. The reprojected raster is missing it's western half.

```

R version 3.3.0 (2016-05-03)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: CentOS release 6.8 (Final)

locale:
  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C LC_TIME=en_US.UTF-8
  [4] LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8 
LC_MESSAGES=en_US.UTF-8
  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 
LC_IDENTIFICATION=C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods base

other attached packages:
[1] rgdal_1.1-8  raster_2.5-2 sp_1.2-3

loaded via a namespace (and not attached):
[1] tools_3.3.0     curl_0.9.7      Rcpp_0.12.4 grid_3.3.0      
lattice_0.20-33

```


On 7/5/2016 8:51 PM, Bacou, Melanie wrote:
> Here is my own sessionInfo (I do have problems using `projectRaster()` 
> across the antemeridian even with R 3.2.1 and rgdal 1.1-10, but 
> `gdalwarp` works for me).
>
> ```
>
> R version 3.2.1 (2015-06-18)
> Platform: x86_64-w64-mingw32/x64 (64-bit)
> Running under: Windows 7 x64 (build 7601) Service Pack 1
>
> locale:
> [1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United 
> States.1252
> [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
> [5] LC_TIME=English_United States.1252
>
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods base
>
> other attached packages:
> [1] rgdal_1.1-10 raster_2.5-8 sp_1.2-3     rj_2.0.5-2
>
> loaded via a namespace (and not attached):
> [1] rj.gd_2.0.0-1   Rcpp_0.12.5     grid_3.2.1 lattice_0.20-33
>
> ```
>
>
> On 7/5/2016 2:25 PM, Roger Bivand wrote:
>> On Tue, 5 Jul 2016, Fischbach, Anthony wrote:
>>
>>> By stepping back to the previous r version (3.2.5) on my workstation
>>> (sessionInfo and rgdal library call printed below),
>>> I am able to obtain the results indicated by Benoit.
>>
>> Could you please repeat the sessionInfo and rgdal startup messages 
>> for the failing R 3.3.1 run? Do you know whether the issue arose 
>> between R 3.2.5 and R 3.3.0, or R 3.3.0 and R 3.3.1 (for the same 
>> version of rgdal)? The C and C++ compilers used < 3.3.0 and >= 3.3.0 
>> changed, so there is a difference in the way that the binaries were 
>> built on Windows. I don't have easy access to Windows, but could 
>> check further if you could provide this information.
>>
>> Roger
>>
>>>
>>> Somewhere in this new r version and rasterProject function, we are 
>>> getting
>>> a problem.
>>>
>>> -Tony
>>> ---
>>> sessionInfo()
>>> R version 3.2.5 (2016-04-14)
>>> Platform: x86_64-w64-mingw32/x64 (64-bit)
>>> Running under: Windows 7 x64 (build 7601) Service Pack 1
>>>
>>> locale:
>>> [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United
>>> States.1252
>>> [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
>>>
>>> [5] LC_TIME=English_United States.1252
>>>
>>> attached base packages:
>>> [1] stats     graphics  grDevices utils     datasets  methods base
>>>
>>> other attached packages:
>>> [1] raster_2.5-2 sp_1.2-3
>>>
>>> loaded via a namespace (and not attached):
>>> [1] rgdal_1.1-10    Rcpp_0.12.5     grid_3.2.5 lattice_0.20-33
>>>> require(rgdal)
>>> Loading required package: rgdal
>>> rgdal: version: 1.1-10, (SVN revision 622)
>>> Geospatial Data Abstraction Library extensions to R successfully loaded
>>> Loaded GDAL runtime: GDAL 2.0.1, released 2015/09/15
>>> Path to GDAL shared files: C:/Program 
>>> Files/R/R-3.2.5/library/rgdal/gdal
>>> GDAL does not use iconv for recoding strings.
>>> Loaded PROJ.4 runtime: Rel. 4.9.1, 04 March 2015, [PJ_VERSION: 491]
>>> Path to PROJ.4 shared files: C:/Program 
>>> Files/R/R-3.2.5/library/rgdal/proj
>>> Linking to sp version: 1.2-3
>>>>
>>> --- END ---
>>>
>>>     [[alternative HTML version deleted]]
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20160705/903d3df3/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nfdkmalcofanekeb.png
Type: image/png
Size: 27097 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20160705/903d3df3/attachment.png>


More information about the R-sig-Geo mailing list