Reputation
2,838
Top tag
Next privilege 3,000 Rep.
Cast close & reopen votes
Badges
9 17
Newest
 Yearling
Impact
~115k people reached

3h
comment Plotting Raster in R give error
Check for sure that the right paired file/s exist, this must have been partnered with a .hdr file - make sure it's there too. Are you sure you didn't copy the .bil from somewhere, without also copying the .hdr? (There may be other auxiliarly files). I'm guessing this driver is the one involved (via rgdal): gdal.org/frmt_various.html#EHdr Also best not to use "setwd", but rather use file.path() to create full paths.
19h
answered Tearing of SpatialLinesDataFrame and problem with gClip function
1d
comment Tearing of SpatialLinesDataFrame and problem with gClip function
Try raster::erase or raster::intersect
Feb
5
comment Spatial index in R - Performance enhancement on geospatial processing
rgeos has some of these facilities, see ?rgeos::poly_findInBoxGEOS - otherwise put your loot in Postgres/PostGIS or similar, do the work there and read from the db back to R.
Feb
5
comment Project NetCDF file using GDAL or QGis
I don't see any error for the process, so did the output not work? (the lookup UOM thing probably doesn't matter). Your "RAIN_MM_1hour_ac_proj.nc" is a GeoTIFF, since you didn't specify the output format with "-of NetCDF" (GDAL doesn't care about extensions, but helpfully reminded you). Why are you trying to warp a file that's pretty clearly in Long-Lat to Long-Lat? I don't understand why you have any work to do here (??). You maybe want to gdal_translate with -a_srs argument to fill in the metadata.
Feb
4
comment Failure of the code calling the function “rasterize” in R: NAs or no slots
r.ext does not have the the same extent as r.ext <- raster(extent(shp.prj)) - so that presumably is not the actual code involved here. (You could just do raster(shp.prj) to avoid the need for projection()<-). Otherwise I would try rasterize(shp.prj, r.ext, field=shp.prj$land_cover_type)
Feb
3
comment create a grid and extract species and climate data
We need much more information in order to help you. With the raster package it's easy to create such a grid, and extract etc. but you'll need some idea of how it works: see the vignette with library(raster); vignette("Raster")
Feb
2
answered How to extract a 5 bands geotiff file into RGB file?
Feb
2
comment What geospatial analysis can be done easily in R but not in Python?
but, you might start with "install.packages('rgdal'); install.packages('rgeos')" and take it from there.
Feb
2
comment What geospatial analysis can be done easily in R but not in Python?
"Easily" is too subjective. :)
Jan
26
comment Spplot SpatialPolygonDataFrame and SpatialGridDataFrame
I think you maybe trying to merge these objects(?) or convert the grid to polygons(?) when all you want is to plot. Does this work? image(SGDF); plot(SPDF, add = TRUE); Otherwise as Spacedman says, please just show us the code
Jan
26
answered High res monthly satellite maps including ice cover
Jan
22
answered Why is geography faster than geometry for distance calculations?
Jan
22
answered How can I convert a large csv of irregular lat/lons to regular raster GRID?
Jan
21
comment R calculate raster cell size in map units
There's no point to convert to polygon after warping since every pixel is the same size in that coordinate system, by definition. It's impossible to answer this generally with a concise solution since there are so many details that matter about where and how extensive etc., but it's more sensible to convert to polygon from the source raster, and use the n option in rasterToPolygons, and then project, to ensure sufficient vertices to carry the polygons into the new coordinate system.
Jan
20
answered R calculate raster cell size in map units
Jan
18
answered Gdal_translate command for SDS having spaces in its name
Jan
15
comment How to compute max and min coordinates for multiple polygons?
Nice! ---------------------
Jan
14
answered How to compute max and min coordinates for multiple polygons?
Jan
10
comment Same coordinate system - shapefiles don't overlap
"Same coordinate system" is the boner here. I don't think it's true.