#include "emdata.h"
#include "ctf.h"
#include "portable_fileio.h"
#include "imageio.h"
#include <cstring>
#include <sstream>
#include <iomanip>
#include <sys/stat.h>
#include <iostream>
Include dependency graph for emdata_metadata.cpp:
Go to the source code of this file.
Functions | |
int | greaterthan (const void *p1, const void *p2) |
vector< Vec3i > | find_region (EMData *image, const vector< Vec3i > &coords, const float value, vector< Vec3i > ®ion) |
|
Definition at line 1335 of file emdata_metadata.cpp. References EMAN::EMData::get_value_at(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), and EMAN::Vec3i. Referenced by EMAN::EMData::mask_contig_region(). 01336 { 01337 static vector<Vec3i> two_six_connected; 01338 if (two_six_connected.size() == 0) { 01339 for(int i = -1; i <= 1; ++i) { 01340 for(int j = -1; j <= 1; ++j) { 01341 for(int k = -1; k <= 1; ++k) { 01342 if ( j != 0 || i != 0 || k != 0) { 01343 two_six_connected.push_back(Vec3i(i,j,k)); 01344 } 01345 } 01346 } 01347 } 01348 } 01349 01350 vector<Vec3i> ret; 01351 for(vector<Vec3i>::const_iterator it = two_six_connected.begin(); it != two_six_connected.end(); ++it ) { 01352 for(vector<Vec3i>::const_iterator it2 = coords.begin(); it2 != coords.end(); ++it2 ) { 01353 if (image->get_value_at((*it2)[0],(*it2)[1],(*it2)[2]) != value) throw; 01354 Vec3i c = (*it)+(*it2); 01355 01356 if ( c[0] < 0 || c[0] >= image->get_xsize()) continue; 01357 if ( c[1] < 0 || c[1] >= image->get_ysize()) continue; 01358 if ( c[2] < 0 || c[2] >= image->get_zsize()) continue; 01359 01360 if( image->get_value_at(c[0],c[1],c[2]) == value ) { 01361 if (find(ret.begin(),ret.end(),c) == ret.end()) { 01362 if (find(region.begin(),region.end(),c) == region.end()) { 01363 region.push_back(c); 01364 ret.push_back(c); 01365 } 01366 } 01367 } 01368 } 01369 } 01370 return ret; 01371 }
|
|
Definition at line 1041 of file emdata_metadata.cpp. Referenced by EMAN::EMData::get_attr(). 01042 { 01043 float* v1 = (float*) p1; 01044 float* v2 = (float*) p2; 01045 01046 if ( *v1 < *v2 ) return 0; 01047 else return 1; 01048 }
|