'Grids selected data columns in selected data file using all 12
'gridding methods and produces contour maps of
'the resulting grids to compare results.

'1. Go to Script | Run
'2. Select data file and click Open.
'3. Enter the column number for the X data and click OK.
'4. Enter the column number for the Y data and click OK.
'5. Enter the column number for the Z data and click OK.

'******************************************

Sub Main
	
	Dim SurferApp As Object
	Set SurferApp = CreateObject("Surfer.Application")
	SurferApp.Visible = True

'Specify the height of the map
	ylength = 2

'Specifies the data file to grid
	Data$ = GetFilePath( , "dat;csv;xls;xlsx",path1, "Data file to grid")
	If Data$ ="" Then End

	Dim Doc As Object
	Set Doc = SurferApp.Documents.Add

'Set data columns
	xcol= Val(InputBox$("Enter X column number. For example, Column A = 1", "X Column", "1"))
	ycol= Val(InputBox$("Enter Y column number. For example, Column B = 2", "Y Column", "2"))
	zcol= Val(InputBox$("Enter Z column number. For example, Column C = 3", "Z Column", "3"))

'1. Inverse Distance
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=1, OutGrid:=Data$+"_InverseDistance.grd")

	Dim Map As Object
	Set Map = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_InverseDistance.grd")

	Map.yLength = ylength
	Map.xMapPerPU = Map.yMapPerPU
	If Map.Width<3.5 Then Map.xMapPerPU=Map.yMapPerPU Else Map.xLength=3.5

	Map.Axes("Bottom Axis").Title = "Inverse Distance"
	Map.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map.Left = 0.25
	Map.Top = 11.0

'2. Kriging
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=2, OutGrid:=Data$+"_Kriging.grd")

	Dim Map2 As Object
	Set Map2 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_Kriging.grd")

	Map2.yLength = ylength
	Map2.xMapPerPU = Map.yMapPerPU
	If Map2.Width<3.5 Then Map2.xMapPerPU=Map2.yMapPerPU Else Map2.xLength=3.5

	Map2.Axes("Bottom Axis").Title = "Kriging"
	Map2.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map2.Left = 2.99
	Map2.Top = 11.0

'3. Minimum Curvature
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=3, OutGrid:=Data$+"_MinimumCurvature.grd")

	Dim Map3 As Object
	Set Map3 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_MinimumCurvature.grd")

	Map3.yLength = ylength
	Map3.xMapPerPU = Map.yMapPerPU
	If Map3.Width<3.5 Then Map3.xMapPerPU=Map3.yMapPerPU Else Map3.xLength=3.5

	Map3.Axes("Bottom Axis").Title = "Minimum Curvature"
	Map3.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map3.Left = 5.74
	Map3.Top = 11.0

'4. Modified Shepard's Method
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=4, OutGrid:=Data$+"_ModifiedShepard.grd")

	Dim Map4 As Object
	Set Map4 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_ModifiedShepard.grd")

	Map4.yLength = ylength
	Map4.xMapPerPU = Map.yMapPerPU
	If Map4.Width<3.5 Then Map4.xMapPerPU=Map4.yMapPerPU Else Map4.xLength=3.5

	Map4.Axes("Bottom Axis").Title = "Modified Shepard's Method"
	Map4.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map4.Left = 0.25
	Map4.Top = 8.23

'5. Natural Neighbor
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=5, OutGrid:=Data$+"_NaturalNeighbor.grd")

	Dim Map5 As Object
	Set Map5 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_NaturalNeighbor.grd")

	Map5.yLength = ylength
	Map5.xMapPerPU = Map.yMapPerPU
	If Map5.Width<3.5 Then Map5.xMapPerPU=Map5.yMapPerPU Else Map5.xLength=3.5

	Map5.Axes("Bottom Axis").Title = "Natural Neighbor"
	Map5.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map5.Left = 2.99
	Map5.Top = 8.23

'6. Nearest Neighbor
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=6, OutGrid:=Data$+"_NearestNeighbor.grd")

	Dim Map6 As Object
	Set Map6 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_NearestNeighbor.grd")

	Map6.yLength = ylength
	Map6.xMapPerPU = Map.yMapPerPU
	If Map6.Width<3.5 Then Map6.xMapPerPU=Map6.yMapPerPU Else Map6.xLength=3.5

	Map6.Axes("Bottom Axis").Title = "Nearest Neighbor"
	Map6.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map6.Left = 5.74
	Map6.Top = 8.23

