gdal获取像元的灰度值java_GDAL利用地理坐标读取图像像元值
最近的一个项目需要在电子海图中下载已知水深点,导出点的地理坐标(经纬度)。然后在arcgis中打开这些地理坐标输出为shp,利用GDAL读取不同波段的点对应的像元值,从而构建水深和像元值的对应关系。
其实想要根据经纬度得到像元值有两种方法。
1.通过arcpy选择点输出像元值,这个方法简单快捷,依赖于arcgis。
2.通过GDAL读取经纬度(度分秒),转换成小数形式。使用GDAL的CoordinateTransformation方法将经纬度转换成投影坐标,再转换成行列坐标,最后读取像元值。这里给出官方文档里面的代码以供参考。
补充一点:这里之所以把经纬度先转成投影坐标是因为遥感影像的坐标系是投影坐标系,而仿射六参的转换方法要求待转换点坐标和影像的投影坐标系一致,所以需要先转换成投影坐标。
如果遥感影像的坐标系是地理坐标系(即没有投影坐标系),那么这里直接对经纬度坐标进行仿射六参转换即可得到行列号。
而按下面的代码此时srs为none,srs与geocs的转换器仍然是对应的,那么经纬度转换后仍然是对应的格式,经过仿射变换最终得到行列号,也同样成立(该步骤可省略)。
其中GetGeoTransform获取获取图像的地理空间范围、分辨率信息。第1、4个值是左上角第一个像元中心的坐标。2、6是xy方向的空间分辨率。3、5是旋转系数和平移系数。
SpatialReference()建立空间参考,GetProjection返回投影信息,以wkt格式输出,最后用ImportFromWkt方法读取字符串并创建地理坐标系、基准面、投影方法、分辨率等。
CloneGeoCS方法创建地理坐标(经纬度),CoordinateTransformation方法建立投影坐标系和地理坐标系的转换关系,TranformPoint方法将地理坐标带入转换关系式,得到投影坐标,最后再计算得到行列数。
Tips
1.list out of range 应该是输入的列表格式有问题。
2.此外注意经度和纬度的输入位置。
gdal获取像元的灰度值java_GDAL利用地理坐标读取图像像元值相关推荐
- Python使用openCV把原始彩色图像转化为灰度图、使用OpenCV把图像二值化(仅仅包含黑色和白色的简化版本)、基于自适应阈值预处理(adaptive thresholding)方法
Python使用openCV把原始彩色图像转化为灰度图.使用OpenCV把图像二值化(仅仅包含黑色和白色的简化版本).基于自适应阈值预处理(adaptive thresholding)方法 目录
- matlab 读取照片imread,利用matlab读取图像
怎样用matlab读取20张图片并依次展示出来 指定路径下 单个文件夹data中所有图像 P = '.\data\';% 图像文件夹路径 img_path_list = dir(strcat(P,'* ...
- 怎么用matlab读一张图像,利用matlab读取图像
怎样用matlab读取20张图片并依次展示出来 指定路径下 单个文件夹data中所有图像 P = '.\data\';% 图像文件夹路径 img_path_list = dir(strcat(P,'* ...
- 图像中值滤波python代码_图像中值滤波FPGA实现
C语言实用数字图像处理.pdf6.34 MB05-11-13|19:30 FPGA实验报告-李炎东.doc633.66 kB16-01-14|10:28 中值滤波在红外成像引信中的应用及硬件实现.ca ...
- matlab 二值化_MATLAB实验,图像二值化处理
防止迷路,记得关注点赞哦.下期见~~ 1.读入一幅彩色图像,变换为灰度图像和二值图像,并在同一个窗口内分成三个子窗口来分别显示RGB图像和灰度图像,注上文字标题. I=imread('001.jpg' ...
- VC++读取图像RGB值
代码: #include <iostream> #include <fstream> #include <string> #include <windows. ...
- python opencv读取图像像素值_python-opencv--图像像素通道读取及修改
data/dtype/size/shape/len ''' import cv2 as cv import numpy as np def access_pixes(image): print(ima ...
- Python-OpenCV 处理图像(六)(七)(八):对象识别 图像灰度化处理 图像二值化处理
为了加快处理速度,在图像处理算法中,往往需要把彩色图像转换为灰度图像. 0x00. 灰度图 灰度数字图像是每个像素只有一个采样颜色的图像,这类图像通常显示为从最暗黑色到最亮的白色的灰度. 灰度图像与黑 ...
- 自适应阈值图像二值化
一.二值化 关于二值化的介绍,以前的博客中有介绍,这里就不再描述了,二值化介绍:二值化分为固定阈值二值化和自适应阈值二值化,固定阈值二值化方式是我们常用的二值化方式,需要自己摸索一个经验阈值,不断调整 ...
最新文章
- excel调用python编程-如何在excel中调用python脚本
- 【问链财经-区块链基础知识系列】 第十九课 区块链应用于电子存证
- 第三章:系统困境之 忽略过去现在未来构成的时间系统
- 初探Golang(2)-常量和命名规范
- linux系统下压缩文件,Linux系统下文件的压缩.打包与解压
- MySQL常用命令大全
- 可视化实战!Python+BI,手把手教你做炫酷的NBA可视化分析
- 前端综合学习笔记---异步、ES6/7、Module、Promise同步 vs 异步
- HTTP、websocket、XMPP、COAP、MQTT和DDS协议对比
- linux子系统安装docker,linux上安装Docker(非常简单的安装方法)
- wsimport的使用
- 微信转发软件后缀_简单修改后缀让微信发送25M以上微信大视频
- 《缠中说禅108课》37:背驰的再分辨
- 【VS Code】Windows10下VS Code配置Ruby语言环境
- activeMQ安装问题wrapper | Launching a JVM... wrapper | Unable to execute Java command.
- windows系统IP地址、localhost、127.0.0.1 、0.0.0.0和 本机IP区别
- 麦田守望者--走出软件作坊:三五个人十来条枪 如何成为开发正规军(四十三)...
- 通用人工智能:我们还有多远?
- JAVA 操作 properties 配置文件
- window10 男装linux,男人帮男装购物app-男人帮手机版下载v2.3.5-Linux公社
热门文章
- 九、华为鸿蒙HarmonyOS应用开发之Java UI框架、常用Image组件使用
- matlab dpsk,2DPSK调制与解调matlab(最新整理)
- 思科isis路由的优先级_华为 路由双点双向引入
- java 反射泛型方法_java基础之反射和泛型以及注解
- python怎么读取石墨表格_每2秒进行一次statsd xaxis石墨统计 - python
- android点击地址调用地图,Android 实现点击按钮 调用手机外部地图导航
- PP视频如何不允许WiFi下自动播放
- win7系统找不到我的电脑图标的应对办法
- win7系统每次开机都需要疑难解答的原因与解决方法
- 手把手教你dns服务器未响应导致无法上网怎么办