Latest Contributions by AnushuyaRamakrishnan

Hi Steve: I have not contacted the ESRI support yet. I do not know how it works. ... View more

Hi Steve: I have not solved the problem that I reported earlier yet. I confirm that there is problem with the python code. When I try to run kriging for certain datasets that already have kriged predictions generated, script is throwing me an error. The error is as follows: Traceback (most recent call last):   File "C:\Users\aramakrishnan\Documents\ArcGIS\scratch\Ordinary Kriging_July2014.py", line 41, in <module>     arcpy.gp.ExtractValuesToPoints_sa(Harris_Lat_Long_csv_Events_lyr, Raster_Path + "OZ8_" + fieldname, Ext_Table_Path +  "OZ8_" + fieldname, "INTERPOLATE", "VALUE_ONLY")   File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\geoprocessing\_base.py", line 498, in <lambda>     return lambda *args: val(*gp_fixargs(args, True)) ExecuteError: ERROR 999999: Error executing function. Failed to execute (ExtractValuesToPoints). Failed to execute (ScriptJul21). Failed at Wed Aug 27 19:35:22 2014 (Elapsed Time: 1 minutes 55 seconds) I tried to run ordinary kriging by manual methods by inputing the lag, sil and nugget values. However, I am not able to even generate the raster to feed it into the script. I attach my complete python code and xml file. Will wait to hear from you Thank you Anushuya ... View more

Hi,       I need some help with the error I get while using spatial analyst tool in a python script (attached).       Due to hard disk space issues, I deleted most of the old data files in the default.gdb folder and tried to run the script, which was working fine before the data clean-up.      Below is the error message I get.      Please help me identify what is missing and how to fix the issue. Thanks, Anushuya Executing: Script Start Time: Sun Aug 10 15:47:46 2014 Running script Script... Traceback (most recent call last):   File "C:\Users\aramakrishnan\Documents\ArcGIS\scratch\Ordinary Kriging_July2014.py", line 41, in <module>     arcpy.gp.ExtractValuesToPoints_sa(Harris_Lat_Long_csv_Events_lyr, Raster_Path + "OZ4" + fieldname, Ext_Table_Path +  "OZ4" + fieldname, "INTERPOLATE", "VALUE_ONLY")   File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\geoprocessing\_base.py", line 498, in <lambda>     return lambda *args: val(*gp_fixargs(args, True)) ExecuteError: ERROR 999999: Error executing function. Failed to execute (ExtractValuesToPoints). Failed to execute (Script). Failed at Sun Aug 10 15:49:38 2014 (Elapsed Time: 1 minutes 51 seconds) ... View more

Hi Eric, I followed your suggestions and ran the loop for 7 days and checked few days manually. I got exactly matching results. Thank you for your help and suggestions. I appreciate it Thanks Anushuya ... View more

Thank you for your suggestion Eric I will try the options that you suggest and let you know if it works. Thanks Anushuya ... View more

Here is the XML input that I used. Thanks Anushuya ... View more

Hi Eric, I generated my ordinary kriging output successfully using the python script for generating geostatistical layer and using it to generate raster and convert it to points and generate the table. However, I see that the output generated by scripting does not exactly match with the output that I generate manually by running ordinary kriging for the same day using the same values of range, sill and nugget. In the script I only set auto=true for lagsize, number of lags, range, sill and nugget values. I attach the python script for your review. Please suggest. Thank you Anushuya ... View more

Hi, Do you input a raster or a layer or a data file like CSV. Please go to add tools and add your specific input. Then go to add xy tool and project it in arcmap. Then perform specific geostatistical operation using the events layer and get the output layer Next right click on layer and convert it to raster. Use this raster to generate your output table. In any module of arcmap the screen should not go blank. Try these steps and check if it works. Also try and reinstall arcmap if the problem persists. Thanks Anushuya ... View more

Thank you Eric... I will do that Anushuya ... View more