'7. Polynomial Regression
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=7, OutGrid:=Data$+"_PolynomialRegression.grd")

	Dim Map7 As Object
	Set Map7 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_PolynomialRegression.grd")

	Map7.yLength = ylength
	Map7.xMapPerPU = Map.yMapPerPU
	If Map7.Width<3.5 Then Map7.xMapPerPU=Map7.yMapPerPU Else Map7.xLength=3.5

	Map7.Axes("Bottom Axis").Title = "Polynomial Regression"
	Map7.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map7.Left = 0.25
	Map7.Top = 5.35

'8. Radial Basis Function
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=8, OutGrid:=Data$+"_RadialBasisFunction.grd")

	Dim Map8 As Object
	Set Map8 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_RadialBasisFunction.grd")

	Map8.yLength = ylength
	Map8.xMapPerPU = Map.yMapPerPU
	If Map8.Width<3.5 Then Map8.xMapPerPU=Map8.yMapPerPU Else Map8.xLength=3.5

	Map8.Axes("Bottom Axis").Title = "Radial Basis Function"
	Map8.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map8.Left = 2.99
	Map8.Top = 5.35

'9. Triangulation with Linear Interpolation
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=9, OutGrid:=Data$+"_TriangulationLinearInterpolation.grd")

	Dim Map9 As Object
	Set Map9 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_TriangulationLinearInterpolation.grd")

	Map9.yLength = ylength
	Map9.xMapPerPU = Map.yMapPerPU
	If Map9.Width<3.5 Then Map9.xMapPerPU=Map9.yMapPerPU Else Map9.xLength=3.5

	Map9.Axes("Bottom Axis").Title = "Triangulation with Linear Interpolation"
	Map9.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map9.Left = 5.74
	Map9.Top = 5.35

'10. Moving Average
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=10, OutGrid:=Data$+"_MovingAverage.grd")

	Dim Map10 As Object
	Set Map10 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_MovingAverage.grd")

	Map10.yLength = ylength
	Map10.xMapPerPU = Map.yMapPerPU
	If Map10.Width<3.5 Then Map10.xMapPerPU=Map10.yMapPerPU Else Map10.xLength=3.5

	Map10.Axes("Bottom Axis").Title = "Moving Average"
	Map10.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map10.Left = 0.25
	Map10.Top = 2.57

'11. Data Metrics
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=11, DataMetric:=23, _
		SearchEnable:=True, SearchNumSectors:=2, SearchRad1:=2000, SearchRad2:=1500, OutGrid:=Data$+"_DataMetrics.grd")

	Dim Map11 As Object
	Set Map11 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_DataMetrics.grd")

	Map11.yLength = ylength
	Map11.xMapPerPU = Map.yMapPerPU
	If Map11.Width<3.5 Then Map11.xMapPerPU=Map11.yMapPerPU Else Map11.xLength=3.5

	Map11.Axes("Bottom Axis").Title = "Data Metrics: Terrain Slope"
	Map11.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map11.Left = 2.99
	Map11.Top = 2.57

'12. Local Polynomial
	SurferApp.GridData(DataFile:=Data$, xcol:=xcol, ycol:=ycol, zcol:=zcol, ShowReport:=False, Algorithm:=12, OutGrid:=Data$+"_LocalPolynomial.grd")

	Dim Map12 As Object
	Set Map12 = Doc.Shapes.AddContourMap(GridFileName:=Data$+"_LocalPolynomial.grd")

	Map12.yLength = ylength
	Map12.xMapPerPU = Map.yMapPerPU
	If Map12.Width<3.5 Then Map12.xMapPerPU=Map12.yMapPerPU Else Map12.xLength=3.5

	Map12.Axes("Bottom Axis").Title = "Local Polynomial"
	Map12.Axes("Bottom Axis").TitleFont.ForeColorRGBA.Color = srfColorRed
	Map12.Left = 5.74
	Map12.Top = 2.57

End Sub