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@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@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 _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo