Skip to content

Commit

Permalink
also added a wkb method to delete vertices
Browse files Browse the repository at this point in the history 
git-svn-id: http://svn.osgeo.org/qgis/trunk@5350 c8812cc2-4d05-0410-92ff-de0c093fc19c
  • Loading branch information
mhugent committed Apr 23, 2006
1 parent 07baabf commit 8577f1d
Show file tree
Hide file tree
Showing 2 changed files with 249 additions and 5 deletions.
249 changes: 248 additions & 1 deletion src/core/qgsgeometry.cpp
View file
Expand Up @@ -524,7 +524,6 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
}

unsigned int wkbType;
double *xPtr, *yPtr;
unsigned char* ptr = mGeometry+1;
memcpy(&wkbType, ptr, sizeof(wkbType));
ptr += sizeof(wkbType);
Expand All @@ -538,6 +537,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr, &x, sizeof(double));
ptr += sizeof(double);
memcpy(ptr, &y, sizeof(double));
mDirtyGeos = true;
return true;
}
else
Expand All @@ -558,6 +558,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr, &x, sizeof(double));
ptr += sizeof(double);
memcpy(ptr, &y, sizeof(double));
mDirtyGeos = true;
return true;
}
case QGis::WKBMultiLineString:
Expand All @@ -575,6 +576,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
ptr += (vertexnr-pointindex)*2*sizeof(double);
memcpy(ptr, &x, sizeof(double));
memcpy(ptr+sizeof(double), &y, sizeof(double));
mDirtyGeos = true;
return true;
}
pointindex += (*nrPoints);
Expand All @@ -599,6 +601,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr+sizeof(double), &y, sizeof(double));
memcpy(ptr+2*sizeof(double)*(*nrPoints-1), &x, sizeof(double));
memcpy(ptr+sizeof(double)+2*sizeof(double)*(*nrPoints-1), &y, sizeof(double));
mDirtyGeos = true;
return true;
}
else if(vertexnr > pointindex && vertexnr < pointindex + (*nrPoints-1))//move only one point
Expand All @@ -607,6 +610,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr, &x, sizeof(double));
ptr += sizeof(double);
memcpy(ptr, &y, sizeof(double));
mDirtyGeos = true;
return true;
}
ptr += 2*sizeof(double)*(*nrPoints);
Expand Down Expand Up @@ -636,6 +640,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr+sizeof(double), &y, sizeof(double));
memcpy(ptr+2*sizeof(double)*(*nrPoints-1), &x, sizeof(double));
memcpy(ptr+sizeof(double)+2*sizeof(double)*(*nrPoints-1), &y, sizeof(double));
mDirtyGeos = true;
return true;
}
else if(vertexnr > pointindex && vertexnr < pointindex + (*nrPoints-1))//move only one point
Expand All @@ -644,6 +649,7 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
memcpy(ptr, &x, sizeof(double));
ptr += sizeof(double);
memcpy(ptr, &y, sizeof(double));
mDirtyGeos = true;
return true;
}
ptr += 2*sizeof(double)*(*nrPoints);
Expand All @@ -655,6 +661,245 @@ bool QgsGeometry::moveVertexAt(double x, double y, QgsGeometryVertexIndex atVert
}
}

