Reputation
3,091
Top tag
Next privilege 5,000 Rep.
Approve tag wiki edits
Badges
9 17
Newest
 Yearling
Impact
~126k people reached

Feb
15
comment Crop raster based on another raster R
Try projectRaster(RasA, RasB)
Feb
13
comment Wrong shape display with North_Pole_Azimuthal_Equidistant projection
Can you check your issue against this fix: hub.qgis.org/issues/7596
Feb
11
answered Plotting netCDF4 data with the right projection
Feb
10
comment How to calculate a curve distance between two xyz points?
Is the curve just the curvature due to depth or also the curvature of the Earth? (You might be surprised how little it matters, but all depends on the extent and the data itself).
Feb
9
answered aggregate raster using the most frequent value
Feb
8
comment Plotting Raster in R gives 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.
Feb
8
answered Tearing of SpatialLinesDataFrame and problem with gClip function
Feb
7
comment Tearing of SpatialLinesDataFrame and problem with gClip function
Try raster::erase or raster::intersect
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
2
answered How to extract a 5 bands geotiff file into RGB file?
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?