gdal(python,C++)
GDAL全称是 Geospatial Data Abstraction Library(地理空间数据抽象库),是一个在XMT许可协议下读写空间数据(包括栅格数据和矢量数据)的开源库,它利用抽象数据模型来表达所支持的各种文件格式,还使用一系列命令行工具来进行数据转换和处理。
陨星落云
这个作者很懒,什么都没留下…
展开
-
OGR2GUI
OGR2GUI 是一个开源项目,为 OGR2OGR 提供图形用户界面 (GUI)。OGR2OGR 是用于转换和操作地理空间数据的命令行实用程序。原创 2024-03-17 13:42:39 · 430 阅读 · 0 评论 -
在Shp文件中添加投影
【代码】在Shp文件中添加投影。原创 2023-10-13 23:19:30 · 296 阅读 · 0 评论 -
vs2019下编译gdal2.4.4
3.第57行,GDAL_HOME = "C:\warmerda\bld"这个路径是用来存储编译后的头文件、静态库、动态库等;路径随你可以修改,这里我未作修改;1:在GDAL目录下的nmake.opt文件中,找到SYM_PREFIX的定义,应该在973行左右。在GDAL目录下的nmake.opt文件中第257行,添加#,注释掉对矢量格式odbc的支持;2.在解压后的文件夹中找到nmake.opt文件,修改几个配置参数;4.第194行,去掉#,用于编译64位的GDAL。就是将后面的@开始,后面的数字删除。原创 2023-08-05 21:01:12 · 419 阅读 · 0 评论 -
GDAL解析xml
CPLXMLNode属性。原创 2023-06-26 17:54:29 · 345 阅读 · 0 评论 -
gdal创建wkt几何对象
【代码】gdal创建wkt几何对象。原创 2023-06-15 14:00:13 · 257 阅读 · 0 评论 -
使用gdal创建圆形、弧线
对原文代码进行部分改写,适用于gdal3。转载 2023-06-14 17:20:33 · 302 阅读 · 0 评论 -
postgis下载地址
postgis下载地址:http://download.osgeo.org/postgis/windows/原创 2023-06-11 18:38:45 · 1779 阅读 · 0 评论 -
GDAL读取属性表值乱码解决方法
不同系统或者服务器之间消息传递经常遇到编码转换问题,这里用C++实现了一个轻量的gbk和utf8互相转换,可跨平台使用。(重量级的可以用libiconv库)在windows下用头文件里的函数进行多字节和宽字符转换,linux下采用头文件里的函数进行编码互相解析。现在乱码的utf-8,原本gbk是正常的,那就将utf-8转换成gbk编码问题就解决了。在网上查看的推荐方法(C++代码),但没解决我遇到的问题;原创 2023-06-10 17:43:05 · 1500 阅读 · 0 评论 -
GDAL图像复制
使用函数指针实现图像复制。原创 2023-05-17 22:35:42 · 373 阅读 · 0 评论 -
GDAL图像重采样
GDAL图像重采样。转载 2022-07-18 21:47:13 · 1443 阅读 · 0 评论 -
GDAL中的地理坐标系、投影坐标系及其相互转换
文章目录地理坐标系例1:国内常用地理坐标系投影坐标系例2:国内常用投影坐标系(不推荐使用)坐标转换例3:地理坐标转投影坐标例4:投影坐标转地理坐标地理坐标系原理参考这篇文章:地理坐标系与投影坐标系区别与联系https://yunxingluoyun.blog.csdn.net/article/details/123970678例1:国内常用地理坐标系#include <cstdio>#include "gdal_priv.h"#include <iostream>原创 2022-04-12 14:46:54 · 5921 阅读 · 9 评论 -
gdal影像镶嵌并裁剪
# -*- coding: utf-8 -*-"""Created on Tue Apr 5 20:18:49 2022@author: 陨星落云"""import os from osgeo import gdaldef mergeAndClip(img_dir,clip_shp,out_dir): files = [os.path.join(img_dir,s)for s in sorted(os.listdir(img_dir))] merged_vrt = os.原创 2022-04-05 23:01:12 · 1466 阅读 · 3 评论 -
地理坐标系与投影坐标系区别与联系
地理坐标系与投影坐标系区别空间参考中主要包含大地水准面、地球椭球体、投影坐标系等几部分内容。地图投影就是把地球表面的任意点,利用一定数学法则,转换到地图平面上的理论和方法,一般有两种坐标系来进行表示,分别是地理坐标系和投影坐标系。如下图所示,描述了地理坐标系与投影坐标系之间的关系。参考椭球体表1 常见的参考椭球体及其参数椭球名称长半轴(米)短半轴(米)扁率的倒数Clarke18666378206.46356583.79999898294.9786982Besse原创 2022-04-05 16:17:55 · 4447 阅读 · 0 评论 -
gdal查看图像信息
#include <iostream>#include "gdal_priv.h"//#include "cpl_conv.h" int main(){ // 注册文件格式 GDALAllRegister(); //影像路径 const char* img_dir = "E:\\beijing_landsat\\LC8_123032_2013244_LGN02\\LC08_L1TP_123032_20130901_20180523_01_T1_B1.TIF"; //使用原创 2022-04-04 20:05:05 · 1529 阅读 · 3 评论 -
编译好的gdal下载
下载编译好的gdalhttps://www.gisinternals.com/release.php转载 2022-03-22 15:48:11 · 1293 阅读 · 1 评论 -
影像中缺失部分,其标签设为背景值0
import numpy as npimport imageioimport cv2from osgeo import osr, gdal def assign_spatial_reference_byfile(src_path, dst_path): src_ds = gdal.Open(src_path, gdal.GA_ReadOnly) sr = osr.SpatialReference() sr.ImportFromWkt(src_ds.GetProjecti原创 2021-11-03 13:40:12 · 247 阅读 · 0 评论 -
GDAL都哪些功能?
GDAL全称是 Geospatial Data Abstraction Library(地理空间数据抽象库),是一个在XMT许可协议下读写空间数据(包括栅格数据和矢量数据)的开源库,它利用抽象数据模型来表达所支持的各种文件格式,还使用一系列命令行工具来进行数据转换和处理。虽然都知道GDAL库功能强大,但是大部分却不知道,有哪些功能。在这里列出GDAL常用的一些命令行的功能,供大家参考:GDAL常用命令描述gdal2tiles.py用于创建png图像金字塔的Python脚本gd原创 2021-10-24 10:19:20 · 2177 阅读 · 0 评论 -
使用GDAL进行影像投影坐标、地理坐标、图上坐标的转换
使用GDAL库写了四个函数分别进行投影坐标与地理坐标(经纬度)之间的转换,投影坐标和图上坐标(行列号)之间的转换。有需要的朋友可以参考。直接上代码吧,因为代码很简单(Python版本)。# -*- encoding: utf-8 -*-from osgeo import gdalfrom osgeo import osrimport numpy as npdef getSRSPair(dataset): ''' 获得给定数据的投影参考系和地理参考系 :param da转载 2021-10-20 13:36:11 · 3208 阅读 · 0 评论 -
gdal栅格数据裁剪
在进行遥感影像处理的时候,我们经常需要进行裁剪的工作,来看看如何使用GDAL工具进行这项操作吧!参考资料:GDAL: gdalwarpGDAL: gdal_translateGDAL/OGR Python API使用GDAL命令GDAL提供了两个命令可以用于影像的裁剪:gdalwarp和gdal_translate,两个命令中我更推荐使用后者。gdalwarp命令可以使用-te制定裁剪范围。默认是在原数据的坐标系下的xmin ymin xmax ymax,当然我们也可以使用-te_srs参数指转载 2021-10-18 13:38:36 · 3535 阅读 · 0 评论 -
gdal将16位的tif转为8位tif
将16位的tif转为8位tifgdal_translate -of GTiff -ot Byte -scale 0 65535 0 255 16bit_img.tif 8bit_img.tif原创 2021-08-26 15:46:40 · 1401 阅读 · 2 评论 -
gdal影像镶嵌
将文件夹下的所有tif文件的文件名保存到filelist.txt中。dir /b *.tif >filelist.txt将文件列表中的tif影像进行合并保存为merged.vrt(虚拟文件格式)gdalbuildvrt -input_file_list filelist.txt merged.vrt文件格式转换将vrt格式文件转换为tif文件,并进行影像压缩gdal_translate -of GTiff merged.vrt merged.tif -co COMPRESS=LZW原创 2021-08-25 22:30:38 · 685 阅读 · 0 评论 -
基于SG滤波的时序序列数据重构(遥感影像)
# -*- coding: utf-8 -*-"""Created on Sun Jul 25 11:00:44 2021@author: xiaohuihui"""import numpy as npfrom osgeo import gdalfrom osgeo.gdalconst import GA_ReadOnlyfrom scipy.signal import savgol_filterimport imageio#import globdef img_to_sigl原创 2021-07-26 16:11:42 · 1476 阅读 · 0 评论 -
16bit转8bit+直方图均衡化
# -*- coding: utf-8 -*-"""Created on Sun Jun 27 22:55:24 2021@author: 陨星落云csdn博客地址: https://blog.csdn.net/qq_28368377"""from osgeo import gdalimport numpy as npimport imageiodef img16to8(imgpath): scale = '-scale min_val max_val' optio原创 2021-06-27 22:58:39 · 777 阅读 · 0 评论 -
img16位转8位
"""将16位遥感图像压缩至8位,并保持色彩一致"""from osgeo import gdalimport osimport globimport numpy as npimport timedef read_tiff(input_file): """ 读取影像 :param input_file:输入影像 :return:波段数据,仿射变换参数,投影信息、行数、列数、波段数 """ dataset = gdal.Open(input_原创 2021-04-19 10:07:13 · 389 阅读 · 1 评论 -
获取带坐标的影像中每个像元的坐标值
from osgeo import gdalds = gdal.Open(r'E:\beijing_landsat\2010-2020\LC81230322015250LGN01\LC08_L1TP_123032_20150907_20170404_01_T1_B1.TIF')# unravel GDAL affine transform parametersc, a, b, f, d, e = ds.GetGeoTransform()def pixel2coord(col, row):原创 2021-01-07 13:52:01 · 645 阅读 · 0 评论 -
使用gdal添加地理坐标
src_path:有地理坐标的影像dst_path:待添加地理坐标的影像from osgeo import osr, gdal def assign_spatial_reference_byfile(src_path, dst_path): src_ds = gdal.Open(src_path, gdal.GA_ReadOnly) sr = osr.SpatialReference() sr.ImportFromWkt(src_ds.GetProjectionRef()原创 2020-12-10 16:32:00 · 1261 阅读 · 0 评论 -
ogrinfo工具
ogrinfoogrinfo工具用于输出OGR支持的矢量数据的信息。查看参数信息(base) E:\遥感图像处理(python)\data>ogrinfoUsage: ogrinfo [--help-general] [-ro] [-q] [-where restricted_where|@filename] [-spat xmin ymin xmax ymax] [-geomfield field] [-fid fid] [-sq原创 2020-07-27 15:49:14 · 630 阅读 · 0 评论 -
从ee中获取Sentinel-1数据
COPERNICUS/S1_GRD数据集Sentinel-1任务从双极化C波段合成孔径雷达(SAR)仪器提供数据。该集合包括使用Sentinel-1工具箱处理的S1地面距离检测(GRD)场景,以生成经过校准的,经过正交校正的产品。该馆藏每天更新。新资产在可用后的两天内会被吸收。该集合包含所有GRD场景。每个场景具有3种分辨率(10、25或40米),4种波段组合(对应于场景极化)和3种乐器模式之一。在镶嵌环境中使用该集合可能需要过滤到一组均匀的波段和参数。请参阅本文 以获取集合使用和预处理的详细信息。每个原创 2020-07-22 16:13:30 · 1474 阅读 · 4 评论 -
国家2000坐标系的理解
文章目录中国大地坐标系2000中国大地坐标系2000 (EPSG:4490)中国大地坐标系2000(EPSG:4479)区别:以中央经度117为例,分析以下三种局部坐标的区别CGCS2000 /高斯-克鲁格CM 117E(EPSG:4509)CGCS2000 / 3度高斯-克鲁格CM 117E(EPSG:4548)CGCS2000 / 3度高斯-克鲁格 区39(**EPSG:4527**)区别:中国大地坐标系2000中国大地坐标系2000 (EPSG:4490)属性单位:度(定义表示形式的供应商)原创 2020-07-01 02:34:15 · 25751 阅读 · 0 评论 -
西安80三度带39带投影坐标系转地理坐标系
# -*- coding: utf-8 -*-"""Created on Sat Jul 18 12:27:15 2020@author: 67231"""try: from osgeo import osrexcept ImportError: import osr def CoordinateTransformation(xy): strXian80 = "+proj=tmerc +lat_0=0 +lon_0=117 +k=1 +x_0=39500000原创 2020-07-20 15:12:15 · 2784 阅读 · 0 评论 -
GDAL支持的矢量数据格式
判断ogr是否安装成功try: from osgeo import ogr print('Import of ogr from osgeo worked. Hurray!')except: print('Import of ogr from osgeo failed')生成ogr帮助文档import osgeo.ogrimport syssys.stdout=open("ogr_help.txt","w")help(osgeo.ogr)sys.stdout.close()原创 2020-07-05 17:58:13 · 1232 阅读 · 0 评论 -
在Python中安装GoogleEarthEngine的错误
GoogleEarthEngine的Python环境安装出错错误1: ee.Initialize()Traceback (most recent call last): File "<stdin>", line 1, in <module> File "D:\Anaconda3\lib\site-packages\ee\__init__.py", line 123, in Initialize http_transport=http_transport) F原创 2020-07-03 18:50:30 · 1476 阅读 · 0 评论 -
GDAL遥感影像读取
读取并显示遥感影像import numpy as npfrom osgeo import gdalfrom osgeo.gdalconst import GA_ReadOnlyimport matplotlib.pyplot as pltdef disp(infile,bandnumber): gdal.AllRegister() inDataset = gdal.Open(infile,GA_ReadOnly) # print(inDataset) cols =原创 2020-07-02 11:31:33 · 1540 阅读 · 1 评论 -
GDAL几何对象
几何对象创建一个点from osgeo import ogrpoint = ogr.Geometry(ogr.wkbPoint)point.AddPoint(1198054.34,648493.09)print(point.ExportToWkt())结果:POINT (1198054.34 648493.09 0)创建一条线from osgeo import ogrline = ogr.Geometry(ogr.wkbLineString)line.AddPoint(1116原创 2020-07-02 11:24:14 · 829 阅读 · 0 评论 -
GDAL / OGR安装测试
GDAL / OGR测试GDAL/OGR是否安装import systry: from osgeo import ogr, osr, gdalexcept: sys.exit('ERROR: cannot find GDAL/OGR modules')检查已安装的GDAL / OGR版本import sysfrom osgeo import gdalversion_num = int(gdal.VersionInfo('VERSION_NUM'))if versio原创 2020-07-01 12:56:25 · 1093 阅读 · 0 评论 -
gdal栅格数据处理(python)
栅格数据处理关闭栅格数据集这里显示如何关闭栅格数据集。在脚本的中间非常有用,它可以恢复访问数据集所拥有的资源,删除文件等。在脚本的末尾没有必要,因为Python垃圾收集器会在脚本运行时自动执行相同的操作。脚本退出。from osgeo import gdal# 打开数据集ds = gdal.Open("jianda.tif")# 关闭数据集ds = None获取栅格数据集的元数据获取栅格元数据以进行快速的分辨率检查from osgeo import gdalgtif = gda原创 2020-06-30 18:08:07 · 2146 阅读 · 0 评论