ich habe folgende Rechnung:

double dU = dPhi / (2.*M_PI);

Wenn ich dann mein Programm durchlafen lasse bekomme ich einen Kleinen Fehler in meinem Bild.

Wenn ich aber M_PI durch 3.2 ersetze dann hab ich ein Perfektes Bild.

Das kann doch eigentlich nur ein Rundungsfehler sein oder??

Und ich wollte fragen wie ich diesen Fehler behebenkann und mir M_PI trotzdem erhalten bleibt.

Danke

Ansonsten ist die Frage schwer zu beantworten , da wir dein "Bild" nicht sehen können und auch nicht erahnen können, was du damit meinst und wie es sich berechnet. Ein paar mehr Infos wären also schon noch nötig.

http://saiyamans-world.pytalhost.de/Rundungsfehler.htm

warum steht da dauernt googel

muss pytalhost heisen

das 1te ist das ohne fehler.

das 2te ist das eingangsbild.

und das dritte ist das mit fehler.

Hoffe das Hilft bei der Lösung.

Danke für die mühe.

for (int y = 0; y < nHeightDst; y++)
    {
      for (int x = 0; x < nWidthDst; x++)
      {
        double dSx = x - dCenterX;

        double dSy = dCenterY - y;

        double dRadius = sqrt(double(dSx*dSx + dSy*dSy));

        double dPhi = atan2(dSy, dSx);

        if (dPhi<0)
        {
          dPhi+= 2.*M_PI;
        }

        if (dPhi>2.*M_PI)
        {
          dPhi-= 2.*M_PI;
        }

        double dU = dPhi / (2.*3.2);

        double dV = dRadius / maxR;

        int i = y*nWidthDst + x;

        if(dRadius > maxR)
        {
          pxMap[i] = -1;

          pyMap[i] = -1;
        }

        pxMap[i] = dU * srcSize.width;

        pyMap[i] = dV * srcSize.height;
      }
    }

Bitteschön

mh, versuch mal das 2. duch 2.00000000 oder so zu ersetzten oder durch static_cast<double>(2.0).

Aber nur bei Vollmond!

Das bringt leider nix.

Danke für die mühe.

Hast Du mal etwas mehr Code? Mit Der Schleife, die zur Generierung durchlaufen wird z.B.?
Wie ist M_PI genau definiert?

Ich hab die math.h includet wo M_PI ja drin ist.

Hmm, Standard ist das aber nicht.

ich definiere Pi immer so:

const double PI = acos(-1.0);

das kann ich ja mal Probieren.

Danke

Edit:

Ohne Veränderung.

if(dRadius > maxR)
{
  pxMap[i] = -1;
  pyMap[i] = -1;
}
pxMap[i] = dU * srcSize.width;
pyMap[i] = dV * srcSize.height;

ist wohl so schon mal richtiger:

if(dRadius > maxR)
{
  pxMap[i] = -1;
  pyMap[i] = -1;
}
else{
  pxMap[i] = dU * srcSize.width;
  pyMap[i] = dV * srcSize.height;
}

Sollte aber dein Problem nicht lösen.
Ansonsten merk dir doch mal eine fehlerhafte Pixelposition im Bild. Deine Schleifen iterieren ja durch dein Zielbild, kannst also beim entsprechenden x,y mit dem Debuger nachschauen, was an der Stelle passiert.
Aber eigentlich bin ich auch der Meinung, dass das Problem nicht in dem Codeteil steckt, den wir hier von dir zu sehen bekommen.

double M_PI = acos(-1);

    MFFTiff source("bild2.tif");

    unsigned char *ucS  = source.read();

    int nHeightDst = 600;
    int nWidthDst = 600;
    int dCenterX = nWidthDst/2;
    int dCenterY = nHeightDst/2;
    int maxR = 300;
    unsigned char *ucT = new unsigned char[nWidthDst*nHeightDst*3];

    for (int y = 0; y < nHeightDst; y++){
      for (int x = 0; x < nWidthDst; x++){
        double dSx = x - dCenterX;
        double dSy = dCenterY - y;
        double dRadius = sqrt(double(dSx*dSx + dSy*dSy));
        double dPhi = atan2(dSy, dSx);

        if (dPhi<0){
          dPhi+= 2.*M_PI;
        }

        if (dPhi>2.*M_PI){
          dPhi-= 2.*M_PI;
        }

        //double dU = dPhi / (2.*3.2);// Edit
        double dU = dPhi / (2.*M_PI);// Edit
        double dV = dRadius / maxR;
        int i = (y*nWidthDst + x)*3;

        if(dRadius > maxR){

        }
	else{
	   //int ix = dU * (source.width()-1);// Edit
	   //int iy = dV * (source.height()-1);// Edit
           int ix = dU * source.width();// Edit
	   int iy = dV * source.height();// Edit

	   ucT[i] = ucS[(ix+iy*source.width())*3];
	   ucT[i+1] = ucS[(ix+iy*source.width())*3+1];
	   ucT[i+2] = ucS[(ix+iy*source.width())*3+2];
	 }
      }
    }

    MFFTiff target("result.tif",nWidthDst,nHeightDst,MRGB);
    target.write(ucT,nWidthDst*nHeightDst*3);

Edit: Hab ein paar Fehler mit übernommen. Aber auch nach Korrektur schliesst sich der Kreis.

Den Fehler hab ich jetzt gefunden:

Ist 100%tig ein Rundungsfehler der dadurch koriegiert werden kann.

else
        {
          pxMap[i] = dU * srcSize.width - 1;

          if(pxMap[i] < 0)
          {
            pxMap[i] = 0;
          }

          pyMap[i] = dV * srcSize.height;
        }

Aber damit du auch mal den Rest vom Code siehst... ist nur ein File Dialog.. und Definition für einige Variablen die ich dann für eine Funktion brauche. poste ich dir mal die Funktion.

ippiRemap_8u_C4R(pSrc, srcSize, srcStep, srcRoi, pxMap, xMapStride, pyMap, yMapStride, pDst, dstStep, dstRoiSize, nInterpolation);

Weis nicht obs dir was bringt aber Was am Code jetzt noch fehlen Würde wäre die Definition der Variablen für dieses Funktion und der FileDialog.

Also vielen Dank an alle und dir MiP ein extra dank dafür das du dich so reihngehängt hast.

Mit was erstellst du eigentlich den Kreis??

Also ich ja mit der ippiRemap weil dein MFFTiff sagt mir reihn garnichts.