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

从R中已有的shapefile转换成crs后,如何保存shapefile?

在R中,可以使用sf包来处理shapefile文件,并将其转换为指定的坐标参考系统(CRS)。要保存shapefile文件,可以使用st_write()函数。

以下是保存shapefile的步骤:

  1. 首先,确保已经安装了sf包。如果没有安装,可以使用以下命令进行安装:
代码语言:txt
复制
install.packages("sf")
  1. 加载sf包:
代码语言:txt
复制
library(sf)
  1. 使用st_read()函数读取shapefile文件,并将其转换为sf对象:
代码语言:txt
复制
shapefile <- st_read("path/to/shapefile.shp")

请将path/to/shapefile.shp替换为实际的shapefile文件路径。

  1. 如果需要转换坐标参考系统(CRS),可以使用st_transform()函数。例如,将shapefile转换为WGS84坐标系(经纬度):
代码语言:txt
复制
shapefile <- st_transform(shapefile, crs = 4326)

请将4326替换为目标CRS的EPSG代码。

  1. 最后,使用st_write()函数保存shapefile文件。指定文件路径和文件格式(例如,"ESRI Shapefile"):
代码语言:txt
复制
st_write(shapefile, "path/to/save.shp", driver = "ESRI Shapefile")

请将path/to/save.shp替换为保存shapefile文件的路径和文件名。

综上所述,以上步骤演示了如何从R中已有的shapefile转换成指定的CRS,并保存shapefile文件。

页面内容是否对你有帮助?
有帮助
没帮助

相关·内容

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

本文将介绍如何使用 geopandas 和 shapely 来读取、处理、可视化和保存地理Shapefile文件。 1....=gdf.crs)) 用户输入了两对坐标之后,我们使用 shapely 库的 box 函数构建了一个矩形框,并使用 difference 函数 GeoDataFrame 清除了该矩形框内的地理要素...# 将删除的图形保存为新的 Shapefile 文件 new_shp_file = "海南省界/海南省界Export_Output.shp" gdf.to_file(new_shp_file) 分析说明...使用 gdf.difference() 函数,原始地理数据删除 bbox 区域内的地理要素,并将结果保存到 gdf 。10. 创建一个新的图形窗口,大小为 10x10 英寸。...定义一个变量 new_shp_file,指定保存Shapefile 文件的路径。 使用 gdf.to_file() 函数将更新的地理数据保存为新的 Shapefile 文件。

1.5K10

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

本文对应代码和数据上传至我的Github仓库https://github.com/CNFeffery/DataScienceStudyNotes 1 简介   在上一篇文章我们对geopandas...图3 缺少投影的shapefile   当shapefile缺失.prj文件时,使用geopandas读入形成的GeoDataFrame会缺失crs属性: ?...图4   如果已经知道数据对应的CRS,可以在读入数据后补充上crs信息以进行其他操作: import pyproj data.crs = pyproj.CRS.from_user_input('EPSG...图6 读取zip压缩包的文件 geopandas通过传入特定语法格式的文件路径信息,以支持直接读取.zip格式压缩包shapefile文件,主要分为两种情况。   ...图21   但我观察到即使出现了上述错误,GeoPackage文件也是成功保存到路径下的且整个程序并未被打断,因此可以无视上述错误: ?

2K31

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

本文对应代码和数据上传至我的Github仓库https://github.com/CNFeffery/DataScienceStudyNotes 1 简介 在上一篇文章我们对geopandas的坐标参考系有了较为深入的学习...crs data.head() # 查看前5行 图3 缺少投影的shapefileshapefile缺失.prj文件时,使用geopandas读入形成的GeoDataFrame会缺失crs属性...:4326') data.crs 图5 直接读取文件夹 当文件夹下只有单个shapefile时,可以直接读取该文件夹: 图6 读取zip压缩包的文件 geopandas通过传入特定语法格式的文件路径信息...,以支持直接读取.zip格式压缩包shapefile文件,主要分为两种情况。...这里需要注意一个bug:在使用geopandas导出GeoPackage文件时,可能会出现图21所示错误: 图21 但我观察到即使出现了上述错误,GeoPackage文件也是成功保存到路径下的且整个程序并未被打断

