Stack *ZFlyEmNeuronImageFactory::createSurfaceImage(const ZObject3dScan &obj) const
{
Stack *stack = NULL;
if (!obj.isEmpty()) {
ZObject3dScan tmpObj = obj;
tmpObj.downsampleMax(m_downsampleInterval[0], m_downsampleInterval[1],
m_downsampleInterval[2]);
int offset[3] = { 0, 0, 0 };
tmpObj.switchYZ();
Stack *objStack = tmpObj.toStack(offset);
offset[1] -= 1500 / (m_downsampleInterval[2] + 1); //hard-coded calibration, need modification later
stack = misc::computeNormal(objStack, NeuTube::Z_AXIS);
C_Stack::kill(objStack);
int height= C_Stack::height(stack);
if (m_sizePolicy[2] == SIZE_SOURCE) {
height = m_sourceDimension[2];
if (m_downsampleInterval[2] > 0) {
height /= m_downsampleInterval[2] + 1;
height += 1;
}
}
if (height > C_Stack::height(stack)) {
Stack *out = C_Stack::crop(stack, 0, -offset[1], 0, C_Stack::width(stack),
height, 1, NULL);
C_Stack::kill(stack);
stack = out;
}
}
return stack;
}
Stack* ZFlyEmNeuronImageFactory::createImage(const ZObject3dScan &obj) const
{
Stack *stack = NULL;
if (!obj.isEmpty()) {
ZObject3dScan objProj = obj.makeYProjection();
objProj.downsampleMax(m_downsampleInterval[0], m_downsampleInterval[2], 0);
Cuboid_I boundBox;
objProj.getBoundBox(&boundBox);
int width = 0;
int height = 0;
//int depth = 0;
if (m_sizePolicy[0] == SIZE_SOURCE) {
width = m_sourceDimension[0] / (m_downsampleInterval[0] + 1) + 1;
}
if (m_sizePolicy[2] == SIZE_SOURCE) {
height = m_sourceDimension[2];
if (m_downsampleInterval[2] > 0) {
height /= m_downsampleInterval[2] + 1;
height += 1;
}
}
/*
if (m_sizePolicy[2] == SIZE_SOURCE) {
depth = m_sourceDimension[2];
}
*/
int offset[3] = { 0, 0, 0 };
if (width == 0) {
width = Cuboid_I_Width(&boundBox);
offset[0] = -boundBox.cb[0];
}
if (height == 0) {
height = Cuboid_I_Height(&boundBox);
offset[1] = -boundBox.cb[1];
}
/*
if (depth == 0) {
depth = Cuboid_I_Depth(&boundBox);
offset[2] = -boundBox.cb[2];
}
*/
if (width > 0 && height > 0/* && depth > 0*/) {
stack = C_Stack::make(GREY, width, height, 1);
C_Stack::setZero(stack);
objProj.drawStack(stack, 255, offset);
}
}
return stack;
}