bool QgsGeometry::deleteVertexAt(QgsGeometryVertexIndex atVertex)
{
int vertexnr = atVertex.back();
bool success = false;

if(mDirtyWkb)
{
exportGeosToWkb();
}

//create a new geometry buffer for the modified geometry
unsigned char* newbuffer = new unsigned char[mGeometrySize-2*sizeof(double)];
memcpy(newbuffer, mGeometry, 1+sizeof(int)); //endian and type are the same

int pointindex = 0;
unsigned int wkbType;
unsigned char* ptr = mGeometry+1;
memcpy(&wkbType, ptr, sizeof(wkbType));
ptr += sizeof(wkbType);
unsigned char* newBufferPtr = newbuffer+1+sizeof(int);

switch(wkbType)
{
case QGis::WKBPoint:
{
break; //cannot remove the only point vertex
}
case QGis::WKBMultiPoint:
case QGis::WKBLineString:
{
int* nPoints = (int*)ptr;
if((*nPoints) < 3 || vertexnr > (*nPoints)-1 || vertexnr < 0)
{
break;
}
int newNPoints = (*nPoints)-1; //new number of points
memcpy(newBufferPtr, &newNPoints, sizeof(int));
ptr += sizeof(int);
newBufferPtr += sizeof(int);
for(int pointnr = 0; pointnr < *nPoints; ++pointnr)
{
if(vertexnr != pointindex)
{
memcpy(newBufferPtr, ptr, sizeof(double));
memcpy(newBufferPtr+sizeof(double), ptr+sizeof(double), sizeof(double));
newBufferPtr += 2*sizeof(double);
}
else
{
success = true;
}
ptr += 2*sizeof(double);
++pointindex;
}
break;
}
case QGis::WKBMultiLineString:
{
int* nLines = (int*)ptr;
memcpy(newBufferPtr, nLines, sizeof(int));
newBufferPtr += sizeof(int);
ptr += sizeof(int);
int* nPoints = 0; //number of points in a line
int pointindex = 0;
for(int linenr = 0; linenr < *nLines; ++linenr)
{
nPoints = (int*)ptr;
ptr += sizeof(int);
int newNPoint;

//find out if the vertex is in this line
if(vertexnr >= pointindex && vertexnr < pointindex + (*nPoints))
{
if(*nPoints < 3)
{
break;
}
newNPoint = (*nPoints)-1;
}
else
{
newNPoint = *nPoints;
}
memcpy(newBufferPtr, &newNPoint, sizeof(int));

for(int pointnr = 0; pointnr < *nPoints; ++pointnr)
{
if(vertexnr != pointindex)
{
memcpy(newBufferPtr, ptr, sizeof(double));//x
memcpy(newBufferPtr+sizeof(double), ptr+sizeof(double), sizeof(double));//y
newBufferPtr += 2*sizeof(double);
}
else
{
success = true;
}
ptr += 2*sizeof(double);
++pointindex;
}
}
break;
}
case QGis::WKBPolygon:
{
int* nRings = (int*)ptr;
memcpy(newBufferPtr, nRings, sizeof(int));
ptr += sizeof(int);
newBufferPtr += sizeof(int);
int* nPoints = 0; //number of points in a ring
int pointindex = 0;

for(int ringnr = 0; ringnr < *nRings; ++ringnr)
{
nPoints = (int*)ptr;
ptr += sizeof(int);
int newNPoints;
if(vertexnr >= pointindex && vertexnr < pointindex + *nPoints)//vertex to delete is in this ring
{
if(*nPoints < 5) //a ring has at least 3 points
{
break;
}
newNPoints = *nPoints - 1;
}
else
{
newNPoints = *nPoints;
}
memcpy(newBufferPtr, &newNPoints, sizeof(int));
newBufferPtr += sizeof(int);

for(int pointnr = 0; pointnr < *nPoints; ++pointnr)
{
if(vertexnr != pointindex)
{
memcpy(newBufferPtr, ptr, sizeof(double));//x
memcpy(newBufferPtr+sizeof(double), ptr+sizeof(double), sizeof(double));//y
newBufferPtr += 2*sizeof(double);
}
else
{
if(pointnr == 0)//adjust the last point of the ring
{
memcpy(ptr+(*nPoints-1)*2*sizeof(double), ptr+2*sizeof(double), sizeof(double));
memcpy(ptr+sizeof(double)+(*nPoints-1)*2*sizeof(double), ptr+3*sizeof(double), sizeof(double));
}
if(pointnr == (*nPoints)-1)
{
memcpy(newBufferPtr-(*nPoints-2)*2*sizeof(double), ptr-2*sizeof(double), sizeof(double));
memcpy(newBufferPtr-(*nPoints-2)*2*sizeof(double)+sizeof(double), ptr-sizeof(double), sizeof(double));
}
success = true;
}
ptr += 2*sizeof(double);
++pointindex;
}
}
break;
}
case QGis::WKBMultiPolygon:
{
int* nPolys = (int*)ptr;
memcpy(newBufferPtr, nPolys, sizeof(int));
newBufferPtr += sizeof(int);
ptr +=sizeof(int);
int* nRings = 0;
int* nPoints = 0;
int pointindex = 0;

for(int polynr = 0; polynr < *nPolys; ++polynr)
{
nRings = (int*)ptr;
memcpy(newBufferPtr, nRings, sizeof(int));
newBufferPtr += sizeof(int);
ptr += sizeof(int);
for(int ringnr = 0; ringnr < *nRings; ++ringnr)
{
nPoints = (int*)ptr;
ptr += sizeof(int);
int newNPoints;
if(vertexnr >= pointindex && vertexnr < pointindex + *nPoints)//vertex to delete is in this ring
{
if(*nPoints < 5) //a ring has at least 3 points
{
break;
}
newNPoints = *nPoints - 1;
}
else
{
newNPoints = *nPoints;
}
memcpy(newBufferPtr, &newNPoints, sizeof(int));
newBufferPtr += sizeof(int);

for(int pointnr = 0; pointnr < *nPoints; ++pointnr)
{
if(vertexnr != pointindex)
{
memcpy(newBufferPtr, ptr, sizeof(double));//x
memcpy(newBufferPtr+sizeof(double), ptr+sizeof(double), sizeof(double));//y
newBufferPtr += 2*sizeof(double);
}
else
{
if(pointnr == 0)//adjust the last point of the ring
{
memcpy(ptr+(*nPoints-1)*2*sizeof(double), ptr+2*sizeof(double), sizeof(double));
memcpy(ptr+sizeof(double)+(*nPoints-1)*2*sizeof(double), ptr+3*sizeof(double), sizeof(double));
}
if(pointnr == (*nPoints)-1)
{
memcpy(newBufferPtr-(*nPoints-2)*2*sizeof(double), ptr-2*sizeof(double), sizeof(double));
memcpy(newBufferPtr-(*nPoints-2)*2*sizeof(double)+sizeof(double), ptr-sizeof(double), sizeof(double));
}
success = true;
}
ptr += 2*sizeof(double);
++pointindex;
}
}
}
break;
}
}
if(success)
{
delete mGeometry;
mGeometry = newbuffer;
mGeometrySize -= (2*sizeof(double));
mDirtyGeos = true;
}
else
{
delete[] newbuffer;
}
}



#if 0
Expand Down Expand Up @@ -697,6 +942,7 @@ bool QgsGeometry::moveVertexAt(double x, double y,
}
#endif

#if 0
bool QgsGeometry::deleteVertexAt(QgsGeometryVertexIndex atVertex)
{

Expand Down Expand Up @@ -789,6 +1035,7 @@ bool QgsGeometry::deleteVertexAt(QgsGeometryVertexIndex atVertex)
return FALSE;

}
#endif


bool QgsGeometry::vertexAt(double &x, double &y,
Expand Down
5 changes: 1 addition & 4 deletions src/core/qgsgeometry.h
View file
Expand Up @@ -104,10 +104,7 @@ class QgsGeometry {
* ring and item (first number is index 0)
* to the given coordinates.
* Returns FALSE if atVertex does not correspond to a valid vertex
* on this geometry (including if this geometry is a Point).
* It is up to the caller to distinguish between
* these error conditions. (Or maybe we add another method to this
* object to help make the distinction?)
* on this geometry
*/
bool moveVertexAt(double x, double y, QgsGeometryVertexIndex atVertex);

Expand Down

0 comments on commit 8577f1d

Please sign in to comment.