1,944 reputation
816
bio website staff.acecrc.org.au/~mdsumner
location
age
visits member for 4 years
seen 9 hours ago

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


Sep
16
comment How to open and plot Globcolour chlorophyll image in R?
Well how about a link to a file or some idea of what you've tried already?
Sep
16
comment Drawing smooth great circle?
Yeah, sorry - just fwiw. If you want to use R I will expand this, just happened to be looking at something similar today.
Sep
16
comment Custom projection problem with proj.4
You'll need to give more information. What steps do you take, what do you get, and how do you determine wrongness?
Sep
16
comment Custom projection problem with proj.4
Eh? The op says "vector"
Sep
16
answered Drawing smooth great circle?
Sep
16
comment Drawing smooth great circle?
FYI: bost.ocks.org/mike/example
Sep
15
comment How to create header files for multiple rasters in R
Other options could be to do esentially the same thing with GDAL's VRT, or raster's .grd format, or write R functions to matrix(readBin(...)) etc - they all amount to the same thing. Note that your ENVI version doesn't have georeferencing or projection, but presumably it could have, and the other pathways could/should too.
Sep
15
revised Avoiding R for windows gui front-end has stopped working error?
raster tag is not relevant
Sep
15
answered Avoiding R for windows gui front-end has stopped working error?
Sep
15
suggested suggested edit on Avoiding R for windows gui front-end has stopped working error?
Sep
15
comment Avoiding R for windows gui front-end has stopped working error?
Try with rgdal instead, but if you don't have the system resources it's not going to work. I.e. forest <- readOGR("F:\\NDVI_MODIS_500\\Mod_ndvi_500_shapes", "mod_ndvi_500_grid_point_03_forest_steppe"); 32-bit R is a real limitation. How many points is it? Can you read a subset first? Also, try with QGIS to do this.
Sep
10
comment Guessing a map projection by its values
Curlew: how is UTM easy to identify? There are 60 (bare minimum) different ones, let alone different hemispheres and different national grids at different times in history. It might be data for a different planet or some other purpose entirely. No answer is possible without some knowledge of the data itself: this has to be heuristic and sociological, the very best way is that the person/s responsible provide the actual required information. (It's mad that it's so often neglected).
Sep
10
comment Is there a limitation on the numbers of layers that can be displayed in QGIS?
Consider Manifold over Mapinfo btw.
Sep
10
comment Is there a limitation on the numbers of layers that can be displayed in QGIS?
Fyi sounds like a "max no. open files limit": hub.qgis.org/issues/4013
Sep
9
comment Polygon Between Contour Line
You could just use Manifold to do this ("areas" is a basic contour option) if you're not wedded to ArcGIS.
Sep
3
comment GDAL for Python: extracting subdomains from NetCDF file?
It is actually really straightforward with GDAL. Start by getting the gdalinfo output of your file, and the variable within that you want. Probably it will present as a subdataset in the file.
Sep
3
comment GDAL for Python: extracting subdomains from NetCDF file?
Why reproject? Especially, why reproject to utm? What have you tried already with GDAL? This reads like "here is my task, please do for me"
Sep
3
comment Reprojecting MODIS Ice Sea Temperature products with gdalwarp
"latitude_of_center" looks wrong, and look at the bbox in longlat from GDAL - you'll need to sort out what it's meant to be
Aug
30
comment Kernel Density Analysis using R and ArcGis
Try library(raster); raster(denRedwood) and see ?writeRaster for tiff etc. You'll need to restore the crs which spatstat does not have. See ?projection
Aug
23
comment How to correctly resize raster GIS images to a given px width?
add -co TILED=YES and -co COMPRESS=LZW for more efficient access and storage, if your target usage can utilize those options