File: occi_featureset.cpp

package info (click to toggle)
mapnik 2.0.0%2Bds1-3
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 35,496 kB
  • sloc: cpp: 91,793; python: 6,051; xml: 3,528; sh: 848; makefile: 70; lisp: 10
file content (490 lines) | stat: -rw-r--r-- 18,072 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
/*****************************************************************************
 * 
 * This file is part of Mapnik (c++ mapping toolkit)
 *
 * Copyright (C) 2007 Artem Pavlenko
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 *
 *****************************************************************************/
//$Id$

// mapnik
#include <mapnik/global.hpp>
#include <mapnik/datasource.hpp>
#include <mapnik/box2d.hpp>
#include <mapnik/geometry.hpp>
#include <mapnik/feature.hpp>
#include <mapnik/feature_layer_desc.hpp>
#include <mapnik/wkb.hpp>
#include <mapnik/unicode.hpp>
#include <mapnik/feature_factory.hpp>

// ogr
#include "occi_featureset.hpp"

using mapnik::query;
using mapnik::box2d;
using mapnik::CoordTransform;
using mapnik::Feature;
using mapnik::feature_ptr;
using mapnik::geometry_type;
using mapnik::geometry_utils;
using mapnik::transcoder;
using mapnik::datasource_exception;
using mapnik::feature_factory;

using oracle::occi::Connection;
using oracle::occi::Statement;
using oracle::occi::ResultSet;
using oracle::occi::StatelessConnectionPool;
using oracle::occi::MetaData;
using oracle::occi::SQLException;
using oracle::occi::Type;
using oracle::occi::Number;
using oracle::occi::Blob;

occi_featureset::occi_featureset(StatelessConnectionPool* pool,
                                 Connection* conn,
                                 std::string const& sqlstring,
                                 std::string const& encoding,
                                 bool multiple_geometries,
                                 bool use_connection_pool,
                                 unsigned prefetch_rows,
                                 unsigned num_attrs)
   : tr_(new transcoder(encoding)),
     multiple_geometries_(multiple_geometries),
     num_attrs_(num_attrs),
     feature_id_(1)
{
    if (use_connection_pool)
        conn_.set_pool(pool);
    else
        conn_.set_connection(conn, false);

    try
    {
        rs_ = conn_.execute_query (sqlstring, prefetch_rows);
    }
    catch (SQLException &ex)
    {
        std::clog << "OCCI Plugin: error processing " << sqlstring << " : " << ex.getMessage() << std::endl;
    }
}

occi_featureset::~occi_featureset()
{
}

