使用Rasterio做投影变换

最新推荐文章于 2024-08-14 11:50:35 发布

卡尔曼和玻尔兹曼谁曼

最新推荐文章于 2024-08-14 11:50:35 发布

阅读量3.2k

点赞数 2

分类专栏： GIS Python # Python空间数据处理实战

本文为博主原创文章，转载请注明原文出处！

本文链接：https://blog.csdn.net/theonegis/article/details/80663218

版权

本文介绍了如何使用Rasterio库进行投影变换。作者阿振详细阐述了思路，指出在Rasterio中可以方便地计算目标空间参考系统的放射变换参数和图像尺寸，从而简化了投影转换的过程。文章提供了一个从计算到更新元信息的完整代码实现。

摘要由CSDN通过智能技术生成

使用Rasterio做投影变换

作者：阿振

邮箱：tanzhenyugis@163.com

博客：https://blog.csdn.net/theonegis/article/details/80089375

修改时间：2018-06-11

声明：本文为博主原创文章，转载请注明原文出处

思路分析

在之前GDAL系列文章中的《栅格数据投影转换》提到过，做投影转换最重要的是计算数据在目标空间参考系统中的放射变换参数（GeoTransform）和图像的尺寸（行数和列数）。而且我们使用GDAL基本库自己写代码进行了计算。

在rasterio中提供了calculate_default_transform，可以直接计算目标系统中的放射变换参数和图像尺寸。

这样我们直接根据计算的结果更新目标文件的元信息即可。

代码实现

import numpy as np
import rasterio
from rasterio.warp import calculate_default_transform, reproject, Resampling

最低0.47元/天解锁文章

确定要放弃本次机会？

福利倒计时

: :

立减 ¥

普通VIP年卡可用

立即使用

卡尔曼和玻尔兹曼谁曼

关注关注

2
点赞
踩
11

收藏

觉得还不错? 一键收藏
1
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

专栏目录

Python 使用rasterio将数组按照指定格式写入内存

糯米丸子TAT的博客

10-28

928

在进行栅格矢量数据处理时，我们常常需要将数组转化为栅格影像，再对栅格数据进行后续处理，输出最终结果。常见的做法是先在本地磁盘生成栅格数据，再进行处理输出。但是，如果数组转化的栅格数据不是最终结果，存储在磁盘中会浪费空间，且整套代码显得不是很优雅。此时，我们可以利用rasterio的MemoryFile将数据写入内存

投影变换

Beware beginning

03-25

228

仿射变换后平行的线还是平行的 投影变换只是保证同一条直线上的点还是在同一条直线上但不再保证平行了因为投影变换是一个二维图像（车牌）经过一个三维变换，然后映射到另外一个二维空间，二维图像的二维空间与映射后的二维空间不一样，如果一样，就是仿射变换 import cv2 import numpy as np #图片的路径 imgname = "img1.jpg" #读取图片 image = cv2.imread(imgname, cv2.IMREAD_COLOR) #图片的高度和宽度 h,w =

1 条评论您还未登录，请先登录后发表或查看评论

[python]rasterio运行代码警告proj_create_from_database: Cannot find proj.db

最新发布

FL1623863129的博客

08-14

881

这个报错要分原因还有rasterio版本讨论，因此官方给出了十分具体回答。

OpenCV : 投影变换

机器会学习的博客

03-21

3449

