随笔分类 -  GDAL篇(遥感图像处理)

• 根据两点经纬度坐标计算距离及方位角
  摘要：# 根据经度、纬度计算两地距离： def get_distance2(lat1, lon1, lat2, lon2): """获取地理坐标系下的两点间距离""" # GetDistanceInGeographyCoordinate, return two point distance radius_ 阅读全文
  posted @ 2023-04-25 20:28 行走的蓑衣客 阅读(481) 评论(0) 推荐(0) 编辑
• 多线程批量解压
  摘要：import threading import time import os import tarfile input_path = r"D:\jieyaqian" out = r'D:\jieyahou' classs = os.listdir(input_path) def repress(fo 阅读全文
  posted @ 2023-04-25 17:27 行走的蓑衣客 阅读(159) 评论(0) 推荐(0) 编辑
• 遥感影像去除nodata值
  摘要：import cv2 import scipy.interpolate import numpy as np from osgeo import gdal def read_img(filename): dataset = gdal.Open(filename) im_width = dataset 阅读全文
  posted @ 2023-04-25 17:22 行走的蓑衣客 阅读(195) 评论(0) 推荐(0) 编辑
• 遥感影像自适应降位拉伸(16_8)
  摘要：import os import sys import glob from osgeo import gdal import numpy as np import cv2 def CalHistogram(img): img_dtype = img.dtype img_hist = img.resh 阅读全文
  posted @ 2022-11-21 19:28 行走的蓑衣客 阅读(207) 评论(0) 推荐(0) 编辑
• 将16位遥感图像压缩至8位，并保持色彩一致
  摘要：""" 将16位遥感图像压缩至8位，并保持色彩一致 """ from osgeo import gdal import os import glob import numpy as np import matplotlib.pyplot as plt # plt 用于显示图片 import matp 阅读全文
  posted @ 2022-11-21 19:22 行走的蓑衣客 阅读(255) 评论(0) 推荐(0) 编辑
• 基于栅格单元与栅格单元的分区统计
  摘要：输入的区域可以是栅格的或者是矢量的，值数据只能是栅格的 import numpy as np import pandas as pd def zonal_statistic(zonal_arr, data_arr, stats_field=['mean'], nodata=255): """ Par 阅读全文
  posted @ 2022-07-23 15:11 行走的蓑衣客 阅读(192) 评论(0) 推荐(0) 编辑
• python 消除二值图像中面积小于某个阈值的连通域(消除孤立点)
  摘要：from skimage import measure import cv2 import numpy as np ## ##image:二值图像 ##threshold_point:符合面积条件大小的阈值 def remove_small_points(image,threshold_point) 阅读全文
  posted @ 2022-05-17 20:59 行走的蓑衣客 阅读(430) 评论(0) 推荐(0) 编辑
• Python 计算八邻域之和及图像连通区域标记
  摘要：#Python计算矩阵8邻域之和 Neighbors = [(1, 1), (1, -1), (1, 0), (-1, 0), (-1, 1), (-1, -1), (0, 1), (0, -1)] def get_neighbors_sum(array_of_arrays, row, col): 阅读全文
  posted @ 2022-05-17 20:58 行走的蓑衣客 阅读(873) 评论(0) 推荐(0) 编辑
• 遥感图像去除NoData及增强处理（笔记）
  摘要：1.去除NoData import cv2 import gdal import scipy.interpolate import numpy as np def read_img(filename): dataset=gdal.Open(filename) im_width = dataset.R 阅读全文
  posted @ 2022-05-17 20:27 行走的蓑衣客 阅读(472) 评论(0) 推荐(0) 编辑
• Python+GDAL 图像坐标、地理坐标、投影坐标的相互转化
  摘要：from osgeo import gdal from osgeo import osr import numpy as np def getSRSPair(dataset): ''' 获得给定数据的投影参考系和地理参考系 :param dataset: GDAL地理数据 :return: 投影参考 阅读全文
  posted @ 2022-04-10 12:53 行走的蓑衣客 阅读(738) 评论(0) 推荐(0) 编辑
• 图像拼接，单应性矩阵计算
  摘要：PCV库链接：https://github.com/Li-Shu14/PCV 将下载的文件解压打开cmd，执行如下指令：（1）执行cd命令，转到你所解压到的PCV的文件夹中。（2）输入python setup.py install # RANSAC算法拼接2张图片 from pylab import 阅读全文
  posted @ 2022-03-17 22:07 行走的蓑衣客 阅读(452) 评论(0) 推荐(0) 编辑
• python在图像上画线
  摘要：python代码 import math import os import sympy as sy import numpy as np import cv2 from osgeo import gdal def Polar2Cartesian(theata): """ 极坐标转直角坐标系 :ret 阅读全文
  posted @ 2021-09-12 22:00 行走的蓑衣客 阅读(1171) 评论(0) 推荐(0) 编辑
• 图像分块和还原
  摘要：在遥感图像处理时,通常因为图像太大,导致计算机内存不够,无法处理.将图像进行分块处理后,再对每一块进行处理将结果进行合并,既能减少计算机内存的负担,又能提高处理速度. python代码 from osgeo import gdal import os import numpy as np def w 阅读全文
  posted @ 2021-09-12 21:56 行走的蓑衣客 阅读(603) 评论(0) 推荐(0) 编辑
• sentinel2波段合成
  摘要：此代码是从sentinel2原始文件包中，选出B、G、R、NIR进行波动合成，保存成tif. python代码 # -*- coding: utf-8 -*- from osgeo import gdal import os import numpy as np from osgeo import 阅读全文
  posted @ 2021-09-12 21:42 行走的蓑衣客 阅读(761) 评论(0) 推荐(0) 编辑
• 遥感图像归一化
  摘要：python代码 import os import numpy as np from osgeo import gdal import glob import datetime # 读图像文件 def read_img(filename): dataset = gdal.Open(filename) 阅读全文
  posted @ 2021-09-12 21:33 行走的蓑衣客 阅读(888) 评论(0) 推荐(0) 编辑
• 批量解压landsat8/sentinel2原始文件包
  摘要：import os import tarfile import zipfile def un_landsat8_file(path,out): classs1 = os.listdir(path) for idx, folder in enumerate(classs1): if folder[-7 阅读全文
  posted @ 2021-09-12 21:10 行走的蓑衣客 阅读(157) 评论(0) 推荐(0) 编辑
• python对landsat8数据进行辐射校正和大气校正
  摘要：Landsat8 L1 T数据是辐射校正数据使用地面控制点和数字高程模型数据进行精确校正后的数据产品，还需要做辐射校正（辐射定标和大气校正）。 一、辐射定标 辐射亮度L=DN*Gain+Bias from osgeo import gdal from osgeo import gdal_array 阅读全文
  posted @ 2021-09-11 12:38 行走的蓑衣客 阅读(4514) 评论(4) 推荐(0) 编辑
• 遥感图像融合
  摘要：一、简介 遥感图像融合的定义是通过将多光谱低分辨率的图像和高分辨率的全色波段进行融合从而得到信息量更丰富的遥感图像。常用的遥感图像融合方法有Brovey\PCA\Gram-Schmidt方法。其中Gram-Schmidt方法效果较好，且应用广泛。 多源图像融合可类比为多源信息融合，即通过获取源数据原 阅读全文
  posted @ 2021-09-09 21:09 行走的蓑衣客 阅读(8054) 评论(0) 推荐(1) 编辑
• 遥感影像重投影
  摘要：一、简介 栅格数据进行重新投影比矢量数据更复杂，对于矢量，你只需要每个顶点的新坐标就可以轻松实现投影转换，但对于栅格，你需要处理像元发生形变和偏移的情况，以及从旧单元格位置到新单元格位置的一对一映射。新单元格位置不存在（如下图）。确定新单元格像元值的最简单方法是使用最接近输出单元格映射的输入单元格中 阅读全文
  posted @ 2021-09-08 21:53 行走的蓑衣客 阅读(1423) 评论(0) 推荐(0) 编辑
• 遥感影像重采样
  摘要：一、简介 图像重采样就是从高分辨率遥感影像中提取出低分辨率影像，或者从低分辨率影像中提取高分辨率影像的过程。常用的方法有最邻近内插法、双线性内插法、三次卷积法等 二、重采样方法 1 使用ReadAsArray函数 def ReadAsArray(self, xoff=0, yoff=0, win_x 阅读全文
  posted @ 2021-09-08 20:37 行走的蓑衣客 阅读(2885) 评论(0) 推荐(3) 编辑