service.c @ 928

Last change on this file since 928 was 917, checked in by djay, 6 years ago
Merge prototype-v0 branch in trunk
Property svn:eol-style set to `native` Property svn:mime-type set to `text/x-csrc`
File size: 47.0 KB

Rev	Line
[1]	1	/******************************************************************************
	2	* $Id: ogr2ogr.cpp 15473 2008-10-07 20:59:24Z warmerdam $
	3	*
	4	* Project: OpenGIS Simple Features Reference Implementation
	5	* Purpose: Simple client for translating between formats.
	6	* Author: Frank Warmerdam, warmerdam@pobox.com
	7	*
	8	******************************************************************************
	9	* Copyright (c) 1999, Frank Warmerdam
	10	*
	11	* Permission is hereby granted, free of charge, to any person obtaining a
	12	* copy of this software and associated documentation files (the "Software"),
	13	* to deal in the Software without restriction, including without limitation
	14	* the rights to use, copy, modify, merge, publish, distribute, sublicense,
	15	* and/or sell copies of the Software, and to permit persons to whom the
	16	* Software is furnished to do so, subject to the following conditions:
	17	*
	18	* The above copyright notice and this permission notice shall be included
	19	* in all copies or substantial portions of the Software.
	20	*
	21	* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
	22	* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
	23	* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
	24	* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
	25	* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
	26	* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
	27	* DEALINGS IN THE SOFTWARE.
	28	****************************************************************************/
	29
	30	#include "ogrsf_frmts.h"
	31	#include "ogr_p.h"
	32	#include "cpl_conv.h"
	33	#include "cpl_string.h"
	34	#include "ogr_api.h"
[766]	35	#if GDAL_VERSION_MAJOR >= 2
	36	#include <gdal_priv.h>
	37	#endif
[1]	38	#ifdef ZOO_SERVICE
	39	#include "service.h"
	40	#endif
	41
	42	CPL_CVSID("$Id: ogr2ogr.cpp 15473 2008-10-07 20:59:24Z warmerdam $");
	43
[395]	44	#ifdef WIN32
	45	#define strcasecmp _stricmp
	46	#define strncasecmp _strnicmp
	47	#endif
	48
[1]	49	#ifdef ZOO_SERVICE
	50	extern "C" {
	51	#endif
	52
	53	static void Usage();
	54
[766]	55	static int TranslateLayer(
	56	#if GDAL_VERSION_MAJOR >= 2
	57	GDALDataset
	58	#else
	59	OGRDataSource
	60	#endif
	61	*poSrcDS,
	62	OGRLayer * poSrcLayer,
	63	#if GDAL_VERSION_MAJOR >= 2
	64	GDALDataset
	65	#else
	66	OGRDataSource
	67	#endif
	68	*poDstDS,
	69	char ** papszLSCO,
	70	const char *pszNewLayerName,
	71	int bTransform,
	72	OGRSpatialReference *poOutputSRS,
	73	OGRSpatialReference *poSourceSRS,
	74	char **papszSelFields,
	75	int bAppend, int eGType,
	76	int bOverwrite,
	77	double dfMaxSegmentLength
	78	);
	79
[1]	80	static int bSkipFailures = FALSE;
	81	static int nGroupTransactions = 200;
[395]	82	static int bPreserveFID = TRUE;
[1]	83	static int nFIDToFetch = OGRNullFID;
	84
	85	/************************************************************************/
	86	/* main() */
	87	/************************************************************************/
	88
	89	#ifdef ZOO_SERVICE
	90	#ifdef WIN32
	91	__declspec(dllexport)
	92	#endif
	93	int Ogr2Ogr(maps& conf,maps& inputs,maps*& outputs)
	94	#else
	95	int main( int nArgc, char ** papszArgv )
	96	#endif
	97	{
[917]	98	char *pszDialect = NULL;
[1]	99	const char *pszFormat = "ESRI Shapefile";
	100	const char *pszDataSource = NULL;
	101	const char *pszDestDataSource = NULL;
[49]	102	const char *pszwebDestData = NULL;
[1]	103	char **papszLayers = NULL;
	104	char papszDSCO = NULL, papszLCO = NULL;
	105	int bTransform = FALSE;
	106	int bAppend = FALSE, bUpdate = FALSE, bOverwrite = FALSE;
	107	const char *pszOutputSRSDef = NULL;
	108	const char *pszSourceSRSDef = NULL;
	109	OGRSpatialReference *poOutputSRS = NULL;
	110	OGRSpatialReference *poSourceSRS = NULL;
	111	const char *pszNewLayerName = NULL;
	112	const char *pszWHERE = NULL;
	113	OGRGeometry *poSpatialFilter = NULL;
	114	const char *pszSelect;
	115	char **papszSelFields = NULL;
	116	const char *pszSQLStatement = NULL;
	117	int eGType = -2;
[252]	118	double dfMaxSegmentLength = 0;
[1]	119
[395]	120	#ifdef ZOO_SERVICE
	121	dumpMaps(inputs);
	122	#endif
[1]	123	/* Check strict compilation and runtime library version as we use C++ API */
	124	if (! GDAL_CHECK_VERSION("ogr2ogr"))
	125	#ifdef ZOO_SERVICE
	126	{
	127	fprintf(stderr,"Not correct version of the gdal library\n");
[38]	128	setMapInMaps(conf,"lenv","message","Unable to check gdal version for ogr2ogr_service.zo");
[1]	129	return SERVICE_FAILED;
	130	}
	131	#else
	132	exit(1);
	133	#endif
	134	/* -------------------------------------------------------------------- */
	135	/* Register format(s). */
	136	/* -------------------------------------------------------------------- */
	137	OGRRegisterAll();
	138
	139	#ifdef ZOO_SERVICE
[40]	140	map *tmpMap=NULL;
[1]	141	char dataPath[1024];
	142	tmpMap=getMapFromMaps(conf,"main","dataPath");
	143	if(tmpMap!=NULL)
	144	sprintf(dataPath,"%s",tmpMap->value);
	145	tmpMap=NULL;
	146	char tempPath[1024];
	147	tmpMap=getMapFromMaps(conf,"main","tmpPath");
	148	if(tmpMap!=NULL){
	149	sprintf(tempPath,"%s",tmpMap->value);
	150	}
[49]	151
	152	tmpMap=NULL;
	153	char serverAddress[1024];
	154	tmpMap=getMapFromMaps(conf,"main","serverAddress");
	155	if(tmpMap!=NULL){
	156	sprintf(serverAddress,"%s",tmpMap->value);
	157	}
[917]	158
[49]	159	tmpMap=NULL;
	160	char tmpurl[1024];
	161	tmpMap=getMapFromMaps(conf,"main","tmpurl");
	162	if(tmpMap!=NULL){
	163	sprintf(tmpurl,"%s",tmpMap->value);
	164	}
[1]	165
[40]	166	tmpMap=NULL;
[1]	167	tmpMap=getMapFromMaps(inputs,"F","value");
[395]	168	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	169	pszFormat=tmpMap->value;
	170	}
	171
[40]	172	tmpMap=NULL;
[1]	173	tmpMap=getMapFromMaps(inputs,"DSCO","value");
[395]	174	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	175	papszDSCO = CSLAddString(papszDSCO, tmpMap->value );
	176	}
	177
[40]	178	tmpMap=NULL;
[1]	179	tmpMap=getMapFromMaps(inputs,"LCO","value");
[395]	180	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	181	papszLCO = CSLAddString(papszLCO, tmpMap->value );
	182	}
	183
[40]	184	tmpMap=NULL;
[1]	185	tmpMap=getMapFromMaps(inputs,"preserve_fid","value");
[395]	186	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	187	bPreserveFID = TRUE;
	188	}
	189