投影变换 物体在三维空间发生了旋转,叫做投影变换.由于可能出现阴影或遮挡,所以变换后较难复原.如果物体时平面的,就可以通过投影变换物体三维变换进行模型化,这叫做专用的二维投影变换. 矩阵表示: [x~y~z~]=[a11a12a13a21a22a23a31a32a33][xyz]\left[\begin{matrix}\tilde{x} \\ \tilde{y}\\ \tilde{z}\end{m...

投影变换[转]

weixin_30868855的博客

03-27

305

4. 投影变换 一个投影变换类似于选择一个透镜的焦距，它是三种变换中最复杂的一个。 4.1 视锥（Viewing Frustum） 4.2 什么是投影变换? 4.3 设置投影矩阵 4.4 一个“W-Friendly”投影矩阵　 4.1 视锥　　视锥就是场景中的一个三维空间，它的位置由视口的摄像机来决定。这个空间的形状决...

使用Rasterio读取栅格数据的实例讲解

09-18

在Rasterio 1.0之后，它不再使用GDAL的传统放射变换方式，而是采用了`affine`库的风格。`affine.Affine(a, b, c, d, e, f)`表示放射变换，与GDAL中的参数顺序(c, a, b, f, d, e)不同。这种新的放射变换模型更直观，...

GIS.rar_gis_python 投影_python 投影变换_坐标系

09-23

Python作为一种强大的编程语言，被广泛应用于GIS领域，尤其是对于数据处理和投影变换。投影是将地球表面的三维空间坐标转换为二维平面坐标的过程，以便于在地图上显示和分析。地球是一个近似的球体，而我们通常在...

地图投影变换代码

12-24

在进行地图投影变换时，我们需要使用编程语言和相关的GIS库。例如，Python中的`geopandas`、`rasterio`或`pyproj`库可以方便地实现这些变换。`geopandas`用于处理地理空间数据，`rasterio`处理栅格数据，`pyproj`则...

使用Rasterio读取栅格数据

给永远比拿愉快

06-08

1万+

使用Rasterio读取栅格数据作者：阿振邮箱：tanzhenyugis@163.com 博客：https://blog.csdn.net/theonegis/article/details/80089375 修改时间：2018-06-06 声明：本文为博主原创文章，转载请注明原文出处 Rasterio简介有没有觉得用GDAL的Python绑定书写的代码很不Pyth...

图像处理中的投影变换（Perspective Transformation）

热门推荐

数据之美的博客

11-27

4万+

透视变换(Perspective Transformation)是将图片投影到一个新的视平面(Viewing Plane)，也称作投影映射(Projective Mapping)。通用的变换公式为： u,v是原始图片左边，对应得到变换后的图片坐标x,y,其中。变换矩阵可以拆成4部分，表示线性变换，比如scaling，shearing和ratotion。用于平移，产生透视变换。

华科_图形学笔记_07_投影变换

精密时空设计局

04-09

1700

计算机图形学_华中科技大学_中国大学MOOC(慕课) 7.1_有趣的投影由于我们的观察设备是二维的，因此，最终都存在一个三维向二维的映射过程。这个映射过程，按照几何意义来说就是投影。而在我们之前讲到的观察变换中，其实还有一个隐含着的概念，观察平面。那么这个平面，也可以叫做投影平面，我们看到的景物，都需要投影到这个平面上。一个是平行投影。工程上我们用到的三视图，学习几何的时候用到的轴测图，都属于平行投影。与平行投影相对而言的是透视投影。最初接触这个概念应该是接

（二十五）ArcGIS空间数据的转换与处理——投影变换

qq_42540466的博客

05-09

1万+

原始数据往往由于在数据结构、数据组织、数据表达等方面与用户需求不一致而要进行转换与处理。本节主要介绍 ArGIS 中数据的投影变换内容。

计算机图形学(五)-MVP变换之投影(Projection)变换

点燃火柴的博客

07-22

3615

MVP变换之投影[Projection]变换1 投影变换介绍2 正交投影2.1 正交投影介绍2.2 正交投影矩阵推导2.3 正交投影矩阵3 透视投影3.1 透视投影介绍3.2 透视投影矩阵推导3.2.1 复习齐次坐标3.2.2 透视投影矩阵推导思路3.2.3 透视投影到正交投影的变换矩阵推导3.3 透视投影矩阵3.4 FOV和Aspect ratio表示透视投影矩阵 1 投影变换介绍之前已经提到MVP变换指的是：模型变换(model)、视图变换(view)、投影变换(projection)。这三大变换好比

Halcon 第七章『图像的几何变换』◆第2节：投影变换

Halcon机器视觉

06-10

5093

投影变换也叫透射变换、投影映射。透射变换是将图像投影到一个新的视平面，是一种二维坐标到三维坐标的变换。透射变换是仿射变换的延续，也可以说仿射变换是透射变换的一种特殊形式。其特殊性在于变换后图像的形状仍然维持原状。投影变换包括的情况很多，有可能变换前后图像的形状发生了很大的改变，如对边不再平行，或者发生了透视畸变，这时可以使用投影变换使其恢复原状。其步骤与仿射变换类似，首先计算投影变换矩阵，然后计算投影变换参数，最后将投影变换矩阵映射到对象上。透射变换是三维空间上的变换，因此，

【计算机图形学理论】三维观察-投影变换

weixin_46013570的博客

04-22

3471

01 三维观察流水线 02 三维观察变换-模视变换 1.如何进行三维场景观察？物体在世界坐标系内的物理位置并没有发生改变。改变观察点的位置。 2.观察变换-建立观察（相机）坐标系

计算机图形学 | 投影变换

ProgramNovice的博客

04-27

2708

计算机图形学 | 投影变化

rasterio库python3.9

09-12

在Python 3.9中使用rasterio库进行地理空间数据处理非常方便。下面是一些常见的用法和示例： 1. 读取地理空间数据：使用rasterio的`open`函数可以打开和读取栅格数据。可以通过指定文件路径或URL来读取数据。例如，使用以下代码读取数据： ```python import rasterio with rasterio.open('path/to/file.tif') as src: data = src.read() ``` 这将返回一个多维数组`data`，其中包含栅格数据的像素值。你还可以使用`read`函数的其他参数来指定读取的窗口、缩放或裁剪等操作。 2. 写入地理空间数据：使用rasterio的`open`函数和`write`方法可以将数据写入新的栅格文件。以下是一个示例： ```python import rasterio with rasterio.open('path/to/output_file.tif', 'w', **src.profile) as dst: dst.write(data) ``` 在这个例子中，`**src.profile`用于复制源文件的元数据（如投影、坐标参考系统等）到输出文件。 3. 栅格数据操作： rasterio提供了许多功能强大的函数和方法来操作栅格数据。例如，你可以使用`numpy`库来计算栅格数据的统计信息、裁剪数据或应用地理变换。以下是一些示例： - 计算栅格数据的统计信息： ```python import rasterio import numpy as np with rasterio.open('path/to/file.tif') as src: data = src.read() statistics = np.mean(data), np.std(data), np.min(data), np.max(data) ``` - 裁剪栅格数据： ```python import rasterio from rasterio import windows with rasterio.open('path/to/file.tif') as src: window = windows.Window(100, 100, 200, 200) # 定义裁剪窗口的位置和大小 data = src.read(window=window) ``` - 应用地理变换： ```python import rasterio import numpy as np with rasterio.open('path/to/file.tif') as src: data = src.read() transform = src.transform new_data = np.dot(data, transform) ``` 希望以上信息对你有所帮助。如果你有其他问题，请随时提问。