首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往

(数据科学学习手札77)基于geopandas空间数据分析——文件IO

layer:str类型,当要读入数据格式为地理数据库.gdb或QGIS.gpkg时,传入对应图层名称   下面结合上述参数,来介绍一下使用geopandas.read_file()不同情况下读取常见格式矢量数据方法...,使用示例数据为中国地图,CRS为EPSG:4326,本文使用所有数据都可以文章开头提及Github仓库对应本文路径下找到: ?...图3 缺少投影shapefile   当shapefile缺失.prj文件时,使用geopandas读入后形成GeoDataFrame会缺失crs属性: ?...图12 2.1.3 GeoJSON   作为web地图中最常使用矢量数据格式,GeoJSON几乎被所有在线地图框架作为数据源格式,geopandas读取GeoJSON非常简单,只需要传入文件路径名称即可...过滤 bbox过滤允许我们read_file()传入一个边界框作为参数bbox,格式为(左下角x, 左下角y, 右上角x, 右上角y),这样在读入过程只会保留几何对象与bbox有相交数据记录

2K31

Python地信专题 | 基于geopandas空间数据分析-文件IO篇

本文使用所有数据都可以文章开头提及Github仓库对应本文路径下找到: 图1 2.1.1 shapefile 作为非常常见一种矢量文件格式,geopandas对shapefile提供了很好读取和写出支持...data.head() # 查看前5行 图3 缺少投影shapefileshapefile缺失.prj文件时,使用geopandas读入后形成GeoDataFrame会缺失crs属性:...当文件压缩包内根目录时,使用下面的语法规则来读取数据: zip://路径/xxx.zip 譬如我们要读取图7所示压缩包内文件: 图7 按照对应语法规则,读取该类型数据方式如下: 图8 而当文件压缩包内文件夹时...bbox有相交数据记录。...不像geopandas.read_file()可以根据传入文件名称信息自动推断类型,我们写出矢量数据时就需要使用driver参数来声明文件类型: ESRI Shapefile 我们将上文最后一次读入

2.1K20
您找到你想要的搜索结果了吗?
是的
没有找到

用Python编写小工具下载OSM路网数据

譬如在geofabrik[3] 中提供了各大洲、国家等大范围数据整合包,可以花费一定时间将其下载下来,再在需要哪些小区域时本地GIS软件或编程工具裁剪下所需范围路网。...具体代码部分本文不做详细说明,我将这部分功能封装到文章开头对应Github仓库下OsmDownloader.py[5],可以自行去下载并使用。 下面我们来学习如何在Python中使用它。...按照上文所述地址下载对应脚本之后,记住其所在路径,接着在其他脚本开头导入模块部分添加: import sys sys.path.append(r'所在路径') from OsmDownloader import...OsmDownloader 这样才能成功导入没有Python中注册过独立模块。...以上就是本文全部内容,对脚本获取或使用有疑问可以留言告诉我。 -END-

1.5K20

(数据科学学习手札80)用Python编写小工具下载OSM路网数据

通常我们可以 https://www.openstreetmap.org/export 中选择矩形区域内路网矢量数据进行下载,但这种方式对选择区域大小有一定限制,想获取较大范围区域路网数据并下载比较费事...GIS软件或编程工具裁剪下所需范围路网,但这种方式一是对电脑资源要求较高,譬如中国范围路网信息shapefile文件大小达到了800多兆,二是OSM路网信息不定期更新之后,要想及时跟上最新数据,就需要人工持续下载数据...)用Python编写小工具下载OSM路网数据/OsmDownloader.py ),可以自行去下载并使用,下面我们来学习如何在Python中使用它。...import sys sys.path.append(r'所在路径') from OsmDownloader import OsmDownloader   这样才能成功导入没有Python中注册过独立模块...如果担心中途网络连接原因导致中断,可以结合Python错误捕捉机制来进行相对应处理,比较简单这里就不再赘述。   以上就是本文全部内容,对脚本获取或使用有疑问可以留言告诉我。

1.6K10

一款具备SAM大模型AI分割,功能强大地理数据生产编辑查看工具 Geobuilding

· 它能加载和维护海量点/线/面数据· 它能设计GIS场景下三维漫游· 它能设计GIS场景下三维模型管理,模型定位/3D点序列· 它能导出支持属性映射第三方数据格式· 它能导出成GeoJSON、ShapeFile...、OSM格式· 它能导出丰富建筑物数据· 它用户主要来自政府、企业、高校、设计院等机构,数据完全本地化· 它同样可作为GIS数据标注来使用出色傻瓜化轮廓绘制能力,绘制精美的轮廓自由绘制、矩形绘制...更直观交互操作方式。轻松创建无缝地理网格,极简图像配准工具将JPG图像矢量化提供方便剪刀裁剪工具,像裁纸一样生产网格数据。...完美的轮廓识别,令人惊叹轮廓直角化设计使用高度框 批量设置建筑物高度支持要素属性扩展,为要素增加更多属性添加更多属性,满足业务系统对数据统计和筛选。...编辑 删除上海市60w+建筑物按轮廓裁剪支持GIS场景下三维漫游设计生产GIS数据有更多玩法,你可以软件自定义漫游路线,漫游镜头设计。

