搜尋
首頁資料庫mysql教程MongoDB中空间数据的存储和操作
MongoDB中空间数据的存储和操作Jun 07, 2016 pm 03:56 PM
dmongodb使用儲存官方操作數據空間

本文使用官方C# Driver,实现在MongoDB中存储,查询空间数据(矢量) 空间数据的存储 本例中,从一个矢量文件(shapefile格式)中读取矢量要素空间信息以及属性表,并写入到MongoDB中去,其中读取shapefile文件以及将空间信息转成json的功能通过Ogr库实现 [csh

本文使用官方C# Driver,实现在MongoDB中存储,查询空间数据(矢量)

空间数据的存储

本例中,从一个矢量文件(shapefile格式)中读取矢量要素空间信息以及属性表,并写入到MongoDB中去,其中读取shapefile文件以及将空间信息转成json的功能通过Ogr库实现
          
[csharp] view plaincopyprint?01.//打开MongoDB的Collection  
02.           MongoDatabase db = server.GetDatabase("aa"); 
03.           MongoCollection colSheng = db.GetCollection("sheng"); 
04.           //使用Ogr库打开Shapefile文件  
05.           DataSource ds = Ogr.Open(@"c:\temp\sheng.shp", 0); 
06.           Layer lyr = ds.GetLayerByIndex(0); 
07.           //读取要素数量和字段数量  
08.           int feaCount = lyr.GetFeatureCount(0); 
09.           int fieldCount = lyr.GetLayerDefn().GetFieldCount(); 
10.           //读取所有字段名  
11.           List fieldNames  =new List(); 
12.           for (int i = 0; i 13.           { 
14.               fieldNames.Add(lyr.GetLayerDefn().GetFieldDefn(i).GetName()); 
15.           } 
16.           //循环将所有要素添加到MongoDB中  
17.           for (int i = 0; i 18.           { 
19.               //使用Ogr库将矢量要素的空间信息转成Json格式  
20.               Feature fea = lyr.GetFeature(i); 
21.               Geometry geo = fea.GetGeometryRef();       
22.               string json = geo.ExportToJson(null); 
23.                
24.               BsonDocument doc = new BsonDocument(); 
25.                
26.               //将Json格式的空间信息存到Collection中  
27.               //BsonValue bs = BsonValue.Create(json);                  //这种方法是不可以的,添加到库里之后无法使用空间查询语句查询  
28.               BsonValue bs2 = BsonDocument.Parse(json);               //这种方法才是正确的  
29.               //doc.Add(new BsonElement("geom", bs));  
30.               doc.Add(new BsonElement("geo",bs2)); 
31.               //通过循环将所有字段的属性信息存入Collection中  
32.               for (int j = 0; j 33.               { 
34.                   string tmpFieldVal = fea.GetFieldAsString(j); 
35.                   doc.Add(new BsonElement(fieldNames[j],tmpFieldVal)); 
36.               } 
37.               var res  = colSheng.Insert(doc);                 
38.           } 
 //打开MongoDB的Collection
            MongoDatabase db = server.GetDatabase("aa");
            MongoCollection colSheng = db.GetCollection("sheng");
            //使用Ogr库打开Shapefile文件
            DataSource ds = Ogr.Open(@"c:\temp\sheng.shp", 0);
            Layer lyr = ds.GetLayerByIndex(0);
            //读取要素数量和字段数量
            int feaCount = lyr.GetFeatureCount(0);
            int fieldCount = lyr.GetLayerDefn().GetFieldCount();
            //读取所有字段名
            List fieldNames  =new List();
            for (int i = 0; i             {
                fieldNames.Add(lyr.GetLayerDefn().GetFieldDefn(i).GetName());
            }
            //循环将所有要素添加到MongoDB中
            for (int i = 0; i             {
                //使用Ogr库将矢量要素的空间信息转成Json格式
                Feature fea = lyr.GetFeature(i);
                Geometry geo = fea.GetGeometryRef();     
                string json = geo.ExportToJson(null);
               
                BsonDocument doc = new BsonDocument();
               
                //将Json格式的空间信息存到Collection中
                //BsonValue bs = BsonValue.Create(json);                  //这种方法是不可以的,添加到库里之后无法使用空间查询语句查询
                BsonValue bs2 = BsonDocument.Parse(json);               //这种方法才是正确的
                //doc.Add(new BsonElement("geom", bs));
                doc.Add(new BsonElement("geo",bs2));
                //通过循环将所有字段的属性信息存入Collection中
                for (int j = 0; j                 {
                    string tmpFieldVal = fea.GetFieldAsString(j);
                    doc.Add(new BsonElement(fieldNames[j],tmpFieldVal));
                }
                var res  = colSheng.Insert(doc);               
            }

然后,可以查看一下存储到MongoDB中的矢量数据是什么样的

在命令行中输入:
[csharp] view plaincopyprint?01.> db.sheng.find().limit(1) 
> db.sheng.find().limit(1)

结果为
[javascript] view plaincopyprint?01.{ "_id" : ObjectId("5371bf4e1dbba31914224563"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 89.8496, 14.093 ], [ 90.3933, 14.004 ], [ 90.2708, 13.4708 ], [ 89.7284, 13.5597 ], [ 89.8496, 14.093 ] ] ] }, "pyname" : "sx", "boxtype" : "inter", "date" : "2012/6/5 12:41:42" } 
{ "_id" : ObjectId("5371bf4e1dbba31914224563"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 89.8496, 14.093 ], [ 90.3933, 14.004 ], [ 90.2708, 13.4708 ], [ 89.7284, 13.5597 ], [ 89.8496, 14.093 ] ] ] }, "pyname" : "sx", "boxtype" : "inter", "date" : "2012/6/5 12:41:42" }

可以看到名称为geo的这个Field,里边存的就是矢量要素的坐标信息

空间查询与空间索引

可用的空间操作包括geointersect,geowithin,near等,参考http://docs.mongodb.org/manual/reference/operator/query-geospatial/
这里使用geointersect为例说明一下:
          
[csharp] view plaincopyprint?01.//获取Collection  
02.           MongoDatabase db = server.GetDatabase("aa"); 
03.           MongoCollection colSheng = db.GetCollection("sheng"); 
04.         
05.           //定义一个查询框或查询多边形  
06.           var poly = GeoJson.Polygon
07.               GeoJson.Position(100, 20), 
08.               GeoJson.Position(110, 20), 
09.               GeoJson.Position(110, 40), 
10.               GeoJson.Position(100, 40), 
11.               GeoJson.Position(100, 20)); 
12.           //以这个查询多边形为条件定义一条查询语句  
13.           var queryFilter2 = Query.GeoIntersects("geo", poly); 
14.           //进行查询,输出MongoCursor  
15.           cur = colSheng.FindAs(queryFilter2).SetFields( "pyname", "date"); 
16.           //获取结果  
17.           var res = cur.ToArray(); 
18.           for (int i = 0; i 19.           { 
20.               BsonDocument tmpDoc = res.ElementAt(i); 
21.               //do something you want  
22.           } 
 //获取Collection
            MongoDatabase db = server.GetDatabase("aa");
            MongoCollection colSheng = db.GetCollection("sheng");
        
            //定义一个查询框或查询多边形
            var poly = GeoJson.Polygon(
                GeoJson.Position(100, 20),
                GeoJson.Position(110, 20),
                GeoJson.Position(110, 40),
                GeoJson.Position(100, 40),
                GeoJson.Position(100, 20));
            //以这个查询多边形为条件定义一条查询语句
            var queryFilter2 = Query.GeoIntersects("geo", poly);
            //进行查询,输出MongoCursor
            cur = colSheng.FindAs(queryFilter2).SetFields( "pyname", "date");
            //获取结果
            var res = cur.ToArray();
            for (int i = 0; i             {
                BsonDocument tmpDoc = res.ElementAt(i);
                //do something you want
            }

关于空间索引,可参考http://docs.mongodb.org/manual/applications/geospatial-indexes/
这里不详细说了

空间查询运算的问题:

在使用GeoIntersect进行空间查询时,遇到了查询结果与ArcGIS不一致的情况,详细看了一下,像是MongoDB的一个BUG(目前使用的是2.6.0版本)

具体信息如下(在命令行中操作):

Collection中的坐标

[csharp] view plaincopyprint?01.> db.test.find() 
02.{ "_id" : ObjectId("535884771dbba31858ad2101"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 96.722, 38.755 ], [ 97.3482, 38.6922 ], [ 97.1674, 38.0752 ], [ 96.5474, 38.1383 ], [ 96.722, 38.755 ] ] ] } } 
> db.test.find()
{ "_id" : ObjectId("535884771dbba31858ad2101"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 96.722, 38.755 ], [ 97.3482, 38.6922 ], [ 97.1674, 38.0752 ], [ 96.5474, 38.1383 ], [ 96.722, 38.755 ] ] ] } }

使用的查询语句

[csharp] view plaincopyprint?01.> db.test.find({ "geo" : { "$geoIntersects" : { "$geometry" : { "type" : "Polygon", "coordinates" : [[[91.0, 33.0], [102.0, 33.0], [102.0, 38.0], [91.0, 38.0], [91.0, 33.0]]] } } } }) 
> db.test.find({ "geo" : { "$geoIntersects" : { "$geometry" : { "type" : "Polygon", "coordinates" : [[[91.0, 33.0], [102.0, 33.0], [102.0, 38.0], [91.0, 38.0], [91.0, 33.0]]] } } } })

查询结果:

[csharp] view plaincopyprint?01.{ "_id" : ObjectId("535884771dbba31858ad2101"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 96.722, 38.755 ], [ 97.3482, 38.6922 ], [ 97.1674, 38.0752 ], [ 96.5474, 38.1383 ], [ 96.722, 38.755 ] ] ] } } 
{ "_id" : ObjectId("535884771dbba31858ad2101"), "geo" : { "type" : "Polygon", "coordinates" : [ [ [ 96.722, 38.755 ], [ 97.3482, 38.6922 ], [ 97.1674, 38.0752 ], [ 96.5474, 38.1383 ], [ 96.722, 38.755 ] ] ] } }

但可以看到,collection中只有一条记录,且该记录所有点的Y坐标均大于38.0,为什么查询结果里,这条记录与语句中的Box相交呢。。。很奇怪

因为有这样的问题,所以还不放心直接将空间查询用于实际应用,而是通过一种变通的方法进行简单的空间查询,测试后发现,可能是由于空间索引的问题,这种方式查询比自带的GeoIntersects方法要快

大致思路为:为每一条记录均生成一个最小外接矩形,得到其xmax,xmin,ymax,ymin四个边界值,用数值的形式保存至Collection中,每次进行空间查询时,首先通过最小外接矩形进行一次筛选,判断这些最小外接矩形与查询语句中多边形的最小外接矩形之间的关系,如果相交,那么进行第二步判断,通过Ogr组件判断实际的多边形是否相交,返回最后结果

首先是生成最小外接矩形的代码:
          
[csharp] view plaincopyprint?01.//获取Collection  
02.            MongoDatabase db = server.GetDatabase("aa"); 
03.            MongoCollection colsheng= db.GetCollection("sheng"); 
04.            //查询所有记录  
05.            var cur = colsheng.FindAllAs(); 
06.            long totalCount = cur.Count(); 
07.            //遍历所有记录  
08.            for (int i = 0; i 09.            { 
10.                if (i * 1000 >= totalCount) continue; 
11.                int skip = i * 1000; 
12.                var cur2 = cur.Clone().SetSkip(skip).SetLimit(1000); 
13.                var lst = cur2.ToArray(); 
14.                for (int j = 0; j 15.                { 
16.                    //获取一条记录对应的BsonDocument  
17.                    BsonDocument doc = lst[j]; 
18.                        var id = doc["_id"];                    //该记录对应的ID  
19.                        BsonDocument geo = doc["geo"].ToBsonDocument();     
20.                        string geostr = geo[1].ToString();        //该记录对应空间信息的Json字符串  
21.                        List coords = GetCoordLstFromString(geostr);        //解析Json串,获得所有点的坐标(这里子函数就省略了)  
22.                        double xmin = 181, xmax = -181, ymin = 91, ymax = -91;        //四个边界值,由于图层为经纬度,所以初值设为这些值  
23.                        //计算最大最小值  
24.                        for (int k = 0; k 25.                        { 
26.                            if (k % 2 == 0) 
27.                            { 
28.                                if (coords[k] 29.                                if (coords[k] > xmax) xmax = coords[k]; 
30.                            } 
31.                            else 
32.                            { 
33.                                if (coords[k] 34.                                if (coords[k] > ymax) ymax = coords[k]; 
35.                            } 
36.                        } 
37.                        //将最大最小值写入Collection  
38.                        var tmpQuery = Query.EQ("_id", id); 
39.                        var tmpUpdate = MongoDB.Driver.Builders.Update.Set("xmax", xmax); 
40.                        var tmpres = col02c.Update(tmpQuery, tmpUpdate); 
41.                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("xmin", xmin); 
42.                        tmpres = col02c.Update(tmpQuery, tmpUpdate); 
43.                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("ymax", ymax); 
44.                        tmpres = col02c.Update(tmpQuery, tmpUpdate); 
45.                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("ymin", ymin); 
46.                        tmpres = col02c.Update(tmpQuery, tmpUpdate); 
47.                    } 
48.                } 
//获取Collection
            MongoDatabase db = server.GetDatabase("aa");
            MongoCollection colsheng= db.GetCollection("sheng");
            //查询所有记录
            var cur = colsheng.FindAllAs();
            long totalCount = cur.Count();
            //遍历所有记录
            for (int i = 0; i             {
                if (i * 1000 >= totalCount) continue;
                int skip = i * 1000;
                var cur2 = cur.Clone().SetSkip(skip).SetLimit(1000);
                var lst = cur2.ToArray();
                for (int j = 0; j                 {
                    //获取一条记录对应的BsonDocument
                    BsonDocument doc = lst[j];
                        var id = doc["_id"];                    //该记录对应的ID
                        BsonDocument geo = doc["geo"].ToBsonDocument();   
                        string geostr = geo[1].ToString();        //该记录对应空间信息的Json字符串
                        List coords = GetCoordLstFromString(geostr);        //解析Json串,获得所有点的坐标(这里子函数就省略了)
                        double xmin = 181, xmax = -181, ymin = 91, ymax = -91;        //四个边界值,由于图层为经纬度,所以初值设为这些值
                        //计算最大最小值
                        for (int k = 0; k                         {
                            if (k % 2 == 0)
                            {
                                if (coords[k]                                 if (coords[k] > xmax) xmax = coords[k];
                            }
                            else
                            {
                                if (coords[k]                                 if (coords[k] > ymax) ymax = coords[k];
                            }
                        }
                        //将最大最小值写入Collection
                        var tmpQuery = Query.EQ("_id", id);
                        var tmpUpdate = MongoDB.Driver.Builders.Update.Set("xmax", xmax);
                        var tmpres = col02c.Update(tmpQuery, tmpUpdate);
                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("xmin", xmin);
                        tmpres = col02c.Update(tmpQuery, tmpUpdate);
                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("ymax", ymax);
                        tmpres = col02c.Update(tmpQuery, tmpUpdate);
                        tmpUpdate = MongoDB.Driver.Builders.Update.Set("ymin", ymin);
                        tmpres = col02c.Update(tmpQuery, tmpUpdate);
                    }
                }

然后是查询的代码:
          
[csharp] view plaincopyprint?01.//获取Collection  
02.            MongoDatabase db = server.GetDatabase("aa"); 
03.            MongoCollection colSheng = db.GetCollection("zy02c"); 
04.            //第一步,通过四边界筛选,  
05.            var query = Query.And(Query.GT("xmax", 91.0), Query.LT("xmin", 102.0), Query.GT("ymax", 33.0), Query.LT("ymin", 38.0)); 
06.            var cur = colSheng.FindAs(query); 
07.            //定义第二空间运算时的条件多边形(Ogr格式的定义)  
08.            Geometry queryGeoLR = new Geometry(wkbGeometryType.wkbLinearRing); 
09.            queryGeoLR.AddPoint(91.0, 33.0,0); 
10.            queryGeoLR.AddPoint(102.0, 33.0,0); 
11.            queryGeoLR.AddPoint(102.0, 38.0,0); 
12.            queryGeoLR.AddPoint(91.0, 38.0,0); 
13.            queryGeoLR.AddPoint(91.0, 33.0,0); 
14.            Geometry queryGeo = new Geometry(wkbGeometryType.wkbPolygon); 
15.            queryGeo.AddGeometry(queryGeoLR); 
16.            //循环查询到的结果  
17.            var lst = cur.ToArray(); 
18.            for (int i = lst.Length-1; i >=0; i--) 
19.            { 
20.                //获取当前记录对应的BsonDocument  
21.                BsonDocument doc = lst[i]; 
22.                var id = doc["_id"];            //当前记录的ID  
23.                BsonDocument geo = doc["geo"].ToBsonDocument(); 
24.                string geostr = geo[1].ToString();        //当前记录对应空间信息的Json字符串  
25.                 
26.                //通过Json串获取坐标值,并生成对应的Geometry对象  
27.                List coords = GetCoordLstFromString(geostr); 
28.                Geometry resGeoLR = new Geometry(wkbGeometryType.wkbLinearRing); 
29.                for (int j = 0; j 30.                { 
31.                    resGeoLR.AddPoint_2D(coords[j], coords[j + 1]); 
32.                } 
33.                resGeoLR.AddPoint_2D(coords[0], coords[1]); 
34.                Geometry resGeo = new Geometry(wkbGeometryType.wkbPolygon); 
35.                resGeo.AddGeometry(resGeoLR); 
36.                //判断是该Geometry与条件多边形是否相交  
37.                if (resGeo.Intersects(queryGeo)) 
38.                { 
39.                        //do something  
40.                } 
41.            } 

陳述
本文內容由網友自願投稿,版權歸原作者所有。本站不承擔相應的法律責任。如發現涉嫌抄襲或侵權的內容,請聯絡admin@php.cn
mongodb php 扩展没有怎么办mongodb php 扩展没有怎么办Nov 06, 2022 am 09:10 AM

mongodb php扩展没有的解决办法:1、在linux中执行“$ sudo pecl install mongo”命令来安装MongoDB的PHP扩展驱动;2、在window中,下载php mongodb驱动二进制包,然后在“php.ini”文件中配置“extension=php_mongo.dll”即可。

Redis和MongoDB的区别与使用场景Redis和MongoDB的区别与使用场景May 11, 2023 am 08:22 AM

Redis和MongoDB都是流行的开源NoSQL数据库,但它们的设计理念和使用场景有所不同。本文将重点介绍Redis和MongoDB的区别和使用场景。Redis和MongoDB简介Redis是一个高性能的数据存储系统,常被用作缓存和消息中间件。Redis以内存为主要存储介质,但它也支持将数据持久化到磁盘上。Redis是一款键值数据库,它支持多种数据结构(例

Go语言中使用MongoDB:完整指南Go语言中使用MongoDB:完整指南Jun 17, 2023 pm 06:14 PM

MongoDB是一种高性能、开源、文档型的NoSQL数据库,被广泛应用于Web应用、大数据以及云计算领域。而Go语言则是一种快速、开发效率高、代码可维护性强的编程语言。本文将为您完整介绍如何在Go语言中使用MongoDB。一、安装MongoDB在使用MongoDB之前,需要先在您的系统中安装MongoDB。在Linux系统下,可以通过如下命令安装:sudo

php7.0怎么安装mongo扩展php7.0怎么安装mongo扩展Nov 21, 2022 am 10:25 AM

php7.0安装mongo扩展的方法:1、创建mongodb用户组和用户;2、下载mongodb源码包,并将源码包放到“/usr/local/src/”目录下;3、进入“src/”目录;4、解压源码包;5、创建mongodb文件目录;6、将文件复制到“mongodb/”目录;7、创建mongodb配置文件并修改配置即可。

php怎么使用mongodb进行增删查改操作php怎么使用mongodb进行增删查改操作Mar 28, 2023 pm 03:00 PM

MongoDB作为一款流行的NoSQL数据库,已经被广泛应用于各种大型Web应用和企业级应用中。而PHP语言也作为一种流行的Web编程语言,与MongoDB的结合也变得越来越重要。在本文中,我们将会学习如何使用PHP语言操作MongoDB数据库进行增删查改的操作。

SpringBoot中logback日志怎么保存到mongoDBSpringBoot中logback日志怎么保存到mongoDBMay 18, 2023 pm 07:01 PM

自定义Appender非常简单,继承一下AppenderBase类即可。可以看到有个AppenderBase,有个UnsynchronizedAppenderBase,还有个AsyncAppenderBase继承了UnsynchronizedAppenderBase。从名字就能看出来区别,异步的、普通的、不加锁的。我们定义一个MongoDBAppender继承UnsynchronizedAppenderBasepublicclassMongoDBAppenderextendsUnsynchron

SpringBoot怎么整合Mongodb实现增删查改SpringBoot怎么整合Mongodb实现增删查改May 13, 2023 pm 02:07 PM

一、什么是MongoDBMongoDB与我们之前熟知的关系型数据库(MySQL、Oracle)不同,MongoDB是一个文档数据库,它具有所需的可伸缩性和灵活性,以及所需的查询和索引。MongoDB将数据存储在灵活的、类似JSON的文档中,这意味着文档的字段可能因文档而异,数据结构也会随着时间的推移而改变。文档模型映射到应用程序代码中的对象,使数据易于处理。MongoDB是一个以分布式数据库为核心的数据库,因此高可用性、横向扩展和地理分布是内置的,并且易于使用。况且,MongoDB是免费的,开源

Swoole与MongoDB的整合:构建高性能的文档数据库系统Swoole与MongoDB的整合:构建高性能的文档数据库系统Jun 14, 2023 am 11:51 AM

在现代企业应用程序开发中,需要处理海量数据和高并发的访问请求。为了满足这些需求,开发人员需要使用高性能的数据库系统,以确保系统的稳定性和可扩展性。本文将介绍如何使用Swoole和MongoDB构建高性能的文档数据库系统。Swoole是一个基于PHP语言开发的异步网络通信框架,它能够大大提高PHP应用程序的性能和并发能力。MongoDB是一种流行的文档数据库,

See all articles

熱AI工具

Undresser.AI Undress

Undresser.AI Undress

人工智慧驅動的應用程序,用於創建逼真的裸體照片

AI Clothes Remover

AI Clothes Remover

用於從照片中去除衣服的線上人工智慧工具。

Undress AI Tool

Undress AI Tool

免費脫衣圖片

Clothoff.io

Clothoff.io

AI脫衣器

AI Hentai Generator

AI Hentai Generator

免費產生 AI 無盡。

熱門文章

R.E.P.O.能量晶體解釋及其做什麼(黃色晶體)
2 週前By尊渡假赌尊渡假赌尊渡假赌
倉庫:如何復興隊友
1 個月前By尊渡假赌尊渡假赌尊渡假赌
Hello Kitty Island冒險:如何獲得巨型種子
4 週前By尊渡假赌尊渡假赌尊渡假赌

熱工具

MantisBT

MantisBT

Mantis是一個易於部署的基於Web的缺陷追蹤工具,用於幫助產品缺陷追蹤。它需要PHP、MySQL和一個Web伺服器。請查看我們的演示和託管服務。

VSCode Windows 64位元 下載

VSCode Windows 64位元 下載

微軟推出的免費、功能強大的一款IDE編輯器

Dreamweaver Mac版

Dreamweaver Mac版

視覺化網頁開發工具

SublimeText3 英文版

SublimeText3 英文版

推薦:為Win版本,支援程式碼提示!

記事本++7.3.1

記事本++7.3.1

好用且免費的程式碼編輯器