protected function createNewLayer($map, $driver, $layerName, $geoType, $bSearchSpatialRef)
{
// if layer contains point objects
if ($geoType == wkbPoint) {
// before creating layer "entities" we create layer "blocks"
$layerBlocks = OGR_DS_CreateLayer($driver, 'blocks', NULL, wkbMultiLineString, NULL);
// in this layer we export representaion of all point objects (it's a square)
$wktLine = 'MULTILINESTRING((-0.5 -0.5, 0.5 -0.5, 0.5 0.5, -0.5 0.5, -0.5 -0.5))';
$geoOGR = ogr_g_createfromwkt($wktLine, NULL, NULL);
$dfnLayerBlocks = OGR_L_GetLayerDefn($layerBlocks);
$dataLine = OGR_F_Create($dfnLayerBlocks);
OGR_F_SetGeometry($dataLine, $geoOGR);
for ($i = 0; $i < OGR_FD_GetFieldCount($dfnLayerBlocks); $i++) {
$fieldDfn = OGR_FD_GetFieldDefn($dfnLayerBlocks, $i);
if (OGR_Fld_GetNameRef($fieldDfn) == 'BlockName') {
OGR_F_SetFieldString($dataLine, $i, 'SQUARE');
}
}
OGR_L_CreateFeature($layerBlocks, $dataLine);
OGR_F_Destroy($dataLine);
OGR_G_DestroyGeometry($geoOGR);
}
return parent::createNewLayer($map, $driver, 'entities', $geoType, false);
}
protected function createNewLayer($map, $driver, $layerName, $geoType, $bSearchSpatialRef)
{
$spatialRef = NULL;
// Looking for spatial référence handle
if ($bSearchSpatialRef === true) {
$shpFile = dirname(__FILE__) . "\\shp\\";
if (!isset($this->driverSHP)) {
// looking for MapInfo driver
for ($iDriver = 0; $iDriver < OGRGetDriverCount(); $iDriver++) {
if (!strcasecmp(OGR_DR_GetName(OGRGetDriver($iDriver)), 'ESRI Shapefile')) {
$this->driverSHP = OGRGetDriver($iDriver);
break;
}
}
}
if (isset($this->driverSHP)) {
// reading layer pojection
if ($map) {
$layerSrc = $map->getLayerByName($layerName);
$proj = $layerSrc->getProjection();
if (empty($proj)) {
$proj = $map->getProjection();
}
if (!empty($proj)) {
$pos = strpos($proj, 'epsg:');
if ($pos !== false) {
$proj = substr($proj, $pos + 5);
$shpFile .= "{$proj}.shp";
}
}
}
// if MapInfo file with same pojection exists
if (file_exists($shpFile) && !is_dir($shpFile)) {
// opening file with OGR and reading OGR projection
$projDataSource = ogr_dr_open($this->driverSHP, $shpFile, 0);
if ($projDataSource && OGR_DS_GetLayerCount($projDataSource) > 0) {
$layer = OGR_DS_GetLayer($projDataSource, 0);
if ($layer) {
$spatialRef = OGR_L_GetSpatialRef($layer);
}
}
}
}
}
$layerDest = OGR_DS_CreateLayer($driver, $layerName, $spatialRef, $geoType, NULL);
return $layerDest;
}
function TranslateLayer($hSrcDS, $hSrcLayer, $hDstDS)
{
/* -------------------------------------------------------------------- */
/* Create the layer. */
/* -------------------------------------------------------------------- */
if (!OGR_DS_TestCapability($hDstDS, ODsCCreateLayer)) {
printf("%s data source does not support layer creation.\n", OGR_DS_GetName($hDstDS));
return OGRERR_FAILURE;
}
$hFDefn = OGR_L_GetLayerDefn($hSrcLayer);
$hDstLayer = OGR_DS_CreateLayer($hDstDS, OGR_FD_GetName($hFDefn), OGR_L_GetSpatialRef($hSrcLayer), OGR_FD_GetGeomType($hFDefn), NULL);
if ($hDstLayer == NULL) {
return FALSE;
}
/* -------------------------------------------------------------------- */
/* Add fields. */
/* -------------------------------------------------------------------- */
for ($iField = 0; $iField < OGR_FD_GetFieldCount($hFDefn); $iField++) {
if (OGR_L_CreateField($hDstLayer, OGR_FD_GetFieldDefn($hFDefn, $iField), 0) != OGRERR_NONE) {
return FALSE;
}
}
/* -------------------------------------------------------------------- */
/* Transfer features. */
/* -------------------------------------------------------------------- */
OGR_L_ResetReading($hSrcLayer);
while (($hFeature = OGR_L_GetNextFeature($hSrcLayer)) != NULL) {
$hDstFeature = OGR_F_Create(OGR_L_GetLayerDefn($hDstLayer));
if (OGR_F_SetFrom($hDstFeature, $hFeature, FALSE) != OGRERR_NONE) {
OGR_F_Destroy($hFeature);
printf("Unable to translate feature %d from layer %s.\n", OGR_F_GetFID($hFeature), OGR_FD_GetName($hFDefn));
return FALSE;
}
OGR_F_Destroy($hFeature);
if (OGR_L_CreateFeature($hDstLayer, $hDstFeature) != OGRERR_NONE) {
OGR_F_Destroy($hDstFeature);
return FALSE;
}
OGR_F_Destroy($hDstFeature);
}
return TRUE;
}
/**
*
* processas as camadas do arquivo importado
*
* @param string $hSrcDS
* @param array $hSrcLayer
* @param string $hDstDS
*/
public function translateLayer($hSrcDS, $hSrcLayer, $hDstDS)
{
// cria layer
if (!OGR_DS_TestCapability($hDstDS, ODsCCreateLayer)) {
printf(STATIC::GEO_UNABLE_CREATE_LAYER, OGR_DS_GetName($hDstDS));
return OGRERR_FAILURE;
}
$hFDefn = OGR_L_GetLayerDefn($hSrcLayer);
$hDstLayer = OGR_DS_CreateLayer($hDstDS, OGR_FD_GetName($hFDefn), OGR_L_GetSpatialRef($hSrcLayer), OGR_FD_GetGeomType($hFDefn), NULL);
if ($hDstLayer == NULL) {
return FALSE;
}
// adiciona campos
for ($iField = 0; $iField < OGR_FD_GetFieldCount($hFDefn); $iField++) {
if (OGR_L_CreateField($hDstLayer, OGR_FD_GetFieldDefn($hFDefn, $iField), 0) != OGRERR_NONE) {
return FALSE;
}
}
// processa os features do layer
OGR_L_ResetReading($hSrcLayer);
while (($hFeature = OGR_L_GetNextFeature($hSrcLayer)) != NULL) {
$hDstFeature = OGR_F_Create(OGR_L_GetLayerDefn($hDstLayer));
if (OGR_F_SetFrom($hDstFeature, $hFeature, FALSE) != OGRERR_NONE) {
OGR_F_Destroy($hFeature);
printf(static::GEO_UNABLE_TRANSLATE_FEATURE_LAYER, OGR_F_GetFID($hFeature), OGR_FD_GetName($hFDefn));
return FALSE;
}
OGR_F_Destroy($hFeature);
if (OGR_L_CreateFeature($hDstLayer, $hDstFeature) != OGRERR_NONE) {
OGR_F_Destroy($hDstFeature);
return FALSE;
}
OGR_F_Destroy($hDstFeature);
}
return TRUE;
}
function ogr2tab_copy_srs()
{
/*Corresponding MapInfo File driver number.*/
$iDriver = 5;
$strFormat = "MapInfo File";
/*Path to existing data source to modify to fit user needs.*/
$strSrcDataSource = "./phpunit_tests/data/tab/mpasabe200/Commune.TAB";
/*Path to new data source to modify to fit user needs.*/
$strDestDataSource = "../ogrtests/tmp/output";
/*New layer name to modify to fit user needs.*/
$strDestLayerName = "MyNewLayer";
/*Layer selected to fetch a feature from. If the data source name in
$strSrcDataSource is a "tab" file, the layer number, specified by
$iLayer, to copy from must be set to zero. Otherwise the layer
number might be any number between 0 and the layer count.*/
// $iLayer = 1;
$iLayer = 0;
/* -------------------------------------------------------------------- */
/* Register format(s). */
/* -------------------------------------------------------------------- */
OGRRegisterAll();
/* -------------------------------------------------------------------- */
/* Open the existing data source. */
/* -------------------------------------------------------------------- */
$hSFDriver = NULL;
$hSrcDS = OGROpen($strSrcDataSource, FALSE, $hSFDriver);
/* -------------------------------------------------------------------- */
/* Report failure */
/* -------------------------------------------------------------------- */
if ($hSrcDS == NULL) {
printf("FAILURE:\nUnable to open datasource `%s'\n", $strSrcDataSource);
return OGRERR_FAILURE;
}
/* -------------------------------------------------------------------- */
/* Find the output driver handle. */
/* -------------------------------------------------------------------- */
if ($hSFDriver == NULL) {
$hSFDriver = OGRGetDriver($iDriver);
}
if ($hSFDriver == NULL) {
printf("Unable to find driver `%s'.\n", $strFormat);
return OGRERR_FAILURE;
}
/* -------------------------------------------------------------------- */
/* Create the destination data source. */
/* -------------------------------------------------------------------- */
$hDestDS = OGR_Dr_CreateDataSource($hSFDriver, $strDestDataSource, NULL);
if ($hDestDS == NULL) {
return OGRERR_FAILURE;
}
/* -------------------------------------------------------------------- */
/* Select an existing layer from the existing data source to get */
/* the projection information from. To modify accordingly to user */
/* needs. */
/* -------------------------------------------------------------------- */
if (OGR_DS_GetLayerCount($hSrcDS) > 0) {
$hSrcLayer = OGR_DS_GetLayer($hSrcDS, $iLayer);
if ($hSrcLayer == NULL) {
printf("FAILURE: Couldn't fetch advertised layer %d!\n", $iLayer);
return OGRERR_FAILURE;
}
$hSrcSRS = OGR_L_GetSpatialRef($hSrcLayer);
}
/* -------------------------------------------------------------------- */
/* Create a new layer based on the existing data source */
/* information. */
/* -------------------------------------------------------------------- */
$hSrcFDefn = OGR_L_GetLayerDefn($hSrcLayer);
$eType = OGR_FD_GetGeomType($hSrcFDefn);
$hDestLayer = OGR_DS_CreateLayer($hDestDS, $strDestLayerName, $hSrcSRS, $eType, NULL);
if ($hDestLayer == NULL) {
return FALSE;
}
/* -------------------------------------------------------------------- */
/* Create the field definitions for the new layer based on the */
/* existing field definitions of the source layer. */
/* -------------------------------------------------------------------- */
for ($iField = 0; $iField < OGR_FD_GetFieldCount($hSrcFDefn); $iField++) {
if (OGR_L_CreateField($hDestLayer, OGR_FD_GetFieldDefn($hSrcFDefn, $iField), 0) != OGRERR_NONE) {
return OGRERR_FAILURE;
}
}
/* -------------------------------------------------------------------- */
/* Get only one feature on the existing layer and copy it to the */
/* new layer. */
/* -------------------------------------------------------------------- */
OGR_L_ResetReading($hSrcLayer);
$hSrcFeature = OGR_L_GetNextFeature($hSrcLayer);
$hDestFeature = OGR_F_Create(OGR_L_GetLayerDefn($hDestLayer));
if (OGR_F_SetFrom($hDestFeature, $hSrcFeature, FALSE) != OGRERR_NONE) {
OGR_F_Destroy($hDestFeature);
printf("Unable to copy feature %d from layer %s.\n", OGR_F_GetFID($hSrcFeature), OGR_FD_GetName($hSrcFDefn));
return OGRERR_FAILURE;
}
OGR_F_Destroy($hSrcFeature);
if (OGR_L_CreateFeature($hDestLayer, $hDestFeature) != OGRERR_NONE) {
OGR_F_Destroy($hDestFeature);
return OGRERR_FAILURE;
}
OGR_F_Destroy($hDestFeature);
/* -------------------------------------------------------------------- */
/* Close down. */
/* -------------------------------------------------------------------- */
OGR_DS_Destroy($hSrcDS);
OGR_DS_Destroy($hDestDS);
return OGRERR_NONE;
}
function TranslateLayer($hSrcDS, $hSrcLayer, $hDstDS, $amAttributeIn)
{
/* -------------------------------------------------------------------- */
/* Create the layer. */
/* -------------------------------------------------------------------- */
if (!OGR_DS_TestCapability($hDstDS, ODsCCreateLayer)) {
printf("%s data source does not support layer creation.\n", OGR_DS_GetName($hDstDS));
return OGRERR_FAILURE;
}
$hFDefn = OGR_L_GetLayerDefn($hSrcLayer);
$hDstLayer = OGR_DS_CreateLayer($hDstDS, OGR_FD_GetName($hFDefn), OGR_L_GetSpatialRef($hSrcLayer), OGR_FD_GetGeomType($hFDefn), NULL);
if ($hDstLayer == NULL) {
return FALSE;
}
$iFieldCount = OGR_FD_GetFieldCount($hFDefn);
printf("field count =%d\n", $iFieldCount);
/* -------------------------------------------------------------------- */
/* Add fields. */
/* -------------------------------------------------------------------- */
for ($iField = 0; $iField < $iFieldCount; $iField++) {
if (OGR_L_CreateField($hDstLayer, OGR_FD_GetFieldDefn($hFDefn, $iField), 0) != OGRERR_NONE) {
return FALSE;
}
}
/* -------------------------------------------------------------------- */
/* Transfer features. */
/* -------------------------------------------------------------------- */
OGR_L_ResetReading($hSrcLayer);
while (($hFeature = OGR_L_GetNextFeature($hSrcLayer)) != NULL) {
$hDstFeature = OGR_F_Create(OGR_L_GetLayerDefn($hDstLayer));
/* Verify if the current feature is corresponding to the feature
to update the attribute to. If yes, the value corresponding to
the attribute name selected is substituted. All other attributes
and geometry are copied as is. If not, the source feature
is copied to the destination feature as is.*/
if ($amAttributeIn && OGR_F_GetFID($hFeature) == $amAttributeIn[0]) {
for ($i = 0; $i < $iFieldCount; $i++) {
$hCurrentFieldDefn = OGR_F_GetFieldDefnRef($hFeature, $i);
$strCurrentFieldName = OGR_FLD_GetNameRef($hCurrentFieldDefn);
if ($strCurrentFieldName == $amAttributeIn[1]) {
$value = $amAttributeIn[2];
} else {
$value = OGR_F_GetFieldAsString($hFeature, $i);
}
OGR_F_SetFieldString($hDstFeature, $i, $value);
}
$hGeometry = OGR_F_GetGeometryRef($hFeature);
OGR_F_SetGeometry($hDstFeature, $hGeometry);
} else {
if (OGR_F_SetFrom($hDstFeature, $hFeature, FALSE) != OGRERR_NONE) {
OGR_F_Destroy($hFeature);
printf("Unable to translate feature %d from layer %s.\n", OGR_F_GetFID($hFeature), OGR_FD_GetName($hFDefn));
return FALSE;
}
}
OGR_F_Destroy($hFeature);
if (OGR_L_CreateFeature($hDstLayer, $hDstFeature) != OGRERR_NONE) {
OGR_F_Destroy($hDstFeature);
return FALSE;
}
OGR_F_Destroy($hDstFeature);
}
return TRUE;
}
