public static int fnFIE_fit_circle(
	DPNT_T_PTR pnts,
	int pnt_num,
	f_fit_mode fitting_mode,
	double param,
	ref DPNT_T center,
	ref double radius
)

Public Shared Function fnFIE_fit_circle ( 
	pnts As DPNT_T_PTR,
	pnt_num As Integer,
	fitting_mode As f_fit_mode,
	param As Double,
	ByRef center As DPNT_T,
	ByRef radius As Double
) As Integer

//    $Revision: 1.1 $

using System;
using System.Collections.Generic;
using System.Text;
using fvalgcli;

namespace TC
{
    public partial class FIE
    {
        [FvPluginExecute]
        public void fnFIE_fit_circle()
        {
            Console.WriteLine("fnFIE_fit_circle");

            int status = (int)f_err.F_ERR_NONE;

            // 定数.
            const int WIDTH = 256;        // 出力画像の幅.
            const int HEIGHT = 256;        // 出力画像の高さ.
            const int PNTS_NUM = 256;    // 入力点数.
            const double PER = 1;        // 円周の何％に点群を配置するか（1で全周）.

            // ルート画像.
            FHANDLE hdst = FHANDLE.Zero;    // 結果画像.

            DOUBLE_PTR val = DOUBLE_PTR.Zero;    // 画像に描画する際の濃度値.

            F_RANDDESC  r = new F_RANDDESC();    // 擬似乱数列生成用データ. 
            uint seed = 1324;                    // 乱数因子（適当な値）.
            double random;                        // 乱数.

            // 近似円生成パラメータ.
            DPNT_T f_center = new DPNT_T();        // 近似円の中心座標.
            double f_radius = 0;                // 近似円の半径.

            // 点列の作成.
            DPNT_T_PTR pnts = DPNT_T_PTR.Zero;

            // 円のノイズ付加点群生成パラメータ.
            double radius = 100;        // 半径.           
            double radian = Math.PI;    // 角度.
            DPNT_T center;                // 点群描画用円中心.
            center.x = WIDTH/2;            // 画像中心X座標.
            center.y = HEIGHT/2;        // 画像中心Y座標.

            try
            {
                // 点列描画用変数初期化.
                val = DOUBLE_PTR.alloc(3);

                // 乱数生成子の初期化.
                api.fnFIE_mtrand_init(seed, ref r);

                // 点列領域確保.
                pnts = DPNT_T_PTR.alloc(PNTS_NUM);

                // 出力画像の領域を確保します.
                hdst = api.fnFIE_img_root_alloc( (int)f_imgtype.F_IMG_UC8, 3, WIDTH, HEIGHT );

                // 出力画像を白色（255)で塗りつぶす.
                api.fnFIE_img_clear(hdst, 255);

                // 濃度値を設定(点群).
                val[0] = 0;        //R.
                val[1] = 0;        //G.
                val[2] = 0;        //B.

                // 入力点列の生成.
                for(int i = 0; i < PNTS_NUM; i++)
                {
                    // 角度.
                    radian += ((2 * Math.PI) * (PER)) / PNTS_NUM;

                    // 初期座標.
                    pnts[i] = DPNT_T.init(center.x + (radius * Math.Cos(radian)), center.y + (radius * Math.Sin(radian)));

                    // 乱数の生成.
                    random = api.fnFIE_mtrand_real2(ref r);

                    // 乱数を付加.
                    pnts[i] = DPNT_T.init( pnts[i].x + 20 * (random - 0.5) , pnts[i].y + 20 * (random - 0.5) );

                    // 点列を描画する（点描画では小さいので，円で描画).
                    api.fnFIE_draw_circle(hdst, val, f_draw_fill_mode.F_DRAW_FILL_IN, pnts[i], 1);
                }

                // 近似円の計算(Ｍ推定法...標準偏差の自動決定).
                status = api.fnFIE_fit_circle(pnts, PNTS_NUM, f_fit_mode.F_FIT_LSM, 0, ref f_center, ref f_radius);
                Assert.IsTrue(status == (int)f_err.F_ERR_NONE, "エラーが発生しました。({0})", (f_err)status);

                // 近似円を描画する.
                val[0] = 255;    //R
                api.fnFIE_draw_circle(hdst, val, f_draw_fill_mode.F_DRAW_LINE, f_center, f_radius);

                // 画像をPNG形式で保存する.
                api.fnFIE_save_png(ResultDir + "/fnFIE_fit_circle.png", hdst, -1);

                // 結果を出力する.
                ConsoleOut.WriteFunctionName(":\t");
                Console.Write(" ...");
                ConsoleOut.IsTrue(hdst != FHANDLE.Zero);
            }
            finally
            {
                // 解放.
                val.Dispose();
                hdst.Dispose();
                pnts.Dispose();
            }
        }
    }
}

