0

ピクセルごとに (Magick++ を使用して) 動作する画像操作用の C++ コードに取り組んでおり、それを OpenMP で使用したいのですが、次の問題があります。

Magick: Semaphore operation failed (unable to destroy semaphore) [Dispositivo o recurso ocupado].
img_test: magick/pixel_cache.c:2765: ModifyCache: La declaración `image->cache != (Cache) ((void *)0)' no se cumple.

また、無限ループに陥り続けます。

コード スニペットは次のとおりです。

int main(int argc,char **argv)
{
    InitializeMagick(*argv);

    Image img1, img2;
    img1.read(argv[1]);
    img2.read(argv[2]);

    int sx = img1.columns();
    int sy = img1.rows();
    Image out;
    out.size(Geometry(sx,sy));

    cout << "Processing pictures..." << endl;

    int iy;
    #pragma omp for private(iy)
    for (iy=0;iy<sy;iy++)
    {
        #pragma omp parallel for
        for (int ix=0;ix<sx;ix++)
        {
            double _r = 0.0, _g = 0.0, _b = 0.0;

            ColorRGB ppix1(img1.pixelColor(ix,iy));
            ColorRGB ppix2(img2.pixelColor(ix,iy));

            // do some image processing...

            ColorRGB opix(_r*MaxRGB,_g*MaxRGB,_b*MaxRGB);
            out.pixelColor(ix,iy,opix);
        }
    }
    out.write("Output.png");
}

これを解決する方法はありますか？

4

1 に答える 1

1

これを解決する方法はありますか？

この例では、 を使用しますschedule ordered。

cout << "Processing pictures..." << endl;

int iy;
#pragma omp for schedule(static) ordered
for (iy=0;iy<sy;iy++)
{
#pragma omp ordered
    for (int ix=0;ix<sx;ix++)
    {
        double _r = 0.0, _g = 0.0, _b = 0.0;

        ColorRGB ppix1(img1.pixelColor(ix,iy));
        ColorRGB ppix2(img2.pixelColor(ix,iy));

        // do some image processing...

        ColorRGB opix(_r*MaxRGB,_g*MaxRGB,_b*MaxRGB);
        out.pixelColor(ix,iy,opix);
    }
}
out.write("Output.png");

編集

並列全体で低レベルのピクセル情報を操作したい場合は、 @NoseKnowsAll が の単一領域に対して正しいですiy。out.pixelColorただし、内部キャッシュが同期しなくなる可能性があるため、呼び出しに問題が発生します。ピクセル データをエクスポートし、作業を並行して実行し、最終結果をインポートすることをお勧めします。

// Allocate three buffers the total size of x * y * RGB
double * buffer1 = new double[sx * sy * 3];
double * buffer2 = new double[sx * sy * 3];
double * buffer3 = new double[sx * sy * 3];

// Write pixel data to first two buffers
img1.write(0,0, sx, sy, "RGB", DoublePixel, buffer1);
img2.write(0,0, sx, sy, "RGB", DoublePixel, buffer2);

cout << "Processing pictures..." << endl;

int iy;
#pragma omp parallel for
for (iy=0;iy<sy;iy++)
{
    for (int ix=0;ix<sx;ix++)
    {
        // Find where in buffer the current pixel is located at
        size_t idx = (iy * sx + ix) * 3;
        // For fun, let's alternate which source to assing to the
        // third buffer.
        if ((iy % 2 && ix % 2) || (!(iy % 2) && !(ix % 2))) {
            buffer3[idx+0] = buffer1[idx+0]; // R
            buffer3[idx+1] = buffer1[idx+1]; // G
            buffer3[idx+2] = buffer1[idx+2]; // B
        } else {
            buffer3[idx+0] = buffer2[idx+0]; // R
            buffer3[idx+1] = buffer2[idx+1]; // G
            buffer3[idx+2] = buffer2[idx+2]; // B
        }
    }
}
// Import the third buffer into out Image
out.read(sx, sy, "RGB", DoublePixel, buffer3);
out.write("Output.png");

YMMV

ピクセルごとの画像処理

于 2016-03-22T02:07:55.343 に答える