2.1K20

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

1简介 在本模块,我们将讨论以下概念: 如何将您自己的数据集引入 GEE。 如何将来自遥感数据的值与您自己的数据相关联。 如何 GEE 导出特征。...3将您自己的数据带入 Earth Engine 在本练习,我们将讨论如何将您自己的数据移动到 GEE、数据集中提取值以及 GEE 中导出这些值。...可以在此处找到有关如何R 中将 csv 文件转换为 shapefile 的详细信息。 代码的一些复杂性来自重命名列以删除“.”。为了符合 GEE 对命名约定的要求,这是必要的。...上传 shapefile:在上面的 R 代码,我们将数据的 csv 文件转换为 shapefile,并定义坐标参考系统 (CRS) 以匹配 GEE 的预期 (WGS 1984)。...运行完此内容,我们建议您加载自己的 shapefile。这可以是您的数据,或者如果您想要快速简便的数据,请尝试使用Natural Earth Datashapefile

74621

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

上一期,对Python绘制气象实用地图做了比较详细的介绍,尽管已经能够满足部分需求了,但是,在实际的应用需求,可能还是别的需求,那么,今天就手把手教大家如何绘制几个省份的白化等值线contour地图。...PathPatch import matplotlib.pyplot as plt from osgeo import gdal import numpy as np import cartopy.crs...此外,建议保存的图片和脚本名称一致,解决方案: #代码头部 import os,sys #代码尾部 (filename, extension) = os.path.splitext(os.path.basename...还是那句话,遇到错误信息了,最值得信赖的还是Google大法,学会如何使用Google,绝对是对debug有极大好处的。...PathPatch import matplotlib.pyplot as plt from osgeo import gdal import numpy as np import cartopy.crs

5.2K65

(数据科学学习手札89)geopandas&geoplot近期重要更新

本文示例代码及数据上传至我的Github仓库https://github.com/CNFeffery/DataScienceStudyNotes 1 简介   最近一段时间(本文写作于2020-07...2 geopandas&geoplot近期重要更新内容 2.1 geopandas近期重要更新 2.1.1 新增高性能文件格式   geopandas0.8.0版本开始,在矢量文件读写方面,新增了.feather...安装完成,我们就来一睹这些新功能的效率如何,首先我们创建一个足够大的虚拟表(200万行11列),并为其新增点要素矢量列: import numpy as np from shapely.geometry...base.progress_apply(lambda row: Point(row['_10'], row['_11']), axis=1) # 添加矢量列 base = gpd.GeoDataFrame(base, crs...()进行过介绍,但在先前的版本只能使用固定的少数几种内置的在线地图,而在最近的版本,webplot()的底图叠加方式进行了非常大的调整,使得可以利用参数provider来像folium那样自由切换底图

83620

geopandas&geoplot近期重要更新

❝本文示例代码及数据上传至我的Github仓库https://github.com/CNFeffery/DataScienceStudyNotes ❞ 1 简介 最近一段时间(本文写作于2020-...2 geopandas&geoplot近期重要更新内容 2.1 geopandas近期重要更新 2.1.1 新增高性能文件格式 geopandas0.8.0版本开始,在矢量文件读写方面,新增了.feather...他们都是Apache Arrow项目下的重要数据格式,提供高性能文件存储服务,使得我们可以既可以快速读写文件,又可以显著减少文件大小,做到了“多快好省”: 图1 在将geopandas更新到0.8.0版本,...安装完成,我们就来一睹这些新功能的效率如何,首先我们创建一个足够大的虚拟表(200万行11列),并为其新增点要素矢量列: import numpy as np from shapely.geometry...base.progress_apply(lambda row: Point(row['_10'], row['_11']), axis=1) # 添加矢量列 base = gpd.GeoDataFrame(base, crs

76630

【GEE】时间序列多源遥感数据随机森林回归预测|反演|验证|散点图|完整代码

