Re: AI-GEOSTATS: moving averages and trend
-1098 (ext. 202) Fax: (734) 913-2201 Courtesy Associate Professor, University of Florida Associate Editor, Mathematical Geosciences Geostatistician, Computer Sciences Corporation President, PGeostat LLC 710 Ridgemont Lane Ann Arbor, MI 48103 Voice: (734) 668-9900 Fax: (734) 668-7788 http://goovaerts.pierre.googlepages.com/ -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: (2) Geostatistics in pain
to a point if a large study area is used. I've read other places that among the commercial stat packages, SAS is best at handling large data sets. Is this true? Also, I've produced my own little routine in IDRISI that can create 'random' samples that are clustered by inverse distance, so that short lags are preferred. Are there any software packages that can create a random sample of points that show a pre-specified clustering pattern in space? JanWMerks wrote: I read with a great deal of interest your emessage about Geostatistics in pain. Read what I have found out. Real statistics turned into surreal geostatistics under the guidance of Professor Dr Georges Matheron... Sebastiano: In general I agree with the comments reported in the preceding replies. Then I would add that the problem, if any one exists, doesn't relies on the lack of a good gui but maybe on the lack of a kind of standard and internationally accepted set of programming routines directed to geostatistical analysis. As a final consideration I think that the world of spatial analysis is more complex than in the past and a kind of holistic view and is needed. For example I'm thinking to other techniques based on statistical learning theory,data mining, etc... Paul wrote: Being an R user myself (gstat, automap) I would like to comment a little on your problem with command line tools. I think what is most important is that each application has its own best tool to use. When a novice user wants to quickly make some maps, a GUI would probably be the preferred tool. But if you, as in my case, want to interpolate thousands of maps, put them on a webservice and allow user to get those maps from the web, a GUI tool such as ArcGIS is probably not the best option. R is great for these kinds of large analysis. In addition, I'm a Linux user and would not trade my command line for any GUI :). Therefore I believe a GUI is not per definition better than command line. It's just that people are used to GUI nowadays, making command line seem old. I agree that it takes quite some time to learn R and that it is not a tool suitable for the casual user. A great combo would be a tool that has the flexibility and power of R and the ease of use of e.g. ArcGIS.Hope you find a solution that suits your particular needs. End of Log == King Regards, Younes __ See what's on at the movies in your area. Find out now: http://au.movies.yahoo.com/session-times/ + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/ -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: (1) Geostatistics in pain
Younes Fadakar wrote: Hi there, This is my first message as a test message checking the usage of the service, working with ai-geostats mailing list. I have many questions to ask you too. To start: The world of Geostatistics seriously needs a tool to present well to novices and professionals. Current availabilities have many of disadvantages. Some are too old, others not user-friendly and the rest more expensive. 1- GsLib seems powerful but too old (DOS-command line in 2010!) 2- WinGsLib is completely confusing despite of logging and automating! no direct input and output!! 3- Variowin is too weak in terms of GUI! 4- mGstat as a Matlab toolbox written too complex not handy program! 5- GeoR as an extention for R makes you to work with R such a command-line environment! what a development rather than GsLib!! Hi Younes, Being an R user myself (gstat, automap) I would like to comment a little on your problem with command line tools. I think what is most important is that each application has its own best tool to use. When a novice user wants to quickly make some maps, a GUI would probably be the preferred tool. But if you, as in my case, want to interpolate thousands of maps, put them on a webservice and allow user to get those maps from the web, a GUI tool such as ArcGIS is probably not the best option. R is great for these kinds of large analysis. In addition, I'm a Linux user and would not trade my command line for any GUI :). Therefore I believe a GUI is not per definition better than command line. It's just that people are used to GUI nowadays, making command line seem old. I agree that it takes quite some time to learn R and that it is not a tool suitable for the casual user. A great combo would be a tool that has the flexibility and power of R and the ease of use of e.g. ArcGIS. Hope you find a solution that suits your particular needs. cheers, Paul 6- Isatis is more expensive; for what?! 7- Gs+ is in pain with weak performance of GUI! 8- Geoeas is something funny just to remember DOS graphics! 9- ... So obviously a serious request remained for more than 20 years without suitable answer! Why? Younes __ See what's on at the movies in your area. Find out now: http://au.movies.yahoo.com/session-times/ + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/ -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: Bringing in ascii DTM file from ArcGIS with -9999 NAN values for KED
Hi Tobin, Tobin Cara wrote: Hello, I am having problems with kriging an ascii file DTM from ArcGIS which has many - (no data) values. Do these normally have to be removed first? I want to make the elevation DTM the predictor base for the kriging with external drift. I would remove them, in this case no prediction would be made at these locations. My code is the following: elev - read.asciigrid(elev.asc, as.image=FALSE, plot.image=TRUE) str(elev) elev_ked - krige(meanRain~Zloc, locations=meanRain, newdata=elev, model=fitted_vario) Error in model.frame.default(terms(formula), as(data, data.frame)) : object is not a matrix Your object is not a matrix. It is hard to give more detail without a piece of R code that reproduces your problem. This particular error message doesn't ring a bell with me. As other people already suggested, use the r-sig-geo mailing list for these kinds of questions. This mailing list is for general discussion, on r-sig-geo you'll get far more response. My data file meanRain has X Y Z coordinates and Precipitation values (4 columns). Thank you so much for your time and help! Cara -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
AI-GEOSTATS: Re: [R-sig-Geo] [R-sig-geo] Interpolation of measures with measurement errors
Hi, Have a look at this: http://books.google.com/books?id=5n_XuL2Wx1ECpg=PA94lpg=PA94dq=interpolating+with+measurement+errorsource=blots=837gfXTqC4sig=qzOg_9cuDDa_Nk5j71b_7P-1iqohl=enei=bcDASriyJ5De-QbJjeHKAQsa=Xoi=book_resultct=resultresnum=3#v=onepageq=interpolating%20with%20measurement%20errorf=false third hit when googling for interpolating with measurement error cheers, Paul note replying from the r-sig-geo list, also cc'ing the ai-geostats list Enrico Guastaldi wrote: Dear list members, I'm looking for some kind of interpolation for values of an environmental variable which has been measured together the measurement errors. I could use some kind of kriging, however I exactly know the magnitude of each error at every sampled location, i.e. the value plus or minus the error gave me by the laboratory. Could anyone tell me what kind of function should I use for handling this problem? It should be nice some R package, of course, but I need to understand the background theory. Thanks in advance, Regards, Enrico Guastaldi [[alternative HTML version deleted]] ___ R-sig-Geo mailing list r-sig-...@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 274 3113 Mon-Tue Phone: +3130 253 5773 Wed-Fri http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: please help me
Hi Sadeghian, You could consider posting this question on the r-help mailing list. (https://stat.ethz.ch/mailman/listinfo/r-help) I think that would give you an answer much sooner. cheers, Paul azadeh sadeghian schreef: dear list I am student M.S. statistics in department statistics . I am working in the function nls in the [R 2.3.1] with 246 data and want to fit the exp model to vectors( v and u ) but I have a problem to use it u 5.00e-13 2.179057e+03 6.537171e+03 1.089529e+04 1.525340e+04 1.961151e+04 2.396963e+04 2.832774e+04 3.268586e+04 3.704397e+04 4.140209e+04 4.576020e+04 5.011831e+04 5.447643e+04 v 8.382562e-01 4.090868e+02 1.311053e+03 2.124143e+03 3.365494e+03 2.138903e+03 7.687774e+03 1.028396e+04 1.004186e+04 2.059798e+04 1.438464e+04 2.861373e+04 2.294919e+04 2.807701e+04 data1-data.frame(u=u ,v=v) nls(v~c0+(ce*(1-exp((-u)/ae))),data=data1,start=list(c0=0,ce=1000,ae=3000)) Error in nls(v ~ c0 + (ce * (1 - exp((-u)/ae))), data = data1, start = list(c0 = 0, : step factor 0.000488281 reduced below 'minFactor' of 0.000976563 i dont know how to solve it, please help me . best regards Sadeghian Building a website is a piece of cake. Yahoo! Small Business gives you all the tools to get online. http://us.rd.yahoo.com/evt=48251/*http://smallbusiness.yahoo.com/webhosting/?p=PASSPORTPLUS -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +31302535773 Fax:+31302531145 http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: difference between correlogram and variogram models
Dear Abani, The following is a quote from a geostatistics book by Chiles and Delfiner (p33): ..gives us two reasons to favor the variogram over the covariance. The first is theoretical: since the class of Intrinsic Random Functions includes the Stationary Random Functions, the variogram is a more general tool than the covariance. ... The second reason is practical: the variogram, unlike the covariance, does not require the knowledge of the mean. Chiles, J., Delfiner, P., 1999. Geostatistics: Modeling Spatial Uncertainty. John Wiley Sons, New York. hth, Paul Abani R Samal schreef: Can any one in the list help me explain the drawbacks of using correlogram for kriging? In other words, what are the advantages of using variogram models instead of correlogram models? why people don't use correlogram models more (as compared to variogram models)? Regards, Abani R Samal *ABANI RANJAN SAMAL **11183 West 17th Avenue, APt 201* *Lakewood, CO 80215* You snooze, you lose. Get messages ASAP with AutoCheck http://us.rd.yahoo.com/evt=47959/*http://advision.webevents.yahoo.com/mailbeta/newmail_html.html in the all-new Yahoo! Mail Beta. -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +31302535773 Fax:+31302531145 http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/
Re: AI-GEOSTATS: multi-point geostatistics
Peter Bossew schreef: Dear list, can somebody recommend an introduction into multiple-point geostatistics ? or a literature review ? Thanks for hints, Peter - Peter Bossew European Commission (EC) Joint Research Centre (JRC) Institute for Environment and Sustainability (IES) TP 441, Via Fermi 1 21020 Ispra (VA) ITALY Tel. +39 0332 78 9109 Fax. +39 0332 78 5466 Email: [EMAIL PROTECTED] WWW: http://rem.jrc.cec.eu.int The views expressed are purely those of the writer and may not in any circumstances be regarded as stating an official position of the European Commission. + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/ Dear Peter, I have some work on multiple point methods for simulating complex geological structures, such as channel belts. These are a few papers I used: Arpat, B.G. (2005), Sequential simulation with patterns, PhD Thesis, Stanford University, Stanford. http://pangea.stanford.edu/~jcaers/theses/thesisBurcArpat.pdf Strebelle, S., (2002), Conditional simulation of complex geological structures using multiple-point statistics. Mathematical Geology, Jan. 2002. Zhang T., Switzer, P., and Journel, A., (2006), Filter-based classification of training image patters for spatial simulation. Mathematical Geology, 38:1, pp. 63-80 Kind regards, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +31302535773 Fax:+31302531145 http://intamap.geo.uu.nl/~paul + + To post a message to the list, send it to ai-geostats@jrc.it + To unsubscribe, send email to majordomo@ jrc.it with no subject and unsubscribe ai-geostats in the message body. DO NOT SEND Subscribe/Unsubscribe requests to the list + As a general service to list users, please remember to post a summary of any useful responses to your questions. + Support to the forum can be found at http://www.ai-geostats.org/