IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> Python知识库 -> 【python + GIS】提取遥感影像经纬度为二维数组 -> 正文阅读

[Python知识库]【python + GIS】提取遥感影像经纬度为二维数组

一、应用场景

最近在求非平坦地区的太阳高度角,其中有一个参数就是纬度,因此就想到了把纬度作为值来创建二维数组的想法。

二、参考文献

《读取遥感影像的信息》
《使用python读取tiff文件中的经纬度,并将数据以excel表的形式输出(详细步骤)》

三、代码展示

import os
import numpy as np
from osgeo import gdal
import glob
import time
import math
import copy


DEM_Path = "D:\\DATA\\SEBAL\\0. InitialData\\DEM1984.tif"  #DEM影像
dataset = gdal.Open(DEM_Path)  # 打开tif
# 获取行数列数和地理信息
# geo_information(0):左上像素左上角的x坐标。    
# geo_information(1):w - e像素分辨率 / 像素宽度。    
# geo_information(2):行旋转(通常为零)。
# geo_information(3):左上像素左上角的y坐标。
# geo_information(4):列旋转(通常为零)。
# geo_information(5):n - s像素分辨率 / 像素高度(北半球上图像为负值)   
width = dataset.RasterXSize  # 宽度
height = dataset.RasterYSize  # 高度
geo_information = dataset.GetGeoTransform()
proj = dataset.GetProjection() #地图投影
filevalue= np.array(dataset.ReadAsArray(0, 0, width, height), dtype=float)  # 将数据写成数组,对应栅格矩阵
LongitudeValue=np.zeros((height,width))  #创建0矩阵
LatitudeValue=np.zeros((height,width))
for x in range(height):  # 行     
   for y in range(width):  # 列
    # 赋值经纬度
       LongitudeValue[x][y]=geo_information[0] + y * geo_information[1] + x * geo_information[2]
       LatitudeValue[x][y]=geo_information[3] + y * geo_information[4] + x * geo_information[5]

最后得到的数组和图像分别是这样的:
经度:
在这里插入图片描述

纬度:
在这里插入图片描述

四、反思与注意事项

1、遥感影像的长度和宽度以及二维数组的行列要对应正确,比如for循环的时候,如果写反了会出现长度不够的错误
2、 创建0矩阵的时候也是要注意参数中的行列对应,如果不行就交换一下
3、这个算法的核心是掌握地理信息的六个参数:
geo_information(0):左上像素左上角的x坐标。
geo_information(1):w - e像素分辨率 / 像素宽度。
geo_information(2):行旋转(通常为零)。
geo_information(3):左上像素左上角的y坐标。
geo_information(4):列旋转(通常为零)。
geo_information(5):n - s像素分辨率 / 像素高度(北半球上图像为负值)

  Python知识库 最新文章
Python中String模块
【Python】 14-CVS文件操作
python的panda库读写文件
使用Nordic的nrf52840实现蓝牙DFU过程
【Python学习记录】numpy数组用法整理
Python学习笔记
python字符串和列表
python如何从txt文件中解析出有效的数据
Python编程从入门到实践自学/3.1-3.2
python变量
上一篇文章      下一篇文章      查看所有文章
加:2022-09-04 01:07:23  更:2022-09-04 01:08:50 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/15 10:13:18-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码