[40]	190	tmpMap=NULL;
[1]	191	tmpMap=getMapFromMaps(inputs,"skipfailure","value");
[395]	192	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	193	bSkipFailures = TRUE;
	194	nGroupTransactions = 1; /* #2409 */
[40]	195	}
[1]	196
[395]	197	/* if exist, overwrite the data with the same name */
[40]	198	tmpMap=NULL;
[395]	199	tmpMap=getMapFromMaps(inputs,"overwrite","value");
	200	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	201	bOverwrite = TRUE;
	202	}
	203
	204	tmpMap=NULL;
[1]	205	tmpMap=getMapFromMaps(inputs,"append","value");
[395]	206	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"TRUE",4)==0){
	207	bAppend = TRUE;
[40]	208	}
[1]	209
[40]	210	tmpMap=NULL;
[1]	211	tmpMap=getMapFromMaps(inputs,"update","value");
[395]	212	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"TRUE",4)==0){
	213	bUpdate = TRUE;
[40]	214	}
[1]	215
[40]	216	tmpMap=NULL;
[1]	217	tmpMap=getMapFromMaps(inputs,"fid","value");
[395]	218	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	219	nFIDToFetch = atoi(tmpMap->value);
[40]	220	}
[1]	221
[40]	222	tmpMap=NULL;
[1]	223	tmpMap=getMapFromMaps(inputs,"sql","value");
[395]	224	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	225	pszSQLStatement = tmpMap->value;
[40]	226	}
[1]	227
[917]	228	tmpMap=getMapFromMaps(inputs,"dialect","value");
	229	if(tmpMap!=NULL){
	230	pszDialect=strdup(tmpMap->value);
	231	}
	232
[40]	233	tmpMap=NULL;
[1]	234	tmpMap=getMapFromMaps(inputs,"nln","value");
[395]	235	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	236	pszNewLayerName = tmpMap->value;
[40]	237	}
[1]	238
[40]	239	tmpMap=NULL;
[1]	240	tmpMap=getMapFromMaps(inputs,"nlt","value");
[395]	241	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	242	pszNewLayerName = tmpMap->value;
	243	if( EQUAL(tmpMap->value,"NONE") )
	244	eGType = wkbNone;
	245	else if( EQUAL(tmpMap->value,"GEOMETRY") )
	246	eGType = wkbUnknown;
	247	else if( EQUAL(tmpMap->value,"POINT") )
	248	eGType = wkbPoint;
	249	else if( EQUAL(tmpMap->value,"LINESTRING") )
	250	eGType = wkbLineString;
	251	else if( EQUAL(tmpMap->value,"POLYGON") )
	252	eGType = wkbPolygon;
	253	else if( EQUAL(tmpMap->value,"GEOMETRYCOLLECTION") )
	254	eGType = wkbGeometryCollection;
	255	else if( EQUAL(tmpMap->value,"MULTIPOINT") )
	256	eGType = wkbMultiPoint;
	257	else if( EQUAL(tmpMap->value,"MULTILINESTRING") )
	258	eGType = wkbMultiLineString;
	259	else if( EQUAL(tmpMap->value,"MULTIPOLYGON") )
	260	eGType = wkbMultiPolygon;
	261	else if( EQUAL(tmpMap->value,"GEOMETRY25D") )
	262	eGType = wkbUnknown \| wkb25DBit;
	263	else if( EQUAL(tmpMap->value,"POINT25D") )
	264	eGType = wkbPoint25D;
	265	else if( EQUAL(tmpMap->value,"LINESTRING25D") )
	266	eGType = wkbLineString25D;
	267	else if( EQUAL(tmpMap->value,"POLYGON25D") )
	268	eGType = wkbPolygon25D;
	269	else if( EQUAL(tmpMap->value,"GEOMETRYCOLLECTION25D") )
	270	eGType = wkbGeometryCollection25D;
	271	else if( EQUAL(tmpMap->value,"MULTIPOINT25D") )
	272	eGType = wkbMultiPoint25D;
	273	else if( EQUAL(tmpMap->value,"MULTILINESTRING25D") )
	274	eGType = wkbMultiLineString25D;
	275	else if( EQUAL(tmpMap->value,"MULTIPOLYGON25D") )
	276	eGType = wkbMultiPolygon25D;
	277	else
	278	{
	279	fprintf( stderr, "-nlt %s: type not recognised.\n",
	280	tmpMap->value );
	281	exit( 1 );
	282	}
[40]	283	}
[1]	284
[40]	285	tmpMap=NULL;
[1]	286	tmpMap=getMapFromMaps(inputs,"tg","value");
[395]	287	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	288	nGroupTransactions = atoi(tmpMap->value);
[40]	289	}
[1]	290
[40]	291	tmpMap=NULL;
[1]	292	tmpMap=getMapFromMaps(inputs,"s_srs","value");
[395]	293	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	294	pszSourceSRSDef = strdup(tmpMap->value);
[40]	295	}
[1]	296
[40]	297	tmpMap=NULL;
[1]	298	tmpMap=getMapFromMaps(inputs,"a_srs","value");
[395]	299	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	300	pszOutputSRSDef = strdup(tmpMap->value);
[40]	301	}
[1]	302
[40]	303	tmpMap=NULL;
[1]	304	tmpMap=getMapFromMaps(inputs,"t_srs","value");
[395]	305	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	306	pszOutputSRSDef = strdup(tmpMap->value);
	307	bTransform = TRUE;
[40]	308	}
[1]	309
[40]	310	tmpMap=NULL;
[1]	311	tmpMap=getMapFromMaps(inputs,"SPAT","value");
[395]	312	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	313	char *tmp=tmpMap->value;
	314	char *t=strtok(tmp,",");
	315	int cnt=0;
	316	double dfULX, dfULY, dfLRX, dfLRY;
	317	while(t!=NULL){
	318	switch(cnt){
	319	case 0:
	320	dfULX = atof(t);
	321	break;
	322	case 1:
	323	dfULY = atof(t);
	324	break;
	325	case 2:
	326	dfLRX = atof(t);
	327	break;
	328	case 3:
	329	dfLRY = atof(t);
	330	break;
	331	}
	332	fprintf(stderr,"%s\n\n",t);
	333	fprintf(stderr,"%f - %f - %f - %f\n\n",dfULX,dfULY,dfLRX,dfLRY);
	334	t=strtok(NULL,",");
	335	cnt++;
	336	}
	337
[40]	338	OGRLinearRing oRing;
	339
	340	oRing.addPoint( dfULX, dfULY );
	341	oRing.addPoint( dfULX, dfLRY );
	342	oRing.addPoint( dfLRX, dfLRY );
	343	oRing.addPoint( dfLRX, dfULY );
	344	oRing.addPoint( dfULX, dfULY );
	345	poSpatialFilter = new OGRPolygon();
	346	((OGRPolygon *) poSpatialFilter)->addRing( &oRing );
	347	}
[1]	348
[40]	349	tmpMap=NULL;
[1]	350	tmpMap=getMapFromMaps(inputs,"where","value");
[395]	351	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	352	pszWHERE = tmpMap->value;
[40]	353	}
[1]	354
[40]	355	tmpMap=NULL;
[1]	356	tmpMap=getMapFromMaps(inputs,"select","value");
[395]	357	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	358	pszSelect = tmpMap->value;
	359	papszSelFields = CSLTokenizeStringComplex(pszSelect, " ,",
	360	FALSE, FALSE );
[40]	361	}
[1]	362
[40]	363	tmpMap=NULL;
[1]	364	tmpMap=getMapFromMaps(inputs,"segmentize","value");
[395]	365	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
[1]	366	dfMaxSegmentLength = atof(tmpMap->value);
[40]	367	}
[917]	368
[1]	369	tmpMap=NULL;
	370	tmpMap=getMapFromMaps(inputs,"InputDSN","value");
