2,233 reputation
817
bio website staff.acecrc.org.au/~mdsumner
location
age
visits member for 4 years, 8 months
seen 3 hours ago

R, Manifold, GDAL, also interested in Python, topojson/D3 and NetCDF/HDF


Jun
30
comment rgdal with multiple Feature Datasets in R
I reckon the versions are the key here, and it may even be different systems. When you start library(rgdal) it prints out what is required, and sessionInfo() completes the story (though devtools::session_info() is much better).
Jun
30
comment rgdal with multiple Feature Datasets in R
When you say "others get" are they using the same R and rgdal versions and are you sure it's the same data source?
Jun
29
comment Getting warning when opening NetCDF file using gdal?
the file is this one: oceandata.sci.gsfc.nasa.gov/cgi/getfile/… available here oceandata.sci.gsfc.nasa.gov/MODISA/L2/2015/069 via navigation from oceandata.sci.gsfc.nasa.gov
Jun
29
comment Getting warning when opening NetCDF file using gdal?
This is one for the gdal-devs, point them to details about the format change to netcdf4: oceancolor.gsfc.nasa.gov/DOCS/FormatChange.html Possibly also this exchange could give you some clues, but I'm not hopeful about it: lists.osgeo.org/pipermail/gdal-dev/2015-February/040973.html
Jun
26
comment What coordinate data could this ASCII file represent?
Always. Reproducibility is the aspirational standard.
Jun
26
comment What coordinate data could this ASCII file represent?
No information means no information I'm afraid. Until you get a nice picture there's no way to be sure.
Jun
26
comment Raster alignment terminology?
I would say that if they share the same resolution, alignment, offset, scale and coordinate system then they are the same raster.
Jun
25
comment Creating a new raster layer based on mean of other raster layers, ignoring NA values in R
Well you must not be using this code I guess. Can you provide useable code? Don't use "mean" as a variable name either, it's already a function.
Jun
25
comment Creating a new raster layer based on mean of other raster layers, ignoring NA values in R
What package provides an overlay method for numeric?
Jun
25
comment Heatmap on globe - distortion problems
(If anyone wants to explore this in R I recreate it here, requires a github package plus CRAN ones: gist.github.com/mdsumner/a3d2277bae7cd8082302)
Jun
21
comment Different kinds of time+spatial data
Pebesma has a paper on this: cran.r-project.org/web/packages/spacetime/vignettes/jss816.pdf It's easy to overcook your own defs though, in ways that don't actually work for others
Jun
19
comment What raster-based GIS systems really work?
R's raster package totally kicks serious raster butt, it's not perfect but there's nothing better.
Jun
14
comment Equal Area Projection for Indonesia - Standard Parallels
The standard parallels are where "true scale" lies, and to the north and south distortion increases. Putting them nearer to the centre means that more of the central region is "optimal". Also, see this for conic discussoin generally en.wikisource.org/wiki/The_American_Practical_Navigator/…
Jun
14
comment Grid-type shapefile from raster GeoTIFF
rasterToContours in the R package raster returns cell edges in a SpatialLinesDataFrame.
Jun
14
comment What projection to use for many parts of the world
Is your poster wrapped around a globe?
Jun
13
revised Where can I get GIS data layers that show sea ice data for the Beaufort and Chukchi Seas?
added 2 characters in body
Jun
13
comment Georeferencing reversal of image
What software are you using? What steps did you take, in detail?
Jun
13
comment Processing large amounts of Rasters in R
I think the call is ok, you just want to assign it somewhere, .i.e. winterSumRaster <- .
Jun
12
answered Where can I get GIS data layers that show sea ice data for the Beaufort and Chukchi Seas?
Jun
12
comment Processing large amounts of Rasters in R
setValues(winterSumRaster,getValues(dailyRaster)) - this is not actually persisted anywhere since there's no assigment