提取刻字中心线源代码

1821117007 · 发表于 2016-10-28 04:35

/*****************************************************************************
**
** 提取刻字中心线【这个程序网上有论文非梅雷原创】【但是大概按照论文写了下效果一般般】
**
** 梅雷QQ1821117007 手机15757753770
** 精力有限，
*****************************************************************************/

/* Include files */
#if ! defined ( __hp9000s800 ) && ! defined ( __sgi ) && ! defined ( __sun )
# include <strstream>
# include <iostream>
using std:

strstream;
using std::endl;
using std::ends;
using std::cerr;
#else
# include <strstream.h>
# include <iostream.h>
#endif
#include <uf.h>
#include <uf_ui.h>
#include <uf_exit.h>
#include <UF_KF.h>
#include <UF_MODL.h>
#include <UF_ASSEM.h>
#include "uf头文件.h"
#include <uf_obj.h>
#include <stdio.h>
#include <math.h>

//设置选择的实体类型
static int M021_init_proc( UF_UI_selection_p_t select, void *user_data )
{
int  errorCode = 0;
int  num_triples = 1; //选择类型数量
UF_UI_mask_t mask_triples[]={22,0,0}; //定义选择类型
errorCode = UF_UI_set_sel_mask(select, UF_UI_SEL_MASK_CLEAR_AND_ENABLE_SPECIFIC,num_triples, mask_triples);
if( errorCode == 0)
{
return UF_UI_SEL_SUCCESS;
}
else
{
return UF_UI_SEL_FAILURE;
}
}

#define UF_CALL(X) (report_error( __FILE__, __LINE__, #X, (X)))

static int report_error( char *file, int line, char *call, int irc)
{
if (irc)
{
      char err[133],
         msg[133];

      sprintf(msg, "*** ERROR code %d at line %d in %s:\n+++ ",
         irc, line, file);
      UF_get_fail_message(irc, err);

      UF_print_syslog(msg, FALSE);
      UF_print_syslog(err, FALSE);
      UF_print_syslog("\n", FALSE);
      UF_print_syslog(call, FALSE);
      UF_print_syslog(";\n", FALSE);

      if (!UF_UI_open_listing_window())
      {
         UF_UI_write_listing_window(msg);
         UF_UI_write_listing_window(err);
         UF_UI_write_listing_window("\n");
         UF_UI_write_listing_window(call);
         UF_UI_write_listing_window(";\n");
      }
}

return(irc);
}

/*****************************************************************************
**  Activation Methods
*****************************************************************************/
/*  Unigraphics Startup
**    This entry point activates the application at Unigraphics startup */
extern DllExport void ufsta( char *param, int *returnCode, int rlen )
{
/* Initialize the API environment */
if( UF_CALL(UF_initialize()) )
{
      /* Failed to initialize */
      return;
}

/* TODO: Add your application code here */

double 建模距离公差=0.001;
double 建模角度公差=0.001;
UF_MODL_ask_distance_tolerance(&建模距离公差);
UF_MODL_ask_angle_tolerance(&建模角度公差);

double  d_tolerance=0.001;
double  a_tolerance=0.001;
int err_flag=0;
UF_MODL_set_distance_tolerance(&d_tolerance,&err_flag);
UF_MODL_set_angle_tolerance(&a_tolerance,&err_flag);

double ZOR[20000];
double ORZJ=0;
double MAXOR=0;
int 点数=0;
double 点[20000][3];
double ijk[3][3];  double  matrix[9];
double p[2000][2][3]; int n=0;
double ap[20000][3]; int an=0;
double uvs[4];
double uv[2];
double point[ 3 ];
double u1[ 3 ];
double v1[ 3 ];
double u2[ 3 ];
double v2[ 3 ];
double radii[ 2 ];

  int err,i,j,k;double u,bb,m;
  double tangent[ 3 ], p_norm[ 3 ],b_norm[ 3 ],torsion,rad_of_cur;
  char message[132]="梅雷提示,请选择刻字平面";
  char title[132]="梅雷-QQ1821117007";
  int  scope=UF_UI_SEL_SCOPE_WORK_PART;
  int  response;
  tag_t object_tag=NULL_TAG;
  double  cursor[3]={0,0,0};
  tag_t view=NULL_TAG;
  UF_UI_lock_ug_access (UF_UI_FROM_CUSTOM);  //加锁
  UF_UI_select_with_single_dialog(message,title,scope,M021_init_proc,NULL,&response,&object_tag,cursor,&view);  //选择框】
  UF_DISP_set_highlight(object_tag,0);    //1高亮显示 0不高亮显示
  UF_UI_unlock_UG_access (UF_UI_FROM_CUSTOM);//解锁
  if (response==UF_UI_OBJECT_SELECTED) //---
  {

   UF_MODL_ask_face_uv_minmax (object_tag, uvs);
   uv[0]=uvs[0]+(uvs[1]-uvs[0])*0.5;
   uv[1]=uvs[2]+(uvs[3]-uvs[2])*0.5;
   UF_MODL_ask_face_props (object_tag,uv,point,u1,v1,u2,v2,ijk[0],radii); //ijk[2] 获得平面方向

   //…………………………………………………………………………………………………………
   tag_t point_tag=NULL_TAG;
   tag_t box_tag=NULL_TAG;
   uf_list_p_t edges_list;
   UF_MODL_create_list(&edges_list);
   UF_MODL_ask_face_edges(object_tag,&edges_list);
   int b=0; err=UF_MODL_ask_list_count(edges_list,&b);
   if (err!=0 || b<1)
   {
      uc1601("无效的表面",1);
   goto LHALT;
   }

   if (b>200)
   {
      uc1601("表面边缘超过200条",1);
   goto LHALT;
   }

   n=0;an=0;
   for (i=0;i<b;i++) //-//
   {
   int edge_type;
      UF_MODL_ask_list_item(edges_list,i,&box_tag);
   UF_MODL_ask_edge_type(box_tag,&edge_type);
   bb=10; //线分成线段
      for (j=1;j<bb+1;j++)
   {
   n++;
   u=(j-1)/bb;
      UF_MODL_ask_curve_props (box_tag,u,p[n-1][0],tangent, p_norm, b_norm, &torsion, &rad_of_cur );
            u=j/bb;
   UF_MODL_ask_curve_props (box_tag,u,p[n-1][1],tangent, p_norm, b_norm, &torsion, &rad_of_cur );
   }
      for (j=0;j<100;j++) //每条直线取点 10*10
   {
   u=j*0.01;
   an++;
       if (u!=1 && u!=0.5) UF_MODL_ask_curve_props (box_tag,u,ap[an-1],tangent, p_norm, b_norm, &torsion, &rad_of_cur );

// if (u!=1) UF_CURVE_create_point(ap[an-1],&point_tag);//创建点
   }


   }//-//
   //…………………………………………………………………………………………………………
   for (i=0;i<n;i++)
   {
   double OR=0.01; //初始化球半径，不能超过文字的最小比划宽度也不能太小
   ijk[1][0]=p[i][1][0]-p[i][0][0];
   ijk[1][1]=p[i][1][1]-p[i][0][1];
   ijk[1][2]=p[i][1][2]-p[i][0][2];
   UF_VEC3_unitize(ijk[1],0.0001,&m,&ijk[1][0]);
   UF_MTX3_initialize (ijk[0],ijk[1],matrix);
   ijk[0][0]=matrix[0];
   ijk[0][1]=matrix[1];
   ijk[0][2]=matrix[2];
   ijk[1][0]=matrix[3];
   ijk[1][1]=matrix[4];
   ijk[1][2]=matrix[5];
   ijk[2][0]=matrix[6];
   ijk[2][1]=matrix[7];
   ijk[2][2]=matrix[8];
   double point1[3],point2[3];int pt_status;
   UF_VEC3_midpt(p[i][0],p[i][1],point1);
        point2[0]=point1[0]+OR*ijk[2][0];
      point2[1]=point1[1]+OR*ijk[2][1];
      point2[2]=point1[2]+OR*ijk[2][2];
   UF_MODL_ask_point_containment(point2,object_tag,&pt_status);
      if (pt_status==2)
   {
      ijk[2][0]=-ijk[2][0];
      ijk[2][1]=-ijk[2][1];
      ijk[2][2]=-ijk[2][2];
        point2[0]=point1[0]+OR*ijk[2][0];
      point2[1]=point1[1]+OR*ijk[2][1];
      point2[2]=point1[2]+OR*ijk[2][2];
      UF_MODL_ask_point_containment(point2,object_tag,&pt_status);
if (pt_status==2)goto LL;
   }
//UF_CURVE_create_point(point2,&point_tag);//创建点


   k=0;
   L10:;
   k++;
   OR=OR*2;
            point2[0]=point1[0]+OR*ijk[2][0];
      point2[1]=point1[1]+OR*ijk[2][1];
      point2[2]=point1[2]+OR*ijk[2][2];
      for (j=0;j<an;j++) //判断所有点是否在圆内
      {
double p[6]={point2[0],point2[1],point2[2],ap[j][0],ap[j][1],ap[j][2]};
double ab=fabs(sqrt((p[0]-p[3])*(p[0]-p[3])+(p[1]-p[4])*(p[1]-p[4])+(p[2]-p[5])*(p[2]-p[5])));
if (ab<OR-0.00005)
{
   OR=OR/2;
            point2[0]=point1[0]+OR*ijk[2][0];
      point2[1]=point1[1]+OR*ijk[2][1];
      point2[2]=point1[2]+OR*ijk[2][2];
   goto  L11;
}
      }
   if (k<10+1) goto L10;
L11:;
   if (k<2) goto LA;
   double RR=OR;
   k=0;
   L20:;
   k++;
   OR=OR+RR*pow(0.5,k);
            point2[0]=point1[0]+OR*ijk[2][0];
      point2[1]=point1[1]+OR*ijk[2][1];
      point2[2]=point1[2]+OR*ijk[2][2];
      for (j=0;j<an;j++) //判断所有点是否在圆外
      {
double p[6]={point2[0],point2[1],point2[2],ap[j][0],ap[j][1],ap[j][2]};
double ab=fabs(sqrt((p[0]-p[3])*(p[0]-p[3])+(p[1]-p[4])*(p[1]-p[4])+(p[2]-p[5])*(p[2]-p[5])));
if (ab<OR-0.00005)
{
        OR=OR-RR*pow(0.5,k);
            point2[0]=point1[0]+OR*ijk[2][0];
         point2[1]=point1[1]+OR*ijk[2][1];
         point2[2]=point1[2]+OR*ijk[2][2];
goto  LB;
      }
      }
   LB:;
   if (k<14+1) //不能超过15
   {
      goto L20;
   }
L21:;
   if (OR<0.028) goto LA;//最小字体间距
      //---取消2个距离小的点
   if (点数>1)
      {
         for (int ii=0;ii<点数;ii++)
         {
         double p[6]={point2[0],point2[1],point2[2],点[ii][0],点[ii][1],点[ii][2]};
         double ab=fabs(sqrt((p[0]-p[3])*(p[0]-p[3])+(p[1]-p[4])*(p[1]-p[4])+(p[2]-p[5])*(p[2]-p[5])));
   if (ab<0.01) goto LA;
      }
      }
   点数++;
   点[点数-1][0]=point2[0];
   点[点数-1][1]=point2[1];
   点[点数-1][2]=point2[2];
   ZOR[点数-1]=OR;
   if (OR>MAXOR) MAXOR=OR;
   ORZJ=ORZJ+OR;
   LA:;
   //

   LL:;
   }
//------------------------------
   // 过滤偏离点
   if (点数>1)
      {
         for (int jj=0;jj<点数;jj++)
{
if (ZOR[jj]>MAXOR/3 && ZOR[jj]>(ORZJ/点数)*(0.75))
{
   UF_CURVE_create_point(点[jj],&point_tag);//创建点
            /*
   if (jj>0)
   tag_t line_tag=NULL_TAG;
   {
          UF_CURVE_line_t line_coords;
          line_coords.start_point[0]=点[jj-1][0];
          line_coords.start_point[1]=点[jj-1][1];
          line_coords.start_point[2]=点[jj-1][2];
          line_coords.end_point[0]=点[jj][0];
          line_coords.end_point[1]=点[jj][1];
          line_coords.end_point[2]=点[jj][2];
          int err=UF_CURVE_create_line(&line_coords,&line_tag);
   }
            */
}
         }
      }

  }//---
  LHALT:;

  UF_MODL_set_distance_tolerance(&建模距离公差,&err_flag);
  UF_MODL_set_angle_tolerance(&建模角度公差,&err_flag);

/* Terminate the API environment */
UF_CALL(UF_terminate());
}

/*****************************************************************************
**  Utilities
*****************************************************************************/

/* Unload Handler
**    This function specifies when to unload your application from Unigraphics.
**    If your application registers a callback (from a MenuScript item or a
**    User Defined Object for example), this function MUST return
**    "UF_UNLOAD_UG_TERMINATE". */
extern int ufusr_ask_unload( void )
{
return (UF_UNLOAD_IMMEDIATELY);
}

gogocat1233 · 发表于 2017-2-5 09:43

我是UG爱好者，看了LZ的帖子，我只想说一句楼上的贴很好！UG论坛好强大！

		自动登录	找回密码
密码			注册

[分享] 提取刻字中心线 源代码

[分享] 提取刻字中心线源代码