Fork me on GitHub

随笔分类 -  GDAL

GDAL几何校正
摘要：GDAL几何校正之Geoloc校正 1. 几何校正原理 ​ 常见的几何校正方法有几何多项式校正、有理函数模型校正、局部区域校正模型、地理查找表校正。GDAL库中可以实现的校正方法有：几何多项式校正、RPC校正（有理函数模型）、TPS(薄板函数模型)校正、Geoloc校正。 2. Geoloc校正（地 阅读全文
posted @ 2020-12-20 23:48 Rser_ljw 阅读(2314) 评论(1) 推荐(0)
使用ogr裁剪矢量数据
摘要：使用ogr裁剪矢量数据 由来： ​ 近期有个需求，内容是这样的：我们有两个矢量数据，现在要求以一个矢量文件为底板，按字段对另一个矢量文件进行分割，生成若干小的shpfile文件 分析： ​ 经过分析之后，将步骤拆解如下： 首先确保两个shpfile投影坐标系统一 ​ 如果出现不统一的情况，那么用Ar 阅读全文
posted @ 2020-12-01 23:57 Rser_ljw 阅读(883) 评论(0) 推荐(0)
gdal笔记--chapter12
摘要：遥感图像分类（监督分类、非监督分类） 1.K均值聚类分析（非监督分类） #非监督分类 #堆叠栅格波段的函数 from osgeo import gdal import numpy as np def stack_bands(filenames): bands = [] for fn in filen 阅读全文
posted @ 2020-04-09 18:45 Rser_ljw 阅读(591) 评论(0) 推荐(0)
gdal笔记--chapter7
摘要：使用GR进行矢量分析 1.计算相交区域面积的两种方法（对要素或图层操作） from osgeo import gdal, ogr water_ds = ogr.Open(r'') #打开矢量数据，沼泽数据 water_lyr = water_ds.GetLayer(0) #打开图层 water_ly 阅读全文
posted @ 2020-04-08 22:15 Rser_ljw 阅读(1205) 评论(0) 推荐(0)
GDAL笔记--chapter11
摘要：本章主要介绍了利用numpy和scipy库做地图运算，并讲述了局部分析、焦点分析、区域分析和全局分析几种不同的方法。 1.保存为新栅格 import gdal import numpy as np #保存为新栅格的函数 def make_raster(in_ds, fn, data, data_ty 阅读全文
posted @ 2020-04-06 19:25 Rser_ljw 阅读(1189) 评论(0) 推荐(0)
GDAL投影转换/坐标偏移笔记
摘要：1.OSR # OSR（矢量数据投影） #作用：投影坐标系之间转换、地理坐标和投影坐标之间转换 #可用于几何对象和点（点属于几何对象） from osgeo import gdal import osr peters_sr = osr.SpatialReference() peters_sr.Imp 阅读全文
posted @ 2020-04-01 22:01 Rser_ljw 阅读(6024) 评论(1) 推荐(2)
GDAL笔记--chapter10
摘要：1.为栅格数据添加地面控制点 import shutil from osgeo import gdal orig_fn = r'' shutil.copy(orig_fn, fn) #因为要更新，所以需要对文件做个备份 ds = gdal.Open(fn, gdal.GA_Update) sr = 阅读全文
posted @ 2020-04-01 21:51 Rser_ljw 阅读(1473) 评论(0) 推荐(1)
GDAL笔记--chapter9
摘要：1.将独立的栅格波段合成为一张图像 #将独立的栅格波段合成一张图像 import os import numpy as np from osgeo import gdal os.chdir(r'E:\桌面文件保存路径\gdal\osgeopy-data\osgeopy-data\Washington 阅读全文
posted @ 2020-03-29 23:42 Rser_ljw 阅读(1084) 评论(0) 推荐(0)
GDAL笔记--chapter8
摘要：#坐标系主要有大地坐标系、地理坐标系和投影坐标系 #其中地理坐标系使用三维球面定义地球表面位置，以实现经纬度对点位的应用。（椭球体参数和基准面）,可以叫做地理投影或者未投影 #投影坐标系是平面坐标，使用X、Y值的坐标系统描述点的位置。它由基准面和投影方法确定。 #等角、等积、等方位投影... fro 阅读全文
posted @ 2020-03-25 20:35 Rser_ljw 阅读(1090) 评论(0) 推荐(0)
GDAL笔记-chapter6
摘要：这一章是利用OGR处理几何要素。 首先介绍一下OGR常量表示的不同几何类型。其中注意线是wkbLineString\多边形环是wkbLinearRing，特别注意多边形环用于多边形中。 点 wkbPoint 多点 wkbMultiPoint 线 wkbLineString 多线 wkbMultiLi 阅读全文
posted @ 2020-03-24 22:27 Rser_ljw 阅读(1119) 评论(0) 推荐(0)
GDAL笔记-chapter5
摘要：1.属性过滤条件（作用于图层lyr） #########属性过滤条件（作用于图层） import sys from osgeo import ogr import ospybook as pb ds = ogr.Open(r'') lyr = ds.GetLayer('') lyr.SetAttri 阅读全文
posted @ 2020-03-23 19:54 Rser_ljw 阅读(918) 评论(0) 推荐(0)
GDAL笔记-chapter4
摘要：ogr支持多图层的数据格式。 1.spatialite数据库 ####支持多图层的数据格式 from osgeo import ogr import ospybook as pb def print_layers(fn): ds = ogr.Open(fn, 0) if ds is None: ra 阅读全文
posted @ 2020-03-23 11:35 Rser_ljw 阅读(545) 评论(0) 推荐(0)
GDAL笔记-chapter3
摘要：1.打印图层中前十个要素的属性字段 方法1 import sys from osgeo import ogr fn = r'E:\桌面文件保存路径\gdal\osgeopy-data\osgeopy-data\global\ne_50m_populated_places.shp' ds = ogr. 阅读全文
posted @ 2020-03-23 11:17 Rser_ljw 阅读(624) 评论(0) 推荐(0)
遥感图像批量下载
摘要：介绍一下遥感图像批量下载的方法，主要是针对Landsat和Sentinel数据进行下载。两种数据在USGS上都能找到，速度比较慢。Sentinel数据还可以在https://scihub.copernicus.eu/上面找到，下载速度很慢，本人下载需要连vpn，否则会报错。今天主要利用国外友人开发的 阅读全文
posted @ 2020-02-29 20:29 Rser_ljw 阅读(1038) 评论(6) 推荐(0)
图像数组运算相关问题（nan/inf）
摘要：新年第一更！祝愿新的一年技术长足进步哈！ 最近在用sklearn的回归分析模型拟合预测遥感图像，遇到了一些问题，好在一一解决，现在总结一下。 1.首先输入sklearn的数据必须reshape(-1,1),只有转换成这种格式才能进行输入。对于多因子的数据，要concatenate((a,b),axi 阅读全文
posted @ 2020-01-01 23:41 Rser_ljw 阅读(434) 评论(0) 推荐(0)