Hi Eric, I could succesfully generated new geostatistical layer from an existing xml input. Can you guide me to the next step of kriging. For kriging using SA I use the following script: # Process: Kriging arcpy.gp.Kriging_sa(merge3_csv_Events_lyr, "longitude", fieldname, "Exponential 0.003512 0.255986 1.092512 0.001093", "3.5120478474024E-03", "VARIABLE 12", Output_variance_of_prediction_raster) 1. I need the script for kriging using ga 2. What is the input for this script? Is it the same as SA. Pls. let me know. Thank you Anushuya ... View more

Hi Eric, I tried to run this in python: # Import arcpy module import arcpy # Check out the ArcGIS Geostatistical Analyst extension license arcpy.CheckOutExtension("GeoStats") merge1.csv.events = arcpy.GeostatisticalDatasets("C:\\temp\\KrigingordinaryJan152.xml") Harris_Lat_Long.csv.events = "C:\\temp\\PM04Jan16.shp" merge1.csv,events.Field = "Y04PMJan16" arcpy.GACreateGeostatisticalLayer_ga("C:\\temp\\KrigingordinaryJan16.xml", newDataset, "PM04Jan16") It is throwing me an error: Traceback (most recent call last):   File "C:\Users\aramakrishnan\Documents\ArcGIS\scratch\Ordinary Kriging_09252013_01.py", line 19, in <module>     merge1.csv.events = arcpy.GeostatisticalDatasets("C:\\temp\\KrigingordinaryJan152.xml")   File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\_ga.py", line 300, in GeostatisticalDatasets     return GADatasets(gapy.GALayerInfo(ga_model_source)) TypeError: bad argument type for built-in operation Am I missing something? >>> ... View more

Thank you Eric!!! I will try this and get back to you Anushuya ... View more

I am using arcmap 10.2 Thanks Anushuya ... View more

HI Eric, I am still not clear about the flow of python script for ordinary kriging. # Process: Kriging arcpy.gp.Kriging_sa(merge1_csv_Events_lyr, "Y04Mar02", Kriging_PM04Mar2, "Spherical 0.003512 0.339420 4.530740 0.517570", "3.5120478474024E-03", "VARIABLE 12",                     Output_variance_of_prediction_raster) (merge_1_csv_Events_lyr contains PM data for all the days of the study; Y04Mar02 is the Z field for one day) I attach the XML file that I generated for one day of observation as ZIP file. Do we built a model in model builder and give that model name as input for xmlpath. What is happening in these steps: xmlPath = "/model[@name = 'Kriging']/model[@name = 'Variogram']/value[@name = 'Range']" newValue = 0 outModel = "C:/gapyexamples/output/outModel.xml" model source="C:\\Users\\aramakrishnan\\Desktop\\Grid combo1\\merge1.csv Events" xmlPath= "/model[@name = 'Kriging']/model[@name = 'Variogram']/model[@name = 'VariogramModel']/value[@name = 'Sill']" newValue = 1 outModel = "C:/gapyexamples/output/outModel.xml" model source="C:\\Users\\aramakrishnan\\Desktop\\Grid combo1\\merge1.csv Events" xmlPath= "/model[@name = 'Kriging']/model[@name = 'Variogram']/value[@name = 'Nugget']/@auto;\ newValue = 2 outModel = "C:/gapyexamples/output/outModel.xml" For this step what is the inlayer? Do we carry out this step after calling in range,sill and nugget values. Execute CreateGeostatisticalLayer arcpy.GACreateGeostatisticalLayer_ga(inLayer, inData, outLayer) Please help me to understand this better. Thank you Anushuya ... View more

Hi Eric, Thank you for the suggestion.  I will certainly try that. EBK is relatively new and has not been reported in scientific literature for prediction of pollutant concentrations in adverse birth outcomes studies that we are involved in. When we send our manuscript for review, there will be questions on the validity of the model and the applicability of it in related studies. As ordinary kriging as an established model reported in a lot of manuscripts, use of it in the study will eliminate question on the validity of it. Thanks Anushuya ... View more