计算仿射变换六参数(Python)
摘要：根据已知点对求得坐标转换的参数是一个值得研究的问题。这里用到的编程技巧不多，关键是要用到线性代数和数值分析的知识。纵观当前地图坐标保密处理或者坐标系转换的实例，其无外乎采用旋转、平移、拉伸等方式，于是数值的计算无外乎于解n个n元一次方程组，最后通过误差分析进行拟合。 下面就是一个形如 x' = Ax
阅读全文
摘要：原文链接：https://blog.csdn.net/qq_27045589/article/details/81062586 一、几何校正方法 图像校正本质是建立一种从原始图像行列号到某种投影的数学关系，即实现图像行列坐标到投影坐标的转换。不同的校正方法利用了不同的方法来表示转换关系，但本质上式相
阅读全文
摘要：1 自定义镶嵌函数 遥感图像的镶嵌，主要分为5大步骤：step1： 1）对于每一幅图像，计算其行与列；2）获取左上角X,Y3）获取像素宽和像素高4）计算max X 和 min Y，切记像素高是负值 maxX1 = minX1 + (cols1 * pixelWidth) minY1 = maxY1
阅读全文
摘要：在同级目录完成tif和jpg的批量转换 import os import cv2 import numpy as np from osgeo import gdal #数据格式转化 def normalization(data): _range = np.max(data) - np.min(dat
阅读全文
摘要：SLIC超像素分割分割后对象掩膜提取 from skimage.segmentation import slic from skimage.segmentation import mark_boundaries from skimage.util import img_as_float import
阅读全文
摘要：本节将介绍如何在Python中用GDAL实现根据矢量边界裁剪栅格数据。 from osgeo import gdal, gdal_array import shapefile import numpy as np import os #批量shp裁剪tiff影像 try: import Image
阅读全文
摘要：在做遥感图像处理时，矢量数据与栅格数据的转换是常有的事，本节l将利用GDAL、ogr来完成矢量与栅格之间的转换. 1.矢量转栅格 import gdal import ogr import gdalconst #from osgeo import ogr, osr, gdal,gdalconst d
阅读全文
摘要：本节主要利用python 和 GDAL来计算一些常用的遥感指数，如，NDVI（归一化植被指数）、RVI（比值植被指数）、NDWI （归一化水体指数）、SAVI（土壤调节植被指数）。 1.读取波段数据 import numpy as np from osgeo import gdal in_ds =
阅读全文
摘要：上一节简单介绍了GDAL，这一节将介绍一些GDAL的基本操作，如影像读写、波段提取、波段合成等。代码均用python编写。 1.遥感影像的读写 众所周知，遥感影像是以栅格形式存储的，GDAL中使用dataset表示一个栅格数据（使用抽象类GDALDataset表示），一个dataset包含了对于栅格
阅读全文
摘要：1.GDAL简介 GDAL(Geospatial Data Abstraction Library)是一个用于栅格数据操作的库，是开源地理空间基金会（Open Source Geospatial Foundation，OSGeo）的一个项目。 GDAL是一个操作各种栅格地理数据格式的库。包括读取、写
阅读全文
