首页 > 编程语言 >Python遥感影像叠加分析:基于一景数据提取另一数据

Python遥感影像叠加分析:基于一景数据提取另一数据

时间:2024-06-17 14:10:32浏览次数:30  
标签:Python 一景 植被 li 栅格数据 vt lcc array 数据

  本文介绍基于PythonGDAL模块,实现基于一景栅格影像,对另一景栅格影像的像元数值加以叠加提取的方法。

  本文期望实现的需求为:现有一景表示6种不同植被类型.tif格式栅格数据,以及另一景与前述栅格数据同区域的、表示植被参数.tif格式栅格数据;我们希望基于前者中的植被类型数据,分别提取6种不同植被类型的植被参数数值。这里需要注意,两景栅格影像的行数、列数也都是一致的。

  了解了具体需求后,我们即可开始代码的实践;本文用到的具体代码如下所示。

# -*- coding: utf-8 -*-
"""
Created on Thu Dec  1 16:56:26 2022

@author: fkxxgis
"""

from osgeo import gdal

vt_file_path = "E:/LC_M/data/LC.tif"
lcc_file_path = "E:/LC_M/data/LC_Clip.tif"

vt_raster = gdal.Open(vt_file_path)
vt_array = vt_raster.ReadAsArray()
lcc_raster = gdal.Open(lcc_file_path)
lcc_array = lcc_raster.ReadAsArray()

raster_row, raster_col = vt_array.shape
li_1, li_2, li_3, li_4, li_5, li_6 = [ [] for i in range(6)]

for i in range(raster_row):
    for j in range(raster_col):
        if vt_array[i][j] == 1 and lcc_array[i][j] != 0:
            li_1.append(lcc_array[i][j])
        elif vt_array[i][j] == 2 and lcc_array[i][j] != 0:
            li_2.append(lcc_array[i][j])
        elif vt_array[i][j] == 3 and lcc_array[i][j] != 0:
            li_3.append(lcc_array[i][j])
        elif vt_array[i][j] == 4 and lcc_array[i][j] != 0:
            li_4.append(lcc_array[i][j])
        elif vt_array[i][j] == 5 and lcc_array[i][j] != 0:
            li_5.append(lcc_array[i][j])
        elif vt_array[i][j] == 6 and lcc_array[i][j] != 0:
            li_6.append(lcc_array[i][j])

  其中,vt_file_path为表示植被类型的栅格数据,lcc_file_path为表示植被参数的栅格数据。

  代码的整体思路其实也非常简单,首先通过gdal.Open()函数与.ReadAsArray()函数,分别读取两个栅格数据,并将两个栅格数据中的像元数值信息转换为数组格式;随后,因为表示不同植被类型.tif格式栅格数据共有6种不同的像元数值,因此我们通过[] for i in range(6)这句代码,批量创建6个空的列表,用于存放6种不同植被类型分别对应的植被参数数值;接下来,同时遍历两个栅格数据,并基于表示不同植被类型.tif格式栅格数据的像元数值,将表示植被参数.tif格式栅格数据的像元数值依次提取、放入不同的列表中。

  这里有一点需要注意,因为在表示植被参数.tif格式栅格数据中0为无效值,因此在提取时,加了一个是否为0的判断;这一点大家在实际应用时结合自己的需求加以修改即可。

  通过上述代码,我们即可将6种不同植被类型分别对应的植被参数数值提取出来,并存放于不同的列表中;随后即可基于不同列表中的数据加以各项空间分析。

标签:Python,一景,植被,li,栅格数据,vt,lcc,array,数据
From: https://www.cnblogs.com/fkxxgis/p/18252259

相关文章

  • 是否可以从存储闪存卡中恢复数据?如何避免存储卡数据丢失?
    天津鸿萌科贸发展有限公司致力于数据恢复业务二十余年,为政府机构、各类大中型企业、医疗教育部门及个人用户等提供了专业的数据恢复服务,避免了重大损失。存储卡无处不在:智能手机、数码相机、平板电脑和摄像机。它们体积小,方便,并以相对较低的成本提供大量的数据存储空间。但是,......
  • 大数据开发语言Scala入门
    大数据开发语言Scala入门指南在当今数据驱动的世界中,Scala语言因其在大数据领域的卓越表现而受到青睐,特别是在ApacheSpark等框架中。本文将为你提供一个Scala语言的入门指南,帮助你理解Scala的基本概念,并指导你如何开始使用Scala进行大数据开发。Scala简介Scala是一种多......
  • 盘点一个Python自动化办公的问题
    大家好,我是Python进阶者。一、前言前几天在Python铂金交流群【逆光】问了一个Python自动化办公的问题,问题如下:问题我现在有两个表a、b,for循环a、b,如果a的条件满足b,则把b的值赋给a,目前a有7万条数据,b有300条。我写的代码20分钟都没跑完。这是代码,请问改怎么解决?二、实现过......
  • python3.10.10安装
    链接:https://www.python.org/选择一个盘建个python文件夹(任意盘,以E盘 python310为例,文件名任意字母数字下划线);安装包可分享路径不要太深E:\python310卸载uninstall 卸载之后可以把之前存储位置的文件夹(E:\python310)删除......
  • Python项目实战:制作一个翻译软件
    大家好,我是你们的老朋友南枫,今天咱们来学一个好玩、有趣、又实用的项目——制作一个翻译软件。我们参考的是有道翻译,首先还是老规矩,咱们先导入所需要使用到的模块:导入进来之后,我们需要把该要的参数,全都给copy下来(有请求头、cookie、防盗链等):既然我们是要做一个翻译的软......
  • 绘制MySQL数据库的实体关系图(ERD)与逻辑模型图
    绘制MySQL的实体关系图(ERD,Entity-RelationshipDiagram)和数据库模型图对于理解和设计数据库结构非常重要。这些图表有助于可视化数据库中的表、列、关系以及约束。以下是一些步骤和工具,你可以使用它们来创建这些图表。步骤确定实体:首先,确定你的数据库中的实体(或称为表......
  • Python中的属性
        Python中的属性主要分为类属性,对象属性。1.类属性    类属性:类所有,所有的实例对象都能够共享,能通过类名和实力对象名访问,当当前的类属性被实例对象通过对象名.属性名的形式调用之后,当前对象就会多出一个实例属性,此后使用对象名.属性名的形式调用的就是对象属......
  • python-不定方程求解
    [题目描述]给定正整数 a,b,c。求不定方程ax+by=c 关于未知数 x 和y 的所有非负整数解组数。输入:一行,包含三个正整数 a,b,c,两个整数之间用单个空格隔开。每个数均不大于 1000。输出:一个整数,即不定方程的非负整数解组数。样例输入12318样例输出14来源/分类(难度系数:......
  • django 接入OIDC认证登录(django admin后台使用OIDC 或github账号登录) django 使用p
    参考文档模块文档:https://python-social-auth.readthedocs.io/en/latest/接入github账号登录参考:https://blog.csdn.net/yannanxiu/article/details/112622781;测试项目地址:https://github.com/AngelLiang/django-social-auth-demo/tree/main--创建githubOAuth应用:https:/......
  • MySQL使用命令导出数据库
    MySQL自带了一个命令行工具mysqldump,可以方便的将数据导出至文件使用案例mysqldump-uroot-p数据库名称>xxx.sql#再输入密码即可#增加过滤条件mysqldump-uroot-p数据库名称表名称--where="id>1000">xxx.sql过滤条件#仅导出数据而不包含表结构--no-crea......