如果你的gdal版本是3.2.2及以上,那么只能:
from osgeo import gdal
ds = gdal.Open(data)
rows = ds.RasterYSize
cols = ds.RasterXSize
bandnum = ds.RasterCount
transform = ds.GetGeoTransform()
ds是一个对象,rows是影像的行,也是Y轴长度。对应的,cols则是影像的列,X轴的长度。bandnum代表影像的波段数。
transform是一个list,存储着栅格数据集的地理坐标信息。
#transform[0] /* top left x 左上角x坐标(经度)*/
#transform[1] /* w--e pixel resolution 东西方向上的像素分辨率*/
#transform[2] /* rotation, 0 if image is "north up" 如果北边朝上,地图的旋转角度*/
#transform[3] /* top left y 左上角y坐标(纬度)*/
#transform[4] ...
rasterio是什么?GDAL (Geospatial Data Abstraction Library) 是地理空间数据处理领域的基石。它是一个用 C++ 编写的开源库
然而,GDAL 的强大也带来了它的“缺点”。它的官方 Python 绑定 (osgeo.gdal) 是从 C++ 代码自动生成的,这导致了它的 API 对于 Python 开发者来说不够友好。
Rasterio 的诞生就是为了解决 GDAL Python 绑定的这些痛点。它把自己定位为“GDAL for Python users”,它不是要取代 GDAL,而是架设在 GDAL 强大的能力之上,提供了一个更加现代、简洁和“Pythonic”的接口。
重要的一点是:Rasterio并没有重新发明轮子,它的核心功能依然依赖于GDAL。
我以前是通过学习rasterio的源代码,进而学习行业内的高手是怎么封装GDAL完成不同的功能。
这个过程的步骤是:
1.安装pycharm
2.安装rasterio
3.打开在pycharm打开某个rasterio的代码。
比如下面的是rasterio封装了matplotlib、nump ...
如果您是地理空间数据处理领域的一员,那么 GDAL/OGR 库想必是您工具箱中不可或缺的瑞士军刀。多年来,我们习惯了使用一系列独立的命令行工具,如 gdalinfo、gdal_translate、gdalwarp 和 ogr2ogr。它们功能强大,但数量繁多,有时我们不得不去翻阅文档,才能记起那个特定任务到底该用哪个命令。
现在,这一切都将成为历史。随着 GDAL 3.11 的发布,一个全新的、革命性的功能登场了——统一命令行界面(Unified CLI)。这不仅仅是一次小修小补,而是一次彻底的重构,旨在为您提供一个更一致、更直观、更强大的 GDAL 体验。
核心理念:从分散到统一过去,每个工具都是一个独立的程序。而现在,所有的功能都被整合到了一个单一的入口点之下:gdal。
新的语法结构清晰明了,采用了分层设计:
gdal [领域] [操作] [选项] [参数]
这里的“领域”通常指代您正在处理的数据类型,如 raster(栅格)、vector(矢量)或 mdim(多维数据)。“操作”则是您想要执行的具体任务,如 info(获取信息)、convert(转换格式)或 reproject ...
GDAL(地理空间数据抽象库,Geospatial Data Abstraction Library)是一个开源库,用于读取、写入和处理多种栅格和矢量地理空间数据格式。它为地理信息系统(GIS)提供了强大而统一的数据访问方式。。简要介绍 GDAL 的架构、关键组件和功能。
GDAL核心架构GDAL 由多个相互连接的子系统组成,它们协同工作,提供全面的地理空间数据处理功能。其核心架构包括栅格 (GDAL) 和矢量 (OGR) 处理组件,以及用于可移植性、算法和虚拟文件系统的支持库。
关键组件GDAL 的架构可分为以下几个核心部分:
核心库(Core GDAL Library)
提供统一的数据模型和 I/O 接口
使用 CPL(Core Portability Library)实现平台兼容性
包含虚拟文件系统(Virtual File Interface)
2. 矢量与栅格支持
OGR(矢量库):支持 GeoJSON、Shapefile、PostGIS 等矢量格式
栅格模块:支持 GeoTIFF、NetCDF、HDF5、JPEG2000 等格式
3. 驱动框架(Driver ...
GDAL_OGR简单的使用教程检测安装from osgeo import gdal
查看版本gdal.VersionInfo('VERSION_NUM')
# '2040100'
开启python异常默认情况下,发生错误时,GDAL/OGR Python绑定不会引发异常。相反,它们返回错误值(例如None),并将错误消息写入sys.stdout。你可以通过调用UseExceptions()函数来开启异常:
from osgeo import gdal
# 开启异常
gdal.UseExceptions()
# 打开不存在的数据集
ds = gdal.Open('test.tif')
# 开启异常前
ERROR 4: test.tif: No such file or directory
# 开启异常后
RuntimeError Traceback (most recent call last)
<ipython-input-5-6ef000fdc647> in <module>
----> ...
各位,今天我们来解构一下“遥感”这个专业。
一、遥感到底是个啥?核心定义: 遥感(Remote Sensing),顾名思义,即“遥远地感知”。它并非超能力,而是一门利用飞机、卫星等平台上的传感器,在不直接接触的情况下,获取地球及其他天体信息的科学与技术。
可以把它想象成给地球做一次全面的“健康体检”。我们使用的工具不是听诊器,而是各种传感器,它们能捕捉到肉眼看不见的光谱信息,比如:
可见光影像: 就是我们常见的“地球写真集”。
红外影像: 可以感知地表温度,监测火灾或农作物长势。
雷达影像: 能穿透云层,全天候无阻碍地“看”地表。
一言以蔽之,遥感的核心流程就两步:获取数据(拍片)和处理分析(修图解读)。
二、听着高大上,学遥感究竟有何用?学习遥感,你将解锁一个独特的视角——“上帝视角”。
当你的同学还在为一张地图上的等高线而苦恼时,你已经可以在电脑上拖动鼠标,以上帝的视角俯瞰山川湖海、城市变迁。这种体验,堪比在玩一款名为“地球Online”的即时战略游戏。
地理相关学科的魅力在于培养人的宏观思维和跨学科整合能力。遥感人,通过分析卫星影像,为国土规划、环境保护、灾害预警、国家安全 ...
原文地址:
https://867810770.github.io/2018/11/09/%E6%B8%B8%E6%88%8F%E5%A4%96%E6%8C%82%E5%88%B6%E4%BD%9C%E5%8E%9F%E7%90%86/
游戏外挂开发原理 外挂现在大体上来讲分为三种,分别是模拟按键,WPE封包和内存挂。比如模拟键盘的,鼠标的,这种就叫做按键模拟,也叫脚本精灵;还有修改数据包的,这种就叫做WPE封包外挂,这种外挂相对而言比较稳定,但是对于技术要求难度也比较高,因为修改WPE封包,你需要知道和了解加密解密算法等等一系列的知识;还有最后一种就是修改本地内存的,这种相对而言比较常见,市场上面比较常见的也是这种游戏外挂,相对而言技术要求也不是太高,资料也比较齐全,比较大众;但好像没有修改服务器内存的哦,呵呵。其实修改服务器也是有办法的,只是技术太高一般人没有办法入手而已。这里就记一下本地修改内存的外挂制作原理。
其实,制作内存挂也不是很难,步骤就这么几步。第一步,找游戏数据内存地址,偏移地址;第二步,修改内存地址的值,达到外挂的目的;第三步,用编程语言写出一个程 ...
正如之前我说的,遥感影像在深度学习的输入形式一般是RGB格式,而遥感影像的原始数据一般由多个波段组合。
所以,我们需要做的一个步骤的:把多波段的遥感影像转换为RGB遥感影像。
RGB,三个字母说明了这个文件由三个矩阵组成,它们的顺序是RED、GREEN、BLUE。也叫做红绿蓝波段。
以国产遥感影像高分一号为例,GF1的原始数据的波段顺序是BLUE、GREEN、RED、NIR。即蓝、绿、红、近红外。
所以把多波段的遥感影像转换为RGB遥感影像的原理很简单,我们只需要保留前三个波段,并且把1、3波段调换顺序即可完成此步骤。
具体的python代码如下:
from osgeo import gdal
import os
def convert_4band_to_rgb(input_tif, output_tif):
"""
Convert a 4-band (RGB+NIR) TIFF to a 3-band RGB TIFF.
Args:
input_tif (str): Path to the input 4-band TIFF file. ...
Bilibili 下载指南本指南提供使用 yt-dlp 下载 Bilibili 视频(ID:BV1n2j2zxEZk)的 ai-zh 字幕(SRT 格式)和 M4A 音频的命令,分别适用于单个视频部分和整个合集。
下载单个视频部分(字幕和音频)要下载某个特定部分(例如第 4 部分)的 ai-zh 字幕和 M4A 音频:
yt-dlp --cookies www.bilibili.com_cookies.txt --write-subs --sub-langs ai-zh -f "bestaudio[ext=m4a]" --restrict-filenames -o "%(title)s_p%(part_number)s.%(ext)s" https://www.bilibili.com/video/BV1n2j2zxEZk?p=4
注意事项:
将 ?p=4 替换为所需的视频部分编号(例如,?p=5 表示第 5 部分)。
确保 www.bilibili.com_cookies.txt 是从浏览器导出的有效 cookies 文件。
输出文件将命名为类似 有声书长生从娶妻开始_p04.m4 ...
landsat9影像下载的渠道有很多个,最常见的是USGS、GEE。
今天是演示怎么在USGS下载Landsat9影像。
下载数据本身是免费的、不收钱。但是最大的阻碍就是科学网上冲浪。
所以也导致了很多小白,跑去国内二手交易平台寻求帮助。
这些帮助不仅限于注册USGS,卖USGS成品号,帮下载数据。
如果你掌握了landsat9的下载全流程,你可以去赚这个辛苦的、繁琐的零花钱。
1. 注册USGS(⊙o⊙)…我在几年前注册过一个,现在为了演示,开始注册第二个USGS账号。
邮箱我使用的是网易的126邮箱。
邮箱是新注册的邮箱,13:55注册邮箱,然后开始在USGS注册,4分钟完成注册USGS。
即,13:59收到USGS注册确认邮箱,4分钟完成USGS账号注册,如下
登录账号地址是:https://ers.cr.usgs.gov/login
成功注册:
不想注册USGS,用我的:账户名:rs4xiaobai12密码:rs4xiaobai12
2. landsat9 数据下载数据下载的页面是https://earthexplorer.usgs.gov/
2.1 设置搜索条件首先要 ...