/////////////////////////////////////////////////////////////////////////////
// Erzeugen eines VT_UI1-SAFEARRAY aus einem Koordinatenfeld
HRESULT CreateSABlobFromVertices (
double_v &rX, double_v &rY, POINTBASE *pNormal, void *pData, REFCLSID rClsId, SAFEARRAY **pSA)
{
if (NULL == pSA)
return E_POINTER;
COM_TRY {
// resultierende Größe feststellen
unsigned long ulSize = 0;
CBlobService MakeBlobs;
_ASSERTE(rX.size() == rY.size());
THROW_FAILED_HRESULT(MakeBlobs.GetBlobSize (rClsId, rX.size(), pData, &ulSize));
// SAFEARRAY anlegen und Daten übernehmen
CSafeArray sa (VT_UI1, ulSize);
if (!sa) _com_issue_error(E_OUTOFMEMORY);
{
CSafeArrayLock<BYTE> lock (sa);
THROW_FAILED_HRESULT(MakeBlobs.MakeBlob (rClsId, rX.size(), rX.begin(), rY.begin(), pNormal, pData, lock));
} // lock goes out of scope
*pSA = sa.Detach(); // Ergebnis liefern
} COM_CATCH;
return S_OK;
}
CArcConversion::Convert1(double StartX, double StartY, double OriginX, double OriginY, double EndX, double EndY, SAFEARRAY ** PolyLineArc)
{
//bool KreisbogenApprox (EFlaeche* pBE, double OriginX, double OriginY, double StartX, double StartY, double EndX,
// double EndY, DoublePair dUF, double dPrec)
// TX_ASSERT (pBE != NULL);
// TX_ASSERT (dPrec > 0);
// if (! pBE || dPrec <= 0) return false;
//double m_dAnisotropy = dUF.Y() / dUF.X(); // Relation bei geodätischen Koordinaten
OriginX *= m_dAnisotropy;
StartX *= m_dAnisotropy;
EndX *= m_dAnisotropy;
// dPrec *= dUF.Y();
double dxa = OriginX - EndX;
double dya = OriginY - EndY;
double dxb = StartX - OriginX;
double dyb = StartY - OriginY;
double dxc = EndX - StartX;
double dyc = EndY - StartY;
double cq = dxc*dxc + dyc*dyc; // c² (c ist die Basis des gleichschenkligen Dreiecks)
// a und b müßten eigentlich beide gleich der Pufferbreite sein; durch die Digitalgeometrie ist
// dies aber nicht garantiert, deshalb nachfolgend die Mittelbildung
double dRadius = (sqrt(dxa*dxa + dya*dya) + sqrt(dxb*dxb + dyb*dyb)) / 2.;
//=== Berechnung des Startwinkels der Strecke (OriginX,OriginY)-(StartX,StartY) ==========================
double dSinStartW = (StartY - OriginY) / dRadius;
if (dSinStartW > 1.) dSinStartW = 1.;
if (dSinStartW < -1.) dSinStartW = -1.;
double dStartW = asin (dSinStartW); // Anstiegswinkels der Strecke (OriginX,OriginY)-(StartX,StartY)
// asin liefert Werte zwischen -90° und +90°, deshalb noch
// evtl. Transformationen entspr. der Quadranten
if (StartX < OriginX)
dStartW = Pi - dStartW;
else if (StartY < OriginY)
dStartW = 2.*Pi + dStartW;
//=== Berechnung des Endwinkels der Strecke (OriginX,OriginY)-(EndX,EndY) ==========================
double dSinEndW = (EndY - OriginY) / dRadius;
if (dSinEndW > 1.) dSinEndW = 1.;
if (dSinEndW < -1.) dSinEndW = -1.;
double dEndW = asin (dSinEndW); // Anstiegswinkels der Strecke (OriginX,OriginY)-(StartX,StartY)
// asin liefert Werte zwischen -90° und +90°, deshalb noch
// evtl. Transformationen entspr. der Quadranten
if (EndX < OriginX)
dEndW = Pi - dEndW;
else if (EndY < OriginY)
dEndW = 2.*Pi + dEndW;
//==================================================================
/*double dCosSekW = 1. - cq/(2.*dRadius*dRadius); // Cosinus des Kreisbogenwinkels (mit Kosinussatz)
if (dCosSekW > 1. )
dCosSekW = 1.;
if (dCosSekW < -1. )
dCosSekW = -1.;
double dSekW = acos(dCosSekW); // Sektorwinkel (acos liefert 0 ... Pi)
*/
//=== Sektorenwinkel ===============================================================================
double dSekW = dEndW - dStartW;
if( dSekW < 0 )
dSekW += 2 * Pi;
double dSegW;
double dBogLaenge;
short iKAnz;
//=== Anzahl der Stützpunkte =======================================================================
if( m_dPrecision < 0 ) { // Angabe in Bogenmaß
iKAnz = (short)ceil(dSekW / ((-m_dPrecision / 360) * 2 * Pi));
}
else {
dBogLaenge = dSekW * dRadius; // Länge des Kreisbogens
iKAnz = (short) floor (dBogLaenge/m_dPrecision); // Anzahl der Kanten
}
if( iKAnz != 0 ) // sonst braucht nichts mehr berechnet zu werden
dSegW = dSekW / iKAnz; // Winkel eines n-Eck-Segmentes
CSafeArray sa (VT_R8, (iKAnz + 1) * 2 );
{
CSafeArrayLock<double> data(sa);
double dWi = dStartW;
*((double*)data) = StartX;
*((double*)data+1) = StartY;
*((double*)data+iKAnz*2) = EndX;
*((double*)data+iKAnz*2+1) = EndY;
for (short i = 1; i < iKAnz; i++) {
dWi += dSegW;
*((double*)data+i*2) = (OriginX + dRadius*cos(dWi))/m_dAnisotropy;
*((double*)data+i*2+1) = OriginY + dRadius*sin(dWi);
}
}
*PolyLineArc = sa.Detach();
return S_OK;
}