feature_ptr occi_featureset::next()
{
    if (rs_ && rs_->next())
    {
        feature_ptr feature(feature_factory::create(feature_id_));
        ++feature_id_;

        boost::scoped_ptr<SDOGeometry> geom (dynamic_cast<SDOGeometry*> (rs_->getObject(1)));
        if (geom.get())
        {
            convert_geometry (geom.get(), feature, multiple_geometries_);
        }

        std::vector<MetaData> listOfColumns = rs_->getColumnListMetaData();

        for (unsigned int i = 1; i < listOfColumns.size(); ++i)
        {
            MetaData columnObj = listOfColumns[i];

            std::string fld_name = columnObj.getString(MetaData::ATTR_NAME);
            int type_oid = columnObj.getInt(MetaData::ATTR_DATA_TYPE);

            /*   
            int type_code = columnObj.getInt(MetaData::ATTR_TYPECODE);
            if (type_code == OCCI_TYPECODE_OBJECT)
            {
                continue;
            }
            */

            switch (type_oid)
            {
                case oracle::occi::OCCIBOOL:
                case oracle::occi::OCCIINT:
                case oracle::occi::OCCIUNSIGNED_INT:
                case oracle::occi::OCCIROWID:
                    boost::put(*feature,fld_name,rs_->getInt (i + 1));
                    break;
                case oracle::occi::OCCIFLOAT:
                case oracle::occi::OCCIBFLOAT:
                case oracle::occi::OCCIDOUBLE:
                case oracle::occi::OCCIBDOUBLE:
                case oracle::occi::OCCINUMBER:
                case oracle::occi::OCCI_SQLT_NUM:
                    boost::put(*feature,fld_name,rs_->getDouble (i + 1));
                    break;
                case oracle::occi::OCCICHAR:
                case oracle::occi::OCCISTRING:
                case oracle::occi::OCCI_SQLT_AFC:
                case oracle::occi::OCCI_SQLT_AVC:
                case oracle::occi::OCCI_SQLT_CHR:
                case oracle::occi::OCCI_SQLT_LVC:
                case oracle::occi::OCCI_SQLT_RDD:
                case oracle::occi::OCCI_SQLT_STR:
                case oracle::occi::OCCI_SQLT_VCS:
                case oracle::occi::OCCI_SQLT_VNU:
                case oracle::occi::OCCI_SQLT_VBI:
                case oracle::occi::OCCI_SQLT_VST:
                    boost::put(*feature,fld_name,(UnicodeString) tr_->transcode (rs_->getString (i + 1).c_str()));
                    break;
                case oracle::occi::OCCIDATE:
                case oracle::occi::OCCITIMESTAMP:
                case oracle::occi::OCCIINTERVALDS:
                case oracle::occi::OCCIINTERVALYM:
                case oracle::occi::OCCI_SQLT_DAT:
                case oracle::occi::OCCI_SQLT_DATE:
                case oracle::occi::OCCI_SQLT_TIME:
                case oracle::occi::OCCI_SQLT_TIME_TZ:
                case oracle::occi::OCCI_SQLT_TIMESTAMP:
                case oracle::occi::OCCI_SQLT_TIMESTAMP_LTZ:
                case oracle::occi::OCCI_SQLT_TIMESTAMP_TZ:
                case oracle::occi::OCCI_SQLT_INTERVAL_YM:
                case oracle::occi::OCCI_SQLT_INTERVAL_DS:
                case oracle::occi::OCCIANYDATA:
                case oracle::occi::OCCIBLOB:
                case oracle::occi::OCCIBFILE:
                case oracle::occi::OCCIBYTES:
                case oracle::occi::OCCICLOB:
                case oracle::occi::OCCIVECTOR:
                case oracle::occi::OCCIMETADATA:
                case oracle::occi::OCCIPOBJECT:
                case oracle::occi::OCCIREF:
                case oracle::occi::OCCIREFANY:
                case oracle::occi::OCCISTREAM:
                case oracle::occi::OCCICURSOR:
                case oracle::occi::OCCI_SQLT_FILE:
                case oracle::occi::OCCI_SQLT_CFILE:
                case oracle::occi::OCCI_SQLT_REF:
                case oracle::occi::OCCI_SQLT_CLOB:
                case oracle::occi::OCCI_SQLT_BLOB:
                case oracle::occi::OCCI_SQLT_RSET:
#ifdef MAPNIK_DEBUG
                    std::clog << "OCCI Plugin: unsupported datatype " << occi_enums::resolve_datatype(type_oid)
                        << " (type_oid=" << type_oid << ")" << std::endl;
#endif
                    break;
                default: // shouldn't get here
#ifdef MAPNIK_DEBUG
                    std::clog << "OCCI Plugin: unknown datatype (type_oid=" << type_oid << ")" << std::endl;
#endif
                    break;
            }    
        }
        
        return feature;
    }

    return feature_ptr();
}


