/**********************************************************************
* DumpCoordsys()
*
* Open a .TAB file and dump coordsys info
**********************************************************************/
static int DumpCoordsys(const char *pszFname)
{
IMapInfoFile *poFile;
double dXMin, dYMin, dXMax, dYMax;
/*---------------------------------------------------------------------
* Try to open source file
*--------------------------------------------------------------------*/
if ((poFile = IMapInfoFile::SmartOpen(pszFname)) == NULL)
{
printf("Failed to open %s\n", pszFname);
return -1;
}
OGRSpatialReference *poSRS = poFile->GetSpatialRef();
char *pszCoordSys = MITABSpatialRef2CoordSys(poSRS);
char *pszProjString=NULL;
printf("CoordSys %s\n", pszCoordSys?pszCoordSys:"(null)");
if (poFile->GetBounds(dXMin, dYMin, dXMax, dYMax) == 0)
{
printf(" Proj. Bounds (%.15g %.15g) (%.15g %.15g)\n", dXMin, dYMin, dXMax, dYMax);
printf(" dX dY = %.15g %.15g\n", dXMax - dXMin, dYMax - dYMin);
}
else
{
printf(" Projection Bounds not available!\n");
}
OGREnvelope oEnv;
if (poFile->GetExtent(&oEnv, TRUE) == 0)
{
printf(" Data Extents (%.15g %.15g) (%.15g %.15g)\n", oEnv.MinX, oEnv.MinY, oEnv.MaxX, oEnv.MaxY);
}
else
{
printf(" Data Extents not available!\n");
}
if (poSRS)
{
// poSRS->exportToWkt(&pszProjString);
// printf(" WKT SRS = %s\n", pszProjString);
// CPLFree(pszProjString);
poSRS->exportToProj4(&pszProjString);
printf(" PROJ4 SRS = %s\n", pszProjString);
// Write bounds to a file and launch 'proj' to convert them to LAT/LON
FILE *fpOut;
if (pszProjString &&
(fpOut = fopen("/tmp/tttbounds.txt", "w")))
{
fprintf(fpOut, "%.15g %.15g\n", dXMin, dYMin);
fprintf(fpOut, "%.15g %.15g\n", dXMax, dYMax);
fclose(fpOut);
fflush(stdout);
system(CPLSPrintf("proj -I %s /tmp/tttbounds.txt", pszProjString));
}
}
/*---------------------------------------------------------------------
* Cleanup and exit.
*--------------------------------------------------------------------*/
CPLFree(pszProjString);
poFile->Close();
delete poFile;
return 0;
}
OGRLayer *
OGRTABDataSource::ICreateLayer( const char *pszLayerName,
OGRSpatialReference *poSRSIn,
OGRwkbGeometryType /* eGeomTypeIn */,
char **papszOptions )
{
if( !m_bUpdate )
{
CPLError( CE_Failure, CPLE_AppDefined,
"Cannot create layer on read-only dataset.");
return nullptr;
}
// If it is a single file mode file, then we may have already
// instantiated the low level layer. We would just need to
// reset the coordinate system and (potentially) bounds.
IMapInfoFile *poFile = nullptr;
char *pszFullFilename = nullptr;
const char *pszEncoding = CSLFetchNameValue( papszOptions, "ENCODING" );
const char *pszCharset( IMapInfoFile::EncodingToCharset( pszEncoding ) );
if( m_bSingleFile )
{
if( m_bSingleLayerAlreadyCreated )
{
CPLError(
CE_Failure, CPLE_AppDefined,
"Unable to create new layers in this single file dataset.");
return nullptr;
}
m_bSingleLayerAlreadyCreated = TRUE;
poFile = (IMapInfoFile *) m_papoLayers[0];
if( pszEncoding )
poFile->SetCharset( pszCharset );
}
else
{
if( m_bCreateMIF )
{
pszFullFilename =
CPLStrdup(CPLFormFilename(m_pszDirectory, pszLayerName, "mif"));
poFile = new MIFFile;
if( poFile->Open(pszFullFilename, TABWrite,
FALSE, pszCharset) != 0 )
{
CPLFree(pszFullFilename);
delete poFile;
return nullptr;
}
}
else
{
pszFullFilename =
CPLStrdup(CPLFormFilename(m_pszDirectory, pszLayerName, "tab"));
TABFile *poTABFile = new TABFile;
if( poTABFile->Open(pszFullFilename, TABWrite, FALSE,
m_nBlockSize, pszCharset) != 0 )
{
CPLFree(pszFullFilename);
delete poTABFile;
return nullptr;
}
poFile = poTABFile;
}
m_nLayerCount++;
m_papoLayers = static_cast<IMapInfoFile **>(
CPLRealloc(m_papoLayers,sizeof(void*)*m_nLayerCount));
m_papoLayers[m_nLayerCount-1] = poFile;
CPLFree(pszFullFilename);
}
poFile->SetDescription(poFile->GetName());
// Assign the coordinate system (if provided) and set
// reasonable bounds.
if( poSRSIn != nullptr )
{
poFile->SetSpatialRef(poSRSIn);
// SetSpatialRef() has cloned the passed geometry
poFile->GetLayerDefn()->GetGeomFieldDefn(0)->SetSpatialRef(
poFile->GetSpatialRef());
}
// Pull out the bounds if supplied
const char *pszOpt = nullptr;
if( (pszOpt = CSLFetchNameValue(papszOptions, "BOUNDS")) != nullptr ) {
double dfBounds[4];
if( CPLsscanf(pszOpt, "%lf,%lf,%lf,%lf", &dfBounds[0],
&dfBounds[1],
&dfBounds[2],
&dfBounds[3]) != 4 )
{
CPLError(
CE_Failure, CPLE_IllegalArg,
"Invalid BOUNDS parameter, expected min_x,min_y,max_x,max_y");
}
else
{
poFile->SetBounds(dfBounds[0], dfBounds[1], dfBounds[2],
dfBounds[3]);
}
}
if( !poFile->IsBoundsSet() && !m_bCreateMIF )
{
if( poSRSIn != nullptr && poSRSIn->GetRoot() != nullptr &&
EQUAL(poSRSIn->GetRoot()->GetValue(),"GEOGCS") )
poFile->SetBounds(-1000, -1000, 1000, 1000);
else
poFile->SetBounds(-30000000, -15000000, 30000000, 15000000);
}
if(m_bQuickSpatialIndexMode == TRUE &&
poFile->SetQuickSpatialIndexMode(TRUE) != 0)
{
CPLError(CE_Warning, CPLE_AppDefined,
"Setting Quick Spatial Index Mode failed.");
}
else if(m_bQuickSpatialIndexMode == FALSE &&
poFile->SetQuickSpatialIndexMode(FALSE) != 0)
{
CPLError(CE_Warning, CPLE_AppDefined,
"Setting Normal Spatial Index Mode failed.");
}
return poFile;
}
/**********************************************************************
* Tab2Tab()
*
* Copy features from source dataset to a new dataset
**********************************************************************/
static int Tab2Tab(const char *pszSrcFname, const char *pszDstFname,
int nMaxFeatures,
GBool bQuickSpatialIndexMode, GBool bOptSpatialIndexMode)
{
IMapInfoFile *poSrcFile = NULL, *poDstFile = NULL;
int nFeatureId, iField, numFeatures=0;
TABFeature *poFeature;
double dXMin, dYMin, dXMax, dYMax;
/*---------------------------------------------------------------------
* If there is a "micdsys.txt" in current directory then load it
*--------------------------------------------------------------------*/
MITABLoadCoordSysTable("micdsys.txt");
/*---------------------------------------------------------------------
* Try to open source file
*--------------------------------------------------------------------*/
if ((poSrcFile = IMapInfoFile::SmartOpen(pszSrcFname)) == NULL)
{
printf("Failed to open %s\n", pszSrcFname);
return -1;
}
OGRFeatureDefn *poDefn = poSrcFile->GetLayerDefn();
/*---------------------------------------------------------------------
* The extension of the output filename tells us if we should create
* a MIF or a TAB file for output.
*--------------------------------------------------------------------*/
if (EQUAL(".mif", pszDstFname + strlen(pszDstFname)-4) ||
EQUAL(".mid", pszDstFname + strlen(pszDstFname)-4) )
{
// Create a MIF file
poDstFile = new MIFFile;
}
else
{
/*-----------------------------------------------------------------
* Create a TAB dataset.
* Find out if the file contains at least 1 unique field... if so we
* will create a TABView instead of a TABFile
*----------------------------------------------------------------*/
GBool bFoundUniqueField = FALSE;
for(iField=0; iField< poDefn->GetFieldCount(); iField++)
{
if (poSrcFile->IsFieldUnique(iField))
bFoundUniqueField = TRUE;
}
if (bFoundUniqueField)
poDstFile = new TABView;
else
poDstFile = new TABFile;
}
/*---------------------------------------------------------------------
* Try to open destination file
*--------------------------------------------------------------------*/
if (poDstFile->Open(pszDstFname, "wb") != 0)
{
printf("Failed to open %s\n", pszDstFname);
poSrcFile->Close();
delete poSrcFile;
delete poDstFile;
return -1;
}
if ( (bQuickSpatialIndexMode &&
poDstFile->SetQuickSpatialIndexMode(TRUE) != 0) ||
(bOptSpatialIndexMode &&
poDstFile->SetQuickSpatialIndexMode(FALSE) != 0) )
{
printf("Failed setting Quick Spatial Index Mode (-q|-o) on %s\n", pszDstFname);
poSrcFile->Close();
delete poSrcFile;
poDstFile->Close();
delete poDstFile;
return -1;
}
// Pass Proj. info directly
// TABProjInfo sProjInfo;
// if (poSrcFile->GetProjInfo(&sProjInfo) == 0)
// poDstFile->SetProjInfo(&sProjInfo);
OGRSpatialReference *poSR;
poSR = poSrcFile->GetSpatialRef();
if( poSR != NULL )
{
poDstFile->SetSpatialRef( poSR );
}
// Set bounds (must be done after setting spatialref)
if (poSrcFile->GetBounds(dXMin, dYMin, dXMax, dYMax) == 0)
poDstFile->SetBounds(dXMin, dYMin, dXMax, dYMax);
/*---------------------------------------------------------------------
* Pass compplete fields information
*--------------------------------------------------------------------*/
for(iField=0; iField< poDefn->GetFieldCount(); iField++)
{
OGRFieldDefn *poFieldDefn = poDefn->GetFieldDefn(iField);
poDstFile->AddFieldNative(poFieldDefn->GetNameRef(),
poSrcFile->GetNativeFieldType(iField),
poFieldDefn->GetWidth(),
poFieldDefn->GetPrecision(),
poSrcFile->IsFieldIndexed(iField),
poSrcFile->IsFieldUnique(iField));
}
/*---------------------------------------------------------------------
* Copy objects until EOF is reached
*--------------------------------------------------------------------*/
nFeatureId = -1;
while ( (nFeatureId = poSrcFile->GetNextFeatureId(nFeatureId)) != -1 &&
(nMaxFeatures < 1 || numFeatures++ < nMaxFeatures ))
{
poFeature = poSrcFile->GetFeatureRef(nFeatureId);
if (poFeature)
{
// poFeature->DumpReadable(stdout);
// poFeature->DumpMIF();
poDstFile->CreateFeature(poFeature);
}
else
{
printf( "Failed to read feature %d.\n",
nFeatureId );
return -1; // GetFeatureRef() failed: Error
}
}
/*---------------------------------------------------------------------
* Cleanup and exit.
*--------------------------------------------------------------------*/
poDstFile->Close();
delete poDstFile;
poSrcFile->Close();
delete poSrcFile;
MITABFreeCoordSysTable();
return 0;
}
