gdal库读取tif影像坐标
# coding=utf-8
from osgeo import gdal
from gdalconst import *def readImage(img_path):data = []# 以只读方式打开遥感影像dataset = gdal.Open(img_path, GA_ReadOnly)if dataset is None:print("Unable to open image file.")return dataelse:print("Open image file success.")geoTransform = dataset.GetGeoTransform()print(geoTransform)im_proj = dataset.GetProjection() # 获取投影信息print (im_proj)bands_num = dataset.RasterCountprint("Image height:" + dataset.RasterYSize.__str__() + " Image width:" + dataset.RasterXSize.__str__())print(bands_num.__str__() + " bands in total.")for i in range(bands_num):# 获取影像的第i+1个波段band_i = dataset.GetRasterBand(i + 1)# 读取第i+1个波段数据band_data = band_i.ReadAsArray(0, 0, band_i.XSize, band_i.YSize)data.append(band_data)print("band " + (i + 1).__str__() + " read success.")return data
print:
Open image file success.
(495341.48890999705, 0.10028999999997644, 0.0, 4243205.078130003, 0.0, -0.10028999999996153)
PROJCS["WGS_1984_Transverse_Mercator",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",114],PARAMETER["scale_factor",1],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]
Image height:4772 Image width:4832
3 bands in total.
band 1 read success.
band 2 read success.
band 3 read success.
输出的六个参数
//adfGeoTransform[0] 左上角x坐标
//adfGeoTransform[1] 东西方向分辨率
//adfGeoTransform[2] 旋转角度, 0表示图像 “北方朝上”
//adfGeoTransform[3] 左上角y坐标
//adfGeoTransform[4] 旋转角度, 0表示图像 “北方朝上”
//adfGeoTransform[5] 南北方向分辨率
gdal库读取tif影像坐标相关推荐
- Python gdal库读取tif文件
from osgeo import gdal # GDAL库主要提供对栅格数据的处理,使用抽象数据模型来解析所支持的数据格式 import filename_cut as fc import matp ...
- c++ 读取 .tfw 文件数据(读取 .tif 影像中的投影信息)
关于读取 .tif 影像的投影信息,下面介绍了两种方法: 1.按读取 .txt 的方法读取 .tfw 文件: 2.使用 GDAL 直接从 .tif 影像 ...
- GDAL库——读取图像并提取基本信息
GDAL库是一个跨平台的栅格地理数据格式库,包括读取.写入.转换.处理各种栅格数据格式(有些特定的格式对一些操作如写入等不支持).它使用了一个单一的抽象数据模型就支持了大多数的栅格数据.这里有GDAL ...
- GDAL库读取Envisat ASAR数据
GDAL库本身就可以读取Envisat的图像数据,具体链接为:http://www.gdal.org/frmt_various.html#Envisat. 但是对于ASAR传感器的数据来说,GDAL在 ...
- python+spyder读取tif影像
第一次遥感作业,老师就让我们用python对遥感图像进行数据处理. import cv2 import numpy as np#读取数据,显示图像 name = 'RGB.tif' img = cv2 ...
- GDAL C#读取shp中文属性值乱码问题
GDAL的C#版本读取shp中,如果属性值中含有中文,读出来有可能是乱码的问题,根据SWIG生成的C#代码调试发现问题所在,在Ogr.cs文件中有这么一个函数,代码如下: internal stati ...
- Linux配置gdal变量环境,基于GDAL库的LINUX环境下的FY-3D和FY-4A遥感数据处理
王燕婷 摘要 基于GDAL库的功能和其开源.可对栅格空间数据操作的特点,将GDAL与C++语言的完美结合,可极大得发挥其在开发遥感图像处理软件中的优势,并可跨平台应用于不同操作系统下的遥感图像数据处理 ...
- python 读取geotiff_科学网—利用python GDAL库读写geotiff格式的遥感影像方法 - 张伟的博文...
(1)利用python GDAL库读写geotiff格式的遥感影像方法,具有很好的参考价值,不错! from osgeo import gdal import numpy as np def read ...
- python读取tiff影像_科学网—利用python GDAL库读写geotiff格式的遥感影像方法 - 张伟的博文...
(1)利用python GDAL库读写geotiff格式的遥感影像方法,具有很好的参考价值,不错! from osgeo import gdal import numpy as np def read ...
最新文章
- TinyOS下TOSSIM仿真
- 有一部分文章设置了私密,论文完成后会公开~
- 有必要服务器虚拟化吗,服务器虚拟化有必要吗
- homebrew长时间停在Updating Homebrew 这个步骤
- java shareable_spring中@Resource和@Autowired理解详解_编程语言_IT虾米网
- java添加锁_java – 如何在这种情况下添加锁?
- asynchttpclient学习
- Jquery中用offset().top和offsetTop的比较
- java se win10_Win10 JAVASE的下载和环境变量设置
- Python爬虫代码示例
- matlab 图像分割并填充颜色
- ubuntu安装anjuta
- Grasshopper 二次开发 (C#) Part 1 - Introductions to the C# Coding in Grasshopper
- 生鲜电商的“7年之谎”
- 简单搜索引擎使用技巧
- 自建网站教程!如何用云服务器搭建个人网站?
- 代理ARP(Proxy Arp)
- html滚动字幕如何调色,如何通过制作滚动字幕的软件实现这种片尾的向上滚动字幕效果...
- 对png空白部分取消button事件教程
- LeetCode1-400题汇总,希望对你有点帮助!
热门文章
- mysql配置方案_MySQL的安装与配置
- elementui :on-remove怎么用_绍兴【办公室绿化】怎么操作
- Javascript第五章获取DOM对象的属性,加游览器搜索框内容点击跳转源码第七课
- 配置 sql server linux,配置SQL Server on Linux(2)
- 语音识别代码_Povey正式出任小米语音首席科学家,小米移动端框架MACE全面支持Kaldi...
- 进程、轻量级进程(LWP)、线程
- C# winform 跨线程修改界面
- 170308、oracle查看被锁的表和解锁
- 【MFC】将当前的日期转化为1970年开始的秒计数
- codeforces 431 B Shower Line【暴力】