python 遥感影像的转投影

最新推荐文章于 2023-04-05 17:16:32 发布
最新推荐文章于 2023-04-05 17:16:32 发布
阅读量523 收藏 2
点赞数
分类专栏： 地理信息系统 文章标签： gis
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42869874/article/details/117030078
版权
1.1 遥感影像Albers等积投影转UTM投影
from osgeo import gdal
path_src = 'H:\\bomi_data\\1KM\\2018\\DAY\\'
dataset = gdal.Open(path_src+'1.tif')  
# 待转换的影像（albers等积投影）
gdal.Warp('F:\\work6\\1_2.tif', dataset, dstSRS='EPSG:32646') 
# 将原影像转换为UTM 46N投影(EPSG:32646), 第一个路径是生成影像的位置
花前墨
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
8影像计算ndvi landsat_基于Python处理遥感影像（1）
weixin_29227425的博客
01-04 1816
最近发现一个遥感数据的模块Pyrsgis，该模块建立在GDAL库的基础上，读取、操作和导出栅格数据比较方便。从地理空间数据云下载Landsat8的Level 1T地形矫正影像的压缩包.tar.gz。pyrsgis可以直接从压缩包中读取相应信息：importosimportnumpyasnpimportpyrsgisasrgimportcv2os.chdir("D:\\la...
遥感影像镶嵌的Python实现
02-25
Python和GDAL的结合使得遥感影像的处理变得更加便捷和高效，无论是小规模的实验还是大规模的生产任务，都能应对自如。通过不断学习和实践，你可以进一步提升遥感影像处理的能力，更好地服务于地理空间数据分析和应用...
pyhton批量处理modis数据最大值合成，转投影，重采样转asc
06-22
在modis网站下载全球数据，为日值数据，每一天的数据大概为500M到1G不等，一年数据365天，大概300多G，数据量算比较大的了。原始数据原来分辨率0.05度，现在需要按照半月合成，即得到1年24期数据，结果转投影成1984，重采样成0.5度。 考虑用python进行批量操作，直接得到最终asc结果，可以在matlab进行批量操作。除了最大值合成，还可以平均值合成，代码并且考虑了闰年语非闰年的差异结果举例，如2008年1月最大值合成文件名：mod09cmg200801a一月上旬;mod09cmg200801b一月下旬。整个程序处理过程中显示每处理完成一个月所用的时间。 对基础数据的处理很方便，不需要打开ArcGIS,安装了GIS，自动有python软件，只需要在开始菜单——ARCGIS——PYHON2.7下打开python，菜单下open打开hdf2asc_15day_max_prj_resample.py文件，修改年份，文件夹路径以及目标投影参考影像路径即可运行。
基于Python依据遥感影像的投影方式对矢量进行重投影
qq_43444176的博客
04-05 251
【代码】基于Python依据遥感影像的投影方式对矢量进行重投影。
Python地学分析 — GDAL对遥感影像重投影
一叶知秋的博客
12-20 9321
欢迎关注博主的微信公众号：“智能遥感”。 该公众号将为您奉上Python地学分析、爬虫、数据分析、Web开发、机器学习、深度学习等热门源代码。 本人的GitHub代码资料主页（持续更新中，多给Star，多Fork）： https://github.com/xbr2017 CSDN也在同步更新： https://blog.csdn.net/XBR_2014 “ 对遥感影像重投影是遥感...
python GDAL遥感影像投影变换
llc的博客
10-02 2424
gdalwarp gdalwarp ttt.tif out.tif -s_srs EPSG:4326 -t_srs EPSG:3857 gdal.Warp from osgeo import gdal options = gdal.WarpOptions(format='GTiff', srcSRS='EPSG:4326', dstSRS='EPSG:3857') gdal.Warp('out.tif', 'ttt.tif', options=options)
基于Python-GDAL的遥感影像镶嵌脚本
03-29
Python作为一门强大的编程语言，结合GDAL库，为处理遥感影像提供了便捷而高效的方法。本文将深入探讨如何利用Python和GDAL进行遥感影像的镶嵌操作。 GDAL，全称Geospatial Data Abstraction Library，是一个开源的...
利用shp矢量数据裁剪遥感影像的Python实现
02-25
本文将深入探讨如何使用Python语言，结合Shp矢量数据对遥感影像进行裁剪操作，以便提取特定区域的信息。Shp文件是ESRI（Environmental Systems Research Institute）公司开发的地理信息系统中广泛使用的矢量数据格式...
Python 实现遥感影像波段组合的示例代码
09-18
Python作为一门功能强大的编程语言，在遥感影像处理领域也得到了广泛的应用。Python提供的第三方库gdal（Geospatial Data Abstraction Library）是专门用于读写栅格和矢量地理空间数据格式的库，非常适合进行遥感...
遥感影像png格式转tif格式
05-18
转换过程中，我们需要注意保持图像的元数据，如地理坐标系、投影信息等，这在遥感影像处理中非常重要。如果PNG文件中包含了这些信息，我们需要在创建TIFF文件时正确设置，以便于后续的分析工作。 通过这个过程，...
python2.7栅格数据批量转换投影
09-15
python2.7栅格数据批量转换投影:ProjectRaster_management (in_raster, out_raster, out_coor_system, {resampling_type}, {cell_size}, {geographic_transform}, {Registration_Point}, {in_coor_system}) in_raster 输入栅格数据集。Mosaic Layer; Raster Layer out_raster 要创建的输出栅格数据集。以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：.bil - Esri BIL, .bip - Esri BIP, .bmp - BMP, .bsq - Esri BSQ, .dat - ENVI DAT,.gif - GIF,.img - ERDAS IMAGINE,.jpg - JPEG,.jp2 - JPEG 2000,.png - PNG,.tif - TIFF,无扩展名 - Esri Grid,以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。 将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。 Raster Dataset out_coor_system 输入栅格待投影到的目标坐标系。默认值将基于“输出坐标系”环境设置进行设定。该参数的有效值是扩展名为 .prj 的文件。现有要素类、要素数据集、栅格目录（基本上包含了与坐标系相关的所有内容）。坐标系的字符串表示。要生成此类较长的字符串，可向模型构建器添加一个坐标系变量，并根据需要设置该变量的值，然后将模型导出到 Python 脚本。 Coordinate System resampling_type (可选) 要使用的重采样算法。默认设置为 NEAREST。 NEAREST —最邻近分配法 BILINEAR —双线性插值法 CUBIC —三次卷积插值法 MAJORITY —众数重采样法 NEAREST 和 MAJORITY 选项用于分类数据，如土地利用分类。NEAREST 选项是默认设置，因为它是最快的插值法，同时也因为它不会更改像元值。请勿对连续数据（如高程表面）使用 NEAREST 或 MAJORITY。BILINEAR 选项和 CUBIC 选项最适用于连续数据。不推荐对分类数据使用 BILINEAR 或者 CUBIC，因为像元值可能被更改。 cell_size (可选) 新栅格数据集的像元大小。默认像元大小为所选栅格数据集的像元大小。 Cell Size XY geographic_transform (可选) 在两个地理坐标系或基准面之间实现变换的方法。当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。 有关各个受支持的地理（基准面）变换的详细信息，请参阅位于 ArcGIS 安装目录的 \Documentation 文件夹下的 geographic_transformations.pdf 文件。 Registration_Point(可选)用于对齐像素的 x 和 y 坐标（位于输出空间中）。配准点的工作原理与捕捉栅格的概念类似。通过配准点可指定用于定位输出像元的原点，而不是仅将输出捕捉到现有栅格像元。所有输出像元与该点之间必须间隔一个像元。该点的坐标不必位于一角，也不必落入栅格数据集中。捕捉栅格环境设置参数将优先于 Registration_Point 参数。因此，如果您要设置配准点，请确保尚未设置捕捉栅格。 in_coor_system (可选) 输入栅格数据集的坐标系。Coordinate System
投影转换工具—Albers等面积.
04-16
投影转换工具—Albers等面积.投影转换工具—Albers等面积.Albers等面积
批量定义投影转换（python）
diz5280的博客
09-01 1585
1.打开ArcMap的python工具 2.输入 # --------------------------------------------------------------------------- # ProjectionConversion.py # Created on: 2019-09-1 # (generated by ArcGIS/ModelBu...
Python读取nc文件转tif
江北20190411的博客
09-10 9514
import netCDF4 as nc data = r'Y:\folder\H08_20200601_0000_1D_ROC010_FLDK.02401_02401.nc' nc_data_obj = nc.Dataset(data) nc_data_obj Out[4]: <class 'netCDF4._netCDF4.Dataset'> root group (NETCDF4 data model, file format HDF5): title: Himawari-
Python rasterio.merge遥感影像融合
weixin_42915923的博客
12-27 3242
最近在尝试用Python做遥感影像融合，用到了rasterio.merge.merge（简称merge）功能，可能还有gdal和arcpy，但是还没有用过。 问题： merge虽然是非常好用的图像融合工具，但是它只提供四种内置融合方式，first,last,min和max，即 MERGE_METHODS = { 'first': copy_first, 'last': copy_last, 'min': copy_min, 'max': copy_max } 因为我想用均值
如何使用Python中的GDAL库对遥感影像进行读取和存储
热门推荐
小E博客
11-17 1万+
与自然影像稍有不同，遥感影像一般有多个波段(波段数大于4)，数据格式多为tiff，比如国产高分一号/二号遥感影像，多光谱影像有4个波段，分别为蓝(1)、绿(2)、红(3)以及近红外(4)波段，真彩色是321显示，假彩色是432显示。 利用遥感或GIS处理软件，如ENVI、ERDAS、ArcGIS等，可以很方便快捷地实现对图像进行读取、处理等操作。 然而，很多时候我们需要对大量的图像数据进行批量处理操作，这时候就需要我们自己编写代码来实现。这篇博客就简单地介绍，如何在Python中实现遥感影像的读取和存储
Python遥感开发之使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
等待着冬天的风的博客
01-04 1926
前言：转载自《使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换》 使用GDAL库写了四个函数分别进行投影坐标与地理坐标（经纬度）之间的转换，投影坐标和图上坐标（行列号）之间的转换。有需要的朋友可以参考。
Python遥感图像处理之图像基本信息读取
bailongyang的博客
03-06 4308
本科大四，小白一枚，编程水平不好，但是现在我有一个大胆的想法，想未来在遥感算法方面有所突破，我要开始好好学了，祝我成功！也祝所有同行节节高升！如有不足，敬请指正！ 工具准备：Vs Code、Landsat 8.tiff 以典型的多光谱数据Landsat8为数据源。 在借助python读取遥感图像时，会用到一些第三方库，最典型的就是GDAL。GDAL是由开源地理空间基金会发布的一个栅格和矢量的地理空间数据转换库，它会创建一个栅格或矢量的抽象数据模型，通过这个抽象模型，可以对栅格或矢量数据进行读取。 1、
基于python与GDAL获取遥感影像
wangzhaoqi_818的博客
11-30 1724
data = band.ReadAsArray( xoff, yoff, xsize, ysize, t_xsize, t_ysize ) //获取影像数据。band_type = fh.GetRasterBand(1).DataType //数据类型。geotransform = fh.GetGeoTransform() //文件几何信息。ysize = ysize = fh.RasterYSize //文件列数。band = fh.GetRasterBand(1) //获取波段。
python 遥感影像校正
最新发布
09-02
遥感影像校正是指对遥感影像进行几何、辐射和大气校正，以提高影像质量和准确性。校正包括获取几何参数、处理错误提示和计算高程。 1. 获取几何参数：对于较大的影像，为了提高校正精度，可以采取影像分块计算模式。可以使用函数`getInfoFromFilename`从影像命名规范中获取几何参数，例如卫星ID、传感器ID和年份等。 2. 处理错误提示：如果在校正过程中出现`UnboundLocalError: local variable 'SR' referenced before assignment`错误提示，可以检查`getInfoFromFilename`函数是否能够处理输入的影像命名规范，并根据需要修改代码。 3. 计算高程：在遥感影像校正中，可以利用全球DEM数据计算影像地区的平均高程。首先可以编写一个函数`geo2imagexy`来将给定的投影或地理坐标转换为影像图上的行列号。然后可以使用函数`extent2xy`来计算影像地理范围对应的影像图上的行列号。 综上所述，遥感影像校正涉及获取几何参数、处理错误提示和计算高程等步骤。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [遥感+python 1.3 大气校正](https://blog.csdn.net/Harry233/article/details/131289833)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值