gvsig-gdal / trunk / org.gvsig.gdal / org.gvsig.gdal.app / org.gvsig.gdal.app.mainplugin / org.gvsig.gdal.app.mainplugin.common / src / main / java / org / gvsig / gdal / app / mainplugin / common / ogr2ogr.java @ 151
History | View | Annotate | Download (70.3 KB)
1 |
/******************************************************************************
|
---|---|
2 |
* $Id$
|
3 |
*
|
4 |
* Project: OpenGIS Simple Features Reference Implementation
|
5 |
* Purpose: Java port of a simple client for translating between formats.
|
6 |
* Author: Even Rouault, <even dot rouault at mines dash paris dot org>
|
7 |
*
|
8 |
* Port from ogr2ogr.cpp by Frank Warmerdam
|
9 |
*
|
10 |
******************************************************************************
|
11 |
* Copyright (c) 2009, Even Rouault
|
12 |
* Copyright (c) 1999, Frank Warmerdam
|
13 |
*
|
14 |
* Permission is hereby granted, free of charge, to any person obtaining a
|
15 |
* copy of this software and associated documentation files (the "Software"),
|
16 |
* to deal in the Software without restriction, including without limitation
|
17 |
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
|
18 |
* and/or sell copies of the Software, and to permit persons to whom the
|
19 |
* Software is furnished to do so, subject to the following conditions:
|
20 |
*
|
21 |
* The above copyright notice and this permission notice shall be included
|
22 |
* in all copies or substantial portions of the Software.
|
23 |
*
|
24 |
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
|
25 |
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
26 |
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
|
27 |
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
28 |
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
|
29 |
* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
|
30 |
* DEALINGS IN THE SOFTWARE.
|
31 |
****************************************************************************/
|
32 |
package org.gvsig.gdal.app.mainplugin.common; |
33 |
|
34 |
import java.io.File; |
35 |
import java.util.Vector; |
36 |
import java.util.Enumeration; |
37 |
import java.util.StringTokenizer; |
38 |
|
39 |
import org.gdal.gdal.gdal; |
40 |
import org.gdal.gdal.ProgressCallback; |
41 |
import org.gdal.gdal.TermProgressCallback; |
42 |
import org.gdal.ogr.ogr; |
43 |
import org.gdal.ogr.ogrConstants; |
44 |
import org.gdal.ogr.Driver; |
45 |
import org.gdal.ogr.DataSource; |
46 |
import org.gdal.ogr.Layer; |
47 |
import org.gdal.ogr.Feature; |
48 |
import org.gdal.ogr.FeatureDefn; |
49 |
import org.gdal.ogr.FieldDefn; |
50 |
import org.gdal.ogr.Geometry; |
51 |
import org.gdal.osr.SpatialReference; |
52 |
import org.gdal.osr.CoordinateTransformation; |
53 |
|
54 |
/* Note : this is the most direct port of ogr2ogr.cpp possible */
|
55 |
/* It could be made much more java'ish ! */
|
56 |
|
57 |
class GDALScaledProgress extends ProgressCallback |
58 |
{ |
59 |
private double pctMin; |
60 |
private double pctMax; |
61 |
private ProgressCallback mainCbk;
|
62 |
|
63 |
public GDALScaledProgress(double pctMin, double pctMax, |
64 |
ProgressCallback mainCbk) |
65 |
{ |
66 |
this.pctMin = pctMin;
|
67 |
this.pctMax = pctMax;
|
68 |
this.mainCbk = mainCbk;
|
69 |
} |
70 |
|
71 |
public int run(double dfComplete, String message) |
72 |
{ |
73 |
return mainCbk.run(pctMin + dfComplete * (pctMax - pctMin), message);
|
74 |
} |
75 |
}; |
76 |
|
77 |
|
78 |
public class ogr2ogr |
79 |
{ |
80 |
static boolean bSkipFailures = false; |
81 |
static int nGroupTransactions = 200; |
82 |
static boolean bPreserveFID = false; |
83 |
static final int OGRNullFID = -1; |
84 |
static int nFIDToFetch = OGRNullFID; |
85 |
|
86 |
static class GeomOperation |
87 |
{ |
88 |
private GeomOperation() {}
|
89 |
public static GeomOperation NONE = new GeomOperation(); |
90 |
public static GeomOperation SEGMENTIZE = new GeomOperation(); |
91 |
public static GeomOperation SIMPLIFY_PRESERVE_TOPOLOGY = new GeomOperation(); |
92 |
} |
93 |
|
94 |
/************************************************************************/
|
95 |
/* main() */
|
96 |
/************************************************************************/
|
97 |
|
98 |
public static void main(String[] args) |
99 |
{ |
100 |
String pszFormat = "ESRI Shapefile"; |
101 |
String pszDataSource = null; |
102 |
String pszDestDataSource = null; |
103 |
Vector papszLayers = new Vector(); |
104 |
Vector papszDSCO = new Vector(), papszLCO = new Vector(); |
105 |
boolean bTransform = false; |
106 |
boolean bAppend = false, bUpdate = false, bOverwrite = false; |
107 |
String pszOutputSRSDef = null; |
108 |
String pszSourceSRSDef = null; |
109 |
SpatialReference poOutputSRS = null;
|
110 |
SpatialReference poSourceSRS = null;
|
111 |
String pszNewLayerName = null; |
112 |
String pszWHERE = null; |
113 |
Geometry poSpatialFilter = null;
|
114 |
String pszSelect;
|
115 |
Vector papszSelFields = null; |
116 |
String pszSQLStatement = null; |
117 |
int eGType = -2; |
118 |
GeomOperation eGeomOp = GeomOperation.NONE; |
119 |
double dfGeomOpParam = 0; |
120 |
Vector papszFieldTypesToString = new Vector(); |
121 |
boolean bDisplayProgress = false; |
122 |
ProgressCallback pfnProgress = null;
|
123 |
boolean bClipSrc = false; |
124 |
Geometry poClipSrc = null;
|
125 |
String pszClipSrcDS = null; |
126 |
String pszClipSrcSQL = null; |
127 |
String pszClipSrcLayer = null; |
128 |
String pszClipSrcWhere = null; |
129 |
Geometry poClipDst = null;
|
130 |
String pszClipDstDS = null; |
131 |
String pszClipDstSQL = null; |
132 |
String pszClipDstLayer = null; |
133 |
String pszClipDstWhere = null; |
134 |
String pszSrcEncoding = null; |
135 |
String pszDstEncoding = null; |
136 |
boolean bExplodeCollections = false; |
137 |
String pszZField = null; |
138 |
|
139 |
ogr.DontUseExceptions(); |
140 |
|
141 |
/* -------------------------------------------------------------------- */
|
142 |
/* Register format(s). */
|
143 |
/* -------------------------------------------------------------------- */
|
144 |
if( ogr.GetDriverCount() == 0 ) |
145 |
ogr.RegisterAll(); |
146 |
|
147 |
/* -------------------------------------------------------------------- */
|
148 |
/* Processing command line arguments. */
|
149 |
/* -------------------------------------------------------------------- */
|
150 |
args = ogr.GeneralCmdLineProcessor( args ); |
151 |
|
152 |
if( args.length < 2 ) |
153 |
{ |
154 |
Usage(); |
155 |
System.exit( -1 ); |
156 |
} |
157 |
|
158 |
for( int iArg = 0; iArg < args.length; iArg++ ) |
159 |
{ |
160 |
if( args[iArg].equalsIgnoreCase("-f") && iArg < args.length-1 ) |
161 |
{ |
162 |
pszFormat = args[++iArg]; |
163 |
} |
164 |
else if( args[iArg].equalsIgnoreCase("-dsco") && iArg < args.length-1 ) |
165 |
{ |
166 |
papszDSCO.addElement(args[++iArg] ); |
167 |
} |
168 |
else if( args[iArg].equalsIgnoreCase("-lco") && iArg < args.length-1 ) |
169 |
{ |
170 |
papszLCO.addElement(args[++iArg] ); |
171 |
} |
172 |
else if( args[iArg].equalsIgnoreCase("-preserve_fid") ) |
173 |
{ |
174 |
bPreserveFID = true;
|
175 |
} |
176 |
else if( args[iArg].length() >= 5 && args[iArg].substring(0, 5).equalsIgnoreCase("-skip") ) |
177 |
{ |
178 |
bSkipFailures = true;
|
179 |
nGroupTransactions = 1; /* #2409 */ |
180 |
} |
181 |
else if( args[iArg].equalsIgnoreCase("-append") ) |
182 |
{ |
183 |
bAppend = true;
|
184 |
bUpdate = true;
|
185 |
} |
186 |
else if( args[iArg].equalsIgnoreCase("-overwrite") ) |
187 |
{ |
188 |
bOverwrite = true;
|
189 |
bUpdate = true;
|
190 |
} |
191 |
else if( args[iArg].equalsIgnoreCase("-update") ) |
192 |
{ |
193 |
bUpdate = true;
|
194 |
} |
195 |
else if( args[iArg].equalsIgnoreCase("-fid") && args[iArg+1] != null ) |
196 |
{ |
197 |
nFIDToFetch = Integer.parseInt(args[++iArg]);
|
198 |
} |
199 |
else if( args[iArg].equalsIgnoreCase("-sql") && args[iArg+1] != null ) |
200 |
{ |
201 |
pszSQLStatement = args[++iArg]; |
202 |
} |
203 |
else if( args[iArg].equalsIgnoreCase("-nln") && iArg < args.length-1 ) |
204 |
{ |
205 |
pszNewLayerName = args[++iArg]; |
206 |
} |
207 |
else if( args[iArg].equalsIgnoreCase("-nlt") && iArg < args.length-1 ) |
208 |
{ |
209 |
if( args[iArg+1].equalsIgnoreCase("NONE") ) |
210 |
eGType = ogr.wkbNone; |
211 |
else if( args[iArg+1].equalsIgnoreCase("GEOMETRY") ) |
212 |
eGType = ogr.wkbUnknown; |
213 |
else if( args[iArg+1].equalsIgnoreCase("POINT") ) |
214 |
eGType = ogr.wkbPoint; |
215 |
else if( args[iArg+1].equalsIgnoreCase("LINESTRING") ) |
216 |
eGType = ogr.wkbLineString; |
217 |
else if( args[iArg+1].equalsIgnoreCase("POLYGON") ) |
218 |
eGType = ogr.wkbPolygon; |
219 |
else if( args[iArg+1].equalsIgnoreCase("GEOMETRYCOLLECTION") ) |
220 |
eGType = ogr.wkbGeometryCollection; |
221 |
else if( args[iArg+1].equalsIgnoreCase("MULTIPOINT") ) |
222 |
eGType = ogr.wkbMultiPoint; |
223 |
else if( args[iArg+1].equalsIgnoreCase("MULTILINESTRING") ) |
224 |
eGType = ogr.wkbMultiLineString; |
225 |
else if( args[iArg+1].equalsIgnoreCase("MULTIPOLYGON") ) |
226 |
eGType = ogr.wkbMultiPolygon; |
227 |
else if( args[iArg+1].equalsIgnoreCase("GEOMETRY25D") ) |
228 |
eGType = ogr.wkbUnknown | ogr.wkb25DBit; |
229 |
else if( args[iArg+1].equalsIgnoreCase("POINT25D") ) |
230 |
eGType = ogr.wkbPoint25D; |
231 |
else if( args[iArg+1].equalsIgnoreCase("LINESTRING25D") ) |
232 |
eGType = ogr.wkbLineString25D; |
233 |
else if( args[iArg+1].equalsIgnoreCase("POLYGON25D") ) |
234 |
eGType = ogr.wkbPolygon25D; |
235 |
else if( args[iArg+1].equalsIgnoreCase("GEOMETRYCOLLECTION25D") ) |
236 |
eGType = ogr.wkbGeometryCollection25D; |
237 |
else if( args[iArg+1].equalsIgnoreCase("MULTIPOINT25D") ) |
238 |
eGType = ogr.wkbMultiPoint25D; |
239 |
else if( args[iArg+1].equalsIgnoreCase("MULTILINESTRING25D") ) |
240 |
eGType = ogr.wkbMultiLineString25D; |
241 |
else if( args[iArg+1].equalsIgnoreCase("MULTIPOLYGON25D") ) |
242 |
eGType = ogr.wkbMultiPolygon25D; |
243 |
else
|
244 |
{ |
245 |
System.err.println("-nlt " + args[iArg+1] + ": type not recognised."); |
246 |
System.exit( 1 ); |
247 |
} |
248 |
iArg++; |
249 |
} |
250 |
else if( (args[iArg].equalsIgnoreCase("-tg") || |
251 |
args[iArg].equalsIgnoreCase("-gt")) && iArg < args.length-1 ) |
252 |
{ |
253 |
nGroupTransactions = Integer.parseInt(args[++iArg]);
|
254 |
} |
255 |
else if( args[iArg].equalsIgnoreCase("-s_srs") && iArg < args.length-1 ) |
256 |
{ |
257 |
pszSourceSRSDef = args[++iArg]; |
258 |
} |
259 |
else if( args[iArg].equalsIgnoreCase("-a_srs") && iArg < args.length-1 ) |
260 |
{ |
261 |
pszOutputSRSDef = args[++iArg]; |
262 |
} |
263 |
else if( args[iArg].equalsIgnoreCase("-t_srs") && iArg < args.length-1 ) |
264 |
{ |
265 |
pszOutputSRSDef = args[++iArg]; |
266 |
bTransform = true;
|
267 |
} |
268 |
else if (args[iArg].equalsIgnoreCase("-spat") && iArg + 4 < args.length) |
269 |
{ |
270 |
Geometry oRing = new Geometry(ogrConstants.wkbLinearRing);
|
271 |
double xmin = new Double(args[++iArg]).doubleValue(); |
272 |
double ymin = new Double(args[++iArg]).doubleValue(); |
273 |
double xmax = new Double(args[++iArg]).doubleValue(); |
274 |
double ymax = new Double(args[++iArg]).doubleValue(); |
275 |
oRing.AddPoint(xmin, ymin); |
276 |
oRing.AddPoint(xmin, ymax); |
277 |
oRing.AddPoint(xmax, ymax); |
278 |
oRing.AddPoint(xmax, ymin); |
279 |
oRing.AddPoint(xmin, ymin); |
280 |
|
281 |
poSpatialFilter = new Geometry(ogrConstants.wkbPolygon);
|
282 |
poSpatialFilter.AddGeometry(oRing); |
283 |
} |
284 |
else if( args[iArg].equalsIgnoreCase("-where") && args[iArg+1] != null ) |
285 |
{ |
286 |
pszWHERE = args[++iArg]; |
287 |
} |
288 |
else if( args[iArg].equalsIgnoreCase("-select") && args[iArg+1] != null) |
289 |
{ |
290 |
pszSelect = args[++iArg]; |
291 |
StringTokenizer tokenizer = new StringTokenizer(pszSelect, " ,"); |
292 |
papszSelFields = new Vector(); |
293 |
while(tokenizer.hasMoreElements())
|
294 |
papszSelFields.addElement(tokenizer.nextToken()); |
295 |
} |
296 |
else if( args[iArg].equalsIgnoreCase("-simplify") && iArg < args.length-1 ) |
297 |
{ |
298 |
eGeomOp = GeomOperation.SIMPLIFY_PRESERVE_TOPOLOGY; |
299 |
dfGeomOpParam = new Double(args[++iArg]).doubleValue(); |
300 |
} |
301 |
else if( args[iArg].equalsIgnoreCase("-segmentize") && iArg < args.length-1 ) |
302 |
{ |
303 |
eGeomOp = GeomOperation.SEGMENTIZE; |
304 |
dfGeomOpParam = new Double(args[++iArg]).doubleValue(); |
305 |
} |
306 |
else if( args[iArg].equalsIgnoreCase("-fieldTypeToString") && iArg < args.length-1 ) |
307 |
{ |
308 |
StringTokenizer tokenizer = new StringTokenizer(args[++iArg], " ,"); |
309 |
while(tokenizer.hasMoreElements())
|
310 |
{ |
311 |
String token = (String)tokenizer.nextToken(); |
312 |
if (token.equalsIgnoreCase("Integer") || |
313 |
token.equalsIgnoreCase("Real") ||
|
314 |
token.equalsIgnoreCase("String") ||
|
315 |
token.equalsIgnoreCase("Date") ||
|
316 |
token.equalsIgnoreCase("Time") ||
|
317 |
token.equalsIgnoreCase("DateTime") ||
|
318 |
token.equalsIgnoreCase("Binary") ||
|
319 |
token.equalsIgnoreCase("IntegerList") ||
|
320 |
token.equalsIgnoreCase("RealList") ||
|
321 |
token.equalsIgnoreCase("StringList"))
|
322 |
{ |
323 |
papszFieldTypesToString.addElement(token); |
324 |
} |
325 |
else if (token.equalsIgnoreCase("All")) |
326 |
{ |
327 |
papszFieldTypesToString = null;
|
328 |
papszFieldTypesToString.addElement("All");
|
329 |
break;
|
330 |
} |
331 |
else
|
332 |
{ |
333 |
System.err.println("Unhandled type for fieldtypeasstring option : " + token); |
334 |
Usage(); |
335 |
} |
336 |
} |
337 |
} |
338 |
else if( args[iArg].equalsIgnoreCase("-progress") ) |
339 |
{ |
340 |
bDisplayProgress = true;
|
341 |
} |
342 |
/*else if( args[iArg].equalsIgnoreCase("-wrapdateline") )
|
343 |
{
|
344 |
bWrapDateline = true;
|
345 |
}
|
346 |
*/
|
347 |
else if( args[iArg].equalsIgnoreCase("-clipsrc") && iArg < args.length-1 ) |
348 |
{ |
349 |
bClipSrc = true;
|
350 |
if ( IsNumber(args[iArg+1]) && iArg < args.length - 4 ) |
351 |
{ |
352 |
Geometry oRing = new Geometry(ogrConstants.wkbLinearRing);
|
353 |
double xmin = new Double(args[++iArg]).doubleValue(); |
354 |
double ymin = new Double(args[++iArg]).doubleValue(); |
355 |
double xmax = new Double(args[++iArg]).doubleValue(); |
356 |
double ymax = new Double(args[++iArg]).doubleValue(); |
357 |
oRing.AddPoint(xmin, ymin); |
358 |
oRing.AddPoint(xmin, ymax); |
359 |
oRing.AddPoint(xmax, ymax); |
360 |
oRing.AddPoint(xmax, ymin); |
361 |
oRing.AddPoint(xmin, ymin); |
362 |
|
363 |
poClipSrc = new Geometry(ogrConstants.wkbPolygon);
|
364 |
poClipSrc.AddGeometry(oRing); |
365 |
} |
366 |
else if( (args[iArg+1].length() >= 7 && args[iArg+1].substring(0, 7).equalsIgnoreCase("POLYGON") ) || |
367 |
(args[iArg+1].length() >= 12 && args[iArg+1].substring(0, 12).equalsIgnoreCase("MULTIPOLYGON") ) ) |
368 |
{ |
369 |
poClipSrc = Geometry.CreateFromWkt(args[iArg+1]);
|
370 |
if (poClipSrc == null) |
371 |
{ |
372 |
System.err.print("FAILURE: Invalid geometry. Must be a valid POLYGON or MULTIPOLYGON WKT\n\n"); |
373 |
Usage(); |
374 |
} |
375 |
iArg ++; |
376 |
} |
377 |
else if (args[iArg+1].equalsIgnoreCase("spat_extent") ) |
378 |
{ |
379 |
iArg ++; |
380 |
} |
381 |
else
|
382 |
{ |
383 |
pszClipSrcDS = args[iArg+1];
|
384 |
iArg ++; |
385 |
} |
386 |
} |
387 |
else if( args[iArg].equalsIgnoreCase("-clipsrcsql") && iArg < args.length-1 ) |
388 |
{ |
389 |
pszClipSrcSQL = args[iArg+1];
|
390 |
iArg ++; |
391 |
} |
392 |
else if( args[iArg].equalsIgnoreCase("-clipsrclayer") && iArg < args.length-1 ) |
393 |
{ |
394 |
pszClipSrcLayer = args[iArg+1];
|
395 |
iArg ++; |
396 |
} |
397 |
else if( args[iArg].equalsIgnoreCase("-clipsrcwhere") && iArg < args.length-1 ) |
398 |
{ |
399 |
pszClipSrcWhere = args[iArg+1];
|
400 |
iArg ++; |
401 |
} |
402 |
else if( args[iArg].equalsIgnoreCase("-clipdst") && iArg < args.length-1 ) |
403 |
{ |
404 |
if ( IsNumber(args[iArg+1]) && iArg < args.length - 4 ) |
405 |
{ |
406 |
Geometry oRing = new Geometry(ogrConstants.wkbLinearRing);
|
407 |
double xmin = new Double(args[++iArg]).doubleValue(); |
408 |
double ymin = new Double(args[++iArg]).doubleValue(); |
409 |
double xmax = new Double(args[++iArg]).doubleValue(); |
410 |
double ymax = new Double(args[++iArg]).doubleValue(); |
411 |
oRing.AddPoint(xmin, ymin); |
412 |
oRing.AddPoint(xmin, ymax); |
413 |
oRing.AddPoint(xmax, ymax); |
414 |
oRing.AddPoint(xmax, ymin); |
415 |
oRing.AddPoint(xmin, ymin); |
416 |
|
417 |
poClipDst = new Geometry(ogrConstants.wkbPolygon);
|
418 |
poClipDst.AddGeometry(oRing); |
419 |
} |
420 |
else if( (args[iArg+1].length() >= 7 && args[iArg+1].substring(0, 7).equalsIgnoreCase("POLYGON") ) || |
421 |
(args[iArg+1].length() >= 12 && args[iArg+1].substring(0, 12).equalsIgnoreCase("MULTIPOLYGON") ) ) |
422 |
{ |
423 |
poClipDst = Geometry.CreateFromWkt(args[iArg+1]);
|
424 |
if (poClipDst == null) |
425 |
{ |
426 |
System.err.print("FAILURE: Invalid geometry. Must be a valid POLYGON or MULTIPOLYGON WKT\n\n"); |
427 |
Usage(); |
428 |
} |
429 |
iArg ++; |
430 |
} |
431 |
else if (args[iArg+1].equalsIgnoreCase("spat_extent") ) |
432 |
{ |
433 |
iArg ++; |
434 |
} |
435 |
else
|
436 |
{ |
437 |
pszClipDstDS = args[iArg+1];
|
438 |
iArg ++; |
439 |
} |
440 |
} |
441 |
else if( args[iArg].equalsIgnoreCase("-clipdstsql") && iArg < args.length-1 ) |
442 |
{ |
443 |
pszClipDstSQL = args[iArg+1];
|
444 |
iArg ++; |
445 |
} |
446 |
else if( args[iArg].equalsIgnoreCase("-clipdstlayer") && iArg < args.length-1 ) |
447 |
{ |
448 |
pszClipDstLayer = args[iArg+1];
|
449 |
iArg ++; |
450 |
} |
451 |
else if( args[iArg].equalsIgnoreCase("-clipdstwhere") && iArg < args.length-1 ) |
452 |
{ |
453 |
pszClipDstWhere = args[iArg+1];
|
454 |
iArg ++; |
455 |
} |
456 |
else if( args[iArg].equalsIgnoreCase("-explodecollections") ) |
457 |
{ |
458 |
bExplodeCollections = true;
|
459 |
} |
460 |
else if( args[iArg].equalsIgnoreCase("-zfield") && iArg < args.length-1 ) |
461 |
{ |
462 |
pszZField = args[iArg+1];
|
463 |
iArg ++; |
464 |
} |
465 |
else if( args[iArg].charAt(0) == '-' ) |
466 |
{ |
467 |
Usage(); |
468 |
} |
469 |
else if( pszDestDataSource == null ) |
470 |
pszDestDataSource = args[iArg]; |
471 |
else if( pszDataSource == null ) |
472 |
pszDataSource = args[iArg]; |
473 |
else
|
474 |
papszLayers.addElement (args[iArg] ); |
475 |
} |
476 |
|
477 |
if( pszDataSource == null ) |
478 |
Usage(); |
479 |
|
480 |
if( bPreserveFID && bExplodeCollections )
|
481 |
{ |
482 |
System.err.print("FAILURE: cannot use -preserve_fid and -explodecollections at the same time\n\n" ); |
483 |
Usage(); |
484 |
} |
485 |
|
486 |
if( bClipSrc && pszClipSrcDS != null) |
487 |
{ |
488 |
poClipSrc = LoadGeometry(pszClipSrcDS, pszClipSrcSQL, pszClipSrcLayer, pszClipSrcWhere); |
489 |
if (poClipSrc == null) |
490 |
{ |
491 |
System.err.print("FAILURE: cannot load source clip geometry\n\n" ); |
492 |
Usage(); |
493 |
} |
494 |
} |
495 |
else if( bClipSrc && poClipSrc == null ) |
496 |
{ |
497 |
if (poSpatialFilter != null) |
498 |
poClipSrc = poSpatialFilter.Clone(); |
499 |
if (poClipSrc == null) |
500 |
{ |
501 |
System.err.print("FAILURE: -clipsrc must be used with -spat option or a\n" + |
502 |
"bounding box, WKT string or datasource must be specified\n\n");
|
503 |
Usage(); |
504 |
} |
505 |
} |
506 |
|
507 |
if( pszClipDstDS != null) |
508 |
{ |
509 |
poClipDst = LoadGeometry(pszClipDstDS, pszClipDstSQL, pszClipDstLayer, pszClipDstWhere); |
510 |
if (poClipDst == null) |
511 |
{ |
512 |
System.err.print("FAILURE: cannot load dest clip geometry\n\n" ); |
513 |
Usage(); |
514 |
} |
515 |
} |
516 |
/* -------------------------------------------------------------------- */
|
517 |
/* Open data source. */
|
518 |
/* -------------------------------------------------------------------- */
|
519 |
DataSource poDS;
|
520 |
|
521 |
poDS = ogr.Open( pszDataSource, false );
|
522 |
|
523 |
/* -------------------------------------------------------------------- */
|
524 |
/* Report failure */
|
525 |
/* -------------------------------------------------------------------- */
|
526 |
if( poDS == null ) |
527 |
{ |
528 |
System.err.println("FAILURE:\n" + |
529 |
"Unable to open datasource ` " + pszDataSource + "' with the following drivers."); |
530 |
|
531 |
for( int iDriver = 0; iDriver < ogr.GetDriverCount(); iDriver++ ) |
532 |
{ |
533 |
System.err.println(" . " + ogr.GetDriver(iDriver).GetName() ); |
534 |
} |
535 |
|
536 |
System.exit( 1 ); |
537 |
} |
538 |
|
539 |
/* -------------------------------------------------------------------- */
|
540 |
/* Try opening the output datasource as an existing, writable */
|
541 |
/* -------------------------------------------------------------------- */
|
542 |
DataSource poODS = null; |
543 |
Driver poDriver = null; |
544 |
|
545 |
if( bUpdate )
|
546 |
{ |
547 |
poODS = ogr.Open( pszDestDataSource, true );
|
548 |
if( poODS == null ) |
549 |
{ |
550 |
if (bOverwrite || bAppend)
|
551 |
{ |
552 |
poODS = ogr.Open( pszDestDataSource, false );
|
553 |
if ( poODS == null ) |
554 |
{ |
555 |
/* ok the datasource doesn't exist at all */
|
556 |
bUpdate = false;
|
557 |
} |
558 |
else
|
559 |
{ |
560 |
poODS.delete(); |
561 |
poODS = null;
|
562 |
} |
563 |
} |
564 |
|
565 |
if (bUpdate)
|
566 |
{ |
567 |
System.err.println("FAILURE:\n" + |
568 |
"Unable to open existing output datasource `" + pszDestDataSource + "'."); |
569 |
System.exit( 1 ); |
570 |
} |
571 |
} |
572 |
|
573 |
else if( papszDSCO.size() > 0 ) |
574 |
{ |
575 |
System.err.println("WARNING: Datasource creation options ignored since an existing datasource\n" + |
576 |
" being updated." );
|
577 |
} |
578 |
|
579 |
if (poODS != null) |
580 |
poDriver = poODS.GetDriver(); |
581 |
} |
582 |
|
583 |
/* -------------------------------------------------------------------- */
|
584 |
/* Find the output driver. */
|
585 |
/* -------------------------------------------------------------------- */
|
586 |
if( !bUpdate )
|
587 |
{ |
588 |
int iDriver;
|
589 |
|
590 |
poDriver = ogr.GetDriverByName(pszFormat); |
591 |
if( poDriver == null ) |
592 |
{ |
593 |
System.err.println("Unable to find driver `" + pszFormat +"'." ); |
594 |
System.err.println( "The following drivers are available:" ); |
595 |
|
596 |
for( iDriver = 0; iDriver < ogr.GetDriverCount(); iDriver++ ) |
597 |
{ |
598 |
System.err.println(" . " + ogr.GetDriver(iDriver).GetName() ); |
599 |
} |
600 |
System.exit( 1 ); |
601 |
} |
602 |
|
603 |
if( poDriver.TestCapability( ogr.ODrCCreateDataSource ) == false ) |
604 |
{ |
605 |
System.err.println( pszFormat + " driver does not support data source creation."); |
606 |
System.exit( 1 ); |
607 |
} |
608 |
|
609 |
/* -------------------------------------------------------------------- */
|
610 |
/* Special case to improve user experience when translating */
|
611 |
/* a datasource with multiple layers into a shapefile. If the */
|
612 |
/* user gives a target datasource with .shp and it does not exist, */
|
613 |
/* the shapefile driver will try to create a file, but this is not */
|
614 |
/* appropriate because here we have several layers, so create */
|
615 |
/* a directory instead. */
|
616 |
/* -------------------------------------------------------------------- */
|
617 |
if (poDriver.GetName().equalsIgnoreCase("ESRI Shapefile") && |
618 |
pszSQLStatement == null &&
|
619 |
(papszLayers.size() > 1 ||
|
620 |
(papszLayers.size() == 0 && poDS.GetLayerCount() > 1)) && |
621 |
pszNewLayerName == null &&
|
622 |
(pszDestDataSource.endsWith(".shp") || pszDestDataSource.endsWith(".SHP"))) |
623 |
{ |
624 |
File f = new File(pszDestDataSource); |
625 |
if (!f.exists())
|
626 |
{ |
627 |
if (!f.mkdir())
|
628 |
{ |
629 |
System.err.println(
|
630 |
"Failed to create directory " + pszDestDataSource + "\n" + |
631 |
"for shapefile datastore.");
|
632 |
System.exit(1); |
633 |
} |
634 |
} |
635 |
} |
636 |
|
637 |
/* -------------------------------------------------------------------- */
|
638 |
/* Create the output data source. */
|
639 |
/* -------------------------------------------------------------------- */
|
640 |
poODS = poDriver.CreateDataSource( pszDestDataSource, papszDSCO ); |
641 |
if( poODS == null ) |
642 |
{ |
643 |
System.err.println( pszFormat + " driver failed to create "+ pszDestDataSource ); |
644 |
System.exit( 1 ); |
645 |
} |
646 |
} |
647 |
|
648 |
/* -------------------------------------------------------------------- */
|
649 |
/* Parse the output SRS definition if possible. */
|
650 |
/* -------------------------------------------------------------------- */
|
651 |
if( pszOutputSRSDef != null ) |
652 |
{ |
653 |
poOutputSRS = new SpatialReference();
|
654 |
if( poOutputSRS.SetFromUserInput( pszOutputSRSDef ) != 0 ) |
655 |
{ |
656 |
System.err.println( "Failed to process SRS definition: " + pszOutputSRSDef ); |
657 |
System.exit( 1 ); |
658 |
} |
659 |
} |
660 |
|
661 |
/* -------------------------------------------------------------------- */
|
662 |
/* Parse the source SRS definition if possible. */
|
663 |
/* -------------------------------------------------------------------- */
|
664 |
if( pszSourceSRSDef != null ) |
665 |
{ |
666 |
poSourceSRS = new SpatialReference();
|
667 |
if( poSourceSRS.SetFromUserInput( pszSourceSRSDef ) != 0 ) |
668 |
{ |
669 |
System.err.println( "Failed to process SRS definition: " + pszSourceSRSDef ); |
670 |
System.exit( 1 ); |
671 |
} |
672 |
} |
673 |
|
674 |
/* -------------------------------------------------------------------- */
|
675 |
/* Special case for -sql clause. No source layers required. */
|
676 |
/* -------------------------------------------------------------------- */
|
677 |
if( pszSQLStatement != null ) |
678 |
{ |
679 |
Layer poResultSet; |
680 |
|
681 |
if( pszWHERE != null ) |
682 |
System.err.println( "-where clause ignored in combination with -sql." ); |
683 |
if( papszLayers.size() > 0 ) |
684 |
System.err.println( "layer names ignored in combination with -sql." ); |
685 |
|
686 |
poResultSet = poDS.ExecuteSQL( pszSQLStatement, poSpatialFilter, |
687 |
null );
|
688 |
|
689 |
if( poResultSet != null ) |
690 |
{ |
691 |
long nCountLayerFeatures = 0; |
692 |
if (bDisplayProgress)
|
693 |
{ |
694 |
if (!poResultSet.TestCapability(ogr.OLCFastFeatureCount))
|
695 |
{ |
696 |
System.err.println( "Progress turned off as fast feature count is not available."); |
697 |
bDisplayProgress = false;
|
698 |
} |
699 |
else
|
700 |
{ |
701 |
nCountLayerFeatures = poResultSet.GetFeatureCount(); |
702 |
pfnProgress = new TermProgressCallback();
|
703 |
} |
704 |
} |
705 |
|
706 |
/* -------------------------------------------------------------------- */
|
707 |
/* Special case to improve user experience when translating into */
|
708 |
/* single file shapefile and source has only one layer, and that */
|
709 |
/* the layer name isn't specified */
|
710 |
/* -------------------------------------------------------------------- */
|
711 |
if (poDriver.GetName().equalsIgnoreCase("ESRI Shapefile") && |
712 |
pszNewLayerName == null)
|
713 |
{ |
714 |
File f = new File(pszDestDataSource); |
715 |
if (f.exists() && f.listFiles() == null) |
716 |
{ |
717 |
pszNewLayerName = f.getName(); |
718 |
int posPoint = pszNewLayerName.lastIndexOf('.'); |
719 |
if (posPoint != -1) |
720 |
pszNewLayerName = pszNewLayerName.substring(0, posPoint);
|
721 |
} |
722 |
} |
723 |
|
724 |
if( !TranslateLayer( poDS, poResultSet, poODS, papszLCO,
|
725 |
pszNewLayerName, bTransform, poOutputSRS, |
726 |
poSourceSRS, papszSelFields, bAppend, eGType, |
727 |
bOverwrite, eGeomOp, dfGeomOpParam, papszFieldTypesToString, |
728 |
nCountLayerFeatures, poClipSrc, poClipDst, bExplodeCollections, |
729 |
pszZField, pszWHERE, pfnProgress )) |
730 |
{ |
731 |
System.err.println(
|
732 |
"Terminating translation prematurely after failed\n" +
|
733 |
"translation from sql statement." );
|
734 |
|
735 |
System.exit( 1 ); |
736 |
} |
737 |
poDS.ReleaseResultSet( poResultSet ); |
738 |
} |
739 |
} |
740 |
else
|
741 |
{ |
742 |
int nLayerCount = 0; |
743 |
Layer[] papoLayers = null; |
744 |
|
745 |
/* -------------------------------------------------------------------- */
|
746 |
/* Process each data source layer. */
|
747 |
/* -------------------------------------------------------------------- */
|
748 |
if ( papszLayers.size() == 0) |
749 |
{ |
750 |
nLayerCount = poDS.GetLayerCount(); |
751 |
papoLayers = new Layer[nLayerCount];
|
752 |
|
753 |
for( int iLayer = 0; |
754 |
iLayer < nLayerCount; |
755 |
iLayer++ ) |
756 |
{ |
757 |
Layer poLayer = poDS.GetLayer(iLayer); |
758 |
|
759 |
if( poLayer == null ) |
760 |
{ |
761 |
System.err.println("FAILURE: Couldn't fetch advertised layer " + iLayer + "!"); |
762 |
System.exit( 1 ); |
763 |
} |
764 |
|
765 |
papoLayers[iLayer] = poLayer; |
766 |
} |
767 |
} |
768 |
/* -------------------------------------------------------------------- */
|
769 |
/* Process specified data source layers. */
|
770 |
/* -------------------------------------------------------------------- */
|
771 |
else
|
772 |
{ |
773 |
nLayerCount = papszLayers.size(); |
774 |
papoLayers = new Layer[nLayerCount];
|
775 |
|
776 |
for( int iLayer = 0; |
777 |
iLayer < papszLayers.size(); |
778 |
iLayer++ ) |
779 |
{ |
780 |
Layer poLayer = poDS.GetLayerByName((String)papszLayers.get(iLayer));
|
781 |
|
782 |
if( poLayer == null ) |
783 |
{ |
784 |
System.err.println("FAILURE: Couldn't fetch advertised layer " + (String)papszLayers.get(iLayer) + "!"); |
785 |
System.exit( 1 ); |
786 |
} |
787 |
|
788 |
papoLayers[iLayer] = poLayer; |
789 |
} |
790 |
} |
791 |
|
792 |
/* -------------------------------------------------------------------- */
|
793 |
/* Special case to improve user experience when translating into */
|
794 |
/* single file shapefile and source has only one layer, and that */
|
795 |
/* the layer name isn't specified */
|
796 |
/* -------------------------------------------------------------------- */
|
797 |
if (poDriver.GetName().equalsIgnoreCase("ESRI Shapefile") && |
798 |
nLayerCount == 1 && pszNewLayerName == null) |
799 |
{ |
800 |
File f = new File(pszDestDataSource); |
801 |
if (f.exists() && f.listFiles() == null) |
802 |
{ |
803 |
pszNewLayerName = f.getName(); |
804 |
int posPoint = pszNewLayerName.lastIndexOf('.'); |
805 |
if (posPoint != -1) |
806 |
pszNewLayerName = pszNewLayerName.substring(0, posPoint);
|
807 |
} |
808 |
} |
809 |
|
810 |
long[] panLayerCountFeatures = new long[nLayerCount]; |
811 |
long nCountLayersFeatures = 0; |
812 |
long nAccCountFeatures = 0; |
813 |
|
814 |
/* First pass to apply filters and count all features if necessary */
|
815 |
for( int iLayer = 0; |
816 |
iLayer < nLayerCount; |
817 |
iLayer++ ) |
818 |
{ |
819 |
Layer poLayer = papoLayers[iLayer]; |
820 |
|
821 |
if( pszWHERE != null ) |
822 |
{ |
823 |
if( poLayer.SetAttributeFilter( pszWHERE ) != ogr.OGRERR_NONE )
|
824 |
{ |
825 |
System.err.println("FAILURE: SetAttributeFilter(" + pszWHERE + ") failed."); |
826 |
if (!bSkipFailures)
|
827 |
System.exit( 1 ); |
828 |
} |
829 |
} |
830 |
|
831 |
if( poSpatialFilter != null ) |
832 |
poLayer.SetSpatialFilter( poSpatialFilter ); |
833 |
|
834 |
if (bDisplayProgress)
|
835 |
{ |
836 |
if (!poLayer.TestCapability(ogr.OLCFastFeatureCount))
|
837 |
{ |
838 |
System.err.println("Progress turned off as fast feature count is not available."); |
839 |
bDisplayProgress = false;
|
840 |
} |
841 |
else
|
842 |
{ |
843 |
panLayerCountFeatures[iLayer] = poLayer.GetFeatureCount(); |
844 |
nCountLayersFeatures += panLayerCountFeatures[iLayer]; |
845 |
} |
846 |
} |
847 |
} |
848 |
|
849 |
/* Second pass to do the real job */
|
850 |
for( int iLayer = 0; |
851 |
iLayer < nLayerCount; |
852 |
iLayer++ ) |
853 |
{ |
854 |
Layer poLayer = papoLayers[iLayer]; |
855 |
|
856 |
if (bDisplayProgress)
|
857 |
{ |
858 |
pfnProgress = new GDALScaledProgress(
|
859 |
nAccCountFeatures * 1.0 / nCountLayersFeatures,
|
860 |
(nAccCountFeatures + panLayerCountFeatures[iLayer]) * 1.0 / nCountLayersFeatures,
|
861 |
new TermProgressCallback());
|
862 |
} |
863 |
|
864 |
nAccCountFeatures += panLayerCountFeatures[iLayer]; |
865 |
|
866 |
if( !TranslateLayer( poDS, poLayer, poODS, papszLCO,
|
867 |
pszNewLayerName, bTransform, poOutputSRS, |
868 |
poSourceSRS, papszSelFields, bAppend, eGType, |
869 |
bOverwrite, eGeomOp, dfGeomOpParam, papszFieldTypesToString, |
870 |
panLayerCountFeatures[iLayer], poClipSrc, poClipDst, bExplodeCollections, |
871 |
pszZField, pszWHERE, pfnProgress) |
872 |
&& !bSkipFailures ) |
873 |
{ |
874 |
System.err.println(
|
875 |
"Terminating translation prematurely after failed\n" +
|
876 |
"translation of layer " + poLayer.GetLayerDefn().GetName() + " (use -skipfailures to skip errors)"); |
877 |
|
878 |
System.exit( 1 ); |
879 |
} |
880 |
} |
881 |
} |
882 |
|
883 |
/* -------------------------------------------------------------------- */
|
884 |
/* Close down. */
|
885 |
/* -------------------------------------------------------------------- */
|
886 |
/* We must explicitly destroy the output dataset in order the file */
|
887 |
/* to be properly closed ! */
|
888 |
poODS.delete(); |
889 |
poDS.delete(); |
890 |
} |
891 |
|
892 |
/************************************************************************/
|
893 |
/* Usage() */
|
894 |
/************************************************************************/
|
895 |
|
896 |
static void Usage() |
897 |
|
898 |
{ |
899 |
System.out.print( "Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update] [-gt n]\n" + |
900 |
" [-select field_list] [-where restricted_where] \n" +
|
901 |
" [-progress] [-sql <sql statement>] \n" +
|
902 |
" [-spat xmin ymin xmax ymax] [-preserve_fid] [-fid FID]\n" +
|
903 |
" [-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def]\n" +
|
904 |
" [-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]\n" +
|
905 |
" [-simplify tolerance]\n" +
|
906 |
// " [-segmentize max_dist] [-fieldTypeToString All|(type1[,type2]*)]\n" +
|
907 |
" [-fieldTypeToString All|(type1[,type2]*)] [-explodecollections]\n" +
|
908 |
" dst_datasource_name src_datasource_name\n" +
|
909 |
" [-lco NAME=VALUE] [-nln name] [-nlt type] [layer [layer ...]]\n" +
|
910 |
"\n" +
|
911 |
" -f format_name: output file format name, possible values are:\n");
|
912 |
|
913 |
for( int iDriver = 0; iDriver < ogr.GetDriverCount(); iDriver++ ) |
914 |
{ |
915 |
Driver poDriver = ogr.GetDriver(iDriver);
|
916 |
|
917 |
if( poDriver.TestCapability( ogr.ODrCCreateDataSource ) )
|
918 |
System.out.print( " -f \"" + poDriver.GetName() + "\"\n" ); |
919 |
} |
920 |
|
921 |
System.out.print( " -append: Append to existing layer instead of creating new if it exists\n" + |
922 |
" -overwrite: delete the output layer and recreate it empty\n" +
|
923 |
" -update: Open existing output datasource in update mode\n" +
|
924 |
" -progress: Display progress on terminal. Only works if input layers have the \"fast feature count\" capability\n" +
|
925 |
" -select field_list: Comma-delimited list of fields from input layer to\n" +
|
926 |
" copy to the new layer (defaults to all)\n" +
|
927 |
" -where restricted_where: Attribute query (like SQL WHERE)\n" +
|
928 |
" -sql statement: Execute given SQL statement and save result.\n" +
|
929 |
" -skipfailures: skip features or layers that fail to convert\n" +
|
930 |
" -gt n: group n features per transaction (default 200)\n" +
|
931 |
" -spat xmin ymin xmax ymax: spatial query extents\n" +
|
932 |
" -simplify tolerance: distance tolerance for simplification.\n" +
|
933 |
//" -segmentize max_dist: maximum distance between 2 nodes.\n" +
|
934 |
//" Used to create intermediate points\n" +
|
935 |
" -dsco NAME=VALUE: Dataset creation option (format specific)\n" +
|
936 |
" -lco NAME=VALUE: Layer creation option (format specific)\n" +
|
937 |
" -nln name: Assign an alternate name to the new layer\n" +
|
938 |
" -nlt type: Force a geometry type for new layer. One of NONE, GEOMETRY,\n" +
|
939 |
" POINT, LINESTRING, POLYGON, GEOMETRYCOLLECTION, MULTIPOINT,\n" +
|
940 |
" MULTIPOLYGON, or MULTILINESTRING. Add \"25D\" for 3D layers.\n" +
|
941 |
" Default is type of source layer.\n" +
|
942 |
" -fieldTypeToString type1,...: Converts fields of specified types to\n" +
|
943 |
" fields of type string in the new layer. Valid types are : \n" +
|
944 |
" Integer, Real, String, Date, Time, DateTime, Binary, IntegerList, RealList,\n" +
|
945 |
" StringList. Special value All can be used to convert all fields to strings.\n");
|
946 |
|
947 |
System.out.print(" -a_srs srs_def: Assign an output SRS\n" + |
948 |
" -t_srs srs_def: Reproject/transform to this SRS on output\n" +
|
949 |
" -s_srs srs_def: Override source SRS\n" +
|
950 |
"\n" +
|
951 |
" Srs_def can be a full WKT definition (hard to escape properly),\n" +
|
952 |
" or a well known definition (i.e. EPSG:4326) or a file with a WKT\n" +
|
953 |
" definition.\n" );
|
954 |
|
955 |
System.exit( 1 ); |
956 |
} |
957 |
|
958 |
static int CSLFindString(Vector v, String str) |
959 |
{ |
960 |
int i = 0; |
961 |
Enumeration e = v.elements();
|
962 |
while(e.hasMoreElements())
|
963 |
{ |
964 |
String strIter = (String)e.nextElement(); |
965 |
if (strIter.equalsIgnoreCase(str))
|
966 |
return i;
|
967 |
i ++; |
968 |
} |
969 |
return -1; |
970 |
} |
971 |
|
972 |
static boolean IsNumber(String pszStr) |
973 |
{ |
974 |
try
|
975 |
{ |
976 |
Double.parseDouble(pszStr);
|
977 |
return true; |
978 |
} |
979 |
catch(Exception ex) |
980 |
{ |
981 |
return false; |
982 |
} |
983 |
} |
984 |
|
985 |
static Geometry LoadGeometry( String pszDS, |
986 |
String pszSQL,
|
987 |
String pszLyr,
|
988 |
String pszWhere)
|
989 |
{ |
990 |
DataSource poDS;
|
991 |
Layer poLyr; |
992 |
Feature poFeat; |
993 |
Geometry poGeom = null;
|
994 |
|
995 |
poDS = ogr.Open( pszDS, false );
|
996 |
if (poDS == null) |
997 |
return null; |
998 |
|
999 |
if (pszSQL != null) |
1000 |
poLyr = poDS.ExecuteSQL( pszSQL, null, null ); |
1001 |
else if (pszLyr != null) |
1002 |
poLyr = poDS.GetLayerByName(pszLyr); |
1003 |
else
|
1004 |
poLyr = poDS.GetLayer(0);
|
1005 |
|
1006 |
if (poLyr == null) |
1007 |
{ |
1008 |
System.err.print("Failed to identify source layer from datasource.\n"); |
1009 |
poDS.delete(); |
1010 |
return null; |
1011 |
} |
1012 |
|
1013 |
if (pszWhere != null) |
1014 |
poLyr.SetAttributeFilter(pszWhere); |
1015 |
|
1016 |
while ((poFeat = poLyr.GetNextFeature()) != null) |
1017 |
{ |
1018 |
Geometry poSrcGeom = poFeat.GetGeometryRef(); |
1019 |
if (poSrcGeom != null) |
1020 |
{ |
1021 |
int eType = wkbFlatten(poSrcGeom.GetGeometryType());
|
1022 |
|
1023 |
if (poGeom == null) |
1024 |
poGeom = new Geometry( ogr.wkbMultiPolygon );
|
1025 |
|
1026 |
if( eType == ogr.wkbPolygon )
|
1027 |
poGeom.AddGeometry( poSrcGeom ); |
1028 |
else if( eType == ogr.wkbMultiPolygon ) |
1029 |
{ |
1030 |
int iGeom;
|
1031 |
int nGeomCount = poSrcGeom.GetGeometryCount();
|
1032 |
|
1033 |
for( iGeom = 0; iGeom < nGeomCount; iGeom++ ) |
1034 |
{ |
1035 |
poGeom.AddGeometry(poSrcGeom.GetGeometryRef(iGeom) ); |
1036 |
} |
1037 |
} |
1038 |
else
|
1039 |
{ |
1040 |
System.err.print("ERROR: Geometry not of polygon type.\n" ); |
1041 |
if( pszSQL != null ) |
1042 |
poDS.ReleaseResultSet( poLyr ); |
1043 |
poDS.delete(); |
1044 |
return null; |
1045 |
} |
1046 |
} |
1047 |
} |
1048 |
|
1049 |
if( pszSQL != null ) |
1050 |
poDS.ReleaseResultSet( poLyr ); |
1051 |
poDS.delete(); |
1052 |
|
1053 |
return poGeom;
|
1054 |
} |
1055 |
|
1056 |
|
1057 |
static int wkbFlatten(int eType) |
1058 |
{ |
1059 |
return eType & (~ogrConstants.wkb25DBit);
|
1060 |
} |
1061 |
|
1062 |
/************************************************************************/
|
1063 |
/* SetZ() */
|
1064 |
/************************************************************************/
|
1065 |
static void SetZ (Geometry poGeom, double dfZ ) |
1066 |
{ |
1067 |
if (poGeom == null) |
1068 |
return;
|
1069 |
switch (wkbFlatten(poGeom.GetGeometryType()))
|
1070 |
{ |
1071 |
case ogr.wkbPoint:
|
1072 |
poGeom.SetPoint(0, poGeom.GetX(), poGeom.GetY(), dfZ);
|
1073 |
break;
|
1074 |
|
1075 |
case ogr.wkbLineString:
|
1076 |
case ogr.wkbLinearRing:
|
1077 |
{ |
1078 |
int i;
|
1079 |
for(i=0;i<poGeom.GetPointCount();i++) |
1080 |
poGeom.SetPoint(i, poGeom.GetX(i), poGeom.GetY(i), dfZ); |
1081 |
break;
|
1082 |
} |
1083 |
|
1084 |
case ogr.wkbPolygon:
|
1085 |
case ogr.wkbMultiPoint:
|
1086 |
case ogr.wkbMultiLineString:
|
1087 |
case ogr.wkbMultiPolygon:
|
1088 |
case ogr.wkbGeometryCollection:
|
1089 |
{ |
1090 |
int i;
|
1091 |
for(i=0;i<poGeom.GetGeometryCount();i++) |
1092 |
SetZ(poGeom.GetGeometryRef(i), dfZ); |
1093 |
break;
|
1094 |
} |
1095 |
|
1096 |
default:
|
1097 |
break;
|
1098 |
} |
1099 |
} |
1100 |
|
1101 |
|
1102 |
/************************************************************************/
|
1103 |
/* TranslateLayer() */
|
1104 |
/************************************************************************/
|
1105 |
|
1106 |
static boolean TranslateLayer( DataSource poSrcDS, |
1107 |
Layer poSrcLayer, |
1108 |
DataSource poDstDS,
|
1109 |
Vector papszLCO,
|
1110 |
String pszNewLayerName,
|
1111 |
boolean bTransform,
|
1112 |
SpatialReference poOutputSRS, |
1113 |
SpatialReference poSourceSRS, |
1114 |
Vector papszSelFields,
|
1115 |
boolean bAppend, int eGType, boolean bOverwrite, |
1116 |
GeomOperation eGeomOp, |
1117 |
double dfGeomOpParam,
|
1118 |
Vector papszFieldTypesToString,
|
1119 |
long nCountLayerFeatures,
|
1120 |
Geometry poClipSrc, |
1121 |
Geometry poClipDst, |
1122 |
boolean bExplodeCollections,
|
1123 |
String pszZField,
|
1124 |
String pszWHERE,
|
1125 |
ProgressCallback pfnProgress) |
1126 |
|
1127 |
{ |
1128 |
Layer poDstLayer; |
1129 |
FeatureDefn poSrcFDefn; |
1130 |
int eErr;
|
1131 |
boolean bForceToPolygon = false; |
1132 |
boolean bForceToMultiPolygon = false; |
1133 |
boolean bForceToMultiLineString = false; |
1134 |
|
1135 |
if( pszNewLayerName == null ) |
1136 |
pszNewLayerName = poSrcLayer.GetLayerDefn().GetName(); |
1137 |
|
1138 |
if( wkbFlatten(eGType) == ogr.wkbPolygon )
|
1139 |
bForceToPolygon = true;
|
1140 |
else if( wkbFlatten(eGType) == ogr.wkbMultiPolygon ) |
1141 |
bForceToMultiPolygon = true;
|
1142 |
else if( wkbFlatten(eGType) == ogr.wkbMultiLineString ) |
1143 |
bForceToMultiLineString = true;
|
1144 |
|
1145 |
/* -------------------------------------------------------------------- */
|
1146 |
/* Setup coordinate transformation if we need it. */
|
1147 |
/* -------------------------------------------------------------------- */
|
1148 |
CoordinateTransformation poCT = null;
|
1149 |
|
1150 |
if( bTransform )
|
1151 |
{ |
1152 |
if( poSourceSRS == null ) |
1153 |
poSourceSRS = poSrcLayer.GetSpatialRef(); |
1154 |
|
1155 |
if( poSourceSRS == null ) |
1156 |
{ |
1157 |
System.err.println("Can't transform coordinates, source layer has no\n" + |
1158 |
"coordinate system. Use -s_srs to set one." );
|
1159 |
System.exit( 1 ); |
1160 |
} |
1161 |
|
1162 |
/*CPLAssert( null != poSourceSRS );
|
1163 |
CPLAssert( null != poOutputSRS );*/
|
1164 |
|
1165 |
/* New in GDAL 1.10. Before was "new CoordinateTransformation(srs,dst)". */
|
1166 |
poCT = CoordinateTransformation.CreateCoordinateTransformation( poSourceSRS, poOutputSRS ); |
1167 |
if( poCT == null ) |
1168 |
{ |
1169 |
String pszWKT = null; |
1170 |
|
1171 |
System.err.println("Failed to create coordinate transformation between the\n" + |
1172 |
"following coordinate systems. This may be because they\n" +
|
1173 |
"are not transformable, or because projection services\n" +
|
1174 |
"(PROJ.4 DLL/.so) could not be loaded." );
|
1175 |
|
1176 |
pszWKT = poSourceSRS.ExportToPrettyWkt( 0 );
|
1177 |
System.err.println( "Source:\n" + pszWKT ); |
1178 |
|
1179 |
pszWKT = poOutputSRS.ExportToPrettyWkt( 0 );
|
1180 |
System.err.println( "Target:\n" + pszWKT ); |
1181 |
System.exit( 1 ); |
1182 |
} |
1183 |
} |
1184 |
|
1185 |
/* -------------------------------------------------------------------- */
|
1186 |
/* Get other info. */
|
1187 |
/* -------------------------------------------------------------------- */
|
1188 |
poSrcFDefn = poSrcLayer.GetLayerDefn(); |
1189 |
|
1190 |
if( poOutputSRS == null ) |
1191 |
poOutputSRS = poSrcLayer.GetSpatialRef(); |
1192 |
|
1193 |
/* -------------------------------------------------------------------- */
|
1194 |
/* Find the layer. */
|
1195 |
/* -------------------------------------------------------------------- */
|
1196 |
|
1197 |
/* GetLayerByName() can instantiate layers that would have been */
|
1198 |
/* 'hidden' otherwise, for example, non-spatial tables in a */
|
1199 |
/* PostGIS-enabled database, so this apparently useless command is */
|
1200 |
/* not useless. (#4012) */
|
1201 |
gdal.PushErrorHandler("CPLQuietErrorHandler");
|
1202 |
poDstLayer = poDstDS.GetLayerByName(pszNewLayerName); |
1203 |
gdal.PopErrorHandler(); |
1204 |
gdal.ErrorReset(); |
1205 |
|
1206 |
int iLayer = -1; |
1207 |
if( poDstLayer != null ) |
1208 |
{ |
1209 |
int nLayerCount = poDstDS.GetLayerCount();
|
1210 |
for( iLayer = 0; iLayer < nLayerCount; iLayer++ ) |
1211 |
{ |
1212 |
Layer poLayer = poDstDS.GetLayer(iLayer); |
1213 |
|
1214 |
if( poLayer != null |
1215 |
&& poLayer.GetName().equals(poDstLayer.GetName()) ) |
1216 |
{ |
1217 |
break;
|
1218 |
} |
1219 |
} |
1220 |
|
1221 |
if (iLayer == nLayerCount)
|
1222 |
/* shouldn't happen with an ideal driver */
|
1223 |
poDstLayer = null;
|
1224 |
} |
1225 |
|
1226 |
/* -------------------------------------------------------------------- */
|
1227 |
/* If the user requested overwrite, and we have the layer in */
|
1228 |
/* question we need to delete it now so it will get recreated */
|
1229 |
/* (overwritten). */
|
1230 |
/* -------------------------------------------------------------------- */
|
1231 |
if( poDstLayer != null && bOverwrite ) |
1232 |
{ |
1233 |
if( poDstDS.DeleteLayer( iLayer ) != 0 ) |
1234 |
{ |
1235 |
System.err.println(
|
1236 |
"DeleteLayer() failed when overwrite requested." );
|
1237 |
return false; |
1238 |
} |
1239 |
poDstLayer = null;
|
1240 |
} |
1241 |
|
1242 |
/* -------------------------------------------------------------------- */
|
1243 |
/* If the layer does not exist, then create it. */
|
1244 |
/* -------------------------------------------------------------------- */
|
1245 |
if( poDstLayer == null ) |
1246 |
{ |
1247 |
if( eGType == -2 ) |
1248 |
{ |
1249 |
eGType = poSrcFDefn.GetGeomType(); |
1250 |
|
1251 |
if ( bExplodeCollections )
|
1252 |
{ |
1253 |
int n25DBit = eGType & ogr.wkb25DBit;
|
1254 |
if (wkbFlatten(eGType) == ogr.wkbMultiPoint)
|
1255 |
{ |
1256 |
eGType = ogr.wkbPoint | n25DBit; |
1257 |
} |
1258 |
else if (wkbFlatten(eGType) == ogr.wkbMultiLineString) |
1259 |
{ |
1260 |
eGType = ogr.wkbLineString | n25DBit; |
1261 |
} |
1262 |
else if (wkbFlatten(eGType) == ogr.wkbMultiPolygon) |
1263 |
{ |
1264 |
eGType = ogr.wkbPolygon | n25DBit; |
1265 |
} |
1266 |
else if (wkbFlatten(eGType) == ogr.wkbGeometryCollection) |
1267 |
{ |
1268 |
eGType = ogr.wkbUnknown | n25DBit; |
1269 |
} |
1270 |
} |
1271 |
|
1272 |
if ( pszZField != null ) |
1273 |
eGType |= ogr.wkb25DBit; |
1274 |
} |
1275 |
|
1276 |
if( poDstDS.TestCapability( ogr.ODsCCreateLayer ) == false) |
1277 |
{ |
1278 |
System.err.println(
|
1279 |
"Layer " + pszNewLayerName + "not found, and CreateLayer not supported by driver."); |
1280 |
return false; |
1281 |
} |
1282 |
|
1283 |
gdal.ErrorReset(); |
1284 |
|
1285 |
poDstLayer = poDstDS.CreateLayer( pszNewLayerName, poOutputSRS, |
1286 |
eGType, papszLCO ); |
1287 |
|
1288 |
if( poDstLayer == null ) |
1289 |
return false; |
1290 |
|
1291 |
bAppend = false;
|
1292 |
} |
1293 |
|
1294 |
/* -------------------------------------------------------------------- */
|
1295 |
/* Otherwise we will append to it, if append was requested. */
|
1296 |
/* -------------------------------------------------------------------- */
|
1297 |
else if( !bAppend ) |
1298 |
{ |
1299 |
System.err.println("FAILED: Layer " + pszNewLayerName + "already exists, and -append not specified.\n" + |
1300 |
" Consider using -append, or -overwrite.");
|
1301 |
return false; |
1302 |
} |
1303 |
else
|
1304 |
{ |
1305 |
if( papszLCO.size() > 0 ) |
1306 |
{ |
1307 |
System.err.println("WARNING: Layer creation options ignored since an existing layer is\n" + |
1308 |
" being appended to." );
|
1309 |
} |
1310 |
} |
1311 |
|
1312 |
/* -------------------------------------------------------------------- */
|
1313 |
/* Add fields. Default to copy all field. */
|
1314 |
/* If only a subset of all fields requested, then output only */
|
1315 |
/* the selected fields, and in the order that they were */
|
1316 |
/* selected. */
|
1317 |
/* -------------------------------------------------------------------- */
|
1318 |
int iField;
|
1319 |
|
1320 |
/* Initialize the index-to-index map to -1's */
|
1321 |
int nSrcFieldCount = poSrcFDefn.GetFieldCount();
|
1322 |
int[] panMap = new int [nSrcFieldCount]; |
1323 |
for( iField=0; iField < nSrcFieldCount; iField++) |
1324 |
panMap[iField] = -1;
|
1325 |
|
1326 |
FeatureDefn poDstFDefn = poDstLayer.GetLayerDefn(); |
1327 |
|
1328 |
if (papszSelFields != null && !bAppend ) |
1329 |
{ |
1330 |
int nDstFieldCount = 0; |
1331 |
if (poDstFDefn != null) |
1332 |
nDstFieldCount = poDstFDefn.GetFieldCount(); |
1333 |
|
1334 |
for( iField=0; iField < papszSelFields.size(); iField++) |
1335 |
{ |
1336 |
int iSrcField = poSrcFDefn.GetFieldIndex((String)papszSelFields.get(iField)); |
1337 |
if (iSrcField >= 0) |
1338 |
{ |
1339 |
FieldDefn poSrcFieldDefn = poSrcFDefn.GetFieldDefn(iSrcField); |
1340 |
FieldDefn oFieldDefn = new FieldDefn( poSrcFieldDefn.GetNameRef(),
|
1341 |
poSrcFieldDefn.GetFieldType() ); |
1342 |
oFieldDefn.SetWidth( poSrcFieldDefn.GetWidth() ); |
1343 |
oFieldDefn.SetPrecision( poSrcFieldDefn.GetPrecision() ); |
1344 |
|
1345 |
if (papszFieldTypesToString != null && |
1346 |
(CSLFindString(papszFieldTypesToString, "All") != -1 || |
1347 |
CSLFindString(papszFieldTypesToString, |
1348 |
ogr.GetFieldTypeName(poSrcFDefn.GetFieldDefn(iSrcField).GetFieldType())) != -1))
|
1349 |
oFieldDefn.SetType(ogr.OFTString); |
1350 |
|
1351 |
/* The field may have been already created at layer creation */
|
1352 |
int iDstField = -1; |
1353 |
if (poDstFDefn != null) |
1354 |
iDstField = poDstFDefn.GetFieldIndex(oFieldDefn.GetNameRef()); |
1355 |
if (iDstField >= 0) |
1356 |
{ |
1357 |
panMap[iSrcField] = iDstField; |
1358 |
} |
1359 |
else if (poDstLayer.CreateField( oFieldDefn ) == 0) |
1360 |
{ |
1361 |
/* now that we've created a field, GetLayerDefn() won't return NULL */
|
1362 |
if (poDstFDefn == null) |
1363 |
poDstFDefn = poDstLayer.GetLayerDefn(); |
1364 |
|
1365 |
/* Sanity check : if it fails, the driver is buggy */
|
1366 |
if (poDstFDefn != null && |
1367 |
poDstFDefn.GetFieldCount() != nDstFieldCount + 1)
|
1368 |
{ |
1369 |
System.err.println(
|
1370 |
"The output driver has claimed to have added the " + oFieldDefn.GetNameRef() + " field, but it did not!"); |
1371 |
} |
1372 |
else
|
1373 |
{ |
1374 |
panMap[iSrcField] = nDstFieldCount; |
1375 |
nDstFieldCount ++; |
1376 |
} |
1377 |
} |
1378 |
|
1379 |
} |
1380 |
else
|
1381 |
{ |
1382 |
System.err.println("Field '" + (String)papszSelFields.get(iField) + "' not found in source layer."); |
1383 |
if( !bSkipFailures )
|
1384 |
return false; |
1385 |
} |
1386 |
} |
1387 |
|
1388 |
/* -------------------------------------------------------------------- */
|
1389 |
/* Use SetIgnoredFields() on source layer if available */
|
1390 |
/* -------------------------------------------------------------------- */
|
1391 |
|
1392 |
/* Here we differ from the ogr2ogr.cpp implementation since the OGRFeatureQuery */
|
1393 |
/* isn't mapped to swig. So in that case just don't use SetIgnoredFields() */
|
1394 |
/* to avoid issue raised in #4015 */
|
1395 |
if (poSrcLayer.TestCapability(ogr.OLCIgnoreFields) && pszWHERE == null) |
1396 |
{ |
1397 |
int iSrcField;
|
1398 |
Vector papszIgnoredFields = new Vector(); |
1399 |
for(iSrcField=0;iSrcField<nSrcFieldCount;iSrcField++) |
1400 |
{ |
1401 |
String pszFieldName =
|
1402 |
poSrcFDefn.GetFieldDefn(iSrcField).GetNameRef(); |
1403 |
boolean bFieldRequested = false; |
1404 |
for( iField=0; iField < papszSelFields.size(); iField++) |
1405 |
{ |
1406 |
if (pszFieldName.equalsIgnoreCase((String)papszSelFields.get(iField))) |
1407 |
{ |
1408 |
bFieldRequested = true;
|
1409 |
break;
|
1410 |
} |
1411 |
} |
1412 |
|
1413 |
if (pszZField != null && pszFieldName.equalsIgnoreCase(pszZField)) |
1414 |
bFieldRequested = true;
|
1415 |
|
1416 |
/* If source field not requested, add it to ignored files list */
|
1417 |
if (!bFieldRequested)
|
1418 |
papszIgnoredFields.addElement(pszFieldName); |
1419 |
} |
1420 |
poSrcLayer.SetIgnoredFields(papszIgnoredFields); |
1421 |
} |
1422 |
} |
1423 |
else if( !bAppend ) |
1424 |
{ |
1425 |
int nDstFieldCount = 0; |
1426 |
if (poDstFDefn != null) |
1427 |
nDstFieldCount = poDstFDefn.GetFieldCount(); |
1428 |
for( iField = 0; iField < nSrcFieldCount; iField++ ) |
1429 |
{ |
1430 |
FieldDefn poSrcFieldDefn = poSrcFDefn.GetFieldDefn(iField); |
1431 |
FieldDefn oFieldDefn = new FieldDefn( poSrcFieldDefn.GetNameRef(),
|
1432 |
poSrcFieldDefn.GetFieldType() ); |
1433 |
oFieldDefn.SetWidth( poSrcFieldDefn.GetWidth() ); |
1434 |
oFieldDefn.SetPrecision( poSrcFieldDefn.GetPrecision() ); |
1435 |
|
1436 |
if (papszFieldTypesToString != null && |
1437 |
(CSLFindString(papszFieldTypesToString, "All") != -1 || |
1438 |
CSLFindString(papszFieldTypesToString, |
1439 |
ogr.GetFieldTypeName(poSrcFDefn.GetFieldDefn(iField).GetFieldType())) != -1))
|
1440 |
oFieldDefn.SetType(ogr.OFTString); |
1441 |
|
1442 |
/* The field may have been already created at layer creation */
|
1443 |
int iDstField = -1; |
1444 |
if (poDstFDefn != null) |
1445 |
iDstField = poDstFDefn.GetFieldIndex(oFieldDefn.GetNameRef()); |
1446 |
if (iDstField >= 0) |
1447 |
{ |
1448 |
panMap[iField] = iDstField; |
1449 |
} |
1450 |
else if (poDstLayer.CreateField( oFieldDefn ) == 0) |
1451 |
{ |
1452 |
/* now that we've created a field, GetLayerDefn() won't return NULL */
|
1453 |
if (poDstFDefn == null) |
1454 |
poDstFDefn = poDstLayer.GetLayerDefn(); |
1455 |
|
1456 |
/* Sanity check : if it fails, the driver is buggy */
|
1457 |
if (poDstFDefn != null && |
1458 |
poDstFDefn.GetFieldCount() != nDstFieldCount + 1)
|
1459 |
{ |
1460 |
System.err.println(
|
1461 |
"The output driver has claimed to have added the " + oFieldDefn.GetNameRef() + " field, but it did not!"); |
1462 |
} |
1463 |
else
|
1464 |
{ |
1465 |
panMap[iField] = nDstFieldCount; |
1466 |
nDstFieldCount ++; |
1467 |
} |
1468 |
} |
1469 |
} |
1470 |
} |
1471 |
else
|
1472 |
{ |
1473 |
/* For an existing layer, build the map by fetching the index in the destination */
|
1474 |
/* layer for each source field */
|
1475 |
if (poDstFDefn == null) |
1476 |
{ |
1477 |
System.err.println("poDstFDefn == NULL.\n" ); |
1478 |
return false; |
1479 |
} |
1480 |
|
1481 |
for( iField = 0; iField < nSrcFieldCount; iField++ ) |
1482 |
{ |
1483 |
FieldDefn poSrcFieldDefn = poSrcFDefn.GetFieldDefn(iField); |
1484 |
int iDstField = poDstFDefn.GetFieldIndex(poSrcFieldDefn.GetNameRef());
|
1485 |
if (iDstField >= 0) |
1486 |
panMap[iField] = iDstField; |
1487 |
} |
1488 |
} |
1489 |
|
1490 |
/* -------------------------------------------------------------------- */
|
1491 |
/* Transfer features. */
|
1492 |
/* -------------------------------------------------------------------- */
|
1493 |
Feature poFeature; |
1494 |
int nFeaturesInTransaction = 0; |
1495 |
long nCount = 0; |
1496 |
|
1497 |
int iSrcZField = -1; |
1498 |
if (pszZField != null) |
1499 |
{ |
1500 |
iSrcZField = poSrcFDefn.GetFieldIndex(pszZField); |
1501 |
} |
1502 |
|
1503 |
poSrcLayer.ResetReading(); |
1504 |
|
1505 |
if( nGroupTransactions > 0) |
1506 |
poDstLayer.StartTransaction(); |
1507 |
|
1508 |
while( true ) |
1509 |
{ |
1510 |
Feature poDstFeature = null;
|
1511 |
|
1512 |
if( nFIDToFetch != OGRNullFID )
|
1513 |
{ |
1514 |
// Only fetch feature on first pass.
|
1515 |
if( nFeaturesInTransaction == 0 ) |
1516 |
poFeature = poSrcLayer.GetFeature(nFIDToFetch); |
1517 |
else
|
1518 |
poFeature = null;
|
1519 |
} |
1520 |
else
|
1521 |
poFeature = poSrcLayer.GetNextFeature(); |
1522 |
|
1523 |
if( poFeature == null ) |
1524 |
break;
|
1525 |
|
1526 |
int nParts = 0; |
1527 |
int nIters = 1; |
1528 |
if (bExplodeCollections)
|
1529 |
{ |
1530 |
Geometry poSrcGeometry = poFeature.GetGeometryRef(); |
1531 |
if (poSrcGeometry != null) |
1532 |
{ |
1533 |
switch (wkbFlatten(poSrcGeometry.GetGeometryType()))
|
1534 |
{ |
1535 |
case ogr.wkbMultiPoint:
|
1536 |
case ogr.wkbMultiLineString:
|
1537 |
case ogr.wkbMultiPolygon:
|
1538 |
case ogr.wkbGeometryCollection:
|
1539 |
nParts = poSrcGeometry.GetGeometryCount(); |
1540 |
nIters = nParts; |
1541 |
if (nIters == 0) |
1542 |
nIters = 1;
|
1543 |
default:
|
1544 |
break;
|
1545 |
} |
1546 |
} |
1547 |
} |
1548 |
|
1549 |
for(int iPart = 0; iPart < nIters; iPart++) |
1550 |
{ |
1551 |
|
1552 |
if( ++nFeaturesInTransaction == nGroupTransactions )
|
1553 |
{ |
1554 |
poDstLayer.CommitTransaction(); |
1555 |
poDstLayer.StartTransaction(); |
1556 |
nFeaturesInTransaction = 0;
|
1557 |
} |
1558 |
|
1559 |
gdal.ErrorReset(); |
1560 |
poDstFeature = new Feature( poDstLayer.GetLayerDefn() );
|
1561 |
|
1562 |
if( poDstFeature.SetFromWithMap( poFeature, 1, panMap ) != 0 ) |
1563 |
{ |
1564 |
if( nGroupTransactions > 0) |
1565 |
poDstLayer.CommitTransaction(); |
1566 |
|
1567 |
System.err.println(
|
1568 |
"Unable to translate feature " + poFeature.GetFID() + " from layer " + |
1569 |
poSrcFDefn.GetName() ); |
1570 |
|
1571 |
poFeature.delete(); |
1572 |
poFeature = null;
|
1573 |
poDstFeature.delete(); |
1574 |
poDstFeature = null;
|
1575 |
return false; |
1576 |
} |
1577 |
|
1578 |
if( bPreserveFID )
|
1579 |
poDstFeature.SetFID( poFeature.GetFID() ); |
1580 |
|
1581 |
Geometry poDstGeometry = poDstFeature.GetGeometryRef(); |
1582 |
if (poDstGeometry != null) |
1583 |
{ |
1584 |
if (nParts > 0) |
1585 |
{ |
1586 |
/* For -explodecollections, extract the iPart(th) of the geometry */
|
1587 |
Geometry poPart = poDstGeometry.GetGeometryRef(iPart).Clone(); |
1588 |
poDstFeature.SetGeometryDirectly(poPart); |
1589 |
poDstGeometry = poPart; |
1590 |
} |
1591 |
|
1592 |
if (iSrcZField != -1) |
1593 |
{ |
1594 |
SetZ(poDstGeometry, poFeature.GetFieldAsDouble(iSrcZField)); |
1595 |
/* This will correct the coordinate dimension to 3 */
|
1596 |
Geometry poDupGeometry = poDstGeometry.Clone(); |
1597 |
poDstFeature.SetGeometryDirectly(poDupGeometry); |
1598 |
poDstGeometry = poDupGeometry; |
1599 |
} |
1600 |
|
1601 |
if (eGeomOp == GeomOperation.SEGMENTIZE)
|
1602 |
{ |
1603 |
/*if (poDstFeature.GetGeometryRef() != null && dfGeomOpParam > 0)
|
1604 |
poDstFeature.GetGeometryRef().segmentize(dfGeomOpParam);*/
|
1605 |
} |
1606 |
else if (eGeomOp == GeomOperation.SIMPLIFY_PRESERVE_TOPOLOGY && dfGeomOpParam > 0) |
1607 |
{ |
1608 |
Geometry poNewGeom = poDstGeometry.SimplifyPreserveTopology(dfGeomOpParam); |
1609 |
if (poNewGeom != null) |
1610 |
{ |
1611 |
poDstFeature.SetGeometryDirectly(poNewGeom); |
1612 |
poDstGeometry = poNewGeom; |
1613 |
} |
1614 |
} |
1615 |
|
1616 |
if (poClipSrc != null) |
1617 |
{ |
1618 |
Geometry poClipped = poDstGeometry.Intersection(poClipSrc); |
1619 |
if (poClipped == null || poClipped.IsEmpty()) |
1620 |
{ |
1621 |
/* Report progress */
|
1622 |
nCount ++; |
1623 |
if (pfnProgress != null) |
1624 |
pfnProgress.run(nCount * 1.0 / nCountLayerFeatures, ""); |
1625 |
poDstFeature.delete(); |
1626 |
continue;
|
1627 |
} |
1628 |
poDstFeature.SetGeometryDirectly(poClipped); |
1629 |
poDstGeometry = poClipped; |
1630 |
} |
1631 |
|
1632 |
if( poCT != null ) |
1633 |
{ |
1634 |
eErr = poDstGeometry.Transform( poCT ); |
1635 |
if( eErr != 0 ) |
1636 |
{ |
1637 |
if( nGroupTransactions > 0) |
1638 |
poDstLayer.CommitTransaction(); |
1639 |
|
1640 |
System.err.println("Failed to reproject feature" + poFeature.GetFID() + " (geometry probably out of source or destination SRS)."); |
1641 |
if( !bSkipFailures )
|
1642 |
{ |
1643 |
poFeature.delete(); |
1644 |
poFeature = null;
|
1645 |
poDstFeature.delete(); |
1646 |
poDstFeature = null;
|
1647 |
return false; |
1648 |
} |
1649 |
} |
1650 |
} |
1651 |
else if (poOutputSRS != null) |
1652 |
{ |
1653 |
poDstGeometry.AssignSpatialReference(poOutputSRS); |
1654 |
} |
1655 |
|
1656 |
if (poClipDst != null) |
1657 |
{ |
1658 |
Geometry poClipped = poDstGeometry.Intersection(poClipDst); |
1659 |
if (poClipped == null || poClipped.IsEmpty()) |
1660 |
{ |
1661 |
/* Report progress */
|
1662 |
nCount ++; |
1663 |
if (pfnProgress != null) |
1664 |
pfnProgress.run(nCount * 1.0 / nCountLayerFeatures, ""); |
1665 |
poDstFeature.delete(); |
1666 |
continue;
|
1667 |
} |
1668 |
poDstFeature.SetGeometryDirectly(poClipped); |
1669 |
poDstGeometry = poClipped; |
1670 |
} |
1671 |
|
1672 |
if( bForceToPolygon )
|
1673 |
{ |
1674 |
poDstFeature.SetGeometryDirectly(ogr.ForceToPolygon(poDstGeometry)); |
1675 |
} |
1676 |
|
1677 |
else if( bForceToMultiPolygon ) |
1678 |
{ |
1679 |
poDstFeature.SetGeometryDirectly(ogr.ForceToMultiPolygon(poDstGeometry)); |
1680 |
} |
1681 |
|
1682 |
else if ( bForceToMultiLineString ) |
1683 |
{ |
1684 |
poDstFeature.SetGeometryDirectly(ogr.ForceToMultiLineString(poDstGeometry)); |
1685 |
} |
1686 |
} |
1687 |
|
1688 |
gdal.ErrorReset(); |
1689 |
if( poDstLayer.CreateFeature( poDstFeature ) != 0 |
1690 |
&& !bSkipFailures ) |
1691 |
{ |
1692 |
if( nGroupTransactions > 0 ) |
1693 |
poDstLayer.RollbackTransaction(); |
1694 |
|
1695 |
poDstFeature.delete(); |
1696 |
poDstFeature = null;
|
1697 |
return false; |
1698 |
} |
1699 |
|
1700 |
poDstFeature.delete(); |
1701 |
poDstFeature = null;
|
1702 |
} |
1703 |
|
1704 |
poFeature.delete(); |
1705 |
poFeature = null;
|
1706 |
|
1707 |
/* Report progress */
|
1708 |
nCount ++; |
1709 |
if (pfnProgress != null) |
1710 |
pfnProgress.run(nCount * 1.0 / nCountLayerFeatures, ""); |
1711 |
|
1712 |
} |
1713 |
|
1714 |
if( nGroupTransactions > 0 ) |
1715 |
poDstLayer.CommitTransaction(); |
1716 |
|
1717 |
return true; |
1718 |
} |
1719 |
} |