Проблемы анимации модели COLLADA
у меня есть некоторые проблемы с анимацией загруженной модели COLLADA. Я написал свой собственный парсер, и теперь я также хочу написать свою собственную процедуру рисования. Проблема в том, что как только я включаю анимацию на моей модели, руки, ноги и голова вытягиваются от начала модели. (Загрузчик реализован на основе учебника здесь:COLLADA учебник)
первое, что я делаю в моей функции рисования модели, - это настройка матриц соединений (не это мир матрицы!) с заданными целями из блоков чтения, Если я, например, читаю канал, Как:
<channel source="#some_sampler" target="some_joint/transform(3)(2)"/>
я изменю компонент матрицы (3) (2) из jointMatrix сустава с помощью sid="transform" на этом первом шаге:
if( mCurrentAnimations_.size() > 0 ) {
unsigned currentFrame = GEAR::Root::getSingleton().getFrameEvent().frame;
bool updateTime = false;
if( currentFrame != mLastFrameUpdate_ ) {
if( timeSinceLastFrame < 1.0f )
updateTime = true;
mLastFrameUpdate_ = currentFrame;
}
/****************************************************
* If we have an active animation, *
* we animate it in each of it's defined channels *
***************************************************/
std::list<DAEAnimation*>::iterator it = mCurrentAnimations_.begin();
while( it != mCurrentAnimations_.end() ) {
for( int c = 0; c < (*it)->animation->channels.size(); ++c ) {
// update the time of the channelanimation if requested
if( updateTime ) {
(*it)->channelStates[c].elapsedTime += timeSinceLastFrame;
}
GEAR::COLLADA::Channel* channel = (*it)->animation->channels[c];
// read the two indices depending on the time we're
int firstKeyframeTimeIndex = 0;
int secondKeyframeTimeIndex = 0;
for( int i = 0; i < channel->sampler->inputSource->mFloatArray_->mCount_; ++i ) {
float time = channel->sampler->inputSource->mFloatArray_->mFloats_[i];
if( firstKeyframeTimeIndex == secondKeyframeTimeIndex && time > (*it)->channelStates[c].elapsedTime && i > 0) {
firstKeyframeTimeIndex = i-1;
secondKeyframeTimeIndex = i;
break;
}
if( firstKeyframeTimeIndex == secondKeyframeTimeIndex && i == channel->sampler->inputSource->mFloatArray_->mCount_-1 ) {
(*it)->channelStates[c].elapsedTime = 0.0f;
firstKeyframeTimeIndex = i;
secondKeyframeTimeIndex = 0;
break;
}
}
// look what kind of TargetAccessor we have
if( channel->targetAccessor != NULL && channel->targetAccessor->type == GEAR::MATRIX_ACCESSOR ) {
// ok we have to read 1 value for first and second index
float firstValue = channel->sampler->outputSource->mFloatArray_->mFloats_[firstKeyframeTimeIndex];
float secondValue = channel->sampler->outputSource->mFloatArray_->mFloats_[secondKeyframeTimeIndex];
float firstTime = channel->sampler->inputSource->mFloatArray_->mFloats_[firstKeyframeTimeIndex];
float secondTime = channel->sampler->inputSource->mFloatArray_->mFloats_[secondKeyframeTimeIndex];
float interpolateValue = 1.0f / (secondTime - firstTime) * (secondTime - (*it)->channelStates[c].elapsedTime);
// now we calculate an linear interpolated value
float value = (secondValue*interpolateValue) + (firstValue*(1.0-interpolateValue));
// now we have to write this value to the Joint's Matrix
int entry = ((COLLADA::MatrixTargetAccessor*)channel->targetAccessor)->firstAccessor*4+((COLLADA::MatrixTargetAccessor*)channel->targetAccessor)->secondAccessor;
channel->targetJoint->matrix->jointSpaceMatrix.entries[entry] = channel->targetJoint->matrix->matrix.entries[entry] + value;
}
}
++it;
}
}
после того, как jointMatrices модифицируются всеми каналами, я пересчитываю worldMatrices сустава, вызывая следующую функцию на корневом суставе:
void
COLLADA::Joint::recalcWorldSpaceTransMat() {
GEAR::Mat4 parentMat;
if( parent != NULL )
parentMat = parent->worldSpaceTransformationMatrix;
// @todo Here we have to test against NULL!
if( matrix != NULL )
this->worldSpaceTransformationMatrix = parentMat * matrix->jointSpaceMatrix;
else {
this->worldSpaceTransformationMatrix = parentMat;
}
//std::cout << "Joint " << sid << " recalculatedn";
for( int i = 0; i < mChildJoints_.size(); ++i )
mChildJoints_[i]->recalcWorldSpaceTransMat();
}
теперь все должно быть готово, чтобы нарисовать мой ширина модели следующая последняя часть моей функции рисования:
for( int i = 0; i < mSubMeshes_.size(); ++i ) {
for( int k = 0; k < mSubMeshes_[i]->mSubMeshes_.size(); ++k ) {
// first we animate it
GEAR::DAESubMesh* submesh = mSubMeshes_[i]->mSubMeshes_[k];
submesh->buffer->lock( true );
{
for( unsigned v = 0; v < submesh->buffer->getNumVertices(); ++v ) {
// get the array of joints, which influence the current vertex
DAEVertexInfo* vertexInfo = submesh->vertexInfo[v];
GEAR::Vec3 vertex; // do not init the vertex with any value!
float totalWeight = 0.0f;
for( int j = 0; j < vertexInfo->joints.size(); ++j ) {
Mat4& invBindPoseMatrix = vertexInfo->joints[j]->joint->invBindPoseMatrix;
Mat4& transMat = vertexInfo->joints[j]->joint->worldSpaceTransformationMatrix;
totalWeight += vertexInfo->joints[j]->weight;
vertex += (transMat*invBindPoseMatrix*(submesh->skin->bindShapeMatrix*vertexInfo->vertex))*vertexInfo->joints[j]->weight;
}
if( totalWeight != 1.0f ) {
float normalizedWeight = 1.0f / totalWeight;
vertex *= normalizedWeight;
}
submesh->buffer->bufferVertexPos( v, vertex );
}
}
submesh->buffer->unlock();
mSubMeshes_[i]->mSubMeshes_[k]->buffer->draw( GEAR::TRIANGLES, 0, mSubMeshes_[i]->mSubMeshes_[k]->buffer->getNumVertices() );
}
}
теперь проблема в том, что вывод выглядит следующим образом:
я уверен, что процедура загрузки данных реализована правильно, потому что общая анимация идущего человека видна, но сетка деформирована:
как я уже сказал, Когда я раскомментируйте строку:
channel->targetJoint->matrix->jointSpaceMatrix.entries[entry] = channel->targetJoint->matrix->matrix.entries[entry] + value;
анимация отключена и модель отображается в стандартной поза:
теперь, кроме того, когда я добавляю нормализацию к первым 3 столбцам jointMatrices, как это, прежде чем я пересчитаю worldMatrix сустава:
GEAR::Vec3 row1( matrix->jointSpaceMatrix.entries[0], matrix->jointSpaceMatrix.entries[1], matrix->jointSpaceMatrix.entries[2] );
row1.normalize();
matrix->jointSpaceMatrix.entries[0] = row1.x;
matrix->jointSpaceMatrix.entries[1] = row1.y;
matrix->jointSpaceMatrix.entries[2] = row1.z;
GEAR::Vec3 row2( matrix->jointSpaceMatrix.entries[4], matrix->jointSpaceMatrix.entries[5], matrix->jointSpaceMatrix.entries[6] );
row2.normalize();
matrix->jointSpaceMatrix.entries[4] = row2.x;
matrix->jointSpaceMatrix.entries[5] = row2.y;
matrix->jointSpaceMatrix.entries[6] = row2.z;
GEAR::Vec3 row3( matrix->jointSpaceMatrix.entries[8], matrix->jointSpaceMatrix.entries[9], matrix->jointSpaceMatrix.entries[10] );
row3.normalize();
matrix->jointSpaceMatrix.entries[8] = row3.x;
matrix->jointSpaceMatrix.entries[9] = row3.y;
matrix->jointSpaceMatrix.entries[10] = row3.z;
проблема все еще существует, но на этот раз в другом выходе. Человек сейчас выглядит как инопланетянин :D, но это уменьшает масштабирование:
я не совсем сейчас, правильно ли я сделал нормализацию. Это normaliation действительно нужен? Это не описано в учебник, и я также не смог найти ничего связанного.
В конце концов, я взглянул на реализацию интерполяции в коде со страницы учебника. И: они не используют никакие кватернионы вообще для того чтобы intrpolate матрица отверстия. То, что они делают, это следующее (что не работает для меня):
Mat4 temp;
for (int i = 0; i < 16; ++i)
temp.entries[i] = interpolatef(matrix->jointSpaceMatrixStart.entries[i],matrix->jointSpaceMatrixFinish.entries[i],matrix->delta);
Vec3 forward,up,right,translation;
forward = Vec3(temp.entries[8], temp.entries[9], temp.entries[10]);
up= Vec3(temp.entries[4], temp.entries[5], temp.entries[6]);
right = Vec3(temp.entries[0], temp.entries[1], temp.entries[2]);
forward.normalize();
up.normalize();
right.normalize();
temp.entries[8] = forward.x; temp.entries[9] = forward.y; temp.entries[10] = forward.z;
temp.entries[4] = up.x; temp.entries[5] = up.y; temp.entries[6] = up.z;
temp.entries[0] = right.x; temp.entries[1] = right.y; temp.entries[2] = right.z;
matrix->jointSpaceMatrix = GEAR::Mat4(temp);
тогда я использую кватернионы в другом подходе (также не работает для меня):
// wat we need for interpolation: rotMatStart, rotMatFinish, delta
// create rotation matrices from our 2 given matrices
GEAR::Mat4 rotMatStart = matrix->jointSpaceMatrixStart;
rotMatStart.setTranslationPart( GEAR::VEC3_ZERO );
GEAR::Mat4 rotMatFinish = matrix->jointSpaceMatrixFinish;
rotMatFinish.setTranslationPart( GEAR::VEC3_ZERO );
rotMatStart.transpose();
rotMatFinish.transpose();
// create Quaternions, which represent these 2 matrices
float w = GEAR::Tools::sqr(1.0 + rotMatStart.entries[0] + rotMatStart.entries[5] + rotMatStart.entries[10]) / 2.0;
float w4 = (4.0 * w);
float x = (rotMatStart.entries[6] - rotMatStart.entries[9]) / w4 ;
float y = (rotMatStart.entries[8] - rotMatStart.entries[2]) / w4 ;
float z = (rotMatStart.entries[1] - rotMatStart.entries[4]) / w4 ;
GEAR::Quaternion rotQuadStart(x, y, z, w);
rotQuadStart.normalize();
w = GEAR::Tools::sqr(1.0 + rotMatFinish.entries[0] + rotMatFinish.entries[5] + rotMatFinish.entries[10]) / 2.0;
w4 = (4.0 * w);
x = (rotMatFinish.entries[6] - rotMatFinish.entries[9]) / w4 ;
y = (rotMatFinish.entries[8] - rotMatFinish.entries[2]) / w4 ;
z = (rotMatFinish.entries[1] - rotMatFinish.entries[4]) / w4 ;
GEAR::Quaternion rotQuadFinish(x, y, z, w);
rotQuadFinish.normalize();
// create the interpolated rotation matrix
GEAR::Quaternion slerpedRotQuat = slerp(rotQuadStart, rotQuadFinish, matrix->delta );
slerpedRotQuat.normalize();
GEAR::Mat4 rotMat;
slerpedRotQuat.createMatrix( rotMat );
// interpolate the translation part
GEAR::Vec3 transVecStart(0.0,0.0,0.0);
matrix->jointSpaceMatrixStart.getTranslatedVector3D( transVecStart );
GEAR::Vec3 transVecFinish(0.0,0.0,0.0);
matrix->jointSpaceMatrixFinish.getTranslatedVector3D( transVecFinish );
GEAR::Mat4 transMat;
transMat.setTranslation( transVecFinish*matrix->delta + (transVecStart*(1.0f-matrix->delta)) );
// now write the resulting Matrix back to the Joint
matrix->jointSpaceMatrix = transMat * rotMat;
он также не будет работать для меня. Ничего похоже, сработало. Я действительно понятия не имею, что происходит.
теперь, через 2 дня, я получил его работу благодаря ответу datenwolf
я хочу сообщить всем, как я заставил его работать. Теперь все кажется ясным, и это был только маленький шаг все время. Теперь начнем с анимационной части. Я перебираю все каналы и сохраняю начальное и конечное значения, а также Значение дельты интерполяции в диапазоне 0.0 1.0 для соединения, канал анимирует:
if( mCurrentAnimations_.size() > 0 ) {
unsigned currentFrame = GEAR::Root::getSingleton().getFrameEvent().frame;
bool updateTime = false;
if( currentFrame != mLastFrameUpdate_ ) {
if( timeSinceLastFrame < 1.0f )
updateTime = true;
mLastFrameUpdate_ = currentFrame;
}
/****************************************************
* If we have an active animation, *
* we animate it in each of it's defined channels *
***************************************************/
std::list<DAEAnimation*>::iterator it = mCurrentAnimations_.begin();
while( it != mCurrentAnimations_.end() ) {
for( int c = 0; c < (*it)->animation->channels.size(); ++c ) {
// update the time of the channelanimation if requested
if( updateTime ) {
(*it)->channelStates[c].elapsedTime += timeSinceLastFrame;
}
GEAR::COLLADA::Channel* channel = (*it)->animation->channels[c];
// read the two indices depending on the time we're
int firstIndex = 0;
int secondIndex = 1;
for( int i = 0; i < channel->sampler->inputSource->mFloatArray_->mCount_; ++i ) {
float time = channel->sampler->inputSource->mFloatArray_->mFloats_[i];
if( time > (*it)->channelStates[c].elapsedTime ) {
firstIndex = i-1;
secondIndex = i;
if( firstIndex == -1 ) // set to last frame
firstIndex = channel->sampler->inputSource->mFloatArray_->mCount_ - 1;
break;
}
else if( i == channel->sampler->inputSource->mFloatArray_->mCount_ - 1 ) {
(*it)->channelStates[c].elapsedTime -= channel->sampler->inputSource->mFloatArray_->mFloats_[i];
firstIndex = 0;
secondIndex = 1;
break;
}
}
// look what kind of TargetAccessor we have
if( channel->targetAccessor != NULL && channel->targetAccessor->type == GEAR::MATRIX_ACCESSOR ) {
/************************************************************************
* Matrix accessors, which are read from a COLLADA <channel> block *
* will always target one matrix component they animate. *
* Such accessors are for example: *
* <channel source"#someSource" target="someJoint/transform(0)(2)"/> *
* *
* @TODO: *
* In a pre processing step, we have to group all channels, which *
* operate on the same joint. In order to accelerate the processing of *
* grouped channels, we have to expand the number of keyframes of all *
* channels to the maximum of all channels. *
************************************************************************/
unsigned entry = ((COLLADA::MatrixTargetAccessor*)channel->targetAccessor)->index;
float firstTime = channel->sampler->inputSource->mFloatArray_->mFloats_[firstIndex];
float secondTime = channel->sampler->inputSource->mFloatArray_->mFloats_[secondIndex];
// in case of matrix accessor, we write the startMatrix and the endMatrix to the Joints accessor, who finally will do the animation interpolation
channel->targetJoint->matrix->interpolationRequired = true;
// write out the start and end value to the jointSpaceMatrix
// this matrix will later be interpolated
channel->targetJoint->matrix->jointSpaceMatrixStart.entries[entry] = channel->sampler->outputSource->mFloatArray_->mFloats_[firstIndex];
channel->targetJoint->matrix->jointSpaceMatrixFinish.entries[entry] = channel->sampler->outputSource->mFloatArray_->mFloats_[secondIndex];
// the delta value is in the range [0.0,1.0]
channel->targetJoint->matrix->delta = 1.0f / (secondTime - firstTime) * (secondTime - (*it)->channelStates[c].elapsedTime);
}
}
++it;
}
}
как вы можете видеть, здесь нет интерполяции на всех. Мы просто кэшируем начальные и конечные значения и дельту для всех анимированных суставов (а также устанавливаем флаг на каждом модифицированном суставе)
теперь, когда все анимации сделаны, мы вызываем функцию interpolateMatrices () на всех корневых суставах:
for( int i = 0; i < mSourceModel_->mVisualSceneLibrary_.mVisualScenes_.size(); ++i ) {
for( int v = 0; v < mSourceModel_->mVisualSceneLibrary_.mVisualScenes_[i]->mSkeleton_.size(); ++v ) {
if( mSourceModel_->mVisualSceneLibrary_.mVisualScenes_[i]->mSkeleton_[v]->mRootJoint_ != NULL ) {
/************************************************************************************
* Now we have constructed all jointSpaceMatrixces for the start and the end and *
* we're ready to interpolate them and to also recalculate the joint's *
* worldSpaceMatrix. *
***********************************************************************************/
mSourceModel_->mVisualSceneLibrary_.mVisualScenes_[i]->mSkeleton_[v]->mRootJoint_->interpolateMatrices();
}
}
}
это не ново, но интересной частью теперь является реализация интерполяции. Ничего не мутит кватернионов в все:
void COLLADA::Joint::interpolateMatrices() {
if( matrix != NULL && matrix->interpolationRequired ) {
for (unsigned i = 0; i < 16; ++i)
matrix->jointSpaceMatrix.entries[i] = interpolatef(matrix->jointSpaceMatrixStart.entries[i],matrix->jointSpaceMatrixFinish.entries[i],matrix->delta);
Vec3 forward,up,right,translation;
forward = Vec3(matrix->jointSpaceMatrix.entries[8], matrix->jointSpaceMatrix.entries[9], matrix->jointSpaceMatrix.entries[10]);
up= Vec3(matrix->jointSpaceMatrix.entries[4], matrix->jointSpaceMatrix.entries[5], matrix->jointSpaceMatrix.entries[6]);
right = Vec3(matrix->jointSpaceMatrix.entries[0], matrix->jointSpaceMatrix.entries[1], matrix->jointSpaceMatrix.entries[2]);
forward.normalize();
up.normalize();
right.normalize();
matrix->jointSpaceMatrix.entries[8] = forward.x; matrix->jointSpaceMatrix.entries[9] = forward.y; matrix->jointSpaceMatrix.entries[10] = forward.z;
matrix->jointSpaceMatrix.entries[4] = up.x; matrix->jointSpaceMatrix.entries[5] = up.y; matrix->jointSpaceMatrix.entries[6] = up.z;
matrix->jointSpaceMatrix.entries[0] = right.x; matrix->jointSpaceMatrix.entries[1] = right.y; matrix->jointSpaceMatrix.entries[2] = right.z;
matrix->jointSpaceMatrix.entries[15] = 1.0f; // this component is always 1.0! In some files, this is exported the wrong way, which causes bugs!
}
/********************************************************
* After the interpolation is finished, *
* we have to recalculate the joint's worldSpaceMatrix. *
********************************************************/
GEAR::Mat4 parentMat;
if( parent != NULL )
parentMat = parent->worldSpaceTransformationMatrix;
if( matrix != NULL )
worldSpaceTransformationMatrix = (parentMat * matrix->jointSpaceMatrix);
else
worldSpaceTransformationMatrix = parentMat;
skinningMatrix = worldSpaceTransformationMatrix*invBindPoseMatrix;
// also interpolate and recalculate all childs
for( unsigned k = 0; k < mChildJoints_.size(); ++k )
mChildJoints_[k]->interpolateMatrices();
}
как вы можете видеть, мы просто интраполируем все значения матрицы и после этого нормализуем верхние 3 столбца матрицы. После этого мы немедленно пересчитываем worldSpaceMatrix для этого соединения, а также полную матрицу снятия кожи для сохранения производительности. Теперь у нас почти все. Последнее, что нужно сделать, это действительно оживить вершины, а затем нарисовать сетку:
for( int i = 0; i < mSubMeshes_.size(); ++i ) {
for( int k = 0; k < mSubMeshes_[i]->mSubMeshes_.size(); ++k ) {
// first we animate it
GEAR::DAESubMesh* submesh = mSubMeshes_[i]->mSubMeshes_[k];
submesh->buffer->lock( true );
{
for( unsigned v = 0; v < submesh->buffer->getNumVertices(); ++v ) {
// get the array of joints, which influence the current vertex
DAEVertexInfo* vertexInfo = submesh->vertexInfo[v];
GEAR::Vec3 vertex; // do not init the vertex with any value!
float totalWeight = 0.0f;
for( int j = 0; j < vertexInfo->joints.size(); ++j ) {
totalWeight += vertexInfo->joints[j]->weight;
vertex += ((vertexInfo->joints[j]->joint->skinningMatrix*(vertexInfo->vertex))*vertexInfo->joints[j]->weight);
}
// since it isn't guaranteed that the total weight is exactly 1.0, we have no normalize it
// @todo this should be moved to the parser
if( totalWeight != 1.0f ) {
float normalizedWeight = 1.0f / totalWeight;
vertex *= normalizedWeight;
}
submesh->buffer->bufferVertexPos( v, vertex );
}
}
submesh->buffer->unlock();
mSubMeshes_[i]->mSubMeshes_[k]->buffer->draw( GEAR::TRIANGLES, 0, mSubMeshes_[i]->mSubMeshes_[k]->buffer->getNumVertices() );
}
}
в целом это было почти то же самое, что код я начал. Но теперь для меня все намного яснее, и я могу начать также поддерживать анимации , и . Не стесняйтесь смотреть на мою реализацию в gear3d.de (скачать магистраль SVN)
я надеюсь, что это поможет некоторым людям реализовать свое собственное решение по этой замечательной теме:)
2 ответов
глядя на эти фотографии, у меня создается впечатление, что ваши совместные матрицы не нормализованы, т. е. верхняя левая часть 3×3 масштабирует вашу сетку. Попробуйте, что произойдет, если вы нормализуете верхние левые векторы столбцов 3.
Если это уменьшает проблему, ее необходимо исследовать, какая часть системы анимации вызывает это.
в моем случае все ссылки на компоненты целевой матрицы. Вот почему я интерполирую только один компонент матрицы.
вы никогда не интерполировать матрицы. когда-нибудь.
как это обычно обрабатывается, при загрузке данных анимации вы разлагаете каждую матрицу на кватернион и позицию (и масштаб, если вы анимируете масштаб). Кватернионы используются, потому что они маленькие, легко интерполируются и легко нормализуются после интерполяция. В отличие от матрицы, которая большая, трудно интерполировать, и трудно потом orthonormalize.
обратите внимание, что вышеизложенное обычно выполняется как шаг предварительной обработки в инструменте. Инструмент загружает анимацию Collada, преобразует в кватернионы и позиции, а затем записывает их в формат файла для последующего чтения.
таким образом, вы затем интерполируете кватернионы (не стесняйтесь использовать Lerp для интерполяции внутри анимации) по мере необходимости, с быстрой нормализацией после этого. Позиции нуждаются в обновлении только в том случае, если позиции фактически изменяются относительно исходного смещения. Вы составляете их обратно в матрицу и продолжаете как обычно.
просто и легко.