void occi_featureset::convert_geometry (SDOGeometry* geom, feature_ptr feature, bool multiple_geometries)
{
    int gtype = (int) geom->getSdo_gtype();
    int dimensions = gtype / 1000;
    int lrsvalue = (gtype - dimensions * 1000) / 100;
    int geomtype = (gtype - dimensions * 1000 - lrsvalue * 100); 

#if 0
    clog << "-----------Geometry Object ------------" << endl;
    clog << "SDO GTYPE = " << gtype << endl;
    clog << "SDO DIMENSIONS = " << dimensions << endl;
    clog << "SDO LRS = " << lrsvalue << endl;
    clog << "SDO GEOMETRY TYPE = " << geomtype << endl;

    Number sdo_srid = geom->getSdo_srid();
    if (sdo_srid.isNull())
        clog << "SDO SRID = " << "Null" << endl;
    else
        clog << "SDO SRID = " << (int) sdo_srid << endl;
#endif

    const std::vector<Number>& elem_info = geom->getSdo_elem_info();
    const std::vector<Number>& ordinates = geom->getSdo_ordinates();
    const int ordinates_size = (int) ordinates.size();

    switch (geomtype)
    {
    case SDO_GTYPE_POINT:
        {
            SDOPointType* sdopoint = geom->getSdo_point();
            if (sdopoint && ! sdopoint->isNull())
            {
                geometry_type* point = new geometry_type(mapnik::Point);
                point->move_to (sdopoint->getX(), sdopoint->getY());
                feature->add_geometry (point);
            }
        }
        break;
    case SDO_GTYPE_LINE:
        {
            if (ordinates_size >= dimensions)
            {
                convert_ordinates (feature,
                                   mapnik::LineString,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   true,   // is_single_geom
                                   false,  // is_point_type
                                   false); // multiple_geometries
            }
        }
        break;
    case SDO_GTYPE_POLYGON:
        {
            if (ordinates_size >= dimensions)
            {
                convert_ordinates (feature,
                                   mapnik::Polygon,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   true,   // is_single_geom
                                   false,  // is_point_type
                                   false); // multiple_geometries
            }
        }
        break;
    case SDO_GTYPE_MULTIPOINT:
        {
            if (ordinates_size >= dimensions)
            {
                const bool is_single_geom = false;
                const bool is_point_type = true;

                // Todo - force using true as multiple_geometries until we have proper multipoint handling
                // http://trac.mapnik.org/ticket/458
            
                convert_ordinates (feature,
                                   mapnik::Point,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   false,  // is_single_geom
                                   true,   // is_point_type
                                   true);  // multiple_geometries
            }
        }
        break;
    case SDO_GTYPE_MULTILINE:
        {
            if (ordinates_size >= dimensions)
            {
                convert_ordinates (feature,
                                   mapnik::LineString,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   false,  // is_single_geom
                                   false,  // is_point_type
                                   multiple_geometries);
            }

        }
        break;
    case SDO_GTYPE_MULTIPOLYGON:
        {
            if (ordinates_size >= dimensions)
            {
                convert_ordinates (feature,
                                   mapnik::Polygon,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   false, // is_single_geom
                                   false, // is_point_type
                                   multiple_geometries);
            }
        
        }
        break;
    case SDO_GTYPE_COLLECTION:
        {
            if (ordinates_size >= dimensions)
            {
                const bool is_single_geom = false;
                const bool is_point_type = false;
            
                convert_ordinates (feature,
                                   mapnik::Polygon,
                                   elem_info,
                                   ordinates,
                                   dimensions,
                                   false, // is_single_geom,
                                   false, // is_point_type
                                   multiple_geometries);
            }
        }
        break;
    case SDO_GTYPE_UNKNOWN:
    default:
#ifdef MAPNIK_DEBUG
        std::clog << "OCCI Plugin: unknown <occi> " << occi_enums::resolve_gtype(geomtype)
            << "(gtype=" << gtype << ")" << std::endl;
#endif
        break;
    }
}

