<div dir="ltr">Nik,<div><br></div><div>One option I frequently use when uploading shapefiles to PostGIS is to set the -nlt option in ogr2ogr. Setting it as multipolygon/multilinestring/multipoint for a shapefile defined as polygon/linestring/point usually overcomes the problem.</div>
<div><br></div><div>Check if Spatialite allows adding single geometries when the datasource is defined as a multigeometry. If not, you can force the geometry to multigeometry before adding.</div><div><br></div></div><div class="gmail_extra">
<br><br><div class="gmail_quote">On Tue, Mar 4, 2014 at 5:39 AM, Nik Sands <span dir="ltr"><<a href="mailto:nixanz@nixanz.com" target="_blank">nixanz@nixanz.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
So now, instead of the 'OGR_DS_CopyLayer', I'm using the code below to run SQL queries to split each layer from the source data set into multiple layers in the destination (Spatialite) dataset based on the geometry types. (The addNewLayerNames: method will not add empty layers).<br>
<br>
It appears to be working OK for now (at least the maps LOOK right on screen).<br>
<br>
Is this a reasonable approach to overcoming the problem?<br>
<br>
Do I need to include more than just the 6 geometry types in the dictionary/hash below?<br>
<br>
Cheers,<br>
Nik.<br>
<br>
<br>
CODE:<br>
-----<br>
NSDictionary *geomTypeExts = @{<br>
@"POINT" : @"pt",<br>
@"MULTIPOINT" : @"mpt",<br>
@"LINESTRING" : @"ls",<br>
@"MULTILINESTRING" : @"mls",<br>
@"POLYGON" : @"pg",<br>
@"MULTIPOLYGON" : @"mpg",<br>
};<br>
<br>
for ( NSString *geomType in [geomTypeExts allKeys] )<br>
{<br>
NSString *sql = [NSString stringWithFormat:@"SELECT * FROM %@ WHERE OGR_GEOMETRY = '%@'", sourceLayerName, geomType];<br>
NSString *destLayerName = [NSString stringWithFormat:@"%@_%@", sourceLayerName, [geomTypeExts valueForKey:geomType]];<br>
<br>
[self addNewLayerNamed:destLayerName withSQL:sql fromDS:ds toDS:destDS];<br>
<div class="HOEnZb"><div class="h5"> }<br>
<br>
<br>
On 4 Mar 2014, at 9:54 am, Nik Sands <<a href="mailto:nixanz@nixanz.com">nixanz@nixanz.com</a>> wrote:<br>
<br>
> Thanks for this tip. I think you're on the right track with mixing single line strings with multi line strings and single polygons with multi polygons.<br>
><br>
> (I don't think the number of features is causing any problems as it happens for some very small groups of features also.)<br>
><br>
> It appears as though the first time it hits a feature type that is different to all of the others so far, it just bails out silently.<br>
><br>
> So this implies that it's OK to mix feature types (multi/single) in Shapefiles (at least for reading) but not in Spatialite (at least for writing). Is this correct?<br>
><br>
> If this is so, then what is the correct way to convert/copy a layer from Shapefiles to Spatialite?<br>
><br>
> (Note that is seems to work OK for plain SQLite but not for Spatialite).<br>
><br>
><br>
> On 3 Mar 2014, at 4:37 pm, Jukka Rahkonen <<a href="mailto:jukka.rahkonen@mmmtike.fi">jukka.rahkonen@mmmtike.fi</a>> wrote:<br>
><br>
>> Nik Sands <nixanz <at> <a href="http://nixanz.com" target="_blank">nixanz.com</a>> writes:<br>
>><br>
>>><br>
>>> Hi all,<br>
>>><br>
>>> I'm planning to switch from using a standard file format of Shapefiles to<br>
>> Spatialite for internal storage<br>
>>> of spatial data.<br>
>>><br>
>>> When importing from other formats, the app attempts to convert these to<br>
>> SQLite (Spatialite) using<br>
>>> 'OGR_DS_CopyLayer' to copy each layer from the source data set to a new<br>
>> layer in the new SQLite dataset.<br>
>>><br>
>>> However, when converting Shapefiles to SQLite, it only copies across about<br>
>> a third of the features for many<br>
>>> of the layers and therefore the new data set is of no use.<br>
>>><br>
>>> The relevant code and output is below.<br>
>>><br>
>>> I'm at a loss to figure out what is going wrong here and I'd be grateful<br>
>> is somebody could assist in getting<br>
>>> this to work as expected (or am I doing the wrong thing completely?).<br>
>>><br>
>>> NB: This is using Spatialite on iOS (simulator).<br>
>>><br>
>>> When I use ogr2ogr on Mac OS X with SQLite (without Spatialite) it works<br>
>> fine and doesn't miss any features.<br>
>>><br>
>>> Cheers,<br>
>>> Nik.<br>
>>><br>
>>> CODE:<br>
>>><br>
>>> OGRLayerH sourceLayer = OGR_DS_GetLayer(ds, i);<br>
>>> int sourceFeatureCount = OGR_L_GetFeatureCount(sourceLayer, YES);<br>
>>><br>
>>> OGR_DS_CopyLayer(destDS, sourceLayer, OGR_L_GetName(sourceLayer), NULL);<br>
>>><br>
>>> OGRLayerH destLayer = OGR_DS_GetLayerByName(destDS,<br>
>> OGR_L_GetName(sourceLayer));<br>
>>> printf("sourceFeatureCount: %d\n", sourceFeatureCount);<br>
>>> printf(" destFeatureCount: %d\n", OGR_L_GetFeatureCount(destLayer, YES));<br>
>>><br>
>>> OUTPUT:<br>
>>><br>
>>> sourceFeatureCount: 3882<br>
>>> destFeatureCount: 1200<br>
>>><br>
>><br>
>> Hi,<br>
>><br>
>> My guess is that CopyLayer is making 200 feature sized transactions and only<br>
>> 6 transactions are accepted by Spatialite which makes 1200 features. The<br>
>> rest transactions are not accepted for some reason. With Spatialite my guess<br>
>> is that it complains about constraints and some of the features do not share<br>
>> the same geometry type and SRID. Perhaps trying to write multipolygons into<br>
>> a layer that accepts only polygons or vice versa.<br>
>><br>
>> -Jukka Rahkonen-<br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> gdal-dev mailing list<br>
>> <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
>> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
><br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Best regards,<br>Chaitanya kumar CH.<br><br>+91-9494447584<br>17.2416N 80.1426E
</div>