void OGRPGResultLayer::SetSpatialFilter( int iGeomField, OGRGeometry * poGeomIn )
{
if( iGeomField < 0 || iGeomField >= GetLayerDefn()->GetGeomFieldCount() ||
GetLayerDefn()->GetGeomFieldDefn(iGeomField)->GetType() == wkbNone )
{
if( iGeomField != 0 )
{
CPLError(CE_Failure, CPLE_AppDefined,
"Invalid geometry field index : %d", iGeomField);
}
return;
}
m_iGeomFieldFilter = iGeomField;
OGRPGGeomFieldDefn* poGeomFieldDefn =
poFeatureDefn->myGetGeomFieldDefn(m_iGeomFieldFilter);
if( InstallFilter( poGeomIn ) )
{
if ( poGeomFieldDefn->ePostgisType == GEOM_TYPE_GEOMETRY ||
poGeomFieldDefn->ePostgisType == GEOM_TYPE_GEOGRAPHY )
{
if( m_poFilterGeom != NULL)
{
char szBox3D_1[128];
char szBox3D_2[128];
OGREnvelope sEnvelope;
m_poFilterGeom->getEnvelope( &sEnvelope );
if( poGeomFieldDefn->ePostgisType == GEOM_TYPE_GEOGRAPHY )
{
if( sEnvelope.MinX < -180.0 )
sEnvelope.MinX = -180.0;
if( sEnvelope.MinY < -90.0 )
sEnvelope.MinY = -90.0;
if( sEnvelope.MaxX > 180.0 )
sEnvelope.MaxX = 180.0;
if( sEnvelope.MaxY > 90.0 )
sEnvelope.MaxY = 90.0;
}
CPLsnprintf(szBox3D_1, sizeof(szBox3D_1), "%.18g %.18g", sEnvelope.MinX, sEnvelope.MinY);
CPLsnprintf(szBox3D_2, sizeof(szBox3D_2), "%.18g %.18g", sEnvelope.MaxX, sEnvelope.MaxY);
osWHERE.Printf("WHERE %s && %s('BOX3D(%s, %s)'::box3d,%d) ",
OGRPGEscapeColumnName(poGeomFieldDefn->GetNameRef()).c_str(),
(poDS->sPostGISVersion.nMajor >= 2) ? "ST_SetSRID" : "SetSRID",
szBox3D_1, szBox3D_2, poGeomFieldDefn->nSRSId );
}
else
{
osWHERE = "";
}
BuildFullQueryStatement();
}
ResetReading();
}
}
void OGRPGResultLayer::ResolveSRID(OGRPGGeomFieldDefn* poGFldDefn)
{
/* We have to get the SRID of the geometry column, so to be able */
/* to do spatial filtering */
int nSRSId = UNDETERMINED_SRID;
if( poGFldDefn->ePostgisType == GEOM_TYPE_GEOMETRY )
{
if (pszGeomTableName != NULL)
{
CPLString osName(pszGeomTableSchemaName);
osName += ".";
osName += pszGeomTableName;
OGRPGLayer* poBaseLayer = (OGRPGLayer*) poDS->GetLayerByName(osName);
if (poBaseLayer)
{
int iBaseIdx = poBaseLayer->GetLayerDefn()->
GetGeomFieldIndex( poGFldDefn->GetNameRef() );
if( iBaseIdx >= 0 )
{
OGRPGGeomFieldDefn* poBaseGFldDefn =
poBaseLayer->myGetLayerDefn()->myGetGeomFieldDefn(iBaseIdx);
poBaseGFldDefn->GetSpatialRef(); /* To make sure nSRSId is resolved */
nSRSId = poBaseGFldDefn->nSRSId;
}
}
}
if( nSRSId == UNDETERMINED_SRID )
{
CPLString osGetSRID;
const char* psGetSRIDFct;
if (poDS->sPostGISVersion.nMajor >= 2)
psGetSRIDFct = "ST_SRID";
else
psGetSRIDFct = "getsrid";
osGetSRID += "SELECT ";
osGetSRID += psGetSRIDFct;
osGetSRID += "(";
osGetSRID += OGRPGEscapeColumnName(poGFldDefn->GetNameRef());
osGetSRID += ") FROM(";
osGetSRID += pszRawStatement;
osGetSRID += ") AS ogrpggetsrid LIMIT 1";
PGresult* hSRSIdResult = OGRPG_PQexec(poDS->GetPGConn(), osGetSRID );
nSRSId = poDS->GetUndefinedSRID();
if( hSRSIdResult && PQresultStatus(hSRSIdResult) == PGRES_TUPLES_OK)
{
if ( PQntuples(hSRSIdResult) > 0 )
nSRSId = atoi(PQgetvalue(hSRSIdResult, 0, 0));
}
else
{
CPLError( CE_Failure, CPLE_AppDefined,
"%s", PQerrorMessage(poDS->GetPGConn()) );
}
OGRPGClearResult(hSRSIdResult);
}
}
else if( poGFldDefn->ePostgisType == GEOM_TYPE_GEOGRAPHY )
{
nSRSId = 4326;
}
poGFldDefn->nSRSId = nSRSId;
}