29910

NCL高效快速精准提取不规则区域内格点数据

而对于不规则范围,数据匹配精度有一定要求,譬如,需要严格按照某个特定区域shapefile文件来截取数据。...虽然,NCL官网提供了可行解决方案,但是 shapefile_mask_data(包含在shapefile_utils.ncl,官网有提供)也仅仅是较好地适用于2维Lat-Lon数据,对于3维或者更高维度数据...所以,针对于这个问题,实际操作我给出了一个快速处理方案,仅供参考: 1 load ".....) 34 f_eccc = addfile(filp1+"ECCC_TP_hr.nc","r") 35 f_ecmwf = addfile(filp1+"ECMWF_TP_hr.nc","r"...mask_io, (/2,3/)), 1) ;处理高维数组进行mask 6 copy_VarCoords(tp_erai,erai_mask) ;复制坐标信息 总结一下:由于使用了自带

6.6K30

Basemap系列教程:使用 shapefiles 文件裁剪栅格

程序 示例使用 pyshp [注3] 读取 shp文件。当然,ogr [注4] 也可以使用,但是 fiona 不行 [注5],因为相同脚本中使用 gdal 时会失败。...一个是包含点数组(脚本 vertices 变量),另一个是作用于每个点函数 此例,仅使用直线,因此 MOVETO 表示多边形开始点, LINETO 表示创建每一段直线,CLOSEPOLY 表示闭合多边形...当然只是用了 Andorra 边界构成多边形 prt 数组可以管理多个多边形,但此例只有一个 polygon(注:因为只有一个多边形),但仍可应用于裁剪多个多边形情况 使用 Path 函数创建...注意 transform = ax.transData 属性,这可以进行多边形坐标类型转换(此例为经纬度坐标) 48-49行进行裁剪操作。...通用白化程序 因为要经常对地图进行白化处理,因此可将上例中使用白化脚本整理为通用函数,这样以后使用中直接调用即可: def basemask(cs, ax, map, shpfile):

1.8K10

arcpy怎么用_python arcpy

r”c:\data” “\”在编码需转义,前面加r,直接按照原文意思输出。...地理处理:地理处理工具函数和非工具函数 工具 地理处理工具是指arcgis软件工具箱工具统称,每个工具都会给出相应图形界面、使用帮助和arcpy语法和使用方式,比如矢量裁剪工具: 打开空间分析裁剪工具图形界面...附加参数常常指处理过程设置条件,如上述裁剪函数附加参数是容差{cluster_tolerance}。裁剪时,我们有可能需要设置裁剪容差。很多时候,设置条件是默认选项,不需要进行设置。...使用列表函数可以轻松地将数据直接导入内存,进而交给工具函数处理。 差异性访问 矢量数据,最出名莫过于shapefile 文件了,其表现方式为: .shp – 存储要素几何主文件;必需文件。...String where_clause 用于限制游标返回可选表达式。有关 where 子句和 SQL 语句详细信息,请参阅查询表达式中使用元素 SQL 参考。

2.3K20

Python-Basemap核密度空间插值可视化绘制

上一篇推文我们使用geopandas+plotnine 完美绘制高斯核密度插值空间可视化结果,并提供了一个简单高效裁剪方法,具体内容点击链接:Python-plotnine 核密度空间插值可视化绘制...从结果我们可以看到,结果是规整网格数据,没有根据目标区域(地图文件) 对结果进行裁剪,接下来我们将使用fiona、shapely包 实现对目标区域裁剪操作。...「接下来就是关键一步操作」:我们根据df_grid数据经纬度信息判断点是否构建面(shp_ploygeon)内,不在点我们赋值为np.nan,点不变,这样即可完成“裁剪”操作,具体操作代码如下...Basemap可视化裁剪结果 可视化之前,我们需要对数据进行reshape操作,代码如下: mask_value_grid = df_grid["mask_value"].values.reshape...本期裁剪操作通用性较大,大家可以好好看下哦!下期我们使用R-ggplot+sf包实现空间插值操作,敬请期待

2.1K20

使用 geopandas 和 shapely(.shp) 进行地理空间数据处理和可视化

通过下面代码,我们可以使用 geopandas 库 read_file 函数来读取 Shapefile 文件,并将其保存为一个 GeoDataFrame 对象: import geopandas...shapely 库 box 函数构建了一个矩形框,并使用 difference 函数从 GeoDataFrame 清除了该矩形框内地理要素。...使用 gdf.plot() 函数将地理数据绘制图形上,ax=ax 参数指定图形窗口。 使用 plt.show() 函数显示图形窗口,展示可视化地图。...使用 gdf.difference() 函数,从原始地理数据删除 bbox 区域内地理要素,并将结果保存到 gdf 。10. 创建一个新图形窗口,大小为 10x10 英寸。...使用 gdf.plot() 函数将更新后地理数据绘制图形上,ax=ax 参数指定图形窗口。 使用 plt.show() 函数显示图形窗口,展示可视化地图。12.