确保将其table重命名为VT_boundary,然后将代码保存到本地存储库。...加载 AOI pedons shapefile 在开始之前,需要将样本数据VT_pedons.shp作为assets加载到 GEE ,并导入到我们的代码,以便接下来在回归中使用。...转到左侧面板的assets选项卡,找到shapefile,然后单击它,此时会弹出预览: 导航到“features”选项卡并浏览 shapefile 的不同属性。...单击右下角的蓝色IMPORT按钮将其添加到我们的代码。 此时会看到一个table变量添加到顶部的导入列表。我们将新 shapefile 的名称默认的table更改为VT_pedons。...: crs, maxPixels: 1e13 }); 运行脚本,窗格右侧的“任务”选项卡将变为橙色,表示有可以运行的导出任务。

1.2K22

GPM 降雨量数据处理 -R(坐标系转换)

如何下载NASA降雨量数据,见此链接。 这里需要强调的一点就是,降雨数据主要在NASA网站主要包括TRMM与GPM项目 下载的数据是HDF5格式,如何R读取HDF与tc文件,戳here。...TRAMM与GRM下载的HDF5格式在R,会出现坐标与我们常用坐标系不一致的情况, 主要投影坐标系不同。 所以这篇文章,这要介绍raster如何转换成常规的4236坐标系。...Raster文件,找到储存在HDF5文件的precipitation位置。...将rasterNoProj 转换成数据库data.frame,包含了x,y坐标信息。 然后我们之前旋转的s2也转换data.frame格式。...但是s2数据太大,转换成sf时间较长, 先喝口水。慢慢等待。 缺点,在制图过程,也需要很长时间才能出图。

1.1K21

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

安装方法如下: conda install -c conda-forge cnmaps 唯一一个小小的问题,就是在使用过程必须保持shapely=1.8.5,不然部分函数在使用时候会出现报错(如下)...具体如何命名,可以根据import salem下面的报错,找到这个~/anaconda3/envs/cnmaps/lib/python3.10/site-packages/salem/utils.py文件...这是修改再次运行import salem,会告诉我们应该修改的具体名字。...使用 先给个亲身体验的总结: cnmaps自带了中国的各级省市县shp,使用起来非常方便,不需要我们额外去下载shp文件了。...也就是说是完全的裁剪:先画好图,然后根据shp文件完整的图中扣出自己感兴趣的区域。注意,用maskout时候无法调整中心经度。

67321

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

图5 Step3: 渲染路网: 图6 了解了上述步骤之后,我们利用requests、geopandas等库仅用不到100行代码就可以参考上述过程,提取所需的shapefile或GeoJSON文件保存到本地...具体的代码部分本文不做详细说明,我将这部分功能封装到文章开头对应Github仓库下的OsmDownloader.py[5],可以自行去下载并使用。 下面我们来学习如何在Python中使用它。...按照上文所述的地址下载对应脚本之后,记住其所在路径,接着在其他脚本开头导入模块部分添加: import sys sys.path.append(r'所在路径') from OsmDownloader import...(path='保存路径') # path参数控制文件保存的路径 程序运行稍等片刻即可完成下载(具体的耗时取决于你的网络状况)。...(path='保存路径') 等待一段时间,我们area_list里多个不同级别行政区的路网数据便下载完成: 图10 如果担心中途网络连接原因导致中断,可以结合Python的错误捕捉机制来进行相对应的处理

1.5K20

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

了解了上述步骤之后,我们利用requests、geopandas等库仅用不到100行代码就可以参考上述过程,提取所需的shapefile或GeoJSON文件保存到本地,具体的代码部分本文不做详细说明,我将这部分功能封装到文章开头对应...可以自行去下载并使用,下面我们来学习如何在Python中使用它。...请提前安装完成requests、geopandas、pandas以及tqdm这几个三方库),记住其所在路径,接着在其他脚本开头导入模块部分添加: import sys sys.path.append(r'...(path='保存路径') # path参数控制文件保存的路径   程序运行稍等片刻即可完成下载(具体的耗时取决于你的网络状况),譬如这里我花了不到20秒就完成成都市路网数据的下载: ?   ...(path='保存路径')   等待一段时间,我们area_list里多个不同级别行政区的路网数据便下载完成: ?

1.6K10
领券