Plik sqlite z NE ma format FDO-OGR, a nie natywną geometrię przestrzenną. Jeśli chcesz wykonać pracę fizyczną, oto sposób na konwersję do bazy danych przestrzennych:
Najpierw utwórz nową, pustą bazę danych przestrzennych (nazywam ją „nev.sqlite”), a następnie w osobnej sesji terminalu otwórz oryginalną natural_earth_vector.sqlite za pomocą spatialite. (Użyłem nowszej wersji 4.1. Nie jestem pewien, czy to zadziała ze starszymi wersjami). Użyj attach
funkcji sqlite, aby połączyć się z nową tabelą nev.sqlite i utworzyć kopie tabel, które chcesz w nowej bazie danych.
Więc:
micha@Wheezy:~$ spatialite natural_earth_vector.sqlite
SpatiaLite version ..: 3.0.0-beta Supported Extensions:
- 'VirtualShape' [direct Shapefile access]
- 'VirtualDbf' [direct DBF access]
- 'VirtualXL' [direct XLS access]
- 'VirtualText' [direct CSV/TXT access]
- 'VirtualNetwork' [Dijkstra shortest path]
- 'RTree' [Spatial Index - R*Tree]
- 'MbrCache' [Spatial Index - MBR cache]
- 'VirtualSpatialIndex' [R*Tree metahandler]
- 'VirtualFDO' [FDO-OGR interoperability]
- 'SpatiaLite' [Spatial SQL - OGC]
PROJ.4 version ......: Rel. 4.7.1, 23 September 2009
GEOS version ........: 3.3.3-CAPI-1.7.4
SQLite version ......: 3.7.13
================ FDO-OGR Spatial Metadata detected ===============
.....
created VirtualFDO table 'fdo_ne_110m_geography_regions_points'
created VirtualFDO table 'fdo_ne_110m_geography_regions_polys'
created VirtualFDO table 'fdo_ne_110m_glaciated_areas'
created VirtualFDO table 'fdo_ne_110m_lakes'
created VirtualFDO table 'fdo_ne_110m_land'
created VirtualFDO table 'fdo_ne_110m_ocean'
created VirtualFDO table 'fdo_ne_110m_rivers_lake_centerlines'
Accessing these fdo_XX tables you can take full advantage of
FDO-OGR auto-wrapping facility
This allows you to access any specific FDO-OGR Geometry as if it
where native SpatiaLite ones in a completely transparent way
==================================================================
Enter ".help" for instructions
spatialite> attach "nev.sqlite" AS nev;
spatialite>
spatialite> CREATE TABLE nev.countries AS SELECT * from fdo_ne_10m_admin_0_countries;
spatialite> CREATE TABLE nev.populated_places AS SELECT * FROM fdo_ne_10m_populated_places;
spatialite> CREATE TABLE nev.railroads AS SELECT * FROM fdo_ne_10m_railroads;
spatialite> .q
*** FDO-OGR auto-wrapping shutdown done ***
Wszystkie wiersze „utworzono VirtualFDO ...” wskazują, że Spatialite rozpoznał dane w formacie FDO i utworzył wirtualne tabele dla każdej z GEOMETRIĄ przekonwertowaną na format przestrzenny. Ja attach
do mojego nowego „Nev” bazy danych i tworzyć nowe tabele dla każdej warstwy Jestem zainteresowany ze CREATE TABLE ... AS SELECT * FROM ...
stwierdzeniami.
Teraz przełączam się z powrotem na nową bazę danych przestrzennych. I uruchom RecoverGeometryColumn()
na każdej tabeli, aby uzyskać odpowiednią bazę danych przestrzennych ze wszystkimi metadanymi itp. Zauważ, że format FDO pozwala na mieszanie typów geometrii WIELU i POJEDYNCZYCH, więc najpierw sprawdzam, które typy geometrii zawiera każda tabela, i upewniam się, że wszystkie funkcje są to samo. Używam CastToMulti()
wszędzie tam, gdzie to konieczne, jak:
micha@Wheezy:~/GIS/World/naturalearthdata.com$ spatialite nev.sqlite
SpatiaLite version ..: 4.1.1 Supported Extensions:
- 'VirtualShape' [direct Shapefile access]
- 'VirtualDbf' [direct DBF access]
- 'VirtualXL' [direct XLS access]
- 'VirtualText' [direct CSV/TXT access]
- 'VirtualNetwork' [Dijkstra shortest path]
- 'RTree' [Spatial Index - R*Tree]
- 'MbrCache' [Spatial Index - MBR cache]
- 'VirtualSpatialIndex' [R*Tree metahandler]
- 'VirtualFDO' [FDO-OGR interoperability]
- 'SpatiaLite' [Spatial SQL - OGC]
PROJ.4 version ......: Rel. 4.7.1, 23 September 2009
GEOS version ........: 3.3.3-CAPI-1.7.4
SQLite version ......: 3.7.13
Enter ".help" for instructions
SQLite version 3.7.13 2012-06-11 02:05:22
Enter ".help" for instructions
Enter SQL statements terminated with a ";"
spatialite> .tables
SpatialIndex geometry_columns_auth spatialite_history
countries populated_places sql_statements_log
geom_cols_ref_sys railroads views_geometry_columns
geometry_columns spatial_ref_sys virts_geometry_columns
spatialite>
spatialite> SELECT GeometryType(GEOMETRY) FROM countries;
POLYGON
POLYGON
MULTIPOLYGON
MULTIPOLYGON
POLYGON
MULTIPOLYGON
POLYGON
MULTIPOLYGON
MULTIPOLYGON
.....
Geometria jest mieszana, więc ustaw wszystko MULTI, a następnie wykonaj RecoverGeometryColumn ():
spatialite> UPDATE countries SET GEOMETRY=CastToMulti(GEOMETRY);
spatialite> SELECT RecoverGeometryColumn('countries','GEOMETRY',4326,'MULTIPOLYGON',2);
1
spatialite>
I tak dalej dla każdego stołu, którego potrzebujesz. Teraz tabele są dostępne w QGIS.
Możesz dodać dane z bazy danych za pomocą
Add vector layer ...
w QGIS 2.0.1.Ale bądź cierpliwy, to dużo danych.
Wtyczka Qspatialite niestety nie jest w stanie poradzić sobie z danymi ani dialogiem Add Spatialite layer.
źródło