nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

automate kriging of multiple (120) layers

LaurenSahagun

Hello,

I am trying to find geographic patterns within pavement distresses at 70+ airfields across the continental United States. I have many layers that I have already completed spatial analysis on (manually). I had to change the values used to Krig from and would like to find a faster way to perform the Krig analysis using Python.

Bottom line:

I need to Krig 120 layers

I would like to use the first layer as a template for the remaining 119 layers because I need to change the parameters of lag size, # of classes, break values, and have it use the maximum of any coincidental points

Once created I would like to rename the Krig layer the same of the inputLayer + "Krig"

Any help would be great!! Thank you very much.

Find more posts tagged with

Accepted answers

DanPatterson_Retired

2017-01-13 Updated script formatting

Alright...an example...just read it carefully

<SPAN class="string token">""</SPAN>"
krigingTest<SPAN class="punctuation token">.</SPAN>py
Dan_Patterson<SPAN class="punctuation token">.</SPAN>carleton<SPAN class="punctuation token">.</SPAN>ca
requires<SPAN class="punctuation token">:</SPAN>&nbsp; read carefully
see http<SPAN class="punctuation token">:</SPAN><SPAN class="operator token">//</SPAN>resources<SPAN class="punctuation token">.</SPAN>arcgis<SPAN class="punctuation token">.</SPAN>com<SPAN class="operator token">/</SPAN>en<SPAN class="operator token">/</SPAN>help<SPAN class="operator token">/</SPAN>main<SPAN class="operator token">/</SPAN><SPAN class="number token">10.2</SPAN><SPAN class="operator token">/</SPAN>index<SPAN class="punctuation token">.</SPAN>html<SPAN class="comment token">#//009z0000006n000000</SPAN>
 <SPAN class="keyword token">as</SPAN> well
<SPAN class="string token">""</SPAN>"

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#get the extension</SPAN>
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#import all the frigging stuff</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#make sure that you overwrite your feable attempts</SPAN>
input_shp <SPAN class="operator token">=</SPAN> <SPAN class="string token">"C:/Temp/krige/pnts.shp"</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#input the shapefile name</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN>input_shp<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"temp"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token">#create a temporary layer</SPAN>
<SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">6</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#run it in a range, make sure it isn't 0 for the start</SPAN>
&nbsp;&nbsp;&nbsp; cellsize <SPAN class="operator token">=</SPAN> i&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#vary the cell size is being modelled *********</SPAN>
&nbsp;&nbsp;&nbsp; outname <SPAN class="operator token">=</SPAN><SPAN class="string token">"C:/Temp/Krige/test"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>cellsize<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN> <SPAN class="string token">"output grid: {0} cellsize {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">,</SPAN>cellsize<SPAN class="punctuation token">)</SPAN> <SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#arcpy.gp.Kriging_sa("pnts","Z","C:/Temp/Krige/output","Spherical 4.110560","2","VARIABLE 12","#") #another example</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#arcpy.gp.Idw_sa("pnts","Z","C:/Temp/Krige/IDW","2","2","VARIABLE 12","#")</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Do the work *****</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>Kriging_sa<SPAN class="punctuation token">(</SPAN>input_shp<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"Z"</SPAN><SPAN class="punctuation token">,</SPAN>outname<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"Spherical 4.110560"</SPAN><SPAN class="punctuation token">,</SPAN>cellsize<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"VARIABLE 12"</SPAN><SPAN class="punctuation token">,</SPAN>"<SPAN class="comment token">#") </SPAN>
&nbsp;&nbsp;&nbsp; outRaster <SPAN class="operator token">=</SPAN> Raster<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; outRaster<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

DanPatterson_Retired

No problem, good luck with your example

LaurenSahagun

Dan- Thank you very much for your help. I really appreciate you taking the time to consider my question.

DanPatterson_Retired

2017-01-13 Updated script formatting

Alright...an example...just read it carefully

<SPAN class="string token">""</SPAN>"
krigingTest<SPAN class="punctuation token">.</SPAN>py
Dan_Patterson<SPAN class="punctuation token">.</SPAN>carleton<SPAN class="punctuation token">.</SPAN>ca
requires<SPAN class="punctuation token">:</SPAN>&nbsp; read carefully
see http<SPAN class="punctuation token">:</SPAN><SPAN class="operator token">//</SPAN>resources<SPAN class="punctuation token">.</SPAN>arcgis<SPAN class="punctuation token">.</SPAN>com<SPAN class="operator token">/</SPAN>en<SPAN class="operator token">/</SPAN>help<SPAN class="operator token">/</SPAN>main<SPAN class="operator token">/</SPAN><SPAN class="number token">10.2</SPAN><SPAN class="operator token">/</SPAN>index<SPAN class="punctuation token">.</SPAN>html<SPAN class="comment token">#//009z0000006n000000</SPAN>
 <SPAN class="keyword token">as</SPAN> well
<SPAN class="string token">""</SPAN>"

<SPAN class="keyword token">import</SPAN> arcpy
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#get the extension</SPAN>
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#import all the frigging stuff</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#make sure that you overwrite your feable attempts</SPAN>
input_shp <SPAN class="operator token">=</SPAN> <SPAN class="string token">"C:/Temp/krige/pnts.shp"</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#input the shapefile name</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN>input_shp<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"temp"</SPAN><SPAN class="punctuation token">)</SPAN>&nbsp; <SPAN class="comment token">#create a temporary layer</SPAN>
<SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">6</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#run it in a range, make sure it isn't 0 for the start</SPAN>
&nbsp;&nbsp;&nbsp; cellsize <SPAN class="operator token">=</SPAN> i&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#vary the cell size is being modelled *********</SPAN>
&nbsp;&nbsp;&nbsp; outname <SPAN class="operator token">=</SPAN><SPAN class="string token">"C:/Temp/Krige/test"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>cellsize<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN> <SPAN class="string token">"output grid: {0} cellsize {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">,</SPAN>cellsize<SPAN class="punctuation token">)</SPAN> <SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#arcpy.gp.Kriging_sa("pnts","Z","C:/Temp/Krige/output","Spherical 4.110560","2","VARIABLE 12","#") #another example</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token">#arcpy.gp.Idw_sa("pnts","Z","C:/Temp/Krige/IDW","2","2","VARIABLE 12","#")</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Do the work *****</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>Kriging_sa<SPAN class="punctuation token">(</SPAN>input_shp<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"Z"</SPAN><SPAN class="punctuation token">,</SPAN>outname<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"Spherical 4.110560"</SPAN><SPAN class="punctuation token">,</SPAN>cellsize<SPAN class="punctuation token">,</SPAN><SPAN class="string token">"VARIABLE 12"</SPAN><SPAN class="punctuation token">,</SPAN>"<SPAN class="comment token">#") </SPAN>
&nbsp;&nbsp;&nbsp; outRaster <SPAN class="operator token">=</SPAN> Raster<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; outRaster<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

DanPatterson_Retired

To get you started...do one manually until you get everything right, then go to the Geoprocessing Results window, right-click on the last process you did and copy it as a python snippet as shown in the figure. From there it is a matter of putting that within a loop fixing your output names as you go along.

# Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script

# The following inputs are layers or table views: "pnts"

arcpy.gp.Kriging_sa("pnts","Z","C:/Temp/Krige/output","Spherical 4.110560","2","VARIABLE 12","#")

is the output from a simple default run with nothing else selected.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.