Option Explicit

Dim StdDev 
Dim Blur

Const mxErosionKF = 9
Const mxUserKF = 10
Const mxGaussianBlur = 11

Dim mxapp	' the MaximDL application object
Dim ImageNum

''''''''''''
'Const imagename = "C:\d\sxv\05-06-26\m16HaRGB.tif"
'Dim imagename
'imagename=inputbox("Image name: ")
StdDev=CDbl(inputbox("Standard Deviation of Background: "))
Blur=CDbl(inputbox("Gaussian Kernel Blur: "))
''''''''''''''

'nKern = int((fwhm-2.184)/0.21)

Dim image
Dim result
Dim image_high_vals
Dim image_low_vals
Dim background
Dim final_result


Set mxapp = CreateObject("Maxim.Application")

Set image = mxapp.CurrentDocument

Set result = image.Duplicate
result.DisplayName = "Result"


	Set background = result.Duplicate
	background.DisplayName = "Background"
	background.KernelFilter mxErosionKF, 7   ' 7x7 erosion
	background.KernelFilter mxGaussianBlur, 3  ' Gaussian blur 3 pixel rad
	background.AddConstant(5*StdDev)

	Set image_high_vals = result.Duplicate
	image_high_vals.DisplayName = "Image high_vals"
	image_high_vals.Subtract(background)

	Set image_low_vals = image_high_vals.Duplicate
	image_low_vals.DisplayName = "Image low_vals"
	image_low_vals.Add(background)
	image_low_vals.Subtract(result)

	image_low_vals.KernelFilter mxGaussianBlur, Blur  

	result.Close
	Set final_result = image_high_vals.Duplicate
	final_result.Add(background)
	final_result.Subtract(image_low_vals)
	final_result.DisplayName = "Neat Image "+ image.DisplayName

	image_low_vals.Close
	image_high_vals.Close
	background.Close

Dim flag
flag = final_result.SetFITSKey("HISTORY", "Neat Image/Vanderbei.")