Выравнивающие биты в BMP

void rangeFilter(BYTE *in,BYTE * out,int w,int h)
{
    
    int k;
    cout<<"  k= ";  
    k=3;//cin>>k; //k параметр рангового фильтра
    cout<<endl;
 
    DWORD matrixDim=3;//getMatrixSize("Enter matrix size ( Number >=3) "); 
    DWORD * matrix = new DWORD [matrixDim * matrixDim]; //матрица апертуры
    DWORD *mPtr=matrix;
    cout<<"Enter aperture elements:"<<endl;
 
    DWORD apertureArraySize=0; //размер апертурного массива для сортировки
    for (DWORD _i =0; _i<matrixDim*matrixDim; )//ввод матрицы
    {
        //cin>>*mPtr;
        *mPtr=1;
        if ((_i++ % matrixDim)==0)       cout<<endl;
 
        if (*mPtr<0)
        {
            cout<<"The number must be >=0"<<endl;
            --_i;
            continue;
        }
        apertureArraySize+=*mPtr;
        mPtr++;     
    }
    if (apertureArraySize==0)
    {
        delete []matrix;
        return;
    }
 
    BYTE * apertureRow=new BYTE[apertureArraySize];
    DWORD halfDim=matrixDim/2;
    DWORD wsize=(halfDim*2+1)*(halfDim*2+1);    
    int nuli=0;
        if ((3*w)%4!=0) nuli=4-(3*w)%4;
 
    for (unsigned y = 0; y < h; ++y)
    {
        for (unsigned x = 0; x < w; ++x)
        {
            if ((y<halfDim)||(y>=(h-halfDim))||(x<halfDim)||(x>=(w-halfDim)))//обработка краев
                continue;
 
            //Обегаем соседние точки и заполняем апертурный ряд
            BYTE * _apertureRow=apertureRow;            
            mPtr=matrix;
            for (int _y=y-halfDim;_y<=y+halfDim;_y++)
            {
                for (int _x=x-halfDim;_x<=x+halfDim;_x++)
                {                   
                    for (unsigned int i=*mPtr;i>0;i--)//заполнение массива без сортировки                 
                        *_apertureRow++=in[(_y) * w + (_x)];
                    
                    mPtr++;
                }
            }
            
            qsort ((void*)(apertureRow), apertureArraySize, sizeof(unsigned char),
                compareUchar); //сортируем апертурный массив
 
            if (abs(apertureRow[k+1]-in[y*w+x]) < (abs(apertureRow[wsize-k]-in[y*w+x])))//
                out[y*(w+nuli)+x]=apertureRow[k+1];
            else
                out[y*(w+nuli)+x]=(apertureRow[wsize-k]);
        }
 
        
    
        for(int i=1;i<=nuli;i++)  out[y*(w+nuli)+w+i]=0;
        
        //for(int i=0;i<=(w+nuli);i++) cout<<i;
        cout<<(w+nuli)<<endl;
        for(int i=0;i<=(w+nuli);i++)  cout<<out[y*w+i];
        cout<<endl;
    }
    delete [] apertureRow;
    delete [] matrix;
}