[R-sig-Geo] Same mask for two rasters, but obtaining different extents in R

[Ben Tupper]

Hi,

Might you show us what bio5 and lg5 look like before you do the masking?

Cheers,
Ben


> On Mar 14, 2017, at 5:31 PM, Vijay Ramesh via R-sig-Geo <r-sig-geo at r-project.org> wrote:
> 
> I loaded a shapefile in R, and cropped and masked it with the same mask,
> but I get different extents. Any suggestions?
> 
> Code below:
> 
>    library(raster)
>    library(rgdal)
>    library(GISTools)
>    library(sp)
>    library(maptools)
> 
>    ##Loading the first mask
>    Mask <- readOGR("C:\\Users\\rameshv\\Downloads\\Climate
> Stability\\SPV_Mask\\mask.shp")
>    proj4string(Mask) <- CRS("+init=epsg:4326")
> 
>    #Loading the Rasters For Present and LGM and clipping them to the right
> extent
>    bio5 <- raster("C:\\Users\\rameshv\\Downloads\\Climate
> Stability\\Data_LGM_Present\\Present\\1_RawData\\bio_5")
>    proj4string(bio5) <- CRS("+init=epsg:4326")
>    lg5 <- raster("C:\\Users\\rameshv\\Downloads\\Climate
> Stability\\Data_LGM_Present\\LGM\\1_RawData\\cclgmbi5.tif")
>    proj4string(lg5) <- CRS("+init=epsg:4326")
> 
>    ##Cropping by using the Crop and Mask functions
>    cr<- crop(bio5, extent(Mask))
>    bio5 <- mask(cr, Mask)
> 
>    cr2<- crop(lg5, extent(Mask))
>    lg5<- mask(cr2, Mask)
> 
>> bio5
>     class       : RasterLayer
>     dimensions  : 339, 246, 83394  (nrow, ncol, ncell)
>     resolution  : 0.04166667, 0.04166667  (x, y)
>     extent      : 72.45835, 82.70835, 8.083337, 22.20834  (xmin, xmax,
> ymin, ymax)
>     coord. ref. : +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs
> +ellps=WGS84 +towgs84=0,0,0
>     data source : in memory
>     names       : bio_5
>     values      : 207, 432  (min, max)
>     attributes  :
>        ID COUNT
>     from: -86     1
>     to  : 489    38
> 
>> lg5
>     class       : RasterLayer
>     dimensions  : 339, 246, 83394  (nrow, ncol, ncell)
>     resolution  : 0.04166667, 0.04166667  (x, y)
>     extent      : 72.45833, 82.70833, 8.083333, 22.20833  (xmin, xmax,
> ymin, ymax)
>     coord. ref. : +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs
> +ellps=WGS84 +towgs84=0,0,0
>     data source : in memory
>     names       : cclgmbi5
>     values      : 187, 392  (min, max)
> 
> 
> -- 
> Data Manager,
> Barbara Han Lab,
> Cary Institute of Ecosystem Studies,
> 2801 Sharon Turnpike, Millbrook, NY 12545
> Lab Site : http://www.hanlab.science/
> Personal Site : http://evolecol.weebly.com/
> Phone : (845)-677-7600 Ext: 241
> 
> 	[[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

Ben Tupper
Bigelow Laboratory for Ocean Sciences
60 Bigelow Drive, P.O. Box 380
East Boothbay, Maine 04544
http://www.bigelow.org