[395]	371	if(tmpMap!=NULL && strncasecmp(tmpMap->value,"NULL",4)!=0){
	372	if(strncasecmp(tmpMap->value,"PG",2)!=0 &&
	373	strncasecmp(tmpMap->value,"My",2)!=0 &&
	374	strncasecmp(tmpMap->value,"OCI",3)!=0 ){
	375	if(strncasecmp(dataPath,tmpMap->value,strlen(dataPath))==0 \|\|
	376	strncasecmp(tempPath,tmpMap->value,strlen(tempPath))==0){
	377	pszDataSource=strdup(tmpMap->value);
	378	}else{
	379	pszDataSource=(char)malloc(sizeof(char)(strlen(dataPath)+strlen(tmpMap->value)+2));
	380	sprintf((char*)pszDataSource,"%s/%s",dataPath,tmpMap->value);
	381	}
	382	}else{
	383	pszDataSource=(char)malloc(sizeof(char)(strlen(tmpMap->value)+1));
	384	sprintf((char*)pszDataSource,"%s",tmpMap->value);
	385	}
[1]	386	}
	387
	388	tmpMap=NULL;
	389	tmpMap=getMapFromMaps(inputs,"OutputDSN","value");
	390	if(tmpMap!=NULL){
[395]	391	if(strncasecmp(tmpMap->value,"PG",2)!=0 &&
	392	strncasecmp(tmpMap->value,"MY",2)!=0 &&
	393	strncasecmp(tmpMap->value,"OCI",3)!=0){
	394	pszDestDataSource=(char)malloc(sizeof(char)(strlen(tempPath)+strlen(tmpMap->value)+4));
	395	sprintf((char)pszDestDataSource,"%s/%s",tempPath,tmpMap->value/,ext*/);
	396	pszwebDestData=(char)malloc(sizeof(char)(strlen(serverAddress)+strlen(tmpurl)+strlen(tmpMap->value)+4));
	397	sprintf((char)pszwebDestData,"%s%s/%s",serverAddress,tmpurl,tmpMap->value/,ext*/);
	398	}
	399	else{
	400	pszDestDataSource=(char)malloc(sizeof(char)(strlen(tmpMap->value)+1));
	401	sprintf((char)pszDestDataSource,"%s",tmpMap->value/,ext*/);
	402	pszwebDestData=(char)malloc(sizeof(char)(strlen(serverAddress)+strlen(tmpurl)+strlen(tmpMap->value)+4));
	403	sprintf((char)pszwebDestData,"%s%s/%s",serverAddress,tmpurl,tmpMap->value/,ext*/);
	404	}
[1]	405	}
	406
[395]	407	fprintf(stderr,"Message %s\n",pszDestDataSource);
[1]	408	#else
	409	/* -------------------------------------------------------------------- */
	410	/* Processing command line arguments. */
	411	/* -------------------------------------------------------------------- */
	412	nArgc = OGRGeneralCmdLineProcessor( nArgc, &papszArgv, 0 );
	413
	414	if( nArgc < 1 )
	415	exit( -nArgc );
	416
	417	for( int iArg = 1; iArg < nArgc; iArg++ )
	418	{
	419	if( EQUAL(papszArgv[iArg], "--utility_version") )
	420	{
	421	printf("%s was compiled against GDAL %s and is running against GDAL %s\n",
	422	papszArgv[0], GDAL_RELEASE_NAME, GDALVersionInfo("RELEASE_NAME"));
	423	return 0;
	424	}
	425	else if( EQUAL(papszArgv[iArg],"-f") && iArg < nArgc-1 )
	426	{
	427	pszFormat = papszArgv[++iArg];
	428	}
	429	else if( EQUAL(papszArgv[iArg],"-dsco") && iArg < nArgc-1 )
	430	{
	431	papszDSCO = CSLAddString(papszDSCO, papszArgv[++iArg] );
	432	}
	433	else if( EQUAL(papszArgv[iArg],"-lco") && iArg < nArgc-1 )
	434	{
	435	papszLCO = CSLAddString(papszLCO, papszArgv[++iArg] );
	436	}
	437	else if( EQUAL(papszArgv[iArg],"-preserve_fid") )
	438	{
	439	bPreserveFID = TRUE;
	440	}
	441	else if( EQUALN(papszArgv[iArg],"-skip",5) )
	442	{
	443	bSkipFailures = TRUE;
	444	nGroupTransactions = 1; /* #2409 */
	445	}
	446	else if( EQUAL(papszArgv[iArg],"-append") )
	447	{
	448	bAppend = TRUE;
	449	}
	450	else if( EQUAL(papszArgv[iArg],"-overwrite") )
	451	{
	452	bOverwrite = TRUE;
	453	}
	454	else if( EQUAL(papszArgv[iArg],"-update") )
	455	{
	456	bUpdate = TRUE;
	457	}
	458	else if( EQUAL(papszArgv[iArg],"-fid") && papszArgv[iArg+1] != NULL )
	459	{
	460	nFIDToFetch = atoi(papszArgv[++iArg]);
	461	}
	462	else if( EQUAL(papszArgv[iArg],"-sql") && papszArgv[iArg+1] != NULL )
	463	{
	464	pszSQLStatement = papszArgv[++iArg];
	465	}
	466	else if( EQUAL(papszArgv[iArg],"-nln") && iArg < nArgc-1 )
	467	{
	468	pszNewLayerName = papszArgv[++iArg];
	469	}
	470	else if( EQUAL(papszArgv[iArg],"-nlt") && iArg < nArgc-1 )
	471	{
	472	if( EQUAL(papszArgv[iArg+1],"NONE") )
	473	eGType = wkbNone;
	474	else if( EQUAL(papszArgv[iArg+1],"GEOMETRY") )
	475	eGType = wkbUnknown;
	476	else if( EQUAL(papszArgv[iArg+1],"POINT") )
	477	eGType = wkbPoint;
	478	else if( EQUAL(papszArgv[iArg+1],"LINESTRING") )
	479	eGType = wkbLineString;
	480	else if( EQUAL(papszArgv[iArg+1],"POLYGON") )
	481	eGType = wkbPolygon;
	482	else if( EQUAL(papszArgv[iArg+1],"GEOMETRYCOLLECTION") )
	483	eGType = wkbGeometryCollection;
	484	else if( EQUAL(papszArgv[iArg+1],"MULTIPOINT") )
	485	eGType = wkbMultiPoint;
	486	else if( EQUAL(papszArgv[iArg+1],"MULTILINESTRING") )
	487	eGType = wkbMultiLineString;
	488	else if( EQUAL(papszArgv[iArg+1],"MULTIPOLYGON") )
	489	eGType = wkbMultiPolygon;
	490	else if( EQUAL(papszArgv[iArg+1],"GEOMETRY25D") )
	491	eGType = wkbUnknown \| wkb25DBit;
	492	else if( EQUAL(papszArgv[iArg+1],"POINT25D") )
	493	eGType = wkbPoint25D;
	494	else if( EQUAL(papszArgv[iArg+1],"LINESTRING25D") )
	495	eGType = wkbLineString25D;
	496	else if( EQUAL(papszArgv[iArg+1],"POLYGON25D") )
	497	eGType = wkbPolygon25D;
	498	else if( EQUAL(papszArgv[iArg+1],"GEOMETRYCOLLECTION25D") )
	499	eGType = wkbGeometryCollection25D;
	500	else if( EQUAL(papszArgv[iArg+1],"MULTIPOINT25D") )
	501	eGType = wkbMultiPoint25D;
	502	else if( EQUAL(papszArgv[iArg+1],"MULTILINESTRING25D") )
	503	eGType = wkbMultiLineString25D;
	504	else if( EQUAL(papszArgv[iArg+1],"MULTIPOLYGON25D") )
	505	eGType = wkbMultiPolygon25D;
	506	else
	507	{
	508	fprintf( stderr, "-nlt %s: type not recognised.\n",
	509	papszArgv[iArg+1] );
	510	exit( 1 );
	511	}
	512	iArg++;
	513	}
	514	else if( (EQUAL(papszArgv[iArg],"-tg") \|\|
	515	EQUAL(papszArgv[iArg],"-gt")) && iArg < nArgc-1 )
	516	{
	517	nGroupTransactions = atoi(papszArgv[++iArg]);
	518	}
	519	else if( EQUAL(papszArgv[iArg],"-s_srs") && iArg < nArgc-1 )
	520	{
	521	pszSourceSRSDef = papszArgv[++iArg];
	522	}
	523	else if( EQUAL(papszArgv[iArg],"-a_srs") && iArg < nArgc-1 )
	524	{
	525	pszOutputSRSDef = papszArgv[++iArg];
	526	}
	527	else if( EQUAL(papszArgv[iArg],"-t_srs") && iArg < nArgc-1 )
	528	{
	529	pszOutputSRSDef = papszArgv[++iArg];
	530	bTransform = TRUE;
	531	}
	532	else if( EQUAL(papszArgv[iArg],"-spat")
	533	&& papszArgv[iArg+1] != NULL
	534	&& papszArgv[iArg+2] != NULL
	535	&& papszArgv[iArg+3] != NULL
	536	&& papszArgv[iArg+4] != NULL )
	537	{
	538	OGRLinearRing oRing;
	539
	540	oRing.addPoint( atof(papszArgv[iArg+1]), atof(papszArgv[iArg+2]) );
	541	oRing.addPoint( atof(papszArgv[iArg+1]), atof(papszArgv[iArg+4]) );
	542	oRing.addPoint( atof(papszArgv[iArg+3]), atof(papszArgv[iArg+4]) );
	543	oRing.addPoint( atof(papszArgv[iArg+3]), atof(papszArgv[iArg+2]) );
	544	oRing.addPoint( atof(papszArgv[iArg+1]), atof(papszArgv[iArg+2]) );
	545
	546	poSpatialFilter = new OGRPolygon();
	547	((OGRPolygon *) poSpatialFilter)->addRing( &oRing );
	548	iArg += 4;
	549	}
	550	else if( EQUAL(papszArgv[iArg],"-where") && papszArgv[iArg+1] != NULL )
	551	{
	552	pszWHERE = papszArgv[++iArg];
	553	}
	554	else if( EQUAL(papszArgv[iArg],"-select") && papszArgv[iArg+1] != NULL)
	555	{
	556	pszSelect = papszArgv[++iArg];
	557	papszSelFields = CSLTokenizeStringComplex(pszSelect, " ,",
	558	FALSE, FALSE );
	559	}
	560	else if( EQUAL(papszArgv[iArg],"-segmentize") && iArg < nArgc-1 )
	561	{
	562	dfMaxSegmentLength = atof(papszArgv[++iArg]);
	563	}
	564	else if( papszArgv[iArg][0] == '-' )
	565	{
	566	Usage();
	567	}
	568	else if( pszDestDataSource == NULL )
	569	pszDestDataSource = papszArgv[iArg];
	570	else if( pszDataSource == NULL )
	571	pszDataSource = papszArgv[iArg];
	572	else
	573	papszLayers = CSLAddString( papszLayers, papszArgv[iArg] );
	574	}
	575	#endif
	576
	577	if( pszDataSource == NULL )
	578	#ifdef ZOO_SERVICE
	579	{
[395]	580	#else
	581	Usage();
[1]	582	#endif
	583	#ifdef ZOO_SERVICE
[39]	584	setMapInMaps(conf,"lenv","message","Wrong parameter");
	585	return SERVICE_FAILED;
[1]	586	}
	587	#endif
	588
	589	/* -------------------------------------------------------------------- */
	590	/* Open data source. */
	591	/* -------------------------------------------------------------------- */
