DBObjRecord *DBVectorIF::Item (DBCoordinate coord) const
{
DBObjRecord *record, *retRecord = (DBObjRecord *) NULL;
switch (DataPTR->Type ())
{
case DBTypeVectorPoint:
{
DBFloat dist, minDist = DBHugeVal;
DBObjTableField *coordFLD = ItemTable->Field (DBrNCoord);
DBMathDistanceFunction distFunc = DBMathGetDistanceFunction (DataPTR);
for (record = FirstItem ();record != (DBObjRecord *) NULL;record = NextItem ())
{
dist = DBMathCoordinateDistance (distFunc,coord,coordFLD->Coordinate (record));
if (dist < minDist) { minDist = dist; retRecord = record; }
}
} break;
case DBTypeVectorLine: break;
case DBTypeVectorPolygon: break;
default: CMmsgPrint (CMmsgAppError, "Unknown Vector Data Type in: %s %d",__FILE__,__LINE__); break;
}
return (retRecord);
}
void DBVLineIF::Vertexes (DBObjRecord *lineRec,DBCoordinate *coord,DBInt vertexNum)
{
DBInt vertex;
DBFloat length = 0;
DBCoordinate *vertexes, prevCoord;
DBObjRecord *dataRec = (DBObjRecord *) NULL;
DBObjectLIST<DBObjRecord> *dataArrays = Data ()->Arrays ();
DBRegion lineExtent, dataExtent = Data ()->Extent ();
DBMathDistanceFunction distFunc = DBMathGetDistanceFunction (Data ());
prevCoord = NodeCoordFLD->Coordinate (FromNodeFLD->Record (lineRec));
lineExtent.Expand (prevCoord);
if (vertexNum > 0)
{
if ((dataRec = dataArrays->Item (lineRec->RowID ())) == (DBObjRecord *) NULL)
{
if ((dataRec = new DBObjRecord (lineRec->Name (),((size_t) vertexNum) * sizeof (DBCoordinate),sizeof (DBFloat))) == (DBObjRecord *) NULL) return;
dataArrays->Add (dataRec);
}
else dataRec->Realloc (vertexNum * sizeof (DBCoordinate));
if ((vertexes = (DBCoordinate *) dataRec->Data ()) == NULL)
{ CMmsgPrint (CMmsgSysError, "Memory Allocation Error in: %s %d",__FILE__,__LINE__); return; }
for (vertex = 0;vertex < vertexNum;vertex++)
{
lineExtent.Expand (vertexes [vertex] = coord [vertex]);
length += DBMathCoordinateDistance (distFunc,prevCoord,coord [vertex]);
prevCoord = coord [vertex];
}
}
VertexesFLD->Record (lineRec,dataRec);
VertexNumFLD->Int (lineRec,vertexNum);
length += DBMathCoordinateDistance (distFunc,prevCoord,NodeCoordFLD->Coordinate (ToNodeFLD->Record (lineRec)));
PerimeterFLD->Float (lineRec,length);
lineExtent.Expand (NodeCoordFLD->Coordinate (ToNodeFLD->Record (lineRec)));
ExtentFLD->Region (lineRec,lineExtent);
dataExtent.Expand (lineExtent);
Data ()->Extent (dataExtent);
}
DBInt DBPointToGrid(DBObjData *pntData, DBObjData *grdData, DBFloat factor) {
DBInt startID, pnt0ID, pntID, id, itemNum;
double dist, minDist, box, box0, bWidth, bHeight;
DBPosition pos;
DBCoordinate gCoord, *pCoord;
DBObjRecord *grdRec, *pntRec, *symRec;
DBObjTable *itemTable, *symTable;
DBObjTableField *valField, *symField;
DBVPointIF *pntIF;
DBGridIF *gridIF;
DBMathDistanceFunction distFunc = DBMathGetDistanceFunction(pntData);
if (distFunc != DBMathGetDistanceFunction(grdData)) {
CMmsgPrint(CMmsgAppError, "Incompatible projections in: %s %d", __FILE__, __LINE__);
return (DBFault);
}
pntIF = new DBVPointIF(pntData);
itemNum = pntIF->ItemNum();
if ((pCoord = (DBCoordinate *) calloc(itemNum, sizeof(DBCoordinate))) == (DBCoordinate *) NULL) {
CMmsgPrint(CMmsgSysError, "Memory allocation Error in: %s %d", __FILE__, __LINE__);
return (DBFault);
}
gridIF = new DBGridIF(grdData);
for (pntID = 0; pntID < itemNum; ++pntID) pCoord[pntID] = pntIF->Coordinate(pntIF->Item(pntID));
if (grdData->Type() == DBTypeGridContinuous)
gridIF->RenameLayer(gridIF->Layer((DBInt) 0), (char *) "Distance to Station");
else {
gridIF->RenameLayer(gridIF->Layer((DBInt) 0), (char *) "Station grid");
itemTable = grdData->Table(DBrNItems);
symTable = grdData->Table(DBrNSymbols);
valField = itemTable->Field(DBrNGridValue);
symField = itemTable->Field(DBrNSymbol);
if ((symRec = symTable->Item(0)) == (DBObjRecord *) NULL)
CMmsgPrint(CMmsgAppError, "Total Metal Gebasz in: %s %d", __FILE__, __LINE__);
for (pntID = 0; pntID < itemNum; ++pntID) {
pntRec = pntIF->Item(pntID);
grdRec = itemTable->Add(pntRec->Name());
valField->Int(grdRec, pntID + 1);
symField->Record(grdRec, symRec);
}
}
startID = 0;
for (pos.Row = 0; pos.Row < gridIF->RowNum(); ++pos.Row) {
DBPause(pos.Row * 100 / gridIF->RowNum());
for (pos.Col = 0; pos.Col < gridIF->ColNum(); ++pos.Col) {
gridIF->Pos2Coord(pos, gCoord);
minDist = box0 = DBHugeVal;
pnt0ID = pntID = startID;
id = DBFault;
do {
bWidth = fabs(gCoord.X - pCoord[pntID].X);
bHeight = fabs(gCoord.Y - pCoord[pntID].Y);
box = bWidth > bHeight ? bWidth : bHeight;
if ((box < box0) && ((dist = DBMathCoordinateDistance(distFunc, gCoord, pCoord[pntID])) < minDist)) {
minDist = dist;
id = startID = pntID;
box *= factor;
if (box0 > box) box0 = box;
}
pntID = pntID + 1 < itemNum ? pntID + 1 : 0;
} while (pntID != pnt0ID);
if (grdData->Type() == DBTypeGridContinuous) gridIF->Value(pos, minDist); else gridIF->Value(pos, id);
}
}
if (grdData->Type() == DBTypeGridContinuous) gridIF->RecalcStats(); else gridIF->DiscreteStats();
delete gridIF;
free(pCoord);
return (DBSuccess);
}
DBMathDistanceFunction DBMathGetDistanceFunction(DBObjData *data) {
return (DBMathGetDistanceFunction(data->Projection()));
}
DBFloat DBMathCoordinateDistance(DBObjData *data, DBCoordinate coord0,
DBCoordinate coord1) { return ((*DBMathGetDistanceFunction(data))(coord0, coord1)); }
DBFloat DBMathCoordinateDistance(DBInt projection, DBCoordinate coord0, DBCoordinate coord1) {
return ((*DBMathGetDistanceFunction(projection))(coord0, coord1));
}
DBInt DBGridIF::Value(DBObjRecord *layerRec, DBCoordinate coord, DBFloat *value) const {
DBInt i, j, pointNum;
DBInt row[9], col[9];
DBPosition pos;
DBCoordinate cellCoord;
DBPosition cellPos;
DBFloat precision, dist, wAvg, sumWeight, retVal;
DBObjRecord *dataRec = LayerFLD->Record(layerRec);
DBFloat missingValue = MissingValueFLD->Float(ItemTable->Item());
DBMathDistanceFunction distFunc = DBMathGetDistanceFunction(DataPTR);
if (DataPTR->Extent().InRegion(coord) == false) return (false);
precision = pow((double) 10.0, (double) DataPTR->Precision());
Coord2Pos(coord, cellPos);
Pos2Coord(cellPos, cellCoord);
if (Flat ||
((fabs(coord.X - cellCoord.X) < precision) &&
(fabs(coord.Y - cellCoord.Y) < precision))) {
j = DimensionVAR.Col * (DimensionVAR.Row - cellPos.Row - 1) + cellPos.Col;
switch (ValueTypeVAR) {
case DBTableFieldFloat:
switch (ValueSizeVAR) {
case sizeof(DBFloat4):
retVal = (DBFloat) ((DBFloat4 *) (dataRec->Data()))[j];
break;
case sizeof(DBFloat):
retVal = (DBFloat) ((DBFloat *) (dataRec->Data()))[j];
break;
}
break;
case DBTableFieldInt:
switch (ValueSizeVAR) {
case sizeof(DBByte):
retVal = (DBFloat) ((DBByte *) (dataRec->Data()))[j];
break;
case sizeof(DBShort):
retVal = (DBFloat) ((DBShort *) (dataRec->Data()))[j];
break;
case sizeof(DBInt):
retVal = (DBFloat) ((DBInt *) (dataRec->Data()))[j];
break;
}
break;
}
if (!CMmathEqualValues(retVal, missingValue)) {
*value = retVal;
return (true);
}
}
col[0] = cellPos.Col;
row[0] = cellPos.Row;
if (coord.X < cellCoord.X) col[0] -= 1;
if (coord.Y < cellCoord.Y) row[0] -= 1;
col[1] = col[0] + 1;
row[1] = row[0];
col[2] = col[0] + 1;
row[2] = row[0] + 1;
col[3] = col[0];
row[3] = row[0] + 1;
pos.Col = col[0];
pos.Row = row[0];
Pos2Coord(pos, cellCoord);
if ((coord.X - cellCoord.X) > (3.0 * CellWidth() / 4.0)) i = 1;
else if ((coord.X - cellCoord.X) > (CellWidth() / 4.0)) i = 0;
else i = -1;
if ((coord.Y - cellCoord.Y) > (3.0 * CellHeight() / 4.0)) j = 1;
else if ((coord.Y - cellCoord.Y) > (CellHeight() / 4.0)) j = 0;
else j = -1;
if ((i != 0) || (j != 0)) {
if (i == 0) {
row[4] = row[5] = j > 0 ? row[2] + 1 : row[0] - 1;
col[4] = col[0];
col[5] = col[2];
pointNum = 6;
}
else if (j == 0) {
row[4] = row[0];
row[5] = row[2];
col[4] = col[5] = i > 0 ? col[2] + 1 : col[0] - 1;
pointNum = 6;
}
else {
row[7] = row[0];
row[8] = row[2];
if (j > 0)
row[4] = row[5] = row[6] = row[2] + 1;
else
row[4] = row[5] = row[6] = row[0] - 1;
if (i > 0) {
col[4] = col[0];
col[5] = col[2];
col[6] = col[7] = col[8] = col[2] + 1;
}
else {
col[5] = col[0];
col[6] = col[2];
col[4] = col[7] = col[8] = col[0] - 1;
}
pointNum = 9;
}
}
else pointNum = 4;
wAvg = sumWeight = 0.0;
for (i = 0; i < pointNum; ++i) {
if (col[i] < 0) continue;
if (row[i] < 0) continue;
if (col[i] >= DimensionVAR.Col) continue;
if (row[i] >= DimensionVAR.Row) continue;
j = DimensionVAR.Col * (DimensionVAR.Row - row[i] - 1) + col[i];
switch (ValueTypeVAR) {
case DBTableFieldFloat:
switch (ValueSizeVAR) {
case sizeof(DBFloat4):
retVal = (DBFloat) ((DBFloat4 *) (dataRec->Data()))[j];
break;
case sizeof(DBFloat):
retVal = (DBFloat) ((DBFloat *) (dataRec->Data()))[j];
break;
}
break;
case DBTableFieldInt:
switch (ValueSizeVAR) {
case sizeof(DBByte):
retVal = (DBFloat) ((DBByte *) (dataRec->Data()))[j];
break;
case sizeof(DBShort):
retVal = (DBFloat) ((DBShort *) (dataRec->Data()))[j];
break;
case sizeof(DBInt):
retVal = (DBFloat) ((DBInt *) (dataRec->Data()))[j];
break;
}
break;
}
if (CMmathEqualValues(retVal, missingValue)) {
if ((col[i] == cellPos.Col) && (row[i] == cellPos.Row)) return (false);
else continue;
}
pos.Row = (DBUShort) row[i];
pos.Col = (DBUShort) col[i];
Pos2Coord(pos, cellCoord);
if (pointNum > 1) {
dist = DBMathCoordinateDistance(distFunc, coord, cellCoord);
dist *= dist;
}
else dist = 1.0;
wAvg = wAvg + retVal / dist;
sumWeight = sumWeight + 1.0 / dist;
}
if (sumWeight > 0) {
*value = wAvg / sumWeight;
return (true);
}
else {
*value = missingValue;
return (false);
}
}
