HouGISer的博客

文章目录背景解决方案背景不同数据源做重采样或者投影变换后，用同一掩膜进行提取，往往会出现范围不一致的情况。这里说的范围不一致是指图像的上下左右平面坐标不完全一致，会出现细微的差别，导致栅格像素无法完全重叠。本文教大家如何解决该问题。解决方案我调研过rasterio，发现并没有解决方案。calculate_default_transform函数无法固定输出图像范围，只能固定输入图像范围，而输入图像因为不同源，往往是没法固定的。我的解决方案是引入arcpy包。arcpy是arcmap中的toolbo

文章目录背景实战背景在测绘领域中提到的坐标系有两种，一种是地理坐标系，用经纬度高程来表达；另一种是投影坐标系，即经过投影变换后的平面坐标系，通常是xy表达。坐标系转换有四种情况：地理坐标系转投影坐标系（也叫大地坐标正算）投影坐标系转地理坐标系（也叫大地坐标反算）一种地理坐标系转另一种地理坐标系一种投影坐标系转另一种投影坐标系本文的坐标系转换理论上涵盖以上四种情况。实战import numpy as npimport rasteriofrom rasterio.warp impor

文章目录背景实战背景图像重采样是指改变图像的分辨率，通常指分辨率变小，即像素的宽高变长，是科研常用的操作。比如我们手头只有300m分辨率的图像，但是为了统一不同数据源的数据精度，需要统一为1000m分辨率，这时候就需要重采样。实战import numpy as npimport rasteriofrom rasterio.warp import calculate_default_transform, reprojectfrom rasterio import crsdataset = r

文章目录背景实战背景掩膜提取是指基于矢量面范围内的栅格像素值，并输出为新图像的操作。因为涉及矢量面数据获取，所以需要引入fiona包，未安装的读者请自行安装。实战import fionaimport rasterioimport rasterio.maskmask = 'D:/A_2021寒假/城市群相关/Data/0Slab中国基础地理数据/China/China_single.shp'src_img = 'input.tif'dst_img = 'output.tif'# 读取

文章目录背景实战注意事项背景图像镶嵌指讲几张栅格图像镶嵌到同一张图像中，也就是拼接图像。栅格图像是有明确的坐标系和空间范围的，因此待镶嵌的图像必须有一样的地理坐标坐标系或者投影坐标系，否则是无法镶嵌在同一坐标系中的。本文使用的是gdal，而不是rasterio，但rasterio是基于gdal的，本文默认你已经安装好了gdal，没有装的请移步我的geopandas安装教程。实战from osgeo import gdal, gdalconsttifPath = './data/' # 待融合的

文章目录背景安装简介基操图像读取图像写入展望背景最近在做一些和图像处理相关的工作，涉及利用rasterio来实现遥感影像的镶嵌、重采样、重投影等，故开一个系列记录下来。rasterio是栅格图像处理的python包，依赖GDAL，如果你想了解如何二安装GDAL，可以转到我的另一篇博文。安装在GDAL安装完成的基础上，到这个网站下载对应版本的rasterio包，然后运行以下命令安装whl文件，即可安装好rasterio。pip install xxx.whl简介rasterio的完整接口可以