Don & Facundo, This is something that I'm interested in as well for coastal marine applications.
The only solution I've found so far is a cumbersome "barrier" function available in the kriging operation of ArcGIS. I would be interested in seeing an (better) implementation of the idea in R. I was able to find a few references on the subject that I pass along to you and anyone else that might find it useful: Krioruchko, K. and A. Gribov. 2002. Geostatistical interpolation and simulation in the presence of barriers. X. Sanchez-Vila et al. (eds.), *geoENV IV -- Geostatistics for Environmental Applications,* pp. 331-342. Kluwer Academic Publishers: the Netherlands. Curriero, F.C. 1997. The use of non-Euclidean distances in geostatistics, PhD thesis. Little, LS, D Edwards and DE Porter. 1997. Kriging in estuaries: as the crow flies or as the fish swims? Journal of experimental marine biology and ecology. 213:1-11. Rathbun, S. 1998. Spatial modeling in irregularly shaped regions: kriging estuaries. Envirometrics. 9:109-129. Jensen, OP, MC Christman and TJ Mller. 2006. Landscape-based geostatistics: a case study of the distribution of blue crab in Chesapeake Bay. Envirometrics. 17:605-621. Regards, Alisa Coffin On Thu, Mar 20, 2008 at 2:42 PM, Don MacQueen <[EMAIL PROTECTED]> wrote: > Facundo, > > I am interested in essentially the same thing, > though my application is inside buildings. I am > wondering what technique(s), software, etc. you > would be using to calculate your distances. > > I asked about this about a month ago, here on > r-sig-geo, and received some helpful replies. The > subject line at that time was > [R-sig-Geo] Shortest path around obstacles (OT) > > -Don > > At 7:05 PM +0100 3/20/08, Facundo Muñoz wrote: > >Hello All, > > > >I'm performing regression kriging (using both gstat and geoR) in an > >urban environment. Thus, i have discontinuities (buildings) in the > >prediction space which make distances to be non-euclidean. > > > >I need to estimate variogram and make predictions using distances > >calculated by myself and passed to the algorithms (say, in a distance > >matrix), instead of letting them calculate distances based on > coordinates. > > > >¿Has anyone faced this issue before? > > > >Thanks in advance, > > Facundo Muñoz > > > >-- > >Departament d'Estadística i Investigació Operativa Universitat de > >València (Estudi General) > >Facultat de Matemàtiques, Dr. Moliner 50, 46100 Burjassot, Spain. > >(+34) 96 354 3987, fax: (+34) 96 354 3238 > >e-mail: [EMAIL PROTECTED], web: > >http://www.uv.es/~famarmu/<http://www.uv.es/%7Efamarmu/> > > > >_______________________________________________ > >R-sig-Geo mailing list > >[email protected] > >https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > > -- > -------------------------------------- > Don MacQueen > Environmental Protection Department > Lawrence Livermore National Laboratory > Livermore, CA, USA > 925-423-1062 > > _______________________________________________ > R-sig-Geo mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo > > -- Alisa Coffin, PhD. Candidate Department of Geography University of Florida Gainesville, FL 32611 [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