'    $Revision: 1.1 $

Imports System.Collections.Generic
Imports System.Text
Imports fvalgcli

Public Partial Class FIE
    <FvPluginExecute> _
    Public Sub fnFIE_fit_circle()
        Dim status As Integer = CInt(f_err.F_ERR_NONE)

        ' 定数.
        Const  WIDTH As Integer = 256
        ' 出力画像の幅.
        Const  HEIGHT As Integer = 256
        ' 出力画像の高さ.
        Const  PNTS_NUM As Integer = 256
        ' 入力点数.
        Const  PER As Double = 1
        ' 円周の何％に点群を配置するか（1で全周）.
        ' ルート画像.
        Dim hdst As FHANDLE = FHANDLE.Zero
        ' 結果画像.
        Dim val As DOUBLE_PTR = DOUBLE_PTR.Zero
        ' 画像に描画する際の濃度値.
        Dim r As New F_RANDDESC()
        ' 擬似乱数列生成用データ. 
        Dim seed As UInteger = 1324
        ' 乱数因子（適当な値）.
        Dim random As Double
        ' 乱数.
        ' 近似円生成パラメータ.
        Dim f_center As New DPNT_T()
        ' 近似円の中心座標.
        Dim f_radius As Double = 0
        ' 近似円の半径.
        ' 点列の作成.
        Dim pnts As DPNT_T_PTR = DPNT_T_PTR.Zero

        ' 円のノイズ付加点群生成パラメータ.
        Dim radius As Double = 100
        ' 半径.           
        Dim radian As Double = Math.PI
        ' 角度.
        Dim center As DPNT_T
        ' 点群描画用円中心.
        center.x = WIDTH \ 2
        ' 画像中心X座標.
        center.y = HEIGHT \ 2
        ' 画像中心Y座標.
        Try
            ' 点列描画用変数初期化.
            val = DOUBLE_PTR.alloc(3)

            ' 乱数生成子の初期化.
            api.fnFIE_mtrand_init(seed, r)

            ' 点列領域確保.
            pnts = DPNT_T_PTR.alloc(PNTS_NUM)

            ' 出力画像の領域を確保します.
            hdst = api.fnFIE_img_root_alloc(CInt(f_imgtype.F_IMG_UC8), 3, WIDTH, HEIGHT)

            ' 出力画像を白色（255)で塗りつぶす.
            api.fnFIE_img_clear(hdst, 255)

            ' 濃度値を設定(点群).
            val(0) = 0
            'R.
            val(1) = 0
            'G.
            val(2) = 0
            'B.
            ' 入力点列の生成.
            For i As Integer = 0 To PNTS_NUM - 1
                ' 角度.
                radian += ((2 * Math.PI) * (PER)) / PNTS_NUM

                ' 初期座標.
                pnts(i) = DPNT_T.init(center.x + (radius * Math.Cos(radian)), center.y + (radius * Math.Sin(radian)))

                ' 乱数の生成.
                random = api.fnFIE_mtrand_real2(r)

                ' 乱数を付加.
                pnts(i) = DPNT_T.init(pnts(i).x + 20 * (random - 0.5), pnts(i).y + 20 * (random - 0.5))

                ' 点列を描画する（点描画では小さいので，円で描画).
                api.fnFIE_draw_circle(hdst, val, f_draw_fill_mode.F_DRAW_FILL_IN, pnts(i), 1)
            Next

            ' 近似円の計算(Ｍ推定法...標準偏差の自動決定).
            status = api.fnFIE_fit_circle(pnts, PNTS_NUM, f_fit_mode.F_FIT_LSM, 0, f_center, f_radius)
            Assert.IsTrue(status = CInt(f_err.F_ERR_NONE), "エラーが発生しました。({0})", CType(status, f_err))

            ' 近似円を描画する.
            val(0) = 255
            'R
            api.fnFIE_draw_circle(hdst, val, f_draw_fill_mode.F_DRAW_LINE, f_center, f_radius)

            ' 画像をPNG形式で保存する.
            api.fnFIE_save_png(ResultDir & "/fnFIE_fit_circle.png", hdst, -1)

            ' 結果を出力する.
            ConsoleOut.WriteFunctionName(":" & vbTab)
            Console.Write(" ...")
            ConsoleOut.IsTrue(hdst <> FHANDLE.Zero)
        Finally
            ' 解放.
            val.Dispose()
            hdst.Dispose()
            pnts.Dispose()
        End Try
    End Sub
End Class