摘要： python代码 import os import numpy as np from osgeo import gdal import glob import datetime # 读图像文件 def read_img(filename): dataset = gdal.Open(filename) 阅读全文

摘要： 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 阅读全文

摘要： Landsat8 L1 T数据是辐射校正数据使用地面控制点和数字高程模型数据进行精确校正后的数据产品，还需要做辐射校正（辐射定标和大气校正）。 一、辐射定标 辐射亮度L=DN*Gain+Bias from osgeo import gdal from osgeo import gdal_array 阅读全文

摘要： 原文链接：https://blog.csdn.net/m0_38106923/article/details/101937827 一张图片可以储存为多种格式，为什么有的几十KB，有的几百MB，有的静止不动，有的是好几个画面循环播放？在项目开发的过程中经常会读取或保存图像文件，不同类型的图像特点不同， 阅读全文

摘要： 一、简介 在对地理信息数据处理时，常常会遇到对线进行平滑和简化的操作。线的平滑可以使用拟合或插值来完成。多段线（Polyline）简化算法可以帮助我们减少Polyline的点数，从而降低输入规模。对多段线简化算法，通常的做法是在一定的近似精度下，删除一些点或者边。 二、线的平滑 插值和拟合有所不同。 阅读全文

摘要： 一维插值 插值不同于拟合。插值函数经过样本点，拟合函数一般基于最小二乘法尽量靠近所有样本点穿过。常见插值方法有拉格朗日插值法、分段插值法、样条插值法。 拉格朗日插值多项式：当节点数n较大时，拉格朗日插值多项式的次数较高，可能出现不一致的收敛情况，而且计算复杂。随着样点增加，高次插值会带来误差的震动现 阅读全文

摘要： 在做图像数据处理时，经常会有栅格数据转矢量数据的操作，转换后的矢量文件会存在锯齿状边缘，不太美观，因此常常需要对矢量（shp）文件做平滑处理。 1 利用arcgis实现shp的平滑和简化 ArcToolbox / Cartography Tool / Generalization / Smooth 阅读全文

摘要： 一、简介 遥感图像融合的定义是通过将多光谱低分辨率的图像和高分辨率的全色波段进行融合从而得到信息量更丰富的遥感图像。常用的遥感图像融合方法有Brovey\PCA\Gram-Schmidt方法。其中Gram-Schmidt方法效果较好，且应用广泛。 多源图像融合可类比为多源信息融合，即通过获取源数据原 阅读全文

摘要： 原文链接：https://blog.csdn.net/lyxleft/article/details/90675666 一、简介 在图像分割领域，一个重要任务便是分割出感兴趣（ROI）区域。如果是简易的矩形ROI区域其实是非常容易分割的，opencv的官方python教程里也有教到最简易的矩形ROI 阅读全文

摘要： 一、简介 栅格数据进行重新投影比矢量数据更复杂，对于矢量，你只需要每个顶点的新坐标就可以轻松实现投影转换，但对于栅格，你需要处理像元发生形变和偏移的情况，以及从旧单元格位置到新单元格位置的一对一映射。新单元格位置不存在（如下图）。确定新单元格像元值的最简单方法是使用最接近输出单元格映射的输入单元格中 阅读全文