Python获取栅格数据的空间坐标系统及EPSG编号

最新推荐文章于 2024-04-25 16:43:36 发布
最新推荐文章于 2024-04-25 16:43:36 发布
阅读量1.4k 收藏 6
点赞数 1
分类专栏： Arcgis and Python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/gengjuan520/article/details/114894535
版权
本文介绍如何使用Python解析栅格数据，获取其空间坐标系统信息和EPSG编号。
摘要由CSDN通过智能技术生成

1、代码如下

from  osgeo import  gdal,osr
#栅格文件存在的地址
f2 = 'C:\\Users\\zijiancao\\Desktop\\osgeopy-data\\osgeopy-data\\osgeopy-data\\Utah\\cache.tif'
ds2 = gdal.Open(f2)
#这里也可以直接用ds2.GetProjection()输出坐标系
sr2 = osr.SpatialReference(wkt=ds2.GetProjection()
最低0.47元/天 解锁文章
曹小小
关注
专栏目录
Python Gdal 栅格数据处理之定义投影系
weixin_43896283的博客
09-29 718
Python Gdal 栅格数据处理之定义投影系
【ArcGIS遇上Python】python批量获取栅格数据四至（top，bottom，left，right）坐标代码
「 刘一哥与GIS的故事」
11-01 3608
ArcGIS python批量获取栅格数据四至（top，bottom，left，right）坐标代码上图所示为ArcGIS自带的影像数据，存放路径为C:\Program Files (x86)\ArcGIS\Desktop10.6\ArcGlobeData\wsiearth.tif"，在源中可以查看该数据的四至坐标，那么，怎样用python批量获取多个栅格数据的四至坐标呢？实现方法：在python中编写如下代码，回车！代码运行结果： import arcpy arcpy.env.workspa
Python-GDAL教程：获取栅格影像的空间参考坐标系编号EPSG
Ethan的博客
08-21 6171
工作中遇到一个问题：需要把栅格影像中的空间坐标系提取出来，然后添加到对应的矢量影像中，而添加只需要知道空间参考坐标的编号就行，所以现在的关键就是提取原始栅格影像的空间参考坐标编号EPSG 1 方法一 安装好gdal后，直接在cmd中输入以下命令： gdalsrsinfo 01.tif # 前面是命令，后面的栅格名称 初步结果如下： 2 方法二 如果代码都是用python写的，所以当然想用py...
使用python的GDAL库将坐标系的wkt转epsg
ylfmsn的博客
06-15 619
坐标系wkt转epsg
geo2epsg:根据地理坐标提供有效的 EPSG 代码 (WGS84 UTM)
06-29
地理2epsg 描述 提供来自地理坐标的有效 EPSG 代码。 目前仅支持 WGS 84和 UTM。 安装 npm install geo2epsg 设置 var Geo2EPSG = require ( 'geo2espg' ) ; var geo2espg = new Geo2EPSG ( ) ; 用法 十进制度 度 分 秒 学分 此 MortonFox 的 Javascript UTM 转换库
基于Python 的简单栅格图像边界提取
兰小莫的博客
05-21 6404
在GIS中，栅格属性里有关于栅格自身的信息，背景（nodata value）对于识别一张图像的边界像元尤为重要，我们目的只要把每行每列中的第一次出现不是nodata的像元和最后一次出现nodata的前一个像元就可以了。对于栅格，可以用ArcPy中的RasterToNumpyArray函数将将栅格转成numpy数组，然后就可以按照所想读取出每行列中首尾像元。以下是部分代码提取边界像元的核心算法，其实
Python Gdal 栅格数据处理之投影与坐标转换
weixin_43896283的博客
06-15 1925
Python Gdal 栅格数据处理之投影与坐标转换
python 矢量数据转栅格数据代码实例
09-18
在Python中，将矢量数据转换为栅格数据是一个常见的地理信息系统（GIS）操作，用于处理地理空间数据。本文将详细介绍如何使用Python库，如`osgeo.osr`、`osgeo.ogr`和`pyproj`来实现这一转换，并提供相关代码实例。 ...
python中读取栅格数据_使用rasterio处理python中的地理空间栅格数据
weixin_26713457的博客
09-01 5983
python中读取栅格数据This article is meant to provide a quick introduction into how to use the Python package Rasterio for common tasks related to geospatial raster data. This is mainly a collection of things...
基于Python和GDAL实现栅格数据中相邻地物边界提取
seibert的博客
06-26 2337
对输入的影像计算指定像元值四邻域内指定像元值的像元数并输出成新文件 此代码块是打包成exe的源文件 1、导入相关python包 在网址 https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal 下载对应python版本的whl文件如GDAL‑3.0.0‑cp38‑cp38m‑win32.whl pip install numpy pip install GD...
基于高德地图api和Python的区县地理边界坐标提取
热门推荐
lzx123338的博客
05-06 1万+
在工作中，经常想用到类似于地热图的方式进行数据展示，奈何工作环境是内网，无法在线进行地图关联，没办法，只好自己想办法上网找边界坐标。 查了很多文档和费心以后，最终发现高德地图api的方法和方式最简单，正好初学Python，分享一下~
python提取二值栅格上边界和中线
alley51775的博客
11-17 2213
比较简单粗暴的方法…缺点就是提取上边界的时候有一些细碎的像元还是没办法删去，如果图案过分扭曲，效果就差很多。（初学python…代码冗长繁琐[扶额]） 适用于值为0和1的栅格图像。 1.提取上边界 原图像： 提取后： 代码： def Get_UpperBoundary(dataarray,rows,cols): boundary=dataarray.copy() boundary[boundary > 0] = 0 for i in range(cols-1):
python安装GDAL，解决ModuleNotFoundError: No module named ‘osgeo‘
记录的博客
01-13 7727
python安装GDAL，解决ModuleNotFoundError: No module named ‘osgeo‘
GDAL库进行坐标转换实例代码
欧特GO
05-19 1471
常用到的地理坐标系有WGS84，beijing54，xian80，CGCS2000这四种。GDAL可以像proj4那样自定义坐标系，也可以仅通过字符串定义一些常用的坐标系，但本文认为最方便的还是通过EPSG数据库定义的代码来定义一个地理坐标系统；
解决import gdal导入失败的问题,但from osgeo import gdal导入成功
qq_42870321的博客
06-20 1万+
解决import gdal导入失败的问题,但from osgeo import gdal导入成功。
python利用gdal读取、写出tif格式的遥感卫星影像，包含超大数据量的读写
weixin_44778507的博客
10-25 733
【代码】python利用gdal读遥感卫星影像。
Python地理数据处理 十四：栅格数据处理（二）
amyniez的博客，欢迎交流讨论
03-14 3403
目录1.地面控制点（GCP）2.多幅图像镶嵌3.颜色表 1.地面控制点（GCP）   如果想将一个旧的航片或者扫描地图，变成地理数据集，可以运用控制点来进行变换。图像转换类型的不同，所需要的地面控制点的数量也不同，常用方法：一阶多项式（至少3个点）、多项式变换（尽可能多的控制点）、样条插值法（对数据的不同部分使用不同的方程，可以精确拟合所提供的控制点，但可能导致图像其他部分扭曲）。   样条插值法：   使用GDAL自带的 gdal warp 程序进行各种插值方法。获取点坐标的过程比较艰难，但是必须人工
windows（python）下osgeo安装
LvManBa的博客
04-25 1902
下载的时候选择与Python版本、操作系统对应的whl，比如我的电脑是Windows10 64位、Python3.11.7，就下载GDAL-3.8.2-cp311-cp311-win_amd64.whl，其中 3.8.2表示的是GDAL的版本，cp311就是Python3.11的意思，win_amd64对应的是操作系统。下载完以后，将文件放到到python安装根路径的下的Scripts. 我的地址是： D:\java_soft\Python311\Scripts。
Python中将nc数据处理为栅格数据并叠加地理网格线
最新发布
08-04
在Python中，处理NetCDF (Network Common Data Form) 数据并将其转换为栅格数据，并叠加地理网格线通常涉及到使用一些专门的数据处理库，比如`xarray`、`rasterio` 和 `geopandas`。 首先，你需要安装这些库： ```bash pip install xarray rasterio geopandas matplotlib cartopy ``` 1. **读取NC文件**: 使用`xarray`库可以方便地加载NetCDF文件，如： ```python import xarray as xr dataset = xr.open_dataset('your_file.nc') ``` 2. **数据转换为栅格**: 如果你想将数据转换为栅格，需要先确定栅格尺寸、坐标系统以及分辨率。例如，如果数据是基于经纬度的，你可以使用`rasterio`创建一个新的GeoTiff文件： ```python from rasterio.crs import CRS from rasterio.plot import show from rasterio.warp import reproject, transform_bounds # 创建目标CRS target_crs = CRS.from_epsg(4326) # 世界经纬度坐标系 # 计算栅格范围和分辨率 min_lon, max_lon, min_lat, max_lat = dataset['lon'].min(), dataset['lon'].max(), dataset['lat'].min(), dataset['lat'].max() width, height = 1000, 1000 # 根据需求调整大小 res_x, res_y = (max_lon - min_lon) / width, (max_lat - min_lat) / height with rasterio.open('output.tif', 'w', driver='GTiff', crs=target_crs, transform=transform_bounds(target_crs, dataset.crs, min_lon, min_lat, width + res_x, height + res_y), width=width, height=height, dtype=data_array.dtype) as dst: reproject(dataset['variable_of_interest'], dst, src_transform=dataset['variable_of_interest'].transform) ``` 3. **叠加地理网格线**: 可以使用`cartopy`库绘制经纬度网格线，结合`matplotlib`画图： ```python import cartopy.crs as ccrs import matplotlib.pyplot as plt fig, ax = plt.subplots(subplot_kw={'projection': ccrs.PlateCarree()}) ax.coastlines(resolution='10m') # 加载海岸线 gl = ax.gridlines(draw_labels=True, linewidth=0.5, color='gray', linestyle='--') # 绘制网格线 gl.xlabels_top = gl.ylabels_right = False # 隐藏顶部和右侧标签 ax.imshow(rasterio.open('output.tif').read(1), transform=rasterio.open('output.tif').transform) plt.show() ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值