1 MAKE_CUBE Output a filled cube from an uncomplete sampled cube MAKE_CUBE This is an image construction task which is able to produce a filled image from one containing many blanked pixels. The reconstructed filled image is not constrained to fit exactly the observed data points. On the opposite, the construction is made by the analogy to a flexible plate attached to fixed points by springs: the plate is the analogous of the surface represented by the image, and the fixed points are the analogous of the observed data points. By adjusting the parameter P P = (plate stiffness) / (springs stiffness) it is possible to control the fidelity to the original data and the amount of smoothing involved in the image reconstruction. Low values of P mean high fidelity to observed data, and negligible amount of smoothing. The original grid is first expanded by a factor EXPANSION\$, new pix- els being attributed the blanking value. Then, the minimization proceeds iteratively to adjust the final image, until convergence is reached. Initially blanked pixels are ignored in the convergence criterium. The algorithm works on cubes, processing each plane independently. It can be used as an alternative to FILL_CUBE for oversampled images, but works also for undersampled data. 2 P\$ TASK\REAL "Smoothing parameter" P\$ This is the value of the smoothing parameters. It should remain be- tween 0.001 and 10. The higher values give stronger smoothing, and should be used for noisy data. 2 NITER\$ TASK\INTEGER "Number of smoothing iterations" NITER\$ This is the number of smoothing iterations. Ten iterations are usu- ally enough to reach convergence, unless the expansion factor is large. 2 EXPANSION\$ TASK\INTEGER "Expansion factor" EXPANSION\$ This is the expansion factor for the output cube, i.e. the number of output pixels for each input pixel along each direction. Minimum value is 2, maximum is 10. Note that the computing time is proportional to the number of output pixels, thus to EXPANSION\$**2. 2 IN\$ TASK\FILE "Input file name" IN\$ This is the name of the input file. 2 OUT\$ TASK\FILE "Output file name" OUT\$ This is the name of the output smoothed image. 2 GUESS\$ TASK\REAL "Initial value for blanked pixels" GUESS\$ This is an initial guess for the original value attributed to blanked pixels. Try to find something reasonable, usually the average of valid pixels. Convergence can be delayed if you choose a crazy value, since the minimization tends to produce a smooth image, but the result doesnot depend on GUESS\$. 1 ENDOFHELP