Re: [R-sig-Geo] Correspondence between gstat and ArcGIS functionality?

/gap_filling_software.html Alternatively GRASS GIS has a function / module called r.fillnulls which uses splines interpolation to fill the gaps. http://grass.itc.it/grass65/manuals/html65_user/r.fillnulls.html Hope this helps, Wesley Wesley Roberts MSc. Researcher: Earth Observation Natural Resources &

Re: [R-sig-Geo] Advice on clustering LiDAR point clouds

ops. > > To get maximum values, maybe sp::overlay can get you started, but I > think it is more to get grid value to point than reverse. I do it in > ArcGIS, but I would be happy to have a solution in R if you find one. > > Etienne Also see r.in.xyz within GRASS GIS. Dylan >

[R-sig-Geo] Advice on clustering LiDAR point clouds

kind regards, Wesley Wesley Roberts MSc. Researcher: Earth Observation Natural Resources & the Environment (NRE) CSIR Tel: +27 (0)21 888-2490 Fax: +27 (0)21 888-2693 "To know the road ahead, ask those coming back." - Chinese proverb -- This message is subject to the CSIR&#x

[R-sig-Geo] writeOGR unknown data type

Both classes seem to be similar and writeOGR should be able to export both data sets. > class(poly) [1] "SpatialPointsDataFrame" attr(,"package") [1] "sp" > class(grd) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" Any he

Re: [R-sig-Geo] how to show polygons from SHP file in Google map?

ve the required packages loaded. Hope this helps Wesley Wesley Roberts MSc. Researcher: Earth Observation Natural Resources & the Environment (NRE) CSIR Tel: +27 (0)21 888-2490 Fax: +27 (0)21 888-2693 "To know the road ahead, ask those coming back." - Chinese proverb >>&g

[R-sig-Geo] calculate minimum radial separation

patial R packages to calculate the minimum maximum distance from a centroid to its polygon boundary? Many thanks for your help, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2693 "

Re: [R-sig-Geo] Calculating area of polygons

Thanks Roger, Works a treat. Wesley >>> Roger Bivand 05/14/09 9:51 AM >>> On Thu, 14 May 2009, Wesley Roberts wrote: > Dear list, > > I would like to calculate the area of a polygon. I have found areapl() in > splancs but am having trouble using it. &g

[R-sig-Geo] Calculating area of polygons

4.760465 square meters (area slot) but cant seem to replicate this using areapl(). Is it possible to access the area slot directly instead of computing the area with areapl()? Many thanks, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environme

Re: [R-sig-Geo] Classification of attribute table

ter of my PhD and would like to get it published. Many thanks to the list for all your assistance. Kind regards, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2693 "To know the road ah

Re: [R-sig-Geo] Classification of attribute table

, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2693 "To know the road ahead, ask those coming back." - Chinese proverb >>> Dan Putler 05/07/09 6:13 PM >&

[R-sig-Geo] Classification of attribute table

kes a really long time. In essence all I want to do is classify the 6 variable data set into 3 classes with the class for each case recorded in the attribute table. Any advice or suggestions would be greatly appreciated. Many thanks and kind regards, Wesley Wesley Roberts MSc. Researcher: Ear

Re: [R-sig-Geo] EPSG codes for LO31 WGS84 projection data

ailing list: http://lists.maptools.org/mailman/listinfo/proj cheers, Paul Wesley Roberts wrote: > Dear Colleagues, > > I have a data set which is projected using the following projection, > Hartebeeshoek94 Projection Gauss Conform LO31. This is a local projection > used in South Africa. So f

[R-sig-Geo] EPSG codes for LO31 WGS84 projection data

;cat", > proj4string=CRS("+init=epsg:2054")) Error in CRS("+init=epsg:2054") : no options found in 'init' file Could someone please provide me with the full proj4 definition as the spatial reference site does not seem to have it. Many thanks and kin

[R-sig-Geo] Time series analysis with irregular time-series

ience with time series analysis and has any suggestions regarding the questions posted, I would really appreciate some help. Many thanks and kind regards, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax

Re: [R-sig-Geo] Search engine for R-sig-Geo archives??

Hi Tyler, I use the following (copy and paste into google search bar) site:https://stat.ethz.ch/pipermail/r-sig-geo/ This way google is restricted to that site only. You can actually use site: for any website. HTH Wesley >>> Paul Hiemstra 03/27/09 12:29 PM >>> Tyler Dean Rudolph schreef:

Re: [R-sig-Geo] Error in if (dxsd > tolerance) { : missing valuewhere TRUE/FALSE needed

the failing function by using traceback() afterwards - here it > is points2grid() that has failed. > > Roger > >> >> -David Wahlund >> >> On Wed, Mar 18, 2009 at 16:01, Wesley Roberts >> wrote: >>> Dear R-sig-Geo'ers >>> >>&g

[R-sig-Geo] Error in if (dxsd > tolerance) { : missing value where TRUE/FALSE needed

-HERE IS THE PROBLEM- a <- mean(as.data.frame(x.clip$band1)) Lst[n] <- list(as.numeric(a)) print(n) x.clip<-0 }

Re: [R-sig-Geo] Masking interpolations

m seems to be an internal coding, not sure if you have seen this before? Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2693 "To know the road ahead, ask those coming back." -

Re: [R-sig-Geo] Masking interpolations

ge <- . and y.range <- . statements but the entire grid (square) is not defined and thus interpolation only occurs where data exists? Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-

Re: [R-sig-Geo] Masking interpolations

e then limiting the interpolation based on kriging variance seems like a good idea? What do you think? Link 3 http://www.flickr.com/photos/35273...@n07/3284446466/ Many thanks, Wesley Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +2

[R-sig-Geo] Masking interpolations

ht = autoKrige(H~1, a, grd, nmax=100) writeGDAL(height$krige_output, fname="test.tiff", drivername ="GTiff", type = "Float32") intensity = autoKrige(I~1, a, grd, nmax=100) writeGDAL(intensity$krige_output, fname="test.tiff", drivername ="GTiff", type

Re: [R-sig-Geo] variogram question

p, by Paul Hiemstra, which > tries to do some of these heuristics--good or bad, who will tell. > -- > Edzer > > Wesley Roberts wrote: >> Dear r-sig-geo users, >> >> I am currently analyzing some Lidar data we have collected over our >> study area. I am interested in

[R-sig-Geo] variogram question

? Many thanks for all your help and suggestions Wesley Wesley Roberts MSc. Researcher Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2490 "To know the road ahead, ask those coming back." - Chinese proverb -- This messag