1.4K10

【GEE】4、 Google 地球引擎数据导入和导出

要将它们带入 GEE,我们需要将它们转换为 shapefile。虽然有很多方法可以将 csv 文件转换为 shapefile,但我们将使用 R。下面的代码包含进行此转换所需所有内容。...上传 shapefile:在上面的 R 代码,我们将数据 csv 文件转换为 shapefile,并定义坐标参考系统 (CRS) 以匹配 GEE 预期 (WGS 1984)。...运行完此内容后,我们建议您加载自己 shapefile。这可以是您数据,或者如果您想要快速简便数据,请尝试使用Natural Earth Data shapefile 。...3.3定义天气变量 本课,我们使用 Google Earth Engine 作为将遥感数据(即我们栅格)与我们点位置相关联方法。...虽然我们可以 GEE 更多地使用这些数据,但很容易将它们导入 R 或 Excel。有几个选项可以定义导出数据最终位置。一般来说,将数据保存到 Google Drive 帐户是一个安全选择。

70521

Arcgis图像裁剪

使用arcmap对数据进行剪裁,Arcgis裁剪分为很多种,有矢量裁剪矢量,矢量裁剪栅格,栅格裁剪栅格。本文主要操作,掩膜裁剪(矢量裁剪栅格)和clip 裁剪。...2.使用这个工具有几点需要注意地方,一是输出范围,这里就选择我们面矢量数据;二是记得勾选使用输入要素裁剪几何选项;三是坐标范围获取;我们影像上移动光标,arcgis面板右下角有XY范围记录即可...3.现在已经建好shapefile文件了,点击上面菜单栏“编辑器”,点击“开始编辑”        4.右边菜单栏点击文件名,例如“测试”,下面会出现一个构造工具菜单栏。        ...shape数据相互转换:我们经常会在Google Earth获取影像数据,要将其Arcgis打开进行使用,经过分析后可能再会回到Google Earth进行对照分析,这就涉及到二者数据格式相互转换问题...二:工具说明 ArcGIS中导入KML(keyhole markup language),arctoolbox,转换工具和KML下都有KML到layer。 三:转换后裁剪,参照掩膜大法

2.1K50

Python绘制气象实用地图(续)

上一期,对Python绘制气象实用地图做了比较详细介绍,尽管已经能够满足部分需求了,但是,实际应用需求,可能还是别的需求,那么,今天就手把手教大家如何绘制几个省份白化等值线contour地图。...最后,对于QGIS强烈安利一波,不光它是免费,而且跨平台,也能够完美的支持Python3.7了,能够替代大部分日常使用ArcGIS功能,用起来不算很笨重!...再选定坐标系方案,最好和原来shp文件一致。我文末会提供相应地图文件!...2.直接在终端使用python xxx.py运行; 需要注意地方:很多人发现输出图片是没有经纬度坐标信息附加在网格线两端,怎么调都还是出不来。...还是那句话,遇到错误信息了,最值得信赖还是Google大法,学会如何使用Google,绝对是对debug有极大好处

5.2K65

【详细总结】cnmaps、maskout、salem正确打开方式

安装方法如下: conda install -c conda-forge cnmaps 唯一一个小小问题,就是使用过程必须保持shapely=1.8.5,不然部分函数使用时候会出现报错(如下)...使用 先给个亲身体验后总结: cnmaps自带了中国各级省市县shp,使用起来非常方便,不需要我们额外去下载shp文件了。...maskout使用方便,需要自己导入shp文件,但不支持掩膜选取数据,且速度稍慢。也就是说是完全裁剪:先画好图,然后根据shp文件从完整图中扣出自己感兴趣区域。...因为裁剪位置整张图片中是固定,但调整中心经度后图片偏移了,而裁剪位置没变,故裁剪结果就不对了。...综上,如果不想麻烦,只需要裁剪自己感兴趣区域,直接上maskout就行,但注意不要使用central_longitude等使图片偏移参数;如果处理省市县级数据,建议使用cnmaps;salem可以作为

65621

打开Shapefile文件正确方式

Shapefile文件简介 Shapefile文件是美国ESRI公司发布文件格式,因其ArcGIS软件推广而得到了普遍使用,是现在GIS领域使用最为广泛矢量数据格式。...shx索引文件用于存储几何数据索引信息,包含对主文件每个记录长度描述(注意不是空间索引) dbf表文件是使用dBase数据库表文件进行空间属性数据存储文件 所以,我们如果要自己完全从底层写代码解析...Shapefile文件的话,需要根据shx文件信息读取shp二进制数据并转化为几何对象,然后再读取dbf表格,将属性添加到几何对象上就完成了对一个Shapefile文件解析....英文好同学,请转移到这里:ESRI Shapefile Technical Desc GDAL矢量数据组织 GDAL栅格数据使用OGRDataSource表示(OGRDataSoruce是抽象类...GDAL空间要素模型是按照OGCSimple Feature规范实现,有兴趣童鞋可以参考官方文档:Simple Feature Access 使用GDAL打开Shapefile文件 下面的例子演示了如何打开

2.6K20
领券