![](https://img-blog.csdnimg.cn/20191127113009487.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
GDAL
llc的足迹
print("hello csdn!")
展开
-
ubuntu从源码编译gdal
【代码】ubuntu从源码编译gdal。原创 2023-11-20 14:55:59 · 400 阅读 · 0 评论 -
ogrmerge 合并矢量数据解决中文乱码问题
lco=["ENCODING=UTF-8"]可以将原始编码转换为UTF8原创 2023-04-03 17:22:10 · 332 阅读 · 0 评论 -
查找遥感影像黑边
其中输入tif数据是墨卡托(3857)投影,所以使用wgs84_to_mercator健左边转换为经纬度(4326)extent 为影像的实际范围(左下角,右上角):【_bottom, _left, _top, _right】,原创 2023-03-21 16:01:49 · 149 阅读 · 0 评论 -
python GDAL遥感影像云优化(COG)
python GDAL遥感影像云优化(COG)原创 2023-03-09 10:02:43 · 521 阅读 · 0 评论 -
python GDAL遥感影像16位转8位——gdal_translate实现
例如,此选项可用于将16位(0至65535)数据缩放为8位(0到255)。前两个描述输入图像的最小值和最大值,后一个描述输出栅格的相同值。(其中bn是频段号),可以多次使用该选项来应用于特定波段。,则值将重新缩放为8位等效值(0到255)。也被省略,该函数将自动计算输入栅格的值范围。文档中提到,通过使用。原创 2023-03-09 09:56:50 · 2325 阅读 · 0 评论 -
python GDAL遥感影像创建缩略图
利用Translate创建缩略图功能,指定缩略图宽度为1024,等比例缩放。原创 2023-03-08 15:45:36 · 890 阅读 · 0 评论 -
GDAL坐标转换——TransformPoint
GDAL坐标转换——TransformPoint原创 2022-09-06 11:38:13 · 1214 阅读 · 0 评论 -
GDAL更新左上角坐标
>>> from osgeo import gdal>>> ds = gdal.Open('tile1.tif', gdal.GA_Update)# 获取旧坐标>>> ds.GetGeoTransform()(12137892.901028791, 0.5971642834779382, 0.0, 2797283.998548957, 0.0, -0.5971642834779382)# 设置新坐标>>> ds.SetGeo原创 2021-09-10 09:59:21 · 234 阅读 · 0 评论 -
python GDAL和PIL图像转换
读取数据先来对比一下GDAL和PIL数据读取的结果GDAL:ds: gdal.Dataset = gdal.Open('test.jpg')gdal_arr = ds.ReadAsArray()print(gdal_arr.shape)# (3, 256, 256)PIL:img = Image.open('test.jpg')pil_arr = np.array(img)print(pil_arr.shape)# (256, 256, 3)转换发现GDAL和PIL的图像存原创 2020-10-26 15:22:28 · 670 阅读 · 0 评论 -
python GDAL可视域分析
gdal_viewshedgdal_viewshed -md 10000 -ox 12509311 -oy 4769484 -a_nodata 0 srtm_59_05_3857.tif viewshed.tifgdal.ViewshedGeneratefrom osgeo import gdalds = gdal.Open('srtm_59_05_3857.tif')band = ds.GetRasterBand(1)gdal.ViewshedGenerate(srcBand=band,原创 2020-10-02 13:46:31 · 1702 阅读 · 2 评论 -
python GDAL遥感影像投影变换
gdalwarpgdalwarp ttt.tif out.tif -s_srs EPSG:4326 -t_srs EPSG:3857gdal.Warpfrom osgeo import gdaloptions = gdal.WarpOptions(format='GTiff', srcSRS='EPSG:4326', dstSRS='EPSG:3857')gdal.Warp('out.tif', 'ttt.tif', options=options)原创 2020-10-02 13:46:17 · 2402 阅读 · 0 评论 -
python GDAL遥感影像格式转换
gdal_translategdal_translate -of PNG geotiff.tif geotiff.pnggdal.Translatefrom osgeo import gdaloptions = gdal.TranslateOptions(format='PNG')gdal.Translate('ttt.png', 'ttt.tif', options=options)原创 2020-10-02 13:45:54 · 2088 阅读 · 0 评论 -
python GDAL遥感影像创建金字塔
gdaladdogdaladdo -ro geotiff.tif 2 4 8 16BuildOverviewsfrom osgeo import gdalds = gdal.Open('ttt.tif')ds.BuildOverviews(overviewlist=[2, 4 ,8, 16])原创 2020-10-02 13:45:40 · 1879 阅读 · 2 评论 -
python GDAL获取栅格数据json格式的元数据
from osgeo import gdal# 打印版本号print(gdal.__version__)# 打开数据ds=gdal.Open('ttt.tif')metadata = gdal.Info(ds, format='json')print(metadata){ "description": "ttt.tif", "driverShortName": "GTiff", "driverLongName": "GeoTIFF", "files": ["ttt.tif"],原创 2020-09-11 17:04:18 · 679 阅读 · 0 评论 -
python+gdal+numpy计算ndvi
from osgeo import gdalimport numpy as npds: gdal.Dataset = gdal.Open('gf.tif')# 红波段red = ds.GetRasterBand(3).ReadAsArray() * 0.0001# 近红外波段nir = ds.GetRasterBand(4).ReadAsArray() * 0.0001np.seterr(all="ignore")ndvi = (nir - red) / (nir + red)#原创 2020-09-05 11:01:09 · 1967 阅读 · 1 评论 -
python+gdal+numpy实现影像uint16转uint8
利用GDAL和numpy实现影像uint16到uint8的转换,支持整图转换和切图转换,高效率转换并保存为tif格式、png格式优点:简单高效,代码逻辑相对简单缺点:使用numpy数组计算,超级耗内存实现逻辑:使用percentile计算每个波段的直方图的2%、98%的值小于2%的值赋值为2%出的值,大于98%的值赋值为98%的值将数组拉伸到0-255from PIL import Imagefrom osgeo import gdalimport numpy as npde原创 2020-08-28 16:18:19 · 2100 阅读 · 0 评论 -
使用python类型注解增强gdal代码自动补全
在pycharm中编写gdal相关代码时,使用dataset=gdal.Open()返回的dataset没法自动补全相关方法,可用下面那方法解决:from osgeo import gdalgdal.UseExceptions()# 这里使用python的”类型注解“语法指定dataset的类型为gdal.Dataset,之后可实现自动补全dataset: gdal.Dataset =...原创 2020-04-20 10:24:10 · 723 阅读 · 3 评论 -
gdal.Translate和gdal.Wrap支持的creationOptions(命令行中的-co参数)
在使用gdal_translate和gdalwarp时,使用-co参数时有时候不知道支持哪些可选项,接下来研究一下:$ gdal_translate.exe --format GTiffFormat Details: Short Name: GTiff Long Name: GeoTIFF Supports: Raster Extensions: tif tiff Mime...原创 2020-04-07 15:38:16 · 3213 阅读 · 1 评论 -
OGR遍历要素及重置
OGR遍历要素及重置feature = layer.GetNextFeature()while feature: # do something ... feature = layer.GetNextFeature()# resetlayer.ResetReading()原创 2020-02-23 15:55:14 · 857 阅读 · 0 评论 -
OGR几何关系与操作
OGR几何关系与操作关系本文脚本需要导入ogr模块from osgeo import ogrEquals两个几何的边界、内部和外部重合a1=ogr.CreateGeometryFromWkt("POINT(1 1)")a2=ogr.CreateGeometryFromWkt("POINT(1 1)")print(a1.Equals(a2)) # TrueContains几...原创 2020-02-04 15:08:57 · 2859 阅读 · 3 评论 -
GDAL设置shapefile编码
GDAL设置shapefile编码设置全局变量# 设置为空时可自动识别编码,GBK和UTF-8都可以识别gdal.SetConfigOption("SHAPE_ENCODING", "")查看options属性字段支持的属性有两种方法可以查看:注:ogrinfo --formats命令可以打印出支持的所有驱动名称1.命令行>ogrinfo --format "ESRI Sh...原创 2020-01-11 09:38:59 · 1570 阅读 · 0 评论 -
GDAL算法进度条使用
GDAL算法进度条使用以Warp算法为例,自定义callback函数,打印出算法处理进度,进度值范围:0-1另:可以通过creationOptions设置相关属性from osgeo import gdalin_ds = gdal.Open('test.tif')out_tif = 'out.tif'def callback(v1, v2, v3): print(v1)...原创 2020-01-11 09:39:39 · 592 阅读 · 4 评论 -
python3安装FileGDB读写驱动
python3安装FileGDB读写驱动环境系统:win10x64python:3.6.5x64安装步骤安装gdal: GDAL‑2.4.1‑cp36‑cp36m‑win_amd64.whlpip install GDAL-2.4.1-cp36-cp36m-win_amd64.whl安装FileGDBAPI.dll:下载地址拷贝FileGDBAPI.dll至C:\Pr...原创 2019-11-27 11:25:17 · 1716 阅读 · 5 评论