[766]	592	#if GDAL_VERSION_MAJOR >= 2
[790]	593	GDALDataset *poDS
	594	= (GDALDataset*) GDALOpenEx( pszDataSource,
	595	GDAL_OF_READONLY \| GDAL_OF_VECTOR,
	596	NULL, NULL, NULL );
[766]	597	GDALDataset *poODS;
	598	GDALDriverManager* poR=GetGDALDriverManager();
	599	GDALDriver *poDriver = NULL;
	600	#else
[790]	601	OGRDataSource* poDS
	602	= OGRSFDriverRegistrar::Open( pszDataSource, FALSE );
[766]	603	OGRDataSource *poODS;
	604	OGRSFDriverRegistrar *poR = OGRSFDriverRegistrar::GetRegistrar();
	605	OGRSFDriver *poDriver = NULL;
	606	#endif
[1]	607
	608	/* -------------------------------------------------------------------- */
	609	/* Report failure */
	610	/* -------------------------------------------------------------------- */
	611	if( poDS == NULL )
	612	{
	613	fprintf( stderr, "FAILURE:\n"
	614	"Unable to open datasource `%s' with the following drivers.\n",
	615	pszDataSource );
	616
	617	for( int iDriver = 0; iDriver < poR->GetDriverCount(); iDriver++ )
	618	{
[766]	619	#if GDAL_VERSION_MAJOR >= 2
	620	fprintf( stderr, " -> %s\n", poR->GetDriver(iDriver)->GetDescription() );
	621	#else
	622	fprintf( stderr, " -> %s\n", poR->GetDriver(iDriver)->GetName() );
	623	#endif
[1]	624	}
	625	#ifdef ZOO_SERVICE
[39]	626	char tmp[1024];
	627	sprintf(tmp,"Unable to open datasource `%s' with the following drivers.",pszDataSource);
	628	setMapInMaps(conf,"lenv","message",tmp);
	629	return SERVICE_FAILED;
[1]	630	#else
	631	exit( 1 );
	632	#endif
	633	}
	634
	635	/* -------------------------------------------------------------------- */
	636	/* Try opening the output datasource as an existing, writable */
	637	/* -------------------------------------------------------------------- */
	638	if( bUpdate )
	639	{
[766]	640	poODS =
	641	#if GDAL_VERSION_MAJOR >= 2
	642	(GDALDataset*) GDALOpenEx( pszDestDataSource,
	643	GDAL_OF_UPDATE \| GDAL_OF_VECTOR,
	644	NULL, NULL, NULL )
	645	#else
	646	OGRSFDriverRegistrar::Open( pszDestDataSource, TRUE )
	647	#endif
	648	;
[1]	649	if( poODS == NULL )
	650	{
	651	fprintf( stderr, "FAILURE:\n"
	652	"Unable to open existing output datasource `%s'.\n",
	653	pszDestDataSource );
	654	#ifdef ZOO_SERVICE
[39]	655	char tmp[1024];
	656	sprintf(tmp,"Unable to open existing output datasource `%s'.",pszDestDataSource);
	657	setMapInMaps(conf,"lenv","message",tmp);
	658	return SERVICE_FAILED;
[1]	659	#else
	660	exit( 1 );
	661	#endif
	662	}
	663
	664	if( CSLCount(papszDSCO) > 0 )
	665	{
	666	fprintf( stderr, "WARNING: Datasource creation options ignored since an existing datasource\n"
	667	" being updated.\n" );
	668	}
	669	}
	670
	671	/* -------------------------------------------------------------------- */
	672	/* Find the output driver. */
	673	/* -------------------------------------------------------------------- */
	674	else
	675	{
	676	int iDriver;
	677
	678	for( iDriver = 0;
	679	iDriver < poR->GetDriverCount() && poDriver == NULL;
	680	iDriver++ )
	681	{
[766]	682	#if GDAL_VERSION_MAJOR >=2
	683	if( EQUAL(poR->GetDriver(iDriver)->GetDescription(),pszFormat) )
	684	#else
	685	if( EQUAL(poR->GetDriver(iDriver)->GetName(),pszFormat) )
	686	#endif
[1]	687	{
	688	poDriver = poR->GetDriver(iDriver);
	689	}
	690	}
	691
	692	if( poDriver == NULL )
	693	{
	694	fprintf( stderr, "Unable to find driver `%s'.\n", pszFormat );
	695	fprintf( stderr, "The following drivers are available:\n" );
	696
	697	for( iDriver = 0; iDriver < poR->GetDriverCount(); iDriver++ )
	698	{
[766]	699	#if GDAL_VERSION_MAJOR >= 2
	700	fprintf( stderr, " -> %s\n", poR->GetDriver(iDriver)->GetDescription() );
	701	#else
	702	fprintf( stderr, " -> %s\n", poR->GetDriver(iDriver)->GetName() );
	703	#endif
[1]	704	}
	705	#ifdef ZOO_SERVICE
[39]	706	char tmp[1024];
	707	sprintf(tmp,"Unable to find driver `%s'.",pszFormat);
	708	setMapInMaps(conf,"lenv","message",tmp);
	709	return SERVICE_FAILED;
[1]	710	#else
	711	exit( 1 );
	712	#endif
	713	}
	714
[766]	715	#if GDAL_VERSION_MAJOR >=2
	716	if( !CPLTestBool( CSLFetchNameValueDef(poDriver->GetMetadata(), GDAL_DCAP_CREATE, "FALSE") ) )
	717	#else
	718	if( !poDriver->TestCapability( ODrCCreateDataSource ) )
	719	#endif
[1]	720	{
	721	fprintf( stderr, "%s driver does not support data source creation.\n",
	722	pszFormat );
	723	#ifdef ZOO_SERVICE
[39]	724	char tmp[1024];
	725	sprintf(tmp,"%s driver does not support data source creation.",pszFormat);
	726	setMapInMaps(conf,"lenv","message",tmp);
	727	return SERVICE_FAILED;
[1]	728	#else
	729	exit( 1 );
	730	#endif
	731	}
	732
	733	/* -------------------------------------------------------------------- */
	734	/* Create the output data source. */
	735	/* -------------------------------------------------------------------- */
[766]	736	#if GDAL_VERSION_MAJOR >=2
	737	poODS = poDriver->Create( pszDestDataSource, 0, 0, 0, GDT_Unknown, papszDSCO );
	738	#else
	739	poODS = poDriver->CreateDataSource( pszDestDataSource, papszDSCO );
	740	#endif
[1]	741	if( poODS == NULL )
	742	{
[790]	743	poODS =
	744	#if GDAL_VERSION_MAJOR >= 2
	745	(GDALDataset*) GDALOpenEx( pszDestDataSource,
	746	GDAL_OF_UPDATE \| GDAL_OF_VECTOR,
	747	NULL, NULL, NULL )
	748	#else
	749	OGRSFDriverRegistrar::Open( pszDestDataSource, TRUE )
	750	#endif
	751	;
	752	if( poODS == NULL )
	753	{
	754	fprintf( stderr, "FAILURE:\n"
	755	"Unable to open existing output datasource `%s'.\n",
	756	pszDestDataSource );
[1]	757	#ifdef ZOO_SERVICE
[790]	758	char tmp[1024];
	759	sprintf(tmp,"Unable to open existing output datasource `%s'.",pszDestDataSource);
	760	setMapInMaps(conf,"lenv","message",tmp);
	761	return SERVICE_FAILED;
[1]	762	#else
[790]	763	exit( 1 );
[1]	764	#endif
[790]	765
	766	fprintf( stderr, "%s driver failed to create %s\n",
	767	pszFormat, pszDestDataSource );
	768	}
[1]	769	}
	770	}
	771
	772	/* -------------------------------------------------------------------- */
	773	/* Parse the output SRS definition if possible. */
	774	/* -------------------------------------------------------------------- */
	775	if( pszOutputSRSDef != NULL )
	776	{
	777	poOutputSRS = new OGRSpatialReference();
	778	if( poOutputSRS->SetFromUserInput( pszOutputSRSDef ) != OGRERR_NONE )
	779	{
	780	fprintf( stderr, "Failed to process SRS definition: %s\n",
	781	pszOutputSRSDef );
	782	#ifdef ZOO_SERVICE
[39]	783	char tmp[1024];
	784	sprintf(tmp,"Failed to process SRS definition: %s",pszOutputSRSDef);
	785	setMapInMaps(conf,"lenv","message",tmp);
	786	return SERVICE_FAILED;
[1]	787	#else
	788	exit( 1 );
	789	#endif
	790	}
	791	}
	792
	793	/* -------------------------------------------------------------------- */
	794	/* Parse the source SRS definition if possible. */
	795	/* -------------------------------------------------------------------- */
	796	if( pszSourceSRSDef != NULL )
	797	{
	798	poSourceSRS = new OGRSpatialReference();
	799	if( poSourceSRS->SetFromUserInput( pszSourceSRSDef ) != OGRERR_NONE )
	800	{
	801	fprintf( stderr, "Failed to process SRS definition: %s\n",
	802	pszSourceSRSDef );
	803	#ifdef ZOO_SERVICE
[39]	804	char tmp[1024];
	805	sprintf(tmp,"Failed to process SRS definition: %s",pszOutputSRSDef);
	806	setMapInMaps(conf,"lenv","message",tmp);
	807	return SERVICE_FAILED;
[1]	808	#else
	809	exit( 1 );
	810	#endif
	811	}
	812	}
	813
	814	/* -------------------------------------------------------------------- */
	815	/* Special case for -sql clause. No source layers required. */
	816	/* -------------------------------------------------------------------- */
	817	if( pszSQLStatement != NULL )
	818	{
	819	OGRLayer *poResultSet;
	820
	821	if( pszWHERE != NULL )
	822	fprintf( stderr, "-where clause ignored in combination with -sql.\n" );
	823	if( CSLCount(papszLayers) > 0 )
[395]	824	fprintf( stderr, "layer names ignored in combination with -sql. (%s)\n", pszSQLStatement);
[1]	825
	826	poResultSet = poDS->ExecuteSQL( pszSQLStatement, poSpatialFilter,
[917]	827	pszDialect );
[1]	828
	829	if( poResultSet != NULL )
	830	{
	831	if( !TranslateLayer( poDS, poResultSet, poODS, papszLCO,
	832	pszNewLayerName, bTransform, poOutputSRS,
	833	poSourceSRS, papszSelFields, bAppend, eGType,
	834	bOverwrite, dfMaxSegmentLength ))
	835	{
	836	CPLError( CE_Failure, CPLE_AppDefined,
	837	"Terminating translation prematurely after failed\n"
	838	"translation from sql statement." );
	839
	840	exit( 1 );
	841	}
	842	poDS->ReleaseResultSet( poResultSet );
[917]	843	}else{
	844	setMapInMaps(conf,"lenv","message","There was an error when running yoru SQL query.");
	845	if(pszDialect!=NULL)
	846	free(pszDialect);
	847	return SERVICE_FAILED;
	848	}
[1]	849	}
[917]	850	if(pszDialect!=NULL)
	851	free(pszDialect);
