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

1d
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.
1d
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.
2d
comment Failure of 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.
Jan
10
comment Identify Hawaiian coordinate system
There are at least 5 State Plane coordinate systems for Hawaii and there's no reason to believe any of them are it. Where did you come across this? The actual source is always the best place to look, the evidence comes in many forms.
Jan
4
answered GDAL “Simple Tiling” (Or yet another - how do I tile question)
Jan
3
comment GDAL “Simple Tiling” (Or yet another - how do I tile question)
I think you could use gdal.org/gdal_retile.html