/*!
\brief Get first found feature based on it's properties (AND)
\param column array of property names
\param value array of property values
\param num number of array items
\param bGeom True to check also geometry != NULL
\return pointer to feature definition or NULL on failure (not found)
*/
VFKFeatureSQLite *VFKDataBlockSQLite::GetFeature(const char **column, GUIntBig *value, int num,
bool bGeom)
{
VFKReaderSQLite *poReader = (VFKReaderSQLite*) m_poReader;
CPLString osSQL;
osSQL.Printf("SELECT %s FROM %s WHERE ", FID_COLUMN, m_pszName);
CPLString osItem;
for( int i = 0; i < num; i++ )
{
if (i > 0)
osItem.Printf(" AND %s = " CPL_FRMT_GUIB, column[i], value[i]);
else
osItem.Printf("%s = " CPL_FRMT_GUIB, column[i], value[i]);
osSQL += osItem;
}
if( bGeom )
{
osItem.Printf(" AND %s IS NOT NULL", GEOM_COLUMN);
osSQL += osItem;
}
sqlite3_stmt *hStmt = poReader->PrepareStatement(osSQL.c_str());
if (poReader->ExecuteSQL(hStmt) != OGRERR_NONE)
return NULL;
int idx = sqlite3_column_int(hStmt, 0) - 1; /* rowid starts at 1 */
sqlite3_finalize(hStmt);
if (idx < 0 || idx >= m_nFeatureCount) // ? assert
return NULL;
return (VFKFeatureSQLite *) GetFeatureByIndex(idx);
}
/*!
\brief Get first found feature based on it's properties (AND)
\param column array of property names
\param value array of property values
\param num number of array items
\param bGeom True to check also geometry != NULL
\return pointer to feature definition or NULL on failure (not found)
*/
VFKFeatureDB *VFKDataBlockDB::GetFeature(const char **column, GUIntBig *value, int num,
bool bGeom)
{
int idx;
CPLString osSQL, osItem;
VFKReaderDB *poReader;
poReader = (VFKReaderDB*) m_poReader;
osSQL.Printf("SELECT %s FROM %s WHERE ", FID_COLUMN, m_pszName);
for (int i = 0; i < num; i++) {
if (i > 0)
osItem.Printf(" AND %s = " CPL_FRMT_GUIB, column[i], value[i]);
else
osItem.Printf("%s = " CPL_FRMT_GUIB, column[i], value[i]);
osSQL += osItem;
}
if (bGeom) {
osItem.Printf(" AND %s IS NOT NULL", GEOM_COLUMN);
osSQL += osItem;
}
if (poReader->ExecuteSQL(osSQL.c_str(), idx) != OGRERR_NONE)
return NULL;
idx -= 1; /* rowid starts at 1 */
if (idx < 0 || idx >= m_nFeatureCount) // ? assert
return NULL;
return (VFKFeatureDB *) GetFeatureByIndex(idx);
}
/*!
\brief Load geometry (point layers)
\return number of invalid features
*/
int VFKDataBlockSQLite::LoadGeometryPoint()
{
int nInvalid, rowId, nGeometries;
bool bSkipInvalid;
/* long iFID; */
double x, y;
CPLString osSQL;
sqlite3_stmt *hStmt;
VFKFeatureSQLite *poFeature;
VFKReaderSQLite *poReader;
nInvalid = nGeometries = 0;
poReader = (VFKReaderSQLite*) m_poReader;
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
bSkipInvalid = EQUAL(m_pszName, "OB") || EQUAL(m_pszName, "OP") || EQUAL(m_pszName, "OBBP");
osSQL.Printf("SELECT SOURADNICE_Y,SOURADNICE_X,%s,rowid FROM %s",
FID_COLUMN, m_pszName);
hStmt = poReader->PrepareStatement(osSQL.c_str());
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
while(poReader->ExecuteSQL(hStmt) == OGRERR_NONE) {
/* read values */
x = -1.0 * sqlite3_column_double(hStmt, 0); /* S-JTSK coordinate system expected */
y = -1.0 * sqlite3_column_double(hStmt, 1);
#ifdef DEBUG
const long iFID = sqlite3_column_double(hStmt, 2);
#endif
rowId = sqlite3_column_int(hStmt, 3);
poFeature = (VFKFeatureSQLite *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
/* create geometry */
OGRPoint pt(x, y);
if (!poFeature->SetGeometry(&pt)) {
nInvalid++;
continue;
}
/* store also geometry in DB */
if (poReader->IsSpatial() &&
SaveGeometryToDB(&pt, rowId) != OGRERR_FAILURE)
nGeometries++;
}
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
return bSkipInvalid ? 0 : nInvalid;
}
/*!
\brief Get feature by FID
Modifies next feature id.
\param nFID feature id
\return pointer to feature definition or NULL on failure (not found)
*/
IVFKFeature *VFKDataBlockDB::GetFeature(GIntBig nFID)
{
int rowId;
CPLString osSQL;
VFKReaderDB *poReader;
if (m_nFeatureCount < 0) {
m_poReader->ReadDataRecords(this);
}
if (nFID < 1 || nFID > m_nFeatureCount)
return NULL;
if (m_bGeometryPerBlock && !m_bGeometry) {
LoadGeometry();
}
poReader = (VFKReaderDB*) m_poReader;
osSQL.Printf("SELECT rowid FROM %s WHERE %s = " CPL_FRMT_GIB,
m_pszName, FID_COLUMN, nFID);
if (EQUAL(m_pszName, "SBP")) {
osSQL += " AND PORADOVE_CISLO_BODU = 1";
}
if (poReader->ExecuteSQL(osSQL.c_str(), rowId) != OGRERR_NONE) {
rowId = -1; // TODO: solve
}
return GetFeatureByIndex(rowId - 1);
}
/*!
\brief Get feature by FID
Modifies next feature id.
\param nFID feature id
\return pointer to feature definition
\return NULL on failure (not found)
*/
IVFKFeature *IVFKDataBlock::GetFeature(long nFID)
{
if (m_nFeatureCount < 0) {
m_poReader->ReadDataRecords(this);
}
if (nFID < 1 || nFID > m_nFeatureCount)
return NULL;
if (m_bGeometryPerBlock && !m_bGeometry) {
LoadGeometry();
}
if (m_nGeometryType == wkbPoint ||
m_nGeometryType == wkbPolygon ||
m_nGeometryType == wkbNone) {
return GetFeatureByIndex(int (nFID) - 1); /* zero-based index */
}
else if (m_nGeometryType == wkbLineString) {
/* line string is built from more data records */
for (int i = 0; i < m_nFeatureCount; i++) {
if (m_papoFeature[i]->GetFID() == nFID) {
m_iNextFeature = i + 1;
return m_papoFeature[i];
}
}
}
return NULL;
}
/*!
\brief Get features based on properties
\param column array of property names
\param value array of property values
\param num number of array items
\return list of features
*/
VFKFeatureSQLiteList VFKDataBlockSQLite::GetFeatures(const char **column, GUIntBig *value, int num)
{
int iRowId;
CPLString osSQL, osItem;
VFKReaderSQLite *poReader;
VFKFeatureSQLiteList fList;
sqlite3_stmt *hStmt;
poReader = (VFKReaderSQLite*) m_poReader;
osSQL.Printf("SELECT rowid from %s WHERE ", m_pszName);
for (int i = 0; i < num; i++) {
if (i > 0)
osItem.Printf(" OR %s = " CPL_FRMT_GUIB, column[i], value[i]);
else
osItem.Printf("%s = " CPL_FRMT_GUIB, column[i], value[i]);
osSQL += osItem;
}
osSQL += " ORDER BY ";
osSQL += FID_COLUMN;
hStmt = poReader->PrepareStatement(osSQL.c_str());
while (poReader->ExecuteSQL(hStmt) == OGRERR_NONE) {
iRowId = sqlite3_column_int(hStmt, 0);
fList.push_back((VFKFeatureSQLite *)GetFeatureByIndex(iRowId - 1));
}
return fList;
}
/*!
\brief Get features based on properties
\param column array of property names
\param value array of property values
\param num number of array items
\return list of features
*/
VFKFeatureDBList VFKDataBlockDB::GetFeatures(const char **column, GUIntBig *value, int num)
{
int iRowId;
CPLString osSQL, osItem;
VFKReaderDB *poReader;
VFKFeatureDBList fList;
poReader = (VFKReaderDB*) m_poReader;
osSQL.Printf("SELECT rowid from %s WHERE ", m_pszName);
for (int i = 0; i < num; i++) {
if (i > 0)
osItem.Printf(" OR %s = " CPL_FRMT_GUIB, column[i], value[i]);
else
osItem.Printf("%s = " CPL_FRMT_GUIB, column[i], value[i]);
osSQL += osItem;
}
osSQL += " ORDER BY ";
osSQL += FID_COLUMN;
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_INT));
poReader->PrepareStatement(osSQL.c_str());
while (poReader->ExecuteSQL(record) == OGRERR_NONE) {
iRowId = static_cast<int> (record[0]);
fList.push_back((VFKFeatureDB *)GetFeatureByIndex(iRowId - 1));
}
return fList;
}
/*!
\brief Get feature by FID
Modifies next feature id.
\param nFID feature id
\return pointer to feature definition or NULL on failure (not found)
*/
IVFKFeature *VFKDataBlockSQLite::GetFeature(GIntBig nFID)
{
if (m_nFeatureCount < 0) {
m_poReader->ReadDataRecords(this);
}
if (nFID < 1 || nFID > m_nFeatureCount)
return NULL;
if( m_bGeometryPerBlock && !m_bGeometry )
{
LoadGeometry();
}
VFKReaderSQLite *poReader = (VFKReaderSQLite*) m_poReader;
CPLString osSQL;
osSQL.Printf("SELECT rowid FROM %s WHERE %s = " CPL_FRMT_GIB,
m_pszName, FID_COLUMN, nFID);
if (EQUAL(m_pszName, "SBP")) {
osSQL += " AND PORADOVE_CISLO_BODU = 1";
}
sqlite3_stmt *hStmt = poReader->PrepareStatement(osSQL.c_str());
int rowId = -1;
if (poReader->ExecuteSQL(hStmt) == OGRERR_NONE) {
rowId = sqlite3_column_int(hStmt, 0);
}
sqlite3_finalize(hStmt);
return GetFeatureByIndex(rowId - 1);
}
/*!
\brief Load geometry (point layers)
\return number of invalid features
*/
int VFKDataBlock::LoadGeometryPoint()
{
/* -> wkbPoint */
long nInvalid;
double x, y;
int i_idxX, i_idxY;
VFKFeature *poFeature;
nInvalid = 0;
i_idxY = GetPropertyIndex("SOURADNICE_Y");
i_idxX = GetPropertyIndex("SOURADNICE_X");
if (i_idxY < 0 || i_idxX < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
for (int j = 0; j < ((IVFKDataBlock *) this)->GetFeatureCount(); j++) {
poFeature = (VFKFeature *) GetFeatureByIndex(j);
x = -1.0 * poFeature->GetProperty(i_idxY)->GetValueD();
y = -1.0 * poFeature->GetProperty(i_idxX)->GetValueD();
OGRPoint pt(x, y);
if (!poFeature->SetGeometry(&pt))
nInvalid++;
}
return nInvalid;
}
/*!
\brief Get first found feature based on it's properties
Note: modifies next feature.
\param idx property index
\param value property value
\param poList list of features (NULL to loop all features)
\return pointer to feature definition
\return NULL on failure (not found)
*/
VFKFeature *VFKDataBlock::GetFeature(int idx, GUIntBig value, VFKFeatureList *poList)
{
GUIntBig iPropertyValue;
VFKFeature *poVfkFeature;
if (poList) {
for (VFKFeatureList::iterator i = poList->begin(), e = poList->end();
i != e; ++i) {
poVfkFeature = *i;
iPropertyValue = strtoul(poVfkFeature->GetProperty(idx)->GetValueS(), NULL, 0);
if (iPropertyValue == value) {
poList->erase(i); /* ??? */
return poVfkFeature;
}
}
}
else {
for (int i = 0; i < m_nFeatureCount; i++) {
poVfkFeature = (VFKFeature *) GetFeatureByIndex(i);
iPropertyValue = strtoul(poVfkFeature->GetProperty(idx)->GetValueS(), NULL, 0);
if (iPropertyValue == value) {
m_iNextFeature = i + 1;
return poVfkFeature;
}
}
}
return NULL;
}
/*!
\brief Get first found feature based on it's property
\param column property name
\param value property value
\param bGeom True to check also geometry != NULL
\return pointer to feature definition or NULL on failure (not found)
*/
VFKFeatureDB *VFKDataBlockDB::GetFeature(const char *column, GUIntBig value,
bool bGeom)
{
int idx;
CPLString osSQL;
VFKReaderDB *poReader;
poReader = (VFKReaderDB*) m_poReader;
osSQL.Printf("SELECT %s from %s WHERE %s = " CPL_FRMT_GUIB,
FID_COLUMN, m_pszName, column, value);
if (bGeom) {
CPLString osColumn;
osColumn.Printf(" AND %s IS NOT NULL", GEOM_COLUMN);
osSQL += osColumn;
}
if (poReader->ExecuteSQL(osSQL.c_str(), idx) != OGRERR_NONE)
return NULL;
idx -= 1;
if (idx < 0 || idx >= m_nFeatureCount) // ? assert
return NULL;
return (VFKFeatureDB *) GetFeatureByIndex(idx);
}
/*!
\brief Get first found feature based on it's property
\param column property name
\param value property value
\param bGeom True to check also geometry != NULL
\return pointer to feature definition or NULL on failure (not found)
*/
VFKFeatureSQLite *VFKDataBlockSQLite::GetFeature(const char *column, GUIntBig value,
bool bGeom)
{
int idx;
CPLString osSQL;
VFKReaderSQLite *poReader;
sqlite3_stmt *hStmt;
poReader = (VFKReaderSQLite*) m_poReader;
osSQL.Printf("SELECT %s from %s WHERE %s = " CPL_FRMT_GUIB,
FID_COLUMN, m_pszName, column, value);
if (bGeom) {
CPLString osColumn;
osColumn.Printf(" AND %s IS NOT NULL", GEOM_COLUMN);
osSQL += osColumn;
}
hStmt = poReader->PrepareStatement(osSQL.c_str());
if (poReader->ExecuteSQL(hStmt) != OGRERR_NONE)
return NULL;
idx = sqlite3_column_int(hStmt, 0) - 1;
sqlite3_finalize(hStmt);
if (idx < 0 || idx >= m_nFeatureCount) // ? assert
return NULL;
return (VFKFeatureSQLite *) GetFeatureByIndex(idx);
}
/*!
\brief Get features based on properties
\param idx property index
\param value property value
\return list of features
*/
VFKFeatureList VFKDataBlock::GetFeatures(int idx, GUIntBig value)
{
GUIntBig iPropertyValue;
VFKFeature *poVfkFeature;
std::vector<VFKFeature *> poResult;
for (int i = 0; i < m_nFeatureCount; i++) {
poVfkFeature = (VFKFeature *) GetFeatureByIndex(i);
iPropertyValue = strtoul(poVfkFeature->GetProperty(idx)->GetValueS(), NULL, 0);
if (iPropertyValue == value) {
poResult.push_back(poVfkFeature);
}
}
return poResult;
}
/*!
\brief Load geometry (linestring HP/DPM layer)
\return number of invalid features
*/
int VFKDataBlock::LoadGeometryLineStringHP()
{
long nInvalid;
int idxId, idxMy_Id, idxPCB;
GUIntBig id;
VFKDataBlock *poDataBlockLines;
VFKFeature *poFeature, *poLine;
VFKFeatureList poLineList;
nInvalid = 0;
poDataBlockLines = (VFKDataBlock *) m_poReader->GetDataBlock("SBP");
if (NULL == poDataBlockLines) {
CPLError(CE_Failure, CPLE_NotSupported,
"Data block %s not found.\n", m_pszName);
return nInvalid;
}
poDataBlockLines->LoadGeometry();
idxId = GetPropertyIndex("ID");
if (EQUAL (m_pszName, "HP"))
idxMy_Id = poDataBlockLines->GetPropertyIndex("HP_ID");
else
idxMy_Id = poDataBlockLines->GetPropertyIndex("DPM_ID");
idxPCB = poDataBlockLines->GetPropertyIndex("PORADOVE_CISLO_BODU");
if (idxId < 0 || idxMy_Id < 0 || idxPCB < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
poLineList = poDataBlockLines->GetFeatures(idxPCB, 1); /* reduce to first segment */
for (int i = 0; i < ((IVFKDataBlock *) this)->GetFeatureCount(); i++) {
poFeature = (VFKFeature *) GetFeatureByIndex(i);
id = strtoul(poFeature->GetProperty(idxId)->GetValueS(), NULL, 0);
poLine = poDataBlockLines->GetFeature(idxMy_Id, id, &poLineList);
if (!poLine || !poLine->GetGeometry())
continue;
if (!poFeature->SetGeometry(poLine->GetGeometry()))
nInvalid++;
}
poDataBlockLines->ResetReading();
return nInvalid;
}
/*!
\brief Load geometry (point layers)
\return number of invalid features
*/
int VFKDataBlockDB::LoadGeometryPoint()
{
int nInvalid, rowId, nGeometries;
bool bSkipInvalid;
/* long iFID; */
double x, y;
CPLString osSQL;
VFKFeatureDB *poFeature;
VFKReaderDB *poReader;
nInvalid = nGeometries = 0;
poReader = (VFKReaderDB*) m_poReader;
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
bSkipInvalid = EQUAL(m_pszName, "OB") || EQUAL(m_pszName, "OP") || EQUAL(m_pszName, "OBBP");
osSQL.Printf("SELECT SOURADNICE_Y,SOURADNICE_X,%s,rowid FROM %s",
FID_COLUMN, m_pszName);
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_DOUBLE));
record.push_back(VFKDbValue(DT_DOUBLE));
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_INT));
poReader->PrepareStatement(osSQL.c_str());
while(poReader->ExecuteSQL(record) == OGRERR_NONE && record.size() > 0) {
/* read values */
x = -1.0 * (double) record[0]; /* S-JTSK coordinate system expected */
y = -1.0 * (double) record[1];
#ifdef DEBUG
const GIntBig iFID = (GIntBig) record[2];
#endif
rowId = (int) record[3];
poFeature = (VFKFeatureDB *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
/* create geometry */
OGRPoint pt(x, y);
if (!poFeature->SetGeometry(&pt)) {
nInvalid++;
continue;
}
/* store also geometry in DB */
if (poReader->IsSpatial() &&
SaveGeometryToDB(&pt, rowId) != OGRERR_FAILURE)
nGeometries++;
}
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
return bSkipInvalid ? 0 : nInvalid;
}
/*!
\brief Load geometry (polygon BUD/PAR layers)
\return number of invalid features
*/
int VFKDataBlockDB::LoadGeometryPolygon()
{
int nInvalidNoLines, nInvalidNoRings, nGeometries, nBridges;
int rowId, nCount, nCountMax;
size_t nLines;
GIntBig iFID;
bool bIsPar, bNewRing, bFound;
CPLString osSQL;
const char *vrColumn[2];
GUIntBig vrValue[2];
GUIntBig id, idOb;
sqlite3_stmt *hStmt;
VFKReaderDB *poReader;
VFKDataBlockDB *poDataBlockLines1, *poDataBlockLines2;
VFKFeatureDB *poFeature;
VFKFeatureDBList poLineList;
/* first is to be considered as exterior */
PointListArray poRingList;
std::vector<OGRLinearRing *> poLinearRingList;
OGRPolygon ogrPolygon;
OGRLinearRing *poOgrRing;
nInvalidNoLines = nInvalidNoRings = nGeometries = 0;
poReader = (VFKReaderDB*) m_poReader;
if (EQUAL (m_pszName, "PAR")) {
poDataBlockLines1 = (VFKDataBlockDB *) m_poReader->GetDataBlock("HP");
poDataBlockLines2 = poDataBlockLines1;
bIsPar = TRUE;
}
else {
poDataBlockLines1 = (VFKDataBlockDB *) m_poReader->GetDataBlock("OB");
poDataBlockLines2 = (VFKDataBlockDB *) m_poReader->GetDataBlock("SBP");
bIsPar = FALSE;
}
if (NULL == poDataBlockLines1) {
CPLError(CE_Warning, CPLE_FileIO,
"Data block %s not found. Unable to build geometry for %s.",
bIsPar ? "HP" : "OB", m_pszName);
return -1;
}
if (NULL == poDataBlockLines2) {
CPLError(CE_Warning, CPLE_FileIO,
"Data block %s not found. Unable to build geometry for %s.",
"SBP", m_pszName);
return -1;
}
poDataBlockLines1->LoadGeometry();
poDataBlockLines2->LoadGeometry();
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
if (bIsPar) {
vrColumn[0] = "PAR_ID_1";
vrColumn[1] = "PAR_ID_2";
}
else {
vrColumn[0] = "OB_ID";
vrColumn[1] = "PORADOVE_CISLO_BODU";
vrValue[1] = 1;
}
osSQL.Printf("SELECT ID,%s,rowid FROM %s", FID_COLUMN, m_pszName);
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_INT));
poReader->PrepareStatement(osSQL.c_str());
while(poReader->ExecuteSQL(record) == OGRERR_NONE) {
nBridges = 0;
/* read values */
id = static_cast<GIntBig> (record[0]);
iFID = static_cast<GIntBig> (record[1]);
rowId = static_cast<int> (record[2]);
poFeature = (VFKFeatureDB *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
if (bIsPar) {
vrValue[0] = vrValue[1] = id;
poLineList = poDataBlockLines1->GetFeatures(vrColumn, vrValue, 2);
}
else {
VFKFeatureDB *poLineSbp;
std::vector<VFKFeatureDB *> poLineListOb;
sqlite3_stmt *hStmtOb;
osSQL.Printf("SELECT ID FROM %s WHERE BUD_ID = " CPL_FRMT_GUIB,
poDataBlockLines1->GetName(), id);
if (poReader->IsSpatial()) {
CPLString osColumn;
osColumn.Printf(" AND %s IS NULL", GEOM_COLUMN);
osSQL += osColumn;
}
std::vector<VFKDbValue> record1;
record1.push_back(VFKDbValue(DT_BIGINT));
poReader->PrepareStatement(osSQL.c_str(), 1);
while(poReader->ExecuteSQL(record, 1) == OGRERR_NONE) {
idOb = static_cast<GIntBig> (record[0]);
vrValue[0] = idOb;
poLineSbp = poDataBlockLines2->GetFeature(vrColumn, vrValue, 2);
if (poLineSbp)
poLineList.push_back(poLineSbp);
}
}
nLines = poLineList.size();
if (nLines < 1) {
CPLDebug("OGR-VFK",
"%s: unable to collect rings for polygon fid = %ld (no lines)",
m_pszName, iFID);
nInvalidNoLines++;
continue;
}
/* clear */
ogrPolygon.empty();
poRingList.clear();
/* collect rings from lines */
bFound = FALSE;
nCount = 0;
nCountMax = static_cast<int>(nLines) * 2;
while (poLineList.size() > 0 && nCount < nCountMax) {
bNewRing = !bFound ? TRUE : FALSE;
bFound = FALSE;
int i = 1;
for (VFKFeatureDBList::iterator iHp = poLineList.begin(), eHp = poLineList.end();
iHp != eHp; ++iHp, ++i) {
const OGRLineString *pLine = (OGRLineString *) (*iHp)->GetGeometry();
if (pLine && AppendLineToRing(&poRingList, pLine, bNewRing)) {
bFound = TRUE;
poLineList.erase(iHp);
break;
}
}
nCount++;
}
CPLDebug("OGR-VFK", "%s: fid = %ld nlines = %d -> nrings = %d", m_pszName,
iFID, (int)nLines, (int)poRingList.size());
if (poLineList.size() > 0) {
CPLDebug("OGR-VFK",
"%s: unable to collect rings for polygon fid = %ld",
m_pszName, iFID);
nInvalidNoRings++;
continue;
}
/* build rings */
poLinearRingList.clear();
int i = 1;
for (PointListArray::const_iterator iRing = poRingList.begin(), eRing = poRingList.end();
iRing != eRing; ++iRing) {
OGRPoint *poPoint;
PointList *poList = *iRing;
poLinearRingList.push_back(new OGRLinearRing());
poOgrRing = poLinearRingList.back();
CPLAssert(NULL != poOgrRing);
for (PointList::iterator iPoint = poList->begin(), ePoint = poList->end();
iPoint != ePoint; ++iPoint) {
poPoint = &(*iPoint);
poOgrRing->addPoint(poPoint);
}
i++;
}
/* find exterior ring */
if (poLinearRingList.size() > 1) {
double dArea, dMaxArea;
std::vector<OGRLinearRing *>::iterator exteriorRing;
exteriorRing = poLinearRingList.begin();
dMaxArea = -1.;
for (std::vector<OGRLinearRing *>::iterator iRing = poLinearRingList.begin(),
eRing = poLinearRingList.end(); iRing != eRing; ++iRing) {
poOgrRing = *iRing;
if (!IsRingClosed(poOgrRing))
continue; /* skip unclosed rings */
dArea = poOgrRing->get_Area();
if (dArea > dMaxArea) {
dMaxArea = dArea;
exteriorRing = iRing;
}
}
if (exteriorRing != poLinearRingList.begin()) {
std::swap(*poLinearRingList.begin(), *exteriorRing);
}
}
/* build polygon from rings */
for (std::vector<OGRLinearRing *>::iterator iRing = poLinearRingList.begin(),
eRing = poLinearRingList.end(); iRing != eRing; ++iRing) {
poOgrRing = *iRing;
/* check if ring is closed */
if (IsRingClosed(poOgrRing)) {
ogrPolygon.addRing(poOgrRing);
}
else {
if (poOgrRing->getNumPoints() == 2) {
CPLDebug("OGR-VFK", "%s: Polygon (fid = %ld) bridge removed",
m_pszName, iFID);
nBridges++;
}
else {
CPLDebug("OGR-VFK",
"%s: Polygon (fid = %ld) unclosed ring skipped",
m_pszName, iFID);
}
}
delete poOgrRing;
*iRing = NULL;
}
/* set polygon */
ogrPolygon.setCoordinateDimension(2); /* force 2D */
if (ogrPolygon.getNumInteriorRings() + nBridges != (int) poLinearRingList.size() - 1 ||
!poFeature->SetGeometry(&ogrPolygon)) {
nInvalidNoRings++;
continue;
}
/* store also geometry in DB */
if (poReader->IsSpatial() &&
SaveGeometryToDB(&ogrPolygon, rowId) != OGRERR_FAILURE)
nGeometries++;
}
/* free ring list */
for (PointListArray::iterator iRing = poRingList.begin(), eRing = poRingList.end();
iRing != eRing; ++iRing) {
delete (*iRing);
*iRing = NULL;
}
CPLDebug("OGR-VFK", "%s: nolines = %d norings = %d",
m_pszName, nInvalidNoLines, nInvalidNoRings);
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
return nInvalidNoLines + nInvalidNoRings;
}
/*!
\brief Load geometry (linestring HP/DPM layer)
\return number of invalid features
*/
int VFKDataBlockDB::LoadGeometryLineStringHP()
{
int nInvalid, nGeometries;
int rowId;
GIntBig iFID;
CPLString osColumn, osSQL;
const char *vrColumn[2];
GUIntBig vrValue[2];
sqlite3_stmt *hStmt;
OGRGeometry *poOgrGeometry;
VFKReaderDB *poReader;
VFKDataBlockDB *poDataBlockLines;
VFKFeatureDB *poFeature, *poLine;
nInvalid = nGeometries = 0;
poReader = (VFKReaderDB*) m_poReader;
poDataBlockLines = (VFKDataBlockDB *) m_poReader->GetDataBlock("SBP");
if (NULL == poDataBlockLines) {
CPLError(CE_Failure, CPLE_FileIO,
"Data block %s not found.", m_pszName);
return nInvalid;
}
poDataBlockLines->LoadGeometry();
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
osColumn.Printf("%s_ID", m_pszName);
vrColumn[0] = osColumn.c_str();
vrColumn[1] = "PORADOVE_CISLO_BODU";
vrValue[1] = 1; /* reduce to first segment */
osSQL.Printf("SELECT ID,%s,rowid FROM %s", FID_COLUMN, m_pszName);
/* TODO: handle points in DPM */
if (EQUAL(m_pszName, "DPM"))
osSQL += " WHERE SOURADNICE_X IS NULL";
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_INT));
poReader->PrepareStatement(osSQL.c_str());
while(poReader->ExecuteSQL(record) == OGRERR_NONE) {
/* read values */
vrValue[0] = static_cast<GIntBig> (record[0]);
iFID = static_cast<GIntBig> (record[1]);
rowId = static_cast<int> (record[2]);
poFeature = (VFKFeatureDB *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
poLine = poDataBlockLines->GetFeature(vrColumn, vrValue, 2, TRUE);
if (!poLine) {
poOgrGeometry = NULL;
}
else {
poOgrGeometry = poLine->GetGeometry();
}
if (!poOgrGeometry || !poFeature->SetGeometry(poOgrGeometry)) {
CPLDebug("OGR-VFK", "VFKDataBlockDB::LoadGeometryLineStringHP(): name=%s fid=%ld "
"id=" CPL_FRMT_GUIB " -> %s geometry", m_pszName, iFID, vrValue[0],
poOgrGeometry ? "invalid" : "empty");
nInvalid++;
continue;
}
/* store also geometry in DB */
if (poReader->IsSpatial() &&
SaveGeometryToDB(poOgrGeometry, rowId) != OGRERR_FAILURE &&
poOgrGeometry)
nGeometries++;
}
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
return nInvalid;
}
/*!
\brief Load geometry (linestring SBP layer)
\return number of invalid features
*/
int VFKDataBlockDB::LoadGeometryLineStringSBP()
{
int nInvalid, nGeometries, rowId, iIdx;
CPLString szFType, szFTypeLine;
GUIntBig id, ipcb;
bool bValid;
std::vector<int> rowIdFeat;
CPLString osSQL;
VFKReaderDB *poReader;
VFKDataBlockDB *poDataBlockPoints;
VFKFeatureDB *poFeature, *poPoint, *poLine;
OGRLineString oOGRLine;
nInvalid = nGeometries = 0;
poReader = (VFKReaderDB*) m_poReader;
poLine = NULL;
poDataBlockPoints = (VFKDataBlockDB *) m_poReader->GetDataBlock("SOBR");
if (NULL == poDataBlockPoints) {
CPLError(CE_Failure, CPLE_FileIO,
"Data block %s not found.\n", m_pszName);
return nInvalid;
}
poDataBlockPoints->LoadGeometry();
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
osSQL.Printf("UPDATE %s SET %s = -1", m_pszName, FID_COLUMN);
poReader->ExecuteSQL(osSQL.c_str());
bValid = TRUE;
iIdx = 0;
for (int i = 0; i < 2; i++) {
/* first collect linestrings related to HP, OB or DPM
then collect rest of linestrings */
if (i == 0)
osSQL.Printf("SELECT BP_ID,PORADOVE_CISLO_BODU,PARAMETRY_SPOJENI,_rowid_ FROM '%s' WHERE "
"HP_ID IS NOT NULL OR OB_ID IS NOT NULL OR DPM_ID IS NOT NULL "
"ORDER BY HP_ID,OB_ID,DPM_ID,PORADOVE_CISLO_BODU", m_pszName);
else
osSQL.Printf("SELECT BP_ID,PORADOVE_CISLO_BODU,PARAMETRY_SPOJENI,_rowid_ FROM '%s' WHERE "
"OB_ID IS NULL AND HP_ID IS NULL AND DPM_ID IS NULL "
"ORDER BY ID,PORADOVE_CISLO_BODU", m_pszName);
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_BIGINT));
record.push_back(VFKDbValue(DT_TEXT));
record.push_back(VFKDbValue(DT_INT));
poReader->PrepareStatement(osSQL.c_str());
while(poReader->ExecuteSQL(record) == OGRERR_NONE) {
// read values
id = (GIntBig) record[0];
ipcb = (GIntBig) record[1];
szFType = (CPLString) record[2];
rowId = (int) record[3];
if (ipcb == 1) {
poFeature = (VFKFeatureDB *) GetFeatureByIndex(iIdx);
if( poFeature == NULL )
{
CPLError(CE_Failure, CPLE_AppDefined, "Cannot retrieve feature %d", iIdx);
break;
}
poFeature->SetRowId(rowId);
/* set geometry & reset */
if (poLine && !SetGeometryLineString(poLine, &oOGRLine,
bValid, szFTypeLine, rowIdFeat, nGeometries)) {
nInvalid++;
}
bValid = TRUE;
poLine = poFeature;
szFTypeLine = szFType;
iIdx++;
}
poPoint = (VFKFeatureDB *) poDataBlockPoints->GetFeature("ID", id);
if (poPoint) {
OGRPoint *pt = (OGRPoint *) poPoint->GetGeometry();
if (pt) {
oOGRLine.addPoint(pt);
}
else {
CPLDebug("OGR-VFK",
"Geometry (point ID = " CPL_FRMT_GUIB ") not valid", id);
bValid = FALSE;
}
}
else {
CPLDebug("OGR-VFK",
"Point ID = " CPL_FRMT_GUIB " not found (rowid = %d)",
id, rowId);
bValid = FALSE;
}
/* add vertex to the linestring */
rowIdFeat.push_back(rowId);
}
/* add last line */
if (poLine && !SetGeometryLineString(poLine, &oOGRLine,
bValid, szFType.c_str(), rowIdFeat, nGeometries)) {
nInvalid++;
}
poLine = NULL;
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
}
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
return nInvalid;
}
/*!
\brief Load geometry (linestring SBP layer)
\return number of invalid features
*/
int VFKDataBlock::LoadGeometryLineStringSBP()
{
int idxId, idxBp_Id, idxPCB;
GUIntBig id, ipcb;
int nInvalid;
VFKDataBlock *poDataBlockPoints;
VFKFeature *poFeature, *poPoint, *poLine;
OGRLineString oOGRLine;
nInvalid = 0;
poLine = NULL;
poDataBlockPoints = (VFKDataBlock *) m_poReader->GetDataBlock("SOBR");
if (NULL == poDataBlockPoints) {
CPLError(CE_Failure, CPLE_NotSupported,
"Data block %s not found.\n", m_pszName);
return nInvalid;
}
poDataBlockPoints->LoadGeometry();
idxId = poDataBlockPoints->GetPropertyIndex("ID");
idxBp_Id = GetPropertyIndex("BP_ID");
idxPCB = GetPropertyIndex("PORADOVE_CISLO_BODU");
if (idxId < 0 || idxBp_Id < 0 || idxPCB < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
for (int j = 0; j < ((IVFKDataBlock *) this)->GetFeatureCount(); j++) {
poFeature = (VFKFeature *) GetFeatureByIndex(j);
poFeature->SetGeometry(NULL);
id = strtoul(poFeature->GetProperty(idxBp_Id)->GetValueS(), NULL, 0);
ipcb = strtoul(poFeature->GetProperty(idxPCB)->GetValueS(), NULL, 0);
if (ipcb == 1) {
if (!oOGRLine.IsEmpty()) {
oOGRLine.setCoordinateDimension(2); /* force 2D */
if (!poLine->SetGeometry(&oOGRLine))
nInvalid++;
oOGRLine.empty(); /* restore line */
}
poLine = poFeature;
}
else {
poFeature->SetGeometryType(wkbUnknown);
}
poPoint = poDataBlockPoints->GetFeature(idxId, id);
if (!poPoint)
continue;
OGRPoint *pt = (OGRPoint *) poPoint->GetGeometry();
oOGRLine.addPoint(pt);
}
/* add last line */
oOGRLine.setCoordinateDimension(2); /* force 2D */
if (poLine) {
if (!poLine->SetGeometry(&oOGRLine))
nInvalid++;
}
poDataBlockPoints->ResetReading();
return nInvalid;
}
/*!
\brief Load geometry (polygon BUD/PAR layers)
\return number of invalid features
*/
int VFKDataBlock::LoadGeometryPolygon()
{
long nInvalid;
bool bIsPar, bNewRing, bFound;
GUIntBig id, idOb;
int nCount, nCountMax;
int idxId, idxPar1, idxPar2, idxBud, idxOb, idxIdOb;
VFKFeature *poFeature;
VFKDataBlock *poDataBlockLines1, *poDataBlockLines2;
VFKFeatureList poLineList;
PointListArray poRingList; /* first is to be considered as exterior */
OGRLinearRing ogrRing;
OGRPolygon ogrPolygon;
idxPar1 = idxPar2 = idxBud = idxOb = idxIdOb = 0;
nInvalid = 0;
if (EQUAL (m_pszName, "PAR")) {
poDataBlockLines1 = (VFKDataBlock *) m_poReader->GetDataBlock("HP");
poDataBlockLines2 = poDataBlockLines1;
bIsPar = TRUE;
}
else {
poDataBlockLines1 = (VFKDataBlock *) m_poReader->GetDataBlock("OB");
poDataBlockLines2 = (VFKDataBlock *) m_poReader->GetDataBlock("SBP");
bIsPar = FALSE;
}
if (NULL == poDataBlockLines1 || NULL == poDataBlockLines2) {
CPLError(CE_Failure, CPLE_NotSupported,
"Data block %s not found.\n", m_pszName);
return nInvalid;
}
poDataBlockLines1->LoadGeometry();
poDataBlockLines2->LoadGeometry();
idxId = GetPropertyIndex("ID");
if (idxId < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
if (bIsPar) {
idxPar1 = poDataBlockLines1->GetPropertyIndex("PAR_ID_1");
idxPar2 = poDataBlockLines1->GetPropertyIndex("PAR_ID_2");
if (idxPar1 < 0 || idxPar2 < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
}
else { /* BUD */
idxIdOb = poDataBlockLines1->GetPropertyIndex("ID");
idxBud = poDataBlockLines1->GetPropertyIndex("BUD_ID");
idxOb = poDataBlockLines2->GetPropertyIndex("OB_ID");
if (idxIdOb < 0 || idxBud < 0 || idxOb < 0) {
CPLError(CE_Failure, CPLE_NotSupported,
"Corrupted data (%s).\n", m_pszName);
return nInvalid;
}
}
for (int i = 0; i < ((IVFKDataBlock *) this)->GetFeatureCount(); i++) {
poFeature = (VFKFeature *) GetFeatureByIndex(i);
id = strtoul(poFeature->GetProperty(idxId)->GetValueS(), NULL, 0);
if (bIsPar) {
poLineList = poDataBlockLines1->GetFeatures(idxPar1, idxPar2, id);
}
else {
VFKFeature *poLineOb, *poLineSbp;
std::vector<VFKFeature *> poLineListOb;
poLineListOb = poDataBlockLines1->GetFeatures(idxBud, id);
for (std::vector<VFKFeature *>::const_iterator iOb = poLineListOb.begin(), eOb = poLineListOb.end();
iOb != eOb; ++iOb) {
poLineOb = (*iOb);
idOb = strtoul(poLineOb->GetProperty(idxIdOb)->GetValueS(), NULL, 0);
poLineSbp = poDataBlockLines2->GetFeature(idxOb, idOb);
if (poLineSbp)
poLineList.push_back(poLineSbp);
}
}
if (poLineList.size() < 1)
continue;
/* clear */
ogrPolygon.empty();
poRingList.clear();
/* collect rings (points) */
bFound = FALSE;
nCount = 0;
nCountMax = poLineList.size() * 2;
while (poLineList.size() > 0 && nCount < nCountMax) {
bNewRing = !bFound ? TRUE : FALSE;
bFound = FALSE;
for (VFKFeatureList::iterator iHp = poLineList.begin(), eHp = poLineList.end();
iHp != eHp; ++iHp) {
const OGRLineString *pLine = (OGRLineString *) (*iHp)->GetGeometry();
if (pLine && AppendLineToRing(&poRingList, pLine, bNewRing)) {
bFound = TRUE;
poLineList.erase(iHp);
break;
}
}
nCount++;
}
/* create rings */
for (PointListArray::const_iterator iRing = poRingList.begin(), eRing = poRingList.end();
iRing != eRing; ++iRing) {
PointList *poList = *iRing;
ogrRing.empty();
for (PointList::iterator iPoint = poList->begin(), ePoint = poList->end();
iPoint != ePoint; ++iPoint) {
ogrRing.addPoint(&(*iPoint));
}
ogrPolygon.addRing(&ogrRing);
}
/* set polygon */
ogrPolygon.setCoordinateDimension(2); /* force 2D */
if (!poFeature->SetGeometry(&ogrPolygon))
nInvalid++;
}
/* free ring list */
for (PointListArray::iterator iRing = poRingList.begin(),
eRing = poRingList.end(); iRing != eRing; ++iRing) {
delete (*iRing);
*iRing = NULL;
}
poDataBlockLines1->ResetReading();
poDataBlockLines2->ResetReading();
return nInvalid;
}
/*!
\brief Load geometry from DB
\return true if geometry successfully loaded otherwise false
*/
bool VFKDataBlockSQLite::LoadGeometryFromDB()
{
VFKReaderSQLite *poReader = (VFKReaderSQLite*) m_poReader;
if (!poReader->IsSpatial()) /* check if DB is spatial */
return false;
CPLString osSQL;
osSQL.Printf("SELECT num_geometries FROM %s WHERE table_name = '%s'",
VFK_DB_TABLE, m_pszName);
sqlite3_stmt *hStmt = poReader->PrepareStatement(osSQL.c_str());
if (poReader->ExecuteSQL(hStmt) != OGRERR_NONE)
return false;
const int nGeometries = sqlite3_column_int(hStmt, 0);
sqlite3_finalize(hStmt);
if( nGeometries < 1 )
return false;
const bool bSkipInvalid =
EQUAL(m_pszName, "OB") ||
EQUAL(m_pszName, "OP") ||
EQUAL(m_pszName, "OBBP");
/* load geometry from DB */
osSQL.Printf("SELECT %s,rowid,%s FROM %s ",
GEOM_COLUMN, FID_COLUMN, m_pszName);
if (EQUAL(m_pszName, "SBP"))
osSQL += "WHERE PORADOVE_CISLO_BODU = 1 ";
osSQL += "ORDER BY ";
osSQL += FID_COLUMN;
hStmt = poReader->PrepareStatement(osSQL.c_str());
int rowId = 0;
int nInvalid = 0;
int nGeometriesCount = 0;
while(poReader->ExecuteSQL(hStmt) == OGRERR_NONE) {
rowId++; // =sqlite3_column_int(hStmt, 1);
#ifdef DEBUG
const GIntBig iFID = sqlite3_column_int64(hStmt, 2);
#endif
VFKFeatureSQLite *poFeature = (VFKFeatureSQLite *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
// read geometry from DB
const int nBytes = sqlite3_column_bytes(hStmt, 0);
OGRGeometry *poGeometry = NULL;
if (nBytes > 0 &&
OGRGeometryFactory::createFromWkb((GByte*) sqlite3_column_blob(hStmt, 0),
NULL, &poGeometry, nBytes) == OGRERR_NONE) {
nGeometriesCount++;
if (!poFeature->SetGeometry(poGeometry)) {
nInvalid++;
}
delete poGeometry;
}
else {
nInvalid++;
}
}
CPLDebug("OGR-VFK", "%s: %d geometries loaded from DB",
m_pszName, nGeometriesCount);
if (nGeometriesCount != nGeometries) {
CPLError(CE_Warning, CPLE_AppDefined,
"%s: %d geometries loaded (should be %d)",
m_pszName, nGeometriesCount, nGeometries);
}
if (nInvalid > 0 && !bSkipInvalid) {
CPLError(CE_Warning, CPLE_AppDefined,
"%s: %d features with invalid or empty geometry",
m_pszName, nInvalid);
}
return true;
}
/*!
\brief Load geometry from DB
\return TRUE geometry successfully loaded otherwise FALSE
*/
bool VFKDataBlockDB::LoadGeometryFromDB()
{
int nInvalid, nGeometries, nGeometriesCount, nBytes, rowId;
#ifdef DEBUG
GIntBig iFID;
#endif
bool bSkipInvalid;
CPLString osSQL;
OGRGeometry *poGeometry;
VFKFeatureDB *poFeature;
VFKReaderDB *poReader;
poReader = (VFKReaderDB*) m_poReader;
if (!poReader->IsSpatial()) /* check if DB is spatial */
return FALSE;
osSQL.Printf("SELECT num_geometries FROM %s WHERE table_name = '%s'",
VFK_DB_TABLE, m_pszName);
if (poReader->ExecuteSQL(osSQL.c_str(), nGeometries) != OGRERR_NONE)
return FALSE;
if (nGeometries < 1)
return FALSE;
bSkipInvalid = EQUAL(m_pszName, "OB") || EQUAL(m_pszName, "OP") || EQUAL(m_pszName, "OBBP");
/* load geometry from DB */
nInvalid = nGeometriesCount = 0;
osSQL.Printf("SELECT %s,rowid,%s FROM %s ",
GEOM_COLUMN, FID_COLUMN, m_pszName);
if (EQUAL(m_pszName, "SBP"))
osSQL += "WHERE PORADOVE_CISLO_BODU = 1 ";
osSQL += "ORDER BY ";
osSQL += FID_COLUMN;
rowId = 0;
std::vector<VFKDbValue> record;
record.push_back(VFKDbValue(DT_INT)); // TODO: blob
record.push_back(VFKDbValue(DT_INT));
record.push_back(VFKDbValue(DT_BIGINT));
poReader->PrepareStatement(osSQL.c_str());
while(poReader->ExecuteSQL(record) == OGRERR_NONE) {
rowId++; // =sqlite3_column_int(hStmt, 1);
#ifdef DEBUG
iFID = static_cast<GIntBig> (record[2]);
#endif
poFeature = (VFKFeatureDB *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
// read geometry from DB
nBytes = 0; // TODO: sqlite3_column_byte(hStmt, 0);
/* TODO: solve
if (nBytes > 0 &&
OGRGeometryFactory::createFromWkb((GByte*) sqlite3_column_blob(hStmt, 0),
NULL, &poGeometry, nBytes) == OGRERR_NONE) {
nGeometriesCount++;
if (!poFeature->SetGeometry(poGeometry)) {
nInvalid++;
}
delete poGeometry;
}
else {
nInvalid++;
}
*/
}
CPLDebug("OGR-VFK", "%s: %d geometries loaded from DB",
m_pszName, nGeometriesCount);
if (nGeometriesCount != nGeometries) {
CPLError(CE_Warning, CPLE_AppDefined,
"%s: %d geometries loaded (should be %d)",
m_pszName, nGeometriesCount, nGeometries);
}
if (nInvalid > 0 && !bSkipInvalid) {
CPLError(CE_Warning, CPLE_AppDefined,
"%s: %d features with invalid or empty geometry",
m_pszName, nInvalid);
}
return TRUE;
}
/*!
\brief Load geometry (linestring HP/DPM layer)
\return number of invalid features
*/
int VFKDataBlockSQLite::LoadGeometryLineStringHP()
{
int nInvalid = 0;
VFKReaderSQLite *poReader = (VFKReaderSQLite*) m_poReader;
VFKDataBlockSQLite *poDataBlockLines =
(VFKDataBlockSQLite *) m_poReader->GetDataBlock("SBP");
if (NULL == poDataBlockLines) {
CPLError(CE_Failure, CPLE_FileIO,
"Data block %s not found.", m_pszName);
return nInvalid;
}
poDataBlockLines->LoadGeometry();
if (LoadGeometryFromDB()) /* try to load geometry from DB */
return 0;
CPLString osColumn;
osColumn.Printf("%s_ID", m_pszName);
const char *vrColumn[2] = {
osColumn.c_str(),
"PORADOVE_CISLO_BODU"
};
GUIntBig vrValue[2] = { 0, 1 }; // Reduce to first segment.
CPLString osSQL;
osSQL.Printf("SELECT ID,%s,rowid FROM %s", FID_COLUMN, m_pszName);
/* TODO: handle points in DPM */
if (EQUAL(m_pszName, "DPM"))
osSQL += " WHERE SOURADNICE_X IS NULL";
sqlite3_stmt *hStmt = poReader->PrepareStatement(osSQL.c_str());
if (poReader->IsSpatial())
poReader->ExecuteSQL("BEGIN");
int nGeometries = 0;
while( poReader->ExecuteSQL(hStmt) == OGRERR_NONE )
{
/* read values */
vrValue[0] = sqlite3_column_int64(hStmt, 0);
const long iFID = static_cast<long>(sqlite3_column_int64(hStmt, 1));
const int rowId = sqlite3_column_int(hStmt, 2);
VFKFeatureSQLite *poFeature =
(VFKFeatureSQLite *) GetFeatureByIndex(rowId - 1);
CPLAssert(NULL != poFeature && poFeature->GetFID() == iFID);
VFKFeatureSQLite *poLine =
poDataBlockLines->GetFeature(vrColumn, vrValue, 2, TRUE);
OGRGeometry *poOgrGeometry = NULL;
if( !poLine )
{
poOgrGeometry = NULL;
}
else
{
poOgrGeometry = poLine->GetGeometry();
}
if (!poOgrGeometry || !poFeature->SetGeometry(poOgrGeometry)) {
CPLDebug("OGR-VFK", "VFKDataBlockSQLite::LoadGeometryLineStringHP(): name=%s fid=%ld "
"id=" CPL_FRMT_GUIB " -> %s geometry", m_pszName, iFID, vrValue[0],
poOgrGeometry ? "invalid" : "empty");
nInvalid++;
continue;
}
/* store also geometry in DB */
if (poReader->IsSpatial() &&
SaveGeometryToDB(poOgrGeometry, rowId) != OGRERR_FAILURE &&
poOgrGeometry)
nGeometries++;
}
/* update number of geometries in VFK_DB_TABLE table */
UpdateVfkBlocks(nGeometries);
if (poReader->IsSpatial())
poReader->ExecuteSQL("COMMIT");
return nInvalid;
}
