繁体   English   中英

P /调用C#与非托管DLL之间的封送和封送2D数组,结构和指针

[英]P/Invoke marshaling and unmarshalling 2D array, structure and pointers between C# and unmanaged DLL

Flann C ++库具有针对CC ++PythonMatlabRuby的包装器,但没有可用的C#包装器。 我试图围绕从此处下载的flann.dll 32位非托管DLL创建一个C#包装器。

作为PInvoke /编组的新手,我敢肯定我没有正确执行对DLL的C#P / Invoke调用。 我基本上是试图镜像C#中可用的Python包装器 混乱的主要方面是:

  • 我不确定如何在C#中的2D托管矩形数组之间进行编组(输入)和解组(输出),其中参数类型为float*即指向以行主要顺序存储的查询集的指针(根据flann.h中的注释 ) 。
  • 我也不确定如何传递对C的结构引用正确,即struct FLANNParameters*
  • IntPtr是否适合引用typedef void*int* indices

非托管C(flann.dll库)

我需要使用的从flann.h公开导出的C ++方法如下:

typedef void* FLANN_INDEX; /* deprecated */
typedef void* flann_index_t;

FLANN_EXPORT extern struct FLANNParameters DEFAULT_FLANN_PARAMETERS;

// dataset = pointer to a query set stored in row major order
FLANN_EXPORT flann_index_t flann_build_index(float* dataset,
                                             int rows,
                                             int cols,
                                             float* speedup,
                                             struct FLANNParameters* flann_params);

FLANN_EXPORT int flann_free_index(flann_index_t index_id,
                                  struct FLANNParameters* flann_params);

FLANN_EXPORT int flann_find_nearest_neighbors(float* dataset,
                                              int rows,
                                              int cols,
                                              float* testset,
                                              int trows,
                                              int* indices,
                                              float* dists,
                                              int nn,
                                              struct FLANNParameters* flann_params);

托管C#包装器(我的实现)

这是基于上述公开方法的C#包装器。

NativeMethods.cs

using System;
using System.Runtime.InteropServices;

namespace FlannWrapper
{
    /// <summary>
    /// Methods to map between native unmanaged C++ DLL and managed C#
    /// Trying to mirror: https://github.com/mariusmuja/flann/blob/master/src/cpp/flann/flann.h
    /// </summary>
    public class NativeMethods
    {
        /// <summary>
        /// 32-bit flann dll obtained from from http://sourceforge.net/projects/pointclouds/files/dependencies/flann-1.7.1-vs2010-x86.exe/download
        /// </summary>
        public const string DllWin32 = @"C:\Program Files (x86)\flann\bin\flann.dll";

        /// <summary>
        /// C++: flann_index_t flann_build_index(float* dataset, int rows, int cols, float* speedup, FLANNParameters* flann_params)
        /// </summary>
        [DllImport(DllWin32, EntryPoint = "flann_build_index", CallingConvention = CallingConvention.Cdecl, SetLastError = true)]
        public static extern IntPtr flannBuildIndex([In] [MarshalAs(UnmanagedType.LPArray, ArraySubType = UnmanagedType.R4)] float[,] dataset,  // ??? [In] IntPtr dataset ???
                                                    int rows, int cols, 
                                                    ref float speedup,      // ???
                                                    [In] ref FlannParameters flannParams);  // ???

        /// <summary>
        /// C++: int flann_free_index(flann_index_t index_ptr, FLANNParameters* flann_params)
        /// </summary>
        [DllImport(DllWin32, EntryPoint = "flann_free_index", CallingConvention = CallingConvention.Cdecl, SetLastError = true)]
        public static extern int flannFreeIndex(IntPtr indexPtr,        // ???
                                                [In] ref FlannParameters flannParams);   // ??? [In, MarshalAs(UnmanagedType.LPStruct)] FlannParameters flannParams);

        /// <summary>
        /// C++: int flann_find_nearest_neighbors_index(flann_index_t index_ptr, float* testset, int tcount, int* result, float* dists, int nn, FLANNParameters* flann_params)
        /// </summary>
        [DllImport(DllWin32, EntryPoint = "flann_find_nearest_neighbors_index", CallingConvention = CallingConvention.Cdecl, SetLastError = true)]
        public static extern int flannFindNearestNeighborsIndex(IntPtr indexPtr,        // ???
                                                                [In] [MarshalAs(UnmanagedType.LPArray, ArraySubType = UnmanagedType.R4)] float[,] testset,  // ??? [In] IntPtr dataset ???
                                                                int tCount,
                                                                [Out] IntPtr result,    // ??? [Out] [MarshalAs(UnmanagedType.LPArray, ArraySubType = UnmanagedType.R4)] int[,] result, 
                                                                [Out] IntPtr dists,     // ???
                                                                int nn,
                                                                [In] ref FlannParameters flannParams);  // ???
    }
}

FlannTest.cs

using System;
using System.IO;
using System.Runtime.InteropServices;
using Microsoft.VisualStudio.TestTools.UnitTesting;

namespace FlannWrapper
{
    [TestClass]
    public class FlannTest : IDisposable
    {
        private IntPtr curIndex; 
        protected FlannParameters flannParams;
        // protected GCHandle gcHandle;

        [TestInitialize]
        public void TestInitialize()
        {
            this.curIndex = IntPtr.Zero;
            // Initialise Flann Parameters
            this.flannParams = new FlannParameters();  // use defaults
            this.flannParams.algorithm = FlannAlgorithmEnum.FLANN_INDEX_KDTREE;
            this.flannParams.trees = 8;
            this.flannParams.logLevel = FlannLogLevelEnum.FLANN_LOG_WARN;
            this.flannParams.checks = 64;
        }

        [TestMethod]
        public void FlannNativeMethodsTestSimple()
        {
            int rows = 3, cols = 5;
            int tCount = 2, nn = 3;

            float[,] dataset2D = { { 1.0f,      1.0f,       1.0f,       2.0f,       3.0f},
                                   { 10.0f,     10.0f,      10.0f,      3.0f,       2.0f},
                                   { 100.0f,    100.0f,     2.0f,       30.0f,      1.0f} };
            //IntPtr dtaasetPtr = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(float)) * dataset2D.Length);

            float[,] testset2D = { { 1.0f,      1.0f,       1.0f,       1.0f,       1.0f},
                                   { 90.0f,     90.0f,      10.0f,      10.0f,      1.0f} };
            //IntPtr testsetPtr = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(float)) * testset2D.Length);

            int outBufferSize = tCount * nn;
            int[] result = new int[outBufferSize];
            int[,] result2D = new int[tCount, nn];
            IntPtr resultPtr = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(int)) * result.Length);

            float[] dists = new float[outBufferSize];
            float[,] dists2D = new float[tCount, nn];
            IntPtr distsPtr = Marshal.AllocHGlobal(Marshal.SizeOf(typeof(float)) * dists.Length);

            try
            {
                // Copy the array to unmanaged memory.
                //Marshal.Copy(testset, 0, testsetPtr, testset.Length);
                //Marshal.Copy(dataset, 0, datasetPtr, dataset.Length);

                if (this.curIndex != IntPtr.Zero)
                {
                    // n - number of bytes which is enough to keep any type used by function
                    NativeMethods.flannFreeIndex(this.curIndex, ref this.flannParams);
                    this.curIndex = IntPtr.Zero;
                }

                //GC.KeepAlive(this.curIndex);    // TODO

                float speedup = 0.0f;  // TODO: ctype float

                Console.WriteLine("Computing index.");
                this.curIndex = NativeMethods.flannBuildIndex(dataset2D, rows, cols, ref speedup, ref this.flannParams);
                NativeMethods.flannFindNearestNeighborsIndex(this.curIndex, testset2D, tCount, resultPtr, distsPtr, nn, ref this.flannParams);

                // Copy unmanaged memory to managed arrays.
                Marshal.Copy(resultPtr, result, 0, result.Length);
                Marshal.Copy(distsPtr, dists, 0, dists.Length);

                // Clutching straws, convert 1D to 2D??
                for(int row=0; row<tCount; row++)
                {
                    for(int col=0; col<nn; col++)
                    {
                        int buffIndex = row*nn + col;
                        result2D[row, col] = result[buffIndex];
                        dists2D[row, col] = dists[buffIndex];
                    }
                }
            }
            finally
            {
                // Free unmanaged memory -- [BREAKPOINT HERE]
                // Free input pointers
                //Marshal.FreeHGlobal(testsetPtr);
                //Marshal.FreeHGlobal(datasetPtr);
                // Free output pointers
                Marshal.FreeHGlobal(resultPtr);
                Marshal.FreeHGlobal(distsPtr);
            }
        }

        [TestCleanup]
        public void TestCleanup()
        {
            if (this.curIndex != IntPtr.Zero)
            {
                NativeMethods.flannFreeIndex(this.curIndex, ref flannParams);
                Marshal.FreeHGlobal(this.curIndex);
                this.curIndex = IntPtr.Zero;
                // gcHandle.Free();
            }
        }
    }
}

FlannParams.cs

尝试镜像Python FLANNParameters类C struct FLANNParameters

using System;
using System.Runtime.InteropServices;

namespace FlannWrapper
{
    // FieldOffsets set based on assumption that C++ equivalent of int, uint, float, enum are all 4 bytes for 32-bit
    [StructLayout(LayoutKind.Explicit)]
    public class FLANNParameters
    {
        [FieldOffset(0)]
        public FlannAlgorithmEnum algorithm;
        [FieldOffset(4)]
        public int checks;
        [FieldOffset(8)]
        public float eps;
        [FieldOffset(12)]
        public int sorted;
        [FieldOffset(16)]
        public int maxNeighbors;
        [FieldOffset(20)]
        public int cores;
        [FieldOffset(24)]
        public int trees;
        [FieldOffset(28)]
        public int leafMaxSize;
        [FieldOffset(32)]
        public int branching;
        [FieldOffset(36)]
        public int iterations;
        [FieldOffset(40)]
        public FlannCentersInitEnum centersInit;
        [FieldOffset(44)]
        public float cbIndex;
        [FieldOffset(48)]
        public float targetPrecision;
        [FieldOffset(52)]
        public float buildWeight;
        [FieldOffset(56)]
        public float memoryWeight;
        [FieldOffset(60)]
        public float sampleFraction;
        [FieldOffset(64)]
        public int tableNumber;
        [FieldOffset(68)]
        public int keySize;
        [FieldOffset(72)]
        public int multiProbeLevel;
        [FieldOffset(76)]
        public FlannLogLevelEnum logLevel;
        [FieldOffset(80)]
        public long randomSeed;

        /// <summary>
        /// Default Constructor
        /// Ref https://github.com/mariusmuja/flann/blob/master/src/python/pyflann/flann_ctypes.py : _defaults
        /// </summary>
        public FlannParameters()
        {
            this.algorithm = FlannAlgorithmEnum.FLANN_INDEX_KDTREE;
            this.checks = 32;
            this.eps = 0.0f;
            this.sorted = 1;
            this.maxNeighbors = -1;
            this.cores = 0;
            this.trees = 1;
            this.leafMaxSize = 4;
            this.branching = 32;
            this.iterations = 5;
            this.centersInit = FlannCentersInitEnum.FLANN_CENTERS_RANDOM;
            this.cbIndex = 0.5f;
            this.targetPrecision = 0.9f;
            this.buildWeight = 0.01f;
            this.memoryWeight = 0.0f;
            this.sampleFraction = 0.1f;
            this.tableNumber = 12;
            this.keySize = 20;
            this.multiProbeLevel = 2;
            this.logLevel = FlannLogLevelEnum.FLANN_LOG_WARN;
            this.randomSeed = -1;
        }
    }
    public enum FlannAlgorithmEnum  : int   
    {
        FLANN_INDEX_KDTREE = 1
    }
    public enum FlannCentersInitEnum : int
    {
        FLANN_CENTERS_RANDOM = 0
    }
    public enum FlannLogLevelEnum : int
    {
        FLANN_LOG_WARN = 3
    }
}

错误的输出-调试模式,立即窗口

?result2D
{int[2, 3]}
    [0, 0]: 7078010
    [0, 1]: 137560165
    [0, 2]: 3014708
    [1, 0]: 3014704
    [1, 1]: 3014704
    [1, 2]: 48
?dists2D
{float[2, 3]}
    [0, 0]: 2.606415E-43
    [0, 1]: 6.06669328E-34
    [0, 2]: 9.275506E-39
    [1, 0]: 1.05612418E-38
    [1, 1]: 1.01938872E-38
    [1, 2]: 1.541428E-43

如您所见,在Debug模式下运行Test时,我没有收到任何错误,但是我知道输出肯定是不正确的-由于内存寻址不正确而导致垃圾值。 我还包括了我尝试过但未成功的其他封送签名(请参阅带有???的评论)。

地面真理Python(调用PyFlann库)

为了找出正确的结果,我使用可用的Python库PyFlann进行了快速测试。

FlannTest.py

import pyflann
import numpy as np

dataset = np.array(
    [[1., 1., 1., 2., 3.],
     [10., 10., 10., 3., 2.],
     [100., 100., 2., 30., 1.] ])
testset = np.array(
    [[1., 1., 1., 1., 1.],
     [90., 90., 10., 10., 1.] ])
flann = pyflann.FLANN()
result, dists = flann.nn(dataset, testset, num_neighbors = 3, 
                         algorithm="kdtree", trees=8, checks=64)  # flann parameters

# Output
print("\nResult:")
print(result)
print("\nDists:")
print(dists)

在幕后,PyFlann.nn()调用了公开的C方法,正如我们从index.py看到的那样

正确的输出

Result:
[[0 1 2]
 [2 1 0]]

Dists:
[[  5.00000000e+00   2.48000000e+02   2.04440000e+04]
 [  6.64000000e+02   1.28500000e+04   1.59910000e+04]]

对此方法的任何正确帮助将不胜感激。 谢谢。

当使用p / invoke时,您必须停止思考“托管”,而应该考虑物理二进制布局,32 vs 64位等。而且,当被调用的本机二进制始终在进程中运行时(例如此处,但是使用COM服务器可能会有所不同),它比进程外处理更容易,因为您不必考虑封送处理/序列化,引用与输出等问题。

另外,您无需告诉.NET它已经知道的信息。 float数组是R4的LPArray,您不必指定它。 越简单越好。

因此,首先是flann_index_t 它在C中定义为void * ,因此它显然必须是IntPtr (“某物”上的不透明指针)。

然后,构造。 如果将其定义为struct则可以将其作为C中的简单指针传递的结构作为ref参数传递。 如果将其定义为class ,请不要使用ref 通常,我更喜欢将struct用于C结构。

您必须确保结构定义正确。 通常,您使用LayoutKind.Sequential因为.NET p / invoke将以与C编译器相同的方式打包参数。 因此,您不必使用显式的,特别是当参数是标准的(不是一点点的东西),例如int,float时,因此,如果所有成员都已正确声明,则可以删除所有FieldOffset并使用LayoutKind.Sequential ...但这不是案件。

对于类型,就像我说的那样,您真的必须考虑二进制,并问自己使用的每种类型,二进制的布局是什么,大小是多少? int是(使用99.9%C编译器)32位。 float和double是IEEE标准,因此永远不会有问题。 枚举通常基于int ,但这可能有所不同(在C和.NET中,为了能够与C匹配)。 long (有99.0%的C编译器)32位,而不是 64位。 因此,.NET等效为Int32(int),而不是Int64(long)。

因此,您应该更正FlannParameters结构,并将long替换为int 如果您确实想确定给定的结构,请使用与用于编译要调用的库相同的C编译器,针对C的sizeof(mystruct)对照Marshal.SizeOf(mystruct)进行检查。 它们应该是相同的。 如果不是,则说明.NET定义中存在错误(包装,成员大小,订单等)。

这是似乎有效的修改后的定义和调用代码。

static void Main(string[] args)
{
    int rows = 3, cols = 5;
    int tCount = 2, nn = 3;

    float[,] dataset2D = { { 1.0f,      1.0f,       1.0f,       2.0f,       3.0f},
                           { 10.0f,     10.0f,      10.0f,      3.0f,       2.0f},
                           { 100.0f,    100.0f,     2.0f,       30.0f,      1.0f} };

    float[,] testset2D = { { 1.0f,      1.0f,       1.0f,       1.0f,       1.0f},
                           { 90.0f,     90.0f,      10.0f,      10.0f,      1.0f} };

    var fparams = new FlannParameters();
    var index = NativeMethods.flannBuildIndex(dataset2D, rows, cols, out float speedup, ref fparams);

    var indices = new int[tCount, nn];
    var idists = new float[tCount, nn];
    NativeMethods.flannFindNearestNeighborsIndex(index, testset2D, tCount, indices, idists, nn, ref fparams);
    NativeMethods.flannFreeIndex(index, ref fparams);
}

[DllImport(DllWin32, EntryPoint = "flann_build_index", CallingConvention = CallingConvention.Cdecl)]
public static extern IntPtr flannBuildIndex(float[,] dataset,
                                            int rows, int cols,
                                            out float speedup, // out because, it's and output parameter, but ref is not a problem
                                            ref FlannParameters flannParams);

[DllImport(DllWin32, EntryPoint = "flann_free_index", CallingConvention = CallingConvention.Cdecl)]
public static extern int flannFreeIndex(IntPtr indexPtr,  ref FlannParameters flannParams);

[DllImport(DllWin32, EntryPoint = "flann_find_nearest_neighbors_index", CallingConvention = CallingConvention.Cdecl, SetLastError = true)]
public static extern int flannFindNearestNeighborsIndex(IntPtr indexPtr,
                                                        float[,] testset,
                                                        int tCount,
                                                        [In, Out] int[,] result, // out because it may be changed by C side
                                                        [In, Out] float[,] dists,// out because it may be changed by C side
                                                        int nn,
                                                        ref FlannParameters flannParams);

[StructLayout(LayoutKind.Sequential)]
public struct FlannParameters
{
    public FlannAlgorithmEnum algorithm;
    public int checks;
    public float eps;
    public int sorted;
    public int maxNeighbors;
    public int cores;
    public int trees;
    public int leafMaxSize;
    public int branching;
    public int iterations;
    public FlannCentersInitEnum centersInit;
    public float cbIndex;
    public float targetPrecision;
    public float buildWeight;
    public float memoryWeight;
    public float sampleFraction;
    public int tableNumber;
    public int keySize;
    public int multiProbeLevel;
    public FlannLogLevelEnum logLevel;
    public int randomSeed;
}

注意:我尝试使用特定的flann参数值,但是在这种情况下库崩溃了,我不知道为什么...

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM