मैं QGIS के साथ प्राकृतिक पृथ्वी SQLite DB का उपयोग कैसे करूँ?


9

मैंने बस http://www.naturalearthdata.com/downloads/ से SQLite प्रारूप में प्राकृतिक पृथ्वी डेटा डाउनलोड किया । मैंने यह मान लिया था कि यह एक SpatiaLite डेटाबेस होगा, लेकिन ऐसा नहीं प्रतीत होता है! QGIS इसे एक स्थानिक डेटाबेस के रूप में पहचानने में विफल हो रहा है। सादे SQLite डेटाबेस में संग्रहीत होने पर भी OGR कथित तौर पर ज्यामिति पढ़ने का समर्थन करता है, लेकिन शायद QGIS OGR का उपयोग SpatiaLite के लिए नहीं करता है?

SQLite डेटाबेस में ज्योमेट्री_ कॉलम और स्थानिक_संग्रह तालिकाएँ हैं। क्या इसे किसी पूर्ण SpatiaLite डेटाबेस में बदलने का कोई तरीका है?

जवाबों:


9

NE से sqlite फ़ाइल FDO-OGR प्रारूप में है, न कि देशी स्थानिक ज्यामिति। यदि आप कुछ मैन्युअल श्रम करने के लिए तैयार हैं, तो यहाँ एक स्थानिक स्थान पर बदलने का तरीका है:

पहले एक नया, खाली स्थानिक डेटाबेस बनाएं (मैं इसे "nev.sqlite" कहता हूं), फिर एक अलग टर्मिनल सत्र में मूल प्राकृतिक_आर्थ_वेक्टर.सक्लाइट को स्थानिक के साथ खोलें। (मैंने नए संस्करण 4.1 का उपयोग किया है। यकीन नहीं है कि यह पुराने संस्करणों के साथ काम करेगा)। attachअपनी नई nev.sqlite तालिका से कनेक्ट करने के लिए साइक्लाइट फ़ंक्शन का उपयोग करें , और नए डेटाबेस में इच्छित तालिकाओं की प्रतियां बनाएं।

इसलिए:

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 ***

सभी लाइनों ने "वर्चुअलाइफ़डीओ बनाया ..." से संकेत मिलता है कि स्पैटियालाइट ने डेटा को एफडीओ के रूप में मान्यता दी, और प्रत्येक के लिए वर्चुअल टेबल बनाई जिसमें जीओएमईटीवाई स्पैटियालाइट प्रारूप में परिवर्तित हो गया। मैं attachअपने नए "नेव" डेटाबेस में आता हूं और प्रत्येक परत के लिए नए टेबल बनाता हूं जिन्हें मैं CREATE TABLE ... AS SELECT * FROM ...बयानों में दिलचस्पी लेता हूं ।

अब मैं नए स्थानिक डेटाबेस पर वापस स्विच करता हूं । और RecoverGeometryColumn()सभी मेटाडेटा, आदि के साथ एक उचित स्थानिक डेटाबेस प्राप्त करने के लिए प्रत्येक तालिका पर चलाएं । ध्यान दें कि FDO प्रारूप मिश्रित MULTI और SINGLE ज्यामिति प्रकारों के लिए अनुमति देता है, इसलिए मैं पहले जाँचता हूं कि प्रत्येक तालिका में किस ज्यामिति के प्रकार हैं, और सुनिश्चित करें कि सभी विशेषताएं हैं वही। मैं CastToMulti()जहाँ भी आवश्यक हो, उपयोग करता हूँ :

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
.....

ज्यामिति मिश्रित होती है, इसलिए सब कुछ MULTI सेट करें, फिर RecoverGeometryColumn () करें:

spatialite> UPDATE countries SET GEOMETRY=CastToMulti(GEOMETRY);
spatialite> SELECT RecoverGeometryColumn('countries','GEOMETRY',4326,'MULTIPOLYGON',2);
1
spatialite> 

और इसलिए आपको प्रत्येक तालिका की आवश्यकता है। अब टेबल QGIS में उपलब्ध हैं।


पूरी तरह से जवाब देने के लिए धन्यवाद। क्या आप जानते हैं कि एफडीओ के साथ केवल कुछ तालिकाओं को "आयात" क्यों किया जाता है? मैं 128 में से 30 तालिकाओं की गिनती करता हूं जिन्हें एफडीओ उपचार मिलता है। मूल और fdo संस्करणों पर PRAGMA समान है, इसलिए मुझे लगता है कि केवल ज्यामिति में ही अंतर है। लेकिन ज्योमेट्री_कैल्यूम सोचता है कि सभी ज्योमेट्रीज WKB हैं।
ली हैचेरडियन

मीका, मैंने इसे सर्वश्रेष्ठ उत्तर के रूप में चुना है। तालिकाओं पर आगे देखते हुए, मुझे यह अजीब लगता है कि डेटाबेस को VirtualFDO के माध्यम से उपयोग किए जाने वाले तालिकाओं में से केवल कुछ के साथ पैक किया गया है, और अन्य को ज्यामितीय लेकिन स्पैटियालाइट के साथ संगत नहीं है। इस तरह से बदलने की कोशिश करने के बजाय, मुझे लगता है कि मैं शेपफाइल्स डाउनलोड करूंगा और स्पैटियालाइट (जो मैंने पहले ही पोस्टजीआईएस के लिए किया है) को आयात करूंगा। लेकिन यह अभी भी बहुत जानकारीपूर्ण उत्तर था, आपने इसे बनाने के लिए जो समय लिया, उसके लिए धन्यवाद।
बजे ली हैचेरियन

यदि आप केवल एक स्थानिक डेटाबेस चाहते हैं, और आपके पास GDAL / OGR है जिसमें संकलित समर्थन के साथ संकलित किया गया है: ogr2ogr -f sqlite -dsco spatialite = yes splite.db nat_thth.db काम करना चाहिए।

3

आप डेटाबेस से डेटा Add vector layer ...QGIS 2.0.1 में जोड़ सकते हैं ।

लेकिन धैर्य रखें, यह बहुत सारा डेटा है।

Qspatialite प्लगइन दुर्भाग्य से डेटा से निपटने में सक्षम नहीं है, और न ही Spatialite परत संवाद जोड़ें।


आंद्रे, यह वास्तव में उपयोगी जानकारी है। मुझे एहसास नहीं था कि क्यूजीआईएस SQLite-not-SpatiaLite डेटाबेस से डेटा जोड़ने के लिए OGR का उपयोग कर सकता है।
ली हैचोरियन

मुझे यह भी समझ में नहीं आता है कि परतों को जोड़ने के दो तरीके क्यों हैं, और केवल एक ही काम करता है ;-)
आंद्रे जे
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.