Is there a way to krige with barriers in ArcGIS 10?

Question

I am trying to do interpolation in a coastal area with data points in several rivers emptying into a bay and on opposite sides of peninsulas. I have created a coastline layer to use it as a barrier but cannot figure out how to do this in ArcGIS 10. I do not see a barriers object anywhere in any of the kriging tools. I tried just doing the kriging then clipping it to the land, but I get unrealistic results where data points close together but separated by land create problems. I see online documentation for doing this in geoprocessing scripts in 9.3 but nothing similar for 10. Thanks so much!

Answer 1

Do the rivers run generally in a specific direction? If so, you could use directional binning along the centerlines of the rivers to constraint the empirical semivariograms only to point pairs running along the length of the rivers. (See here for more details: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Exploring_directional_autocorrelation/0031000000m2000000/)

Another option would by to partition the point samples into separate datasets and run your kriging on each set separately. This might be more appropriate anyway, since each river might have different statistical models to them; but would not work if your sample size is small for each river.

Do you need to use a geostatistical method, or would a deterministic method work? I think you might be able to achieve what you want using Diffusion Interpolation with Barriers. You could use your land as an extremely high cost barrier. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00310000002s000000.htm

Also, here is Diffusion Interpolation with absolute barriers as a geoprocessing script tool. Your coastline would make a perfect absolute barrier. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Diffusion_Interpolation_With_Barriers/003000000005000000/

Answer 2

Have you tried the Topo to Raster tool? You can specify the type of the input features (Point Elevation, Stream, sink, contour, boundary, lake.) You could try setting your coast line layer features to "TopoLake"

Also, the result of the Topo to Raster tools is a "hydologically correct" raster

Quoting from Help: TopoLake A polygon feature class that specifies the location of lakes. All output raster cells within a lake will be assigned to the minimum elevation value of all cells along the shoreline.

Answer 3

I'm a late adopter so I can only tell you what works in 9.3 with the Spatial Analyst extension -- posting just in case the method still works in 10. First, create a closed linestring representing your area of interest (sounds like the inverse of your coastline here). Specify your "container" in the Spatial Analyst options dialog on the extent tab. Perform your analysis by selecting Kriging from the Interpolate to Raster in the Spatial Analyst menu.

Good luck!

Answer 4

i think that you can use topo to raster tool and then use output extent for interpolating or use kriging then you use extract by mask tool for barrier . you should create a barrier layer (polygon layer) and then use extract by mask tool

sincerely