void occi_featureset::convert_ordinates (mapnik::feature_ptr feature,
                                         const mapnik::eGeomType& geom_type,
                                         const std::vector<Number>& elem_info,
                                         const std::vector<Number>& ordinates,
                                         const int dimensions,
                                         const bool is_single_geom,
                                         const bool is_point_geom,
                                         const bool multiple_geometries)
{
    const int elem_size = elem_info.size();
    const int ord_size = ordinates.size();
    
    if (elem_size >= 0)
    {
        int offset = elem_info [0];
        int etype = elem_info [1];
        int interp = elem_info [2];

        if (! is_single_geom && elem_size > SDO_ELEM_INFO_SIZE)
        {
            geometry_type* geom = multiple_geometries ? 0 : new geometry_type(geom_type);
            if (geom) geom->set_capacity (ord_size);
        
            for (int i = SDO_ELEM_INFO_SIZE; i < elem_size; i+=3)
            {
                int next_offset = elem_info [i];
                int next_etype = elem_info [i + 1];
                int next_interp = elem_info [i + 2];
                bool is_linear_element = true;
                bool is_unknown_etype = false;
                mapnik::eGeomType gtype = mapnik::Point;

                switch (etype)
                {
                case SDO_ETYPE_POINT:
                    if (interp == SDO_INTERPRETATION_POINT) {}
                    if (interp > SDO_INTERPRETATION_POINT)  {}
                    gtype = mapnik::Point;
                    break;
                
                case SDO_ETYPE_LINESTRING:
                    if (interp == SDO_INTERPRETATION_STRAIGHT) {}
                    if (interp == SDO_INTERPRETATION_CIRCULAR) {}
                    gtype = mapnik::LineString;
                    break;

                case SDO_ETYPE_POLYGON:
                case SDO_ETYPE_POLYGON_INTERIOR:
                    if (interp == SDO_INTERPRETATION_STRAIGHT)  {}
                    if (interp == SDO_INTERPRETATION_CIRCULAR)  {}
                    if (interp == SDO_INTERPRETATION_RECTANGLE) {}
                    if (interp == SDO_INTERPRETATION_CIRCLE)    {}
                    gtype = mapnik::Polygon;
                    break;

                case SDO_ETYPE_COMPOUND_LINESTRING:
                case SDO_ETYPE_COMPOUND_POLYGON:
                case SDO_ETYPE_COMPOUND_POLYGON_INTERIOR:
                    // interp = next ETYPE to consider
                    is_linear_element = false;
                    gtype = mapnik::Polygon;
                    break;

                case SDO_ETYPE_UNKNOWN:    // unknown
                default:
                    is_unknown_etype = true;
                    break;
                }

                if (is_unknown_etype)
                    break;

                if (is_linear_element)
                {
                    if (multiple_geometries)
                    {
                        if (geom)
                            feature->add_geometry (geom);
                            
                        geom = new geometry_type(gtype);
                        geom->set_capacity ((next_offset - 1) - (offset - 1 - dimensions));
                    }

                    fill_geometry_type (geom,
                                        offset - 1,
                                        next_offset - 1,
                                        ordinates,
                                        dimensions,
                                        is_point_geom);
                }

                offset = next_offset;
                etype = next_etype;
                interp = next_interp;
            }
            
            if (geom)
            {
                feature->add_geometry (geom);
                geom = 0;
            }
        }
        else
        {
            geometry_type * geom = new geometry_type(geom_type);
            geom->set_capacity (ord_size);
        
            fill_geometry_type (geom,
                                offset - 1,
                                ord_size,
                                ordinates,
                                dimensions,
                                is_point_geom);

            feature->add_geometry (geom);
        }
    }
}

void occi_featureset::fill_geometry_type (geometry_type * geom,
                                          const int real_offset,
                                          const int next_offset,
                                          const std::vector<Number>& ordinates,
                                          const int dimensions,
                                          const bool is_point_geom)
{
    geom->move_to ((double) ordinates[real_offset], (double) ordinates[real_offset + 1]);

    if (is_point_geom)
    {
        for (int p = real_offset + dimensions; p < next_offset; p += dimensions)
            geom->move_to ((double) ordinates[p], (double) ordinates[p + 1]);
    }
    else
    {
        for (int p = real_offset + dimensions; p < next_offset; p += dimensions)
            geom->line_to ((double) ordinates[p], (double) ordinates[p + 1]);
    }
}