[1]	852
	853	/* -------------------------------------------------------------------- */
	854	/* Process each data source layer. */
	855	/* -------------------------------------------------------------------- */
	856	for( int iLayer = 0;
	857	pszSQLStatement == NULL && iLayer < poDS->GetLayerCount();
	858	iLayer++ )
	859	{
	860	OGRLayer *poLayer = poDS->GetLayer(iLayer);
	861
	862	if( poLayer == NULL )
	863	{
	864	fprintf( stderr, "FAILURE: Couldn't fetch advertised layer %d!\n",
	865	iLayer );
	866	#ifdef ZOO_SERVICE
[39]	867	char tmp[1024];
	868	sprintf(tmp,"Couldn't fetch advertised layer %d!",iLayer);
	869	setMapInMaps(conf,"lenv","message",tmp);
	870	return SERVICE_FAILED;
[1]	871	#else
[39]	872	exit( 1 );
[1]	873	#endif
	874	}
	875
	876	if( CSLCount(papszLayers) == 0
	877	\|\| CSLFindString( papszLayers,
	878	poLayer->GetLayerDefn()->GetName() ) != -1 )
	879	{
	880	if( pszWHERE != NULL )
	881	poLayer->SetAttributeFilter( pszWHERE );
	882
	883	if( poSpatialFilter != NULL )
	884	poLayer->SetSpatialFilter( poSpatialFilter );
	885
	886	if( !TranslateLayer( poDS, poLayer, poODS, papszLCO,
	887	pszNewLayerName, bTransform, poOutputSRS,
	888	poSourceSRS, papszSelFields, bAppend, eGType,
	889	bOverwrite, dfMaxSegmentLength )
	890	&& !bSkipFailures )
	891	{
	892	#ifdef ZOO_SERVICE
[39]	893	char tmp[1024];
	894	sprintf(tmp,"Terminating translation prematurely after failed of layer %s",poLayer->GetLayerDefn()->GetName() );
	895	setMapInMaps(conf,"lenv","message",tmp);
	896	return SERVICE_FAILED;
[1]	897	#else
[395]	898	CPLError( CE_Failure, CPLE_AppDefined,
	899	"Terminating translation prematurely after failed\n"
	900	"translation of layer %s (use -skipfailures to skip errors)\n",
	901	poLayer->GetLayerDefn()->GetName() );
	902
[1]	903	exit( 1 );
	904	#endif
	905	}
	906	}
	907	}
	908
[395]	909	#ifdef ZOO_SERVICE
[917]	910	setMapInMaps(outputs,"OutputedDataSourceName","value",(char*)pszwebDestData);
	911	//outputs->content=createMap("value",(char*)pszwebDestData);
[395]	912	#endif
	913
[1]	914	/* -------------------------------------------------------------------- */
	915	/* Close down. */
	916	/* -------------------------------------------------------------------- */
[395]	917	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	918	delete poOutputSRS;
[395]	919	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	920	delete poSourceSRS;
[395]	921	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	922	delete poODS;
[395]	923	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	924	delete poDS;
[395]	925	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	926
[395]	927	#ifndef ZOO_SERVICE
[1]	928	CSLDestroy(papszSelFields);
[395]	929	CSLDestroy( papszArgv );
[1]	930	#endif
[395]	931	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	932	CSLDestroy( papszLayers );
[395]	933	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
	934	#ifndef ZOO_SERVICE
[1]	935	CSLDestroy( papszDSCO );
	936	CSLDestroy( papszLCO );
[395]	937	#endif
	938	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	939
[917]	940	//OGRCleanupAll();
[395]	941	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	942
	943	#ifdef DBMALLOC
	944	malloc_dump(1);
	945	#endif
[395]	946	fprintf(stderr,"%s %d\n",__FILE__,__LINE__);
[1]	947
	948	#ifdef ZOO_SERVICE
[917]	949	//sleep(10);
[1]	950	return SERVICE_SUCCEEDED;
	951	#else
[395]	952	return 0;
[1]	953	#endif
	954	}
	955
	956	/************************************************************************/
	957	/* Usage() */
	958	/************************************************************************/
	959
	960	static void Usage()
	961
	962	{
[766]	963	#if GDAL_VERSION_MAJOR >= 2
	964	GDALDataset *poDS;
	965	GDALDataset *poODS;
	966	GDALDriverManager* poR=GetGDALDriverManager();
	967	GDALDriver *poDriver = NULL;
	968	#else
	969	OGRDataSource* poDS;
	970	OGRDataSource *poODS;
	971	OGRSFDriverRegistrar *poR = OGRSFDriverRegistrar::GetRegistrar();
	972	OGRSFDriver *poDriver = NULL;
	973	#endif
[1]	974
	975	#ifdef ZOO_SERVICE
	976	fprintf(stderr,
	977	#else
	978	printf(
	979	#endif
	980	"Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update] [-gt n]\n"
	981	" [-select field_list] [-where restricted_where] \n"
	982	" [-sql <sql statement>] \n"
	983	" [-spat xmin ymin xmax ymax] [-preserve_fid] [-fid FID]\n"
	984	" [-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def]\n"
	985	" [-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]\n"
	986	" [-segmentize max_dist]\n"
	987	" dst_datasource_name src_datasource_name\n"
	988	" [-lco NAME=VALUE] [-nln name] [-nlt type] [layer [layer ...]]\n"
	989	"\n"
	990	" -f format_name: output file format name, possible values are:\n");
	991
	992	for( int iDriver = 0; iDriver < poR->GetDriverCount(); iDriver++ )
	993	{
[766]	994	poDriver = poR->GetDriver(iDriver);
[1]	995
[766]	996	#if GDAL_VERSION_MAJOR >=2
	997	if( !CPLTestBool( CSLFetchNameValueDef(poDriver->GetMetadata(), GDAL_DCAP_CREATE, "FALSE") ) )
	998	#else
	999	if( !poDriver->TestCapability( ODrCCreateDataSource ) )
	1000	#endif
	1001	#if GDAL_VERSION_MAJOR >= 2
	1002	fprintf( stderr, " -> %s\n", poDriver->GetDescription() );
	1003	#else
	1004	fprintf( stderr, " -> %s\n", poDriver->GetName() );
	1005	#endif
[1]	1006	}
	1007
	1008	#ifdef ZOO_SERVICE
	1009	fprintf(stderr,
	1010	#else
	1011	printf(
	1012	#endif
	1013	" -append: Append to existing layer instead of creating new if it exists\n"
	1014	" -overwrite: delete the output layer and recreate it empty\n"
	1015	" -update: Open existing output datasource in update mode\n"
	1016	" -select field_list: Comma-delimited list of fields from input layer to\n"
	1017	" copy to the new layer (defaults to all)\n"
	1018	" -where restricted_where: Attribute query (like SQL WHERE)\n"
	1019	" -sql statement: Execute given SQL statement and save result.\n"
	1020	" -skipfailures: skip features or layers that fail to convert\n"
	1021	" -gt n: group n features per transaction (default 200)\n"
	1022	" -spat xmin ymin xmax ymax: spatial query extents\n"
	1023	" -segmentize max_dist: maximum distance between 2 nodes.\n"
	1024	" Used to create intermediate points\n"
	1025	" -dsco NAME=VALUE: Dataset creation option (format specific)\n"
	1026	" -lco NAME=VALUE: Layer creation option (format specific)\n"
	1027	" -nln name: Assign an alternate name to the new layer\n"
	1028	" -nlt type: Force a geometry type for new layer. One of NONE, GEOMETRY,\n"
	1029	" POINT, LINESTRING, POLYGON, GEOMETRYCOLLECTION, MULTIPOINT,\n"
	1030	" MULTIPOLYGON, or MULTILINESTRING. Add \"25D\" for 3D layers.\n"
	1031	" Default is type of source layer.\n" );
	1032
	1033	#ifdef ZOO_SERVICE
	1034	fprintf(stderr,
	1035	#else
	1036	printf(
	1037	#endif
	1038	" -a_srs srs_def: Assign an output SRS\n"
	1039	" -t_srs srs_def: Reproject/transform to this SRS on output\n"
	1040	" -s_srs srs_def: Override source SRS\n"
	1041	"\n"
	1042	" Srs_def can be a full WKT definition (hard to escape properly),\n"
	1043	" or a well known definition (ie. EPSG:4326) or a file with a WKT\n"
	1044	" definition.\n" );
	1045
	1046
	1047	#ifndef ZOO_SERVICE
	1048	exit( 1 );
	1049	#endif
	1050	}
	1051
	1052	/************************************************************************/
	1053	/* TranslateLayer() */
	1054	/************************************************************************/
	1055
[766]	1056	static int TranslateLayer(
	1057	#if GDAL_VERSION_MAJOR >= 2
	1058	GDALDataset
	1059	#else
	1060	OGRDataSource
	1061	#endif
	1062	*poSrcDS,
	1063	OGRLayer * poSrcLayer,
	1064	#if GDAL_VERSION_MAJOR >= 2
	1065	GDALDataset
	1066	#else
	1067	OGRDataSource
	1068	#endif
	1069	*poDstDS,
	1070	char **papszLCO,
	1071	const char *pszNewLayerName,
	1072	int bTransform,
	1073	OGRSpatialReference *poOutputSRS,
	1074	OGRSpatialReference *poSourceSRS,
	1075	char **papszSelFields,
	1076	int bAppend, int eGType, int bOverwrite,
	1077	double dfMaxSegmentLength)
[40]	1078
[1]	1079	{
	1080	OGRLayer *poDstLayer;
	1081	OGRFeatureDefn *poFDefn;
	1082	OGRErr eErr;
	1083	int bForceToPolygon = FALSE;
	1084	int bForceToMultiPolygon = FALSE;
	1085
	1086	if( pszNewLayerName == NULL )
	1087	pszNewLayerName = poSrcLayer->GetLayerDefn()->GetName();
	1088
	1089	if( wkbFlatten(eGType) == wkbPolygon )
	1090	bForceToPolygon = TRUE;
	1091	else if( wkbFlatten(eGType) == wkbMultiPolygon )
	1092	bForceToMultiPolygon = TRUE;
	1093
	1094	/* -------------------------------------------------------------------- */
	1095	/* Setup coordinate transformation if we need it. */
	1096	/* -------------------------------------------------------------------- */
	1097	OGRCoordinateTransformation *poCT = NULL;
	1098
	1099	if( bTransform )
	1100	{
	1101	if( poSourceSRS == NULL )
	1102	poSourceSRS = poSrcLayer->GetSpatialRef();
	1103
	1104	if( poSourceSRS == NULL )
	1105	{
	1106	fprintf( stderr, "Can't transform coordinates, source layer has no\n"
	1107	"coordinate system. Use -s_srs to set one.\n" );
	1108	#ifdef ZOO_SERVICE
	1109	return SERVICE_FAILED;
	1110	#else
	1111	exit( 1 );
	1112	#endif
	1113	}
	1114
	1115	CPLAssert( NULL != poSourceSRS );
	1116	CPLAssert( NULL != poOutputSRS );
	1117
	1118	poCT = OGRCreateCoordinateTransformation( poSourceSRS, poOutputSRS );
	1119	if( poCT == NULL )
	1120	{
	1121	char *pszWKT = NULL;
	1122
	1123	fprintf( stderr, "Failed to create coordinate transformation between the\n"
	1124	"following coordinate systems. This may be because they\n"
	1125	"are not transformable, or because projection services\n"
	1126	"(PROJ.4 DLL/.so) could not be loaded.\n" );
	1127
	1128	poSourceSRS->exportToPrettyWkt( &pszWKT, FALSE );
	1129	fprintf( stderr, "Source:\n%s\n", pszWKT );
	1130
	1131	poOutputSRS->exportToPrettyWkt( &pszWKT, FALSE );
	1132	fprintf( stderr, "Target:\n%s\n", pszWKT );
	1133	#ifdef ZOO_SERVICE
	1134	return SERVICE_FAILED;
	1135	#else
	1136	exit( 1 );
	1137	#endif
	1138	}
	1139	}
	1140
	1141	/* -------------------------------------------------------------------- */
	1142	/* Get other info. */
	1143	/* -------------------------------------------------------------------- */
	1144	poFDefn = poSrcLayer->GetLayerDefn();
	1145
	1146	if( poOutputSRS == NULL )
	1147	poOutputSRS = poSrcLayer->GetSpatialRef();
	1148
	1149	/* -------------------------------------------------------------------- */
	1150	/* Find the layer. */
	1151	/* -------------------------------------------------------------------- */
	1152	int iLayer = -1;
	1153	poDstLayer = NULL;
	1154
	1155	for( iLayer = 0; iLayer < poDstDS->GetLayerCount(); iLayer++ )
	1156	{
	1157	OGRLayer *poLayer = poDstDS->GetLayer(iLayer);
	1158
	1159	if( poLayer != NULL
	1160	&& EQUAL(poLayer->GetLayerDefn()->GetName(),pszNewLayerName) )
	1161	{
	1162	poDstLayer = poLayer;
	1163	break;
	1164	}
	1165	}
	1166
	1167	/* -------------------------------------------------------------------- */
	1168	/* If the user requested overwrite, and we have the layer in */
	1169	/* question we need to delete it now so it will get recreated */
	1170	/* (overwritten). */
	1171	/* -------------------------------------------------------------------- */
	1172	if( poDstLayer != NULL && bOverwrite )
	1173	{
	1174	if( poDstDS->DeleteLayer( iLayer ) != OGRERR_NONE )
	1175	{
	1176	fprintf( stderr,
	1177	"DeleteLayer() failed when overwrite requested.\n" );
	1178	return FALSE;
	1179	}
	1180	poDstLayer = NULL;
	1181	}
	1182
	1183	/* -------------------------------------------------------------------- */
	1184	/* If the layer does not exist, then create it. */
	1185	/* -------------------------------------------------------------------- */
	1186	if( poDstLayer == NULL )
	1187	{
	1188	if( eGType == -2 )
	1189	eGType = poFDefn->GetGeomType();
	1190
	1191	if( !poDstDS->TestCapability( ODsCCreateLayer ) )
	1192	{
	1193	fprintf( stderr,
	1194	"Layer %s not found, and CreateLayer not supported by driver.",
	1195	pszNewLayerName );
	1196	return FALSE;
	1197	}
	1198
	1199	CPLErrorReset();
	1200
	1201	poDstLayer = poDstDS->CreateLayer( pszNewLayerName, poOutputSRS,
	1202	(OGRwkbGeometryType) eGType,
	1203	papszLCO );
	1204
	1205	if( poDstLayer == NULL )
	1206	return FALSE;
	1207
	1208	bAppend = FALSE;
	1209	}
	1210
	1211	/* -------------------------------------------------------------------- */
	1212	/* Otherwise we will append to it, if append was requested. */
	1213	/* -------------------------------------------------------------------- */
	1214	else if( !bAppend )
	1215	{
	1216	fprintf( stderr, "FAILED: Layer %s already exists, and -append not specified.\n"
	1217	" Consider using -append, or -overwrite.\n",
	1218	pszNewLayerName );
	1219	return FALSE;
	1220	}
	1221	else
	1222	{
	1223	if( CSLCount(papszLCO) > 0 )
	1224	{
	1225	fprintf( stderr, "WARNING: Layer creation options ignored since an existing layer is\n"
	1226	" being appended to.\n" );
	1227	}
	1228	}
	1229
	1230	/* -------------------------------------------------------------------- */
	1231	/* Add fields. Default to copy all field. */
	1232	/* If only a subset of all fields requested, then output only */
	1233	/* the selected fields, and in the order that they were */
	1234	/* selected. */
	1235	/* -------------------------------------------------------------------- */
	1236	int iField;
[395]	1237	int hasMmField=-1;
[1]	1238
	1239	if (papszSelFields && !bAppend )
	1240	{
	1241	for( iField=0; papszSelFields[iField] != NULL; iField++)
	1242	{
	1243	int iSrcField = poFDefn->GetFieldIndex(papszSelFields[iField]);
[395]	1244	OGRFieldDefn *tmp=poFDefn->GetFieldDefn(iSrcField);
	1245	fprintf(stderr,"NAME: %s\n",tmp->GetNameRef());
	1246	fflush(stderr);
	1247	if(tmp!=NULL && strncasecmp(tmp->GetNameRef(),"MMID",4)==0)
	1248	hasMmField=1;
[1]	1249	if (iSrcField >= 0)
	1250	poDstLayer->CreateField( poFDefn->GetFieldDefn(iSrcField) );
	1251	else
	1252	{
	1253	fprintf( stderr, "Field '%s' not found in source layer.\n",
	1254	papszSelFields[iField] );
	1255	if( !bSkipFailures )
	1256	return FALSE;
	1257	}
	1258	}
[395]	1259
[1]	1260	}
	1261	else if( !bAppend )
	1262	{
[395]	1263	for( iField = 0; iField < poFDefn->GetFieldCount(); iField++ ){
[1]	1264	poDstLayer->CreateField( poFDefn->GetFieldDefn(iField) );
[395]	1265	OGRFieldDefn *tmp=poFDefn->GetFieldDefn(iField);
	1266	if(tmp!=NULL && strncasecmp(tmp->GetNameRef(),"MMID",4)==0)
	1267	hasMmField=1;
	1268	fprintf(stderr,"NAME: %s\n",tmp->GetNameRef());
	1269	}
[1]	1270	}
[395]	1271	/*if(hasMmField<0){
	1272	OGRFieldDefn oField( "MMID", OFTInteger );
	1273	poDstLayer->CreateField( &oField );
	1274	}*/
[1]	1275
	1276	/* -------------------------------------------------------------------- */
	1277	/* Transfer features. */
	1278	/* -------------------------------------------------------------------- */
	1279	OGRFeature *poFeature;
	1280	int nFeaturesInTransaction = 0;
[395]	1281	int fCount=0;
[1]	1282	poSrcLayer->ResetReading();
	1283
	1284	if( nGroupTransactions )
	1285	poDstLayer->StartTransaction();
	1286
	1287	while( TRUE )
	1288	{
	1289	OGRFeature *poDstFeature = NULL;
	1290
	1291	if( nFIDToFetch != OGRNullFID )
	1292	{
	1293	// Only fetch feature on first pass.
	1294	if( nFeaturesInTransaction == 0 )
	1295	poFeature = poSrcLayer->GetFeature(nFIDToFetch);
	1296	else
	1297	poFeature = NULL;
	1298	}
	1299	else
	1300	poFeature = poSrcLayer->GetNextFeature();
[395]	1301
[1]	1302	if( poFeature == NULL )
	1303	break;
	1304
[395]	1305	//poFeature->SetField((poFeature->GetFieldCount()-1),fCount);
	1306
[1]	1307	if( ++nFeaturesInTransaction == nGroupTransactions )
	1308	{
	1309	poDstLayer->CommitTransaction();
	1310	poDstLayer->StartTransaction();
	1311	nFeaturesInTransaction = 0;
	1312	}
	1313
	1314	CPLErrorReset();
	1315	poDstFeature = OGRFeature::CreateFeature( poDstLayer->GetLayerDefn() );
	1316
	1317	if( poDstFeature->SetFrom( poFeature, TRUE ) != OGRERR_NONE )
	1318	{
	1319	if( nGroupTransactions )
	1320	poDstLayer->CommitTransaction();
	1321
	1322	CPLError( CE_Failure, CPLE_AppDefined,
	1323	"Unable to translate feature %ld from layer %s.\n",
	1324	poFeature->GetFID(), poFDefn->GetName() );
	1325
	1326	OGRFeature::DestroyFeature( poFeature );
	1327	OGRFeature::DestroyFeature( poDstFeature );
	1328	return FALSE;
	1329	}
	1330
	1331	if( bPreserveFID )
	1332	poDstFeature->SetFID( poFeature->GetFID() );
	1333
[395]	1334	/*if(hasMmField<0){
	1335	poDstFeature->SetField((poDstFeature->GetFieldCount()-1),fCount);
	1336	fCount++;
	1337	}*/
	1338
[1]	1339	#ifndef GDAL_1_5_0
	1340	if (poDstFeature->GetGeometryRef() != NULL && dfMaxSegmentLength > 0)
	1341	poDstFeature->GetGeometryRef()->segmentize(dfMaxSegmentLength);
	1342	#endif
	1343
	1344	if( poCT && poDstFeature->GetGeometryRef() != NULL )
	1345	{
	1346	eErr = poDstFeature->GetGeometryRef()->transform( poCT );
	1347	if( eErr != OGRERR_NONE )
	1348	{
	1349	if( nGroupTransactions )
	1350	poDstLayer->CommitTransaction();
	1351
	1352	fprintf( stderr, "Failed to reproject feature %d (geometry probably out of source or destination SRS).\n",
	1353	(int) poFeature->GetFID() );
	1354	if( !bSkipFailures )
	1355	{
	1356	OGRFeature::DestroyFeature( poFeature );
	1357	OGRFeature::DestroyFeature( poDstFeature );
	1358	return FALSE;
	1359	}
	1360	}
	1361	}
	1362
	1363	if( poDstFeature->GetGeometryRef() != NULL && bForceToPolygon )
	1364	{
	1365	poDstFeature->SetGeometryDirectly(
	1366	OGRGeometryFactory::forceToPolygon(
	1367	poDstFeature->StealGeometry() ) );
	1368	}
	1369
	1370	if( poDstFeature->GetGeometryRef() != NULL && bForceToMultiPolygon )
	1371	{
	1372	poDstFeature->SetGeometryDirectly(
	1373	OGRGeometryFactory::forceToMultiPolygon(
	1374	poDstFeature->StealGeometry() ) );
	1375	}
	1376
	1377	OGRFeature::DestroyFeature( poFeature );
	1378
	1379	CPLErrorReset();
	1380	if( poDstLayer->CreateFeature( poDstFeature ) != OGRERR_NONE
	1381	&& !bSkipFailures )
	1382	{
	1383	if( nGroupTransactions )
	1384	poDstLayer->RollbackTransaction();
	1385
	1386	OGRFeature::DestroyFeature( poDstFeature );
	1387	return FALSE;
	1388	}
	1389
	1390	OGRFeature::DestroyFeature( poDstFeature );
	1391	}
	1392
	1393	if( nGroupTransactions )
	1394	poDstLayer->CommitTransaction();
	1395
	1396	/* -------------------------------------------------------------------- */
	1397	/* Cleaning */
	1398	/* -------------------------------------------------------------------- */
	1399	delete poCT;
	1400
	1401	return TRUE;
	1402	}
	1403
	1404	#ifdef ZOO_SERVICE
	1405	}
	1406	#endif

Note: See TracBrowser for help on using the repository browser.

Open WPS platform

source: trunk/zoo-project/zoo-services/ogr/ogr2ogr/service.c @ 928

Download in other formats:

Download now !

Support ZOO-Project

Search

Context Navigation

Main Navigation

ZOO Sponsors

Knowledge partners

Related links