遥感数据
浅野川茶
CV & RS
展开
-
基于Arcpy的批量二值图转矢量文件并选取特定值
最近导师要求实现一个批量的二值图转矢量文件,发现用python调用arcgis的arcpy库非常方便,记录一下。我的数据都是0,1分布的二值图,实际使用看需求来修改吧,我要提取值为1的区域,所以代码中gridcode值为1。# -*- coding: UTF-8 -*-#导入arcgis的pythonIDE使用或使用arcgis自带的python.exeimport arcpyfrom arcpy import envimport osimport timeenv.workspac原创 2020-12-29 15:00:06 · 589 阅读 · 0 评论 -
Python+OpenCV对灰度图像进行伪彩色处理
伪彩色处理:根据一定准则给灰度值赋予彩色值的处理。宏观来说就是将黑白图像转化为彩色图像,或者是将单色图像变换成给定彩色分布的图像。博主在学习遥感图像处理时接触到NDVI图像的计算,对生成的灰度图像需要进行伪彩色处理以增强对比。代码如下:(python+OpenCV)import cv2im_gray = cv2.imread("H:/gdal/rs-data/result/ndvi.jpg", cv2.IMREAD_GRAYSCALE)im_color = cv2.applyColorM原创 2020-07-12 18:07:45 · 8522 阅读 · 7 评论 -
Python批量计算遥感图像NDVI(归一化植被指数)
归一化植被指数:遥感影像中,近红外波段的反射值与红光波段的反射值之差比上两者之和。计算公式:NDVI = (NIR-R)/(NIR+R)注:NIR为近红外波段的反射值 R为红光波段的反射值。由于学习需要,从网上寻找方法试图批量计算该指数,防止以后忘记。附代码如下:(主要使用gdal库)import osfrom PIL import Imageimport numpy as npfrom osgeo import gdalimport globimport ...原创 2020-07-11 13:47:39 · 10015 阅读 · 5 评论 -
Python+GDAL实现批量进行遥感图像的拼接
主要的方法就是GDAL作者提供的gdal_merge.py,通过python调用控制台实现。#如果图片大于4G,需要在输入图片路径前加入 -co BIGTIFF=YESBIGTIFF=YESimport subprocessimport oscreate_slope = '''D:/python/python.exe D:/python/Scripts/gdal_merge.py -of GTiff -o '''list=[]fileDir = 'E:/湖北影像/武汉/2018'file原创 2020-09-04 13:22:41 · 5724 阅读 · 35 评论