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

Oct
22
comment Improves performance on R working with raster files
fwiw, try using basename(all.files[i]) rather than substr.
Oct
22
comment How to make a wedge shapefile in R?
Very nice indeed!
Oct
5
answered How to make a wedge shapefile in R?
Sep
17
comment How do you geocode addresses in R?
Also dismo::geocode - nice because dismo also has gmap which gives raster objects.
Sep
17
comment Can you check if a gps coordinate is inside a gps polygon in R?
The high-level ?over or the low-level ?point.in.polygon is what you are after.
Sep
17
comment Why does gdalwarp's cubicspline method expand into nodata areas?
Also, I'm actually looking at some benthic data from this area with a colleague, so I'm genuinely interested.
Sep
17
comment Why does gdalwarp's cubicspline method expand into nodata areas?
Not necessarily, I'd be happy to have a closer look at your workflow - which btw I think your question really needs specifics so we can reproduce the issue. I don't think there's any easy answer! Is this a pretty localized part of the grid, or is it most of it? (Ah I see Freycinet now).
Sep
17
comment Why does gdalwarp's cubicspline method expand into nodata areas?
Warping data like this is treacherous. Do you really need to warp grids? I work with these kinds of data a lot and I always resample/extract rather than warp - and just leave grid reprojection for occasional nice pictures.
Sep
15
comment Loading *.adf files into QGIS using Python
You need to find out what the layer name is, if it is a vector "coverage". Try PAL, ARC, or CNT rather than "layer_name". gdal.org/drv_avcbin.html (Also, for sanity purposes please confirm that it's not a raster ADF according to the notes on that page).
Sep
15
answered Selecting features within polygon from another layer using R or Python?
Sep
12
comment NASA world view (EOSDIS) coverage delay per day - scanning footprint issue?
Search for details on satellite remote sensing coverage. The Earth is rotating, and it's large. High resolution imagery requires that you be relatively close and for a satellite to be close it must be moving (there's only one orbit distance that is "geo stationary", and that's not "close"). There's only so much a satellite can "see", and there's only so many we can launch. It takes time to orbit the Earth.
Sep
12
comment NASA world view (EOSDIS) coverage delay per day - scanning footprint issue?
Because the viewing platform is a finite number of satellites in orbit . . .
Sep
11
comment Contours from delaunay triangles
There is this in R: cran.r-project.org/web/packages/contoureR/index.html
Sep
7
comment How to export a dataset with “SpatialPolygonsDataFrame” as a shapefile
What does Forced mean here?
Sep
3
comment Programs/ languages to process netcdf files “fast”
The main thing (referring to your example code) is to do the extract step only once (outside the i loop) as 'extract(brick,shape, cellnumbers = TRUE)', and apply that each time. Since polygons can be multiple, the cell numbers are organized in a list of matrices, so it needs some extra handling to apply. The other thing is not to rbind() in the loop - preallocate the matrix, or use a list and build the matrix later.
Sep
3
comment Descrambling scrambled coordinates?
There's no coherent question here. It seems like you've figured out the problem, perhaps you just need a way to implement it to these data to check your interpration, and that using these data lon/lat pairs translate to the decimal ones?
Sep
3
answered Programs/ languages to process netcdf files “fast”
Jul
22
comment Merge multiple SpatialPolygonDataFrame into 1 SPDF in R (analogue to “Merge” Tool in ArcGIS)
Super easy solution provided by Robert Hijmans, the author of raster :)
Jul
21
comment Merge multiple SpatialPolygonDataFrame into 1 SPDF in R (analogue to “Merge” Tool in ArcGIS)
See ?rgeos::gUnion and / or ?raster::union
Jul
15
revised How to find the distance from a point to a polygon, via another point, in R
need raster for projection, data( wrld_simpl) too