1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
/* ==================================================================== */
/* Author :Priyanka Hiranandani NIT Surat */
/* ==================================================================== */
/* Syntax : [a,b,c,e,f]=fitellipse((InputArray points);
a,b,c,d,e are parameters of rotated rect returned
/* ==================================================================== */
#include <numeric>
#include "opencv2/core/core.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/opencv.hpp"
#include <iostream>
#include <sciprint.h>
using namespace cv;
using namespace std;
extern "C"
{
#include "api_scilab.h"
#include "Scierror.h"
#include "BOOL.h"
#include <localization.h>
#include "../common.h"
int fitellipse(char *fname, unsigned long fname_len)
{
// Error management variable
SciErr sciErr;
//variable info
int iRows = 0;
int iCols = 0;
int piRows = 0;
int piCols = 0;
int* piAddr = NULL;
int* piAddr2 = NULL;
int* piAddr3 = NULL;
int* piAddr4 = NULL;
int* piAddr5 = NULL;
int* piAddr6 = NULL;
int* piAddrNew = NULL;
int* piLen = NULL;
double* pstData = NULL;
int* pstData2 = NULL;
double *rrows;
double *rcols;
double sigmaSpace;
int error;
////checking input argument
CheckInputArgument(pvApiCtx,1,1);
//checking output argument
CheckOutputArgument(pvApiCtx,5,5);
//for first argument
// get Address of first input
sciErr = getVarAddressFromPosition(pvApiCtx,1,&piAddr);
//check for any error
if(sciErr.iErr)
{
printError(&sciErr, 0);
return 0;
}
//retrieve input array
//SciErr getMatrixOfDouble(void* _pvCtx, int* _piAddress, int* _piRows, int* _piCols, double** _pdblReal)
sciErr = getMatrixOfDouble(pvApiCtx, piAddr, &iRows, &iCols,&pstData);
if(sciErr.iErr)
{
printError(&sciErr, 0);
return 0;
}
int k=0;
vector<Point> contours;
for(int i=0;i<iCols;i++)
{
contours.push_back(Point2f(pstData[i],pstData[i+1])) ;
i++;
}
RotatedRect rect;
rect=fitEllipse(contours);
double size1,size2,center1,center2,angle;
//from rotated rect datatype we can fetch required size,point,angle
Size2f sz; //size2f function for floating point values
Point2f cen;
sz=rect.size;
cen=rect.center;
angle=rect.angle;
size1=sz.width;
size2=sz.height;
center1=cen.x;
center2=cen.y;
error=createScalarDouble(pvApiCtx,nbInputArgument(pvApiCtx)+1,size1);
if(error!=0)
{
sciprint("error occurred");
return 0;
}
error=createScalarDouble(pvApiCtx,nbInputArgument(pvApiCtx)+2,size2);
if(error!=0)
{
sciprint("error occurred");
return 0;
}
error=createScalarDouble(pvApiCtx,nbInputArgument(pvApiCtx)+3,center1);
if(error!=0)
{
sciprint("error occurred");
return 0;
}
error=createScalarDouble(pvApiCtx,nbInputArgument(pvApiCtx)+4,center2);
if(error!=0)
{
sciprint("error occurred");
return 0;
}
error=createScalarDouble(pvApiCtx,nbInputArgument(pvApiCtx)+5,angle);
if(error!=0)
{
sciprint("error occurred");
return 0;
}
AssignOutputVariable(pvApiCtx, 1) = nbInputArgument(pvApiCtx) + 1;
AssignOutputVariable(pvApiCtx, 2) = nbInputArgument(pvApiCtx) + 2;
AssignOutputVariable(pvApiCtx, 3) = nbInputArgument(pvApiCtx) + 3;
AssignOutputVariable(pvApiCtx, 4) = nbInputArgument(pvApiCtx) + 4;
AssignOutputVariable(pvApiCtx, 5) = nbInputArgument(pvApiCtx) + 5;
ReturnArguments(pvApiCtx);
}
}
|
