EMAN::TomoObject Class Reference

#include <tomoseg.h>

Public Member Functions
	TomoObject (vector< Vec3i > allpt, float md=1, int slice=0)
void	write_imod (FILE *fp)
int	get_size ()
int	bend ()
void	enlong (EMData *bwmap, EMData *skelmap)

Public Attributes
vector< int >	ptid
vector< int >	segid
vector< Vec3i >	allpoints
float	maxdist
int	nowslice

Detailed Description

Definition at line 35 of file tomoseg.h.

Constructor & Destructor Documentation

TomoObject()

TomoObject::TomoObject	(	vector< Vec3i >	allpt,
		float	md = 1,
		int	slice = 0
	)

Definition at line 36 of file tomoseg.cpp.

                                                              {
        nowslice=slice;
        maxdist=md;
        // find one end point
        for (int i=0; i<(int)allpt.size(); i++){
                if (allpt[i][2]==1){
                        allpoints.push_back(allpt[i]);
                        allpt[i][2]=-1;
                        break;
                }
        }
        
        // sort points from one endpoint to the other
        for (int i=0; i<(int)allpt.size()-1; i++){
                float mindst=1e10;
                int nxti=-1;
                Vec3i nowp=allpoints.back();
                for (int j=0; j<(int)allpt.size(); j++){
                        if (allpt[j][2]>0){
                                float dst=(allpt[j][0]-nowp[0])*(allpt[j][0]-nowp[0])+(allpt[j][1]-nowp[1])*(allpt[j][1]-nowp[1]);
                                if (dst<mindst){
                                        mindst=dst;
                                        nxti=j;
                                }
                        }
                }
                if (nxti>=0){
                        allpoints.push_back(allpt[nxti]);
                        allpt[nxti][2]=-1;
                }
                else
                        break;
        }
//      for (int i=0; i<(int)allpoints.size(); i++)
//              printf("(%d,%d)",allpoints[i][0],allpoints[i][1]);
//      printf("\n");
        ptid.push_back(0);
        ptid.push_back((int)allpoints.size()-1);
        
        
        // break into line segments
        while(1){
                int newpt=bend();
                if (newpt>0){
                        ptid.push_back(newpt);
                        sort(ptid.begin(),ptid.end());
//                      for (int i=0; i<ptid.size(); i++) printf("%d\t",ptid[i]);
//                      printf("\n");
                        
                }
                else
                        break;
        }
        
}

References allpoints, bend(), maxdist, nowslice, and ptid.

Member Function Documentation

bend()

int TomoObject::bend

(

)

Definition at line 92 of file tomoseg.cpp.

                    {
        // distance to current line segments
        float maxd=maxdist;
        int newpt=-1;
        for (int i=1; i<(int)ptid.size(); i++){
                int x1=allpoints[ptid[i-1]][0],y1=allpoints[ptid[i-1]][1];
                int x2=allpoints[ptid[i]][0],y2=allpoints[ptid[i]][1];
                for (int j=ptid[i-1]; j<ptid[i]; j++){
                        int x0=allpoints[j][0],y0=allpoints[j][1];
                        float dst=DistToLine(x0,y0,x1,y1,x2,y2);
                        if (dst>maxd){
                                maxd=dst;
                                newpt=j;
                        }
                }
                
        }
        return newpt;
}

References allpoints, DistToLine(), maxdist, and ptid.

Referenced by TomoObject().

enlong()

void TomoObject::enlong	(	EMData *	bwmap,
		EMData *	skelmap
	)

Definition at line 123 of file tomoseg.cpp.

                                                    {
        // first push the current segment in.
        int myid=skelmap->get_value_at(allpoints[0][0],allpoints[0][1]);
        segid.push_back(myid);
        
        //
        
}

References allpoints, and segid.

get_size()

int TomoObject::get_size

(

)

Definition at line 119 of file tomoseg.cpp.

                        {
        return (int)ptid.size();
}

References ptid.

write_imod()

void TomoObject::write_imod

(

FILE *

fp

)

Definition at line 112 of file tomoseg.cpp.

                                   {
        fprintf(fp,"color 0 1 0 0\nopen\ncontour 0 0 %d\n",(int)ptid.size());
        for (int i=0; i<(int)ptid.size(); i++){
                fprintf(fp,"%d %d %d\n",allpoints[ptid[i]][0],allpoints[ptid[i]][1],nowslice);
        }
}

References allpoints, nowslice, and ptid.

Member Data Documentation

allpoints

vector<Vec3i> EMAN::TomoObject::allpoints

Definition at line 51 of file tomoseg.h.

Referenced by bend(), enlong(), TomoObject(), and write_imod().

maxdist

float EMAN::TomoObject::maxdist

Definition at line 52 of file tomoseg.h.

Referenced by bend(), and TomoObject().

nowslice

int EMAN::TomoObject::nowslice

Definition at line 53 of file tomoseg.h.

Referenced by TomoObject(), and write_imod().

ptid

vector<int> EMAN::TomoObject::ptid

Definition at line 49 of file tomoseg.h.

Referenced by bend(), get_size(), TomoObject(), and write_imod().

segid

vector<int> EMAN::TomoObject::segid

Definition at line 50 of file tomoseg.h.

Referenced by enlong().

---

The documentation for this class was generated from the following files:

• tomoseg.h
• tomoseg.cpp