相关文章推荐
多情的剪刀  ·  政务公开_上海杨浦_关于印发《杨浦区工业通信 ...·  9 月前    · 
淡定的乒乓球  ·  CBETA 線上閱讀·  1 年前    · 
爱喝酒的楼房  ·  不二倾谈——没有陌生人的世界 | ...·  1 年前    · 
打盹的水煮肉  ·  【宝马2系多功能旅行车】宝马(进口)_宝马2 ...·  2 年前    · 
爽快的生姜  ·  估值500亿！中国数据库第一股准备上市！！- ...·  2 年前    · 
Code  ›  geopandas 笔记：geometry上的操作汇总
osg
https://blog.csdn.net/qq_40206371/article/details/132767108
跑龙套的开水瓶
2 年前

如无特殊说明，数据主要来自： GeoDataFrame 应用：公园分布映射至subzone_UQI-LIUWJ的博客-CSDN博客

0 读入数据 

subzone = gpd.read_file('ura-mp19-subzone-no-sea-pl.geojson')
subzone
subzone_tst=subzone[0:5]
subzone_tst
subzone_tst.plot()

1 area 计算面积 

subzone_tst.area
0    0.000036
1    0.000033
2    0.000048
3    0.000003
4    0.000015
dtype: float64

2 boundary 边界
 

subzone_tst.boundary
0    LINESTRING Z (103.81454 1.28239 0.00000, 103.8...
1    LINESTRING Z (103.82209 1.28049 0.00000, 103.8...
2    LINESTRING Z (103.84375 1.28508 0.00000, 103.8...
3    LINESTRING Z (103.84962 1.28412 0.00000, 103.8...
4    LINESTRING Z (103.85253 1.28617 0.00000, 103.8...
dtype: geometry
subzone_tst.boundary.plot()
 

3 centroid 中心点
 

subzone_tst.centroid
0    POINT (103.80856 1.28222)
1    POINT (103.81859 1.28201)
2    POINT (103.84369 1.27997)
3    POINT (103.84865 1.28528)
4    POINT (103.85101 1.28372)
dtype: geometry

ax=subzone_tst.plot()
subzone_tst.centroid.plot(ax=ax,color='red')
 

4 convex_hull 凸包
 

subzone_tst.convex_hull
0    POLYGON Z ((103.81749 1.28004 0.00000, 103.813...
1    POLYGON Z ((103.82153 1.27882 0.00000, 103.821...
2    POLYGON Z ((103.84137 1.27415 0.00000, 103.841...
3    POLYGON Z ((103.84955 1.28391 0.00000, 103.847...
4    POLYGON Z ((103.85108 1.28077 0.00000, 103.849...
dtype: geometry
ax=subzone_tst.convex_hull.plot()
subzone_tst.plot(ax=ax,color='red',alpha=0.2)
 

5 envelope 最小旋转矩形
 

subzone_tst.envelope
0    POLYGON ((103.80126 1.28004, 103.81774 1.28004...
1    POLYGON ((103.81306 1.27882, 103.82372 1.27882...
2    POLYGON ((103.83960 1.27415, 103.84919 1.27415...
3    POLYGON ((103.84718 1.28391, 103.84995 1.28391...
4    POLYGON ((103.84852 1.28077, 103.85321 1.28077...
dtype: geometry
ax=subzone_tst.envelope.plot()
subzone_tst.plot(ax=ax,color='red',alpha=0.2)
 

6 buffer 缓冲区
 

ax=subzone_tst.buffer(0.01,0.05).plot()
subzone_tst.plot(ax=ax,color='red',alpha=0.2)
 

7 intersects 是否与另一个几何形状相交
 

subzone_tst.intersects(subzone.at[5,'geometry'])
0    False
1    False
2     True
3     True
4     True
dtype: bool

ax=subzone_tst.plot(figsize=(10,7))
subzone_tst.boundary.plot(ax=ax,color='black')
gpd.GeoSeries(subzone.loc[5,'geometry']).plot(ax=ax,color='red',alpha=0.4)
 

8 contains  是否包含另一个几何形状
 

subzone_tst.contains(subzone.at[5,'geometry'])
0    False
1    False
2    False
3    False
4    False
dtype: bool

9 crosses 是否与另一个几何形状交叉
 

subzone_tst.crosses(subzone.at[5,'geometry'])
0    False
1    False
2    False
3    False
4    False
dtype: bool

9.1 crosses和intersects的区别
 

  •  intersects 
  
    • 如果两个几何形状共享任何点、边或面，则它们相交
  • crosses 
  
    • crosses 描述了一个几何形状与另一个不同维度的几何形状之间的关系
    • 如果一个几何形状（线或多边形）穿越另一个几何形状，但不完全包含在其中，那么它们就交叉
 

10 distance 距离
 

subzone_tst.distance(subzone.at[5,'geometry'])
0    0.026645
1    0.020509
2    0.000000
3    0.000000
4    0.000000
dtype: float64

11 union 并集
 

subzone_tst.union(subzone)
0      POLYGON Z ((103.81774 1.28043 0.00000, 103.817...
1      POLYGON Z ((103.82210 1.28011 0.00000, 103.822...
2      POLYGON Z ((103.84400 1.28491 0.00000, 103.844...
3      POLYGON Z ((103.84955 1.28391 0.00000, 103.849...
4      POLYGON Z ((103.85253 1.28615 0.00000, 103.852...
327                                                 None
328                                                 None
329                                                 None
330                                                 None
331                                                 None
Length: 332, dtype: geometry

12 unary_union
 

合并所有几何形状为一个单一的几何形状
 

gpd.GeoSeries(subzone_tst.unary_union)
0    MULTIPOLYGON Z (((103.84418 1.28480 0.00000, 1...
dtype: geometry
Geometry Dash专用服务器
基本上是Geometry Dash服务器仿真器
受支持的Geometry Dash版本：1.0-2.11（因此，撰写此[2020年2月2日]时，任何版本的Geometry Dash均可使用）
 所需PHP版本：5.4+（经过7.3.11测试）
将文件上传到网络服务器上
将database.sql导入MySQL / MariaDB数据库
编辑GeometryDash.exe中的链接（自2.1开始，其中一些是base64编码的，请记住）
someguy28的帐户设置和私人消息传递系统的基础
使用这种用于XOR加密-  - （包括/ lib目录/ XORCipher.php）
 使用这种云保存加密-  - （包括/ lib目录/化解-crypto.phar）
 Jscolor（在packCreate.php颜色选
				
python-geopandas读取、创建shapefile文件、geopandas学习教程
shapefile是GIS中非常重要的一种数据类型，在ArcGIS中被称为要素类(Feature Class)，主要包括点(point)、线(polyline)和多边形(polygon)。作为一种十分常见的矢量文件格式，geopandas对shapefile提供了很好的读取和写出支持 。
geopandas库允许对几何类型进行空间操作，其DataFrame结构相当于GIS数据中的一张属性表，使得可以直接操作矢量数据

file_name = "****.shp"
data = gpd.read_file(file_name, encoding="utf-8")
data.crs # {'init': 'epsg:4326'}
reprojeted_data = data.to_crs(crs="+proj =utm +zone=48 +...
gdf = gs.read_file("F:\\Zhujie\\JY_Data\\test_blocks.shp")
polygon = gdf.geometry.to_json()
polygon_dict = json.loads(polygon)
print(polygon_dict["features"][0]["geometry"]["c
一、geopandas介绍
GeoPandas是一个开源项目，主要为了简化python中地理空间数据处理的过程。GeoPandas扩展了pandas使用的数据类型，允许对几何类型进行空间操作。
GeoPandas几何运算由shape执行。Geopandas进一步依赖于fiona的文件存取和matplotlib的绘图。
GeoPandas结合了pandas和shaely的功能，提供了pandas中的地理空间操作和多个几何
GeoPandas是一个开源项目，可以更轻松地使用python处理地理空间数据。
GeoPandas扩展了Pandas中使用的数据类型DataFrame，允许对几何类型进行空间操作。
GeoPandas的目标是使在python中使用地理空间数据更容易。它结合了Pandas和Shapely的能力，提供了Pandas的地理空间操作和多种Shapely的高级接口。GeoPandas可以让您轻松地在python中进行操作，否则将需要...
				
将一个 `osg::Geometry` 对象拆分成多个 `osg::Geometry` 对象，可以通过以下步骤实现：
1. 遍历原始 `osg::Geometry` 对象的 `osg::PrimitiveSet`，并将其中不同类型的 `PrimitiveSet` 分别存储到不同的 `osg::PrimitiveSet` 容器中（如 `osg::DrawArrays` 存储到 `osg::DrawArrays` 容器中，`osg::DrawElementsUInt` 存储到 `osg::DrawElementsUInt` 容器中等等）。
2. 遍历存储不同类型的 `PrimitiveSet` 的容器，根据每个容器中存储的 `PrimitiveSet` 创建对应的 `osg::Geometry` 对象，并将这些 `osg::Geometry` 对象添加到场景图中。
以下是一个示例代码，可以将一个 `osg::Geometry` 对象拆分成多个 `osg::Geometry` 对象：
```cpp
osg::ref_ptr<osg::Geometry> originalGeometry = ...;
// 定义不同类型的 PrimitiveSet 容器
std::vector<osg::ref_ptr<osg::DrawArrays>> drawArraysList;
std::vector<osg::ref_ptr<osg::DrawElementsUInt>> drawElementsList;
// 遍历原始 Geometry 对象的 PrimitiveSet，将不同类型的 PrimitiveSet 存储到不同容器中
for (unsigned int i = 0; i < originalGeometry->getNumPrimitiveSets(); ++i)
    osg::PrimitiveSet* primitiveSet = originalGeometry->getPrimitiveSet(i);
    if (primitiveSet->getType() == osg::PrimitiveSet::DrawArraysPrimitiveType)
        drawArraysList.push_back(static_cast<osg::DrawArrays*>(primitiveSet));
    else if (primitiveSet->getType() == osg::PrimitiveSet::DrawElementsUIntPrimitiveType)
        drawElementsList.push_back(static_cast<osg::DrawElementsUInt*>(primitiveSet));
// 创建多个 Geometry 对象，并将它们添加到场景图中
osg::ref_ptr<osg::Group> group = new osg::Group;
for (unsigned int i = 0; i < drawArraysList.size(); ++i)
    osg::ref_ptr<osg::Geometry> newGeometry = new osg::Geometry;
    newGeometry->addPrimitiveSet(drawArraysList[i].get());
    group->addChild(newGeometry);
for (unsigned int i = 0; i < drawElementsList.size(); ++i)
    osg::ref_ptr<osg::Geometry> newGeometry = new osg::Geometry;
    newGeometry->addPrimitiveSet(drawElementsList[i].get());
    group->addChild(newGeometry);
请注意，这只是一个简单的示例代码，实际应用中可能需要更多的逻辑处理。

                论文笔记：TILDE-Q: A TRANSFORMATION INVARIANT LOSS FUNCTION FOR TIME-SERIES FORECASTING
                    faith~~: 
                    能找到相关代码吗，实验过这个方法对时间序列有效果吗
                PYTORCH笔记 actor-critic （A2C）
                    2301_79774295: 
                    博主，计算td-error 减号两边反了吧，公式是v-tdtarget
 
推荐文章
多情的剪刀  ·  政务公开_上海杨浦_关于印发《杨浦区工业通信业服务业节能减排和合同能源管理专项扶持办法》的通知
9 月前
淡定的乒乓球  ·  CBETA 線上閱讀
1 年前
爱喝酒的楼房  ·  不二倾谈——没有陌生人的世界 | 香港中文大学（深圳）
1 年前
打盹的水煮肉  ·  【宝马2系多功能旅行车】宝马(进口)_宝马2系多功能旅行车_参数_图片_投诉_口碑_答疑_资讯-车质网
2 年前
爽快的生姜  ·  估值500亿！中国数据库第一股准备上市！！-腾讯云开发者社区-腾讯云
2 年前
今天看啥   ·   Py中国   ·   codingpro   ·   小百科   ·   link之家   ·   卧龙AI搜索
删除内容请联系邮箱 2879853325@qq.com
Code - 代码工具平台
© 2024 ~ 沪ICP备11025650号