![](https://img-blog.csdnimg.cn/20201014180756919.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Gdal
John H.
一直在学习,但还是很菜...
展开
-
基于OSGEO的矢量文件多部件拆分(不使用arcpy)
矢量文件多部件拆分原创 2022-09-23 11:42:49 · 191 阅读 · 0 评论 -
基于OSGEO的矢量文件的空间融合(不使用Arcpy)
矢量文件的空间融合,可以消除重叠图斑原创 2022-09-23 11:41:35 · 166 阅读 · 0 评论 -
基于OSGEO对单个矢量文件的多个feature(图斑)分别导出(不使用arcpy)
导出矢量文件中包含的多个图斑原创 2022-09-23 11:39:34 · 177 阅读 · 0 评论 -
基于OSGEO的遥感影像窗口滑动裁剪
对遥感影像滑动窗口裁剪.原创 2022-09-23 11:37:47 · 358 阅读 · 0 评论 -
基于osgeo创建渔网(不使用arcpy)
def shapefile_fishnet_creat( _input_boundary_shapefile, _output_grid_fishnet, flag=1, _rows=5, _cols=5, _gridHeight=100, _gridWidth=100): """ 根据矢量文件边界创建渔网 _input_boundary_shapefile需要创建渔网的矢量边界,坐标系一定要是投影坐标系 _output_gr.原创 2022-05-03 11:58:22 · 433 阅读 · 1 评论 -
矢量面文件指定距离做缓冲(不用arcpy函数)
代码如下.......def shapefile_create_buffer( _input_shapefile, buffer_result_shapefile, buffer_distance): """ _input_shapefil需要缓冲的矢量文件 buffer_result_shapefile缓冲后的矢量文件 buffer_distance缓冲距离,单位:米 """ input_ds = ogr.Open(_input_sh原创 2022-05-03 11:53:39 · 435 阅读 · 0 评论 -
基于osgeo判断矢量文件是点、线还是面文件
# -*- coding: utf-8 -*-"""Created on Tue Jun 22 09:02:18 2021@author: DeepLearning"""from osgeo import ogrdef judge_geom_type(_input_shapefile): '''判断一个矢量文件是面、线还得点, 如果是面返回值为:Polygon 如果是线返回时为:Line String: 如果是点返回值为:Point''' _ds = .原创 2021-08-12 16:39:52 · 382 阅读 · 0 评论 -
基于GDAL矢量文件边界平滑
见如下函数# -*- coding: utf-8 -*-"""Created on Tue Jun 22 09:25:41 2021@author: DeepLearning"""from osgeo import ogrfrom pathlib import Pathdef shapefile_edge_smooth(_input_shapefile, _smoothed_shapefile_output_path, _buffer_distance=2): '''平滑矢量文原创 2021-08-12 15:55:46 · 995 阅读 · 3 评论 -
多个矢量文件合并为一个矢量文件
见如下代码。这个代码是根据我自己业务需要写的,整体是可以根据自己的需要进行适当调整和修改。许多坑已经都跨过去了。def mutil_shapefiles_merge(input_folder_path, output_shpfile_abspath, delete_area_threshold=0.03): '''本函数用于多个矢量文件的合并 第一个参数为需要合并矢量文件的文件夹路径 第二个是合并后矢量文件的完整路径''' ptdriver = ogr.GetDri原创 2021-08-10 17:29:18 · 2582 阅读 · 2 评论 -
从矢量文件导出指定的一条记录,基于GDAL的矢量数据操作
本函数用于从shpfile文件从导出一条shp格式记录(图斑)def extract_single_feature(input_shpfile, output_shpfile, specific_field_name, specific_value): '''本函数用于从shpfile文件从导出一条shp格式记录(图斑) 第一个参数是输入shpfile的完整路径和文件名 第二个参数是输出shpfile的完整路径和文件名 第三个参数是根据哪个字段名称导出记录 第四原创 2021-08-10 17:25:12 · 255 阅读 · 0 评论 -
基于矢量文件的栅格数据批量裁剪与导出
见如下代码。def clip_raster_by_shpfile_save_into_tif_formate_jpg_formate(input_raster, input_shpfile, output_raster, output_jpg): '''这个函数是用shpfile裁剪raster,同时保存两个结果到指定文件夹中 一个是tif格式带坐标系的tif栅格数原创 2021-08-10 17:22:13 · 818 阅读 · 0 评论 -
将多个shapfile矢量文件合并为一整个shapfile矢量文件
用于多个矢量文件的合并。还有一些关于坐标系的问题尚未解决,不过整体没有问题可以使用,后续填补这些小问题。import osfrom osgeo import ogr, osrimport globdef mutil_layers_union(input_folder_path, output_shpfile_abspath): '''本函数用于多个矢量文件的合并 第一个参数为需要合并矢量文件的文件夹路径 第二个是合并后矢量文件的绝对路径''' ptdriv原创 2021-07-23 09:44:05 · 2444 阅读 · 0 评论 -
栅格数据波段计算与合并
本函数是利用rgb三波段影像,计算RGBGDSI植被指数,并将该植被指数作为1个波段与原始影像合并生成一个4波段影像数据,用以增加影像数据信息量。由于影像数据量可能很大,一次将1个波段读入内存可能会发生内存溢出的情况 本函数使用了分块读取和写入,block_size是分块大小,可根据实际情况调整,默认值为一次读入、写入5000行列数据import sysimport mathimport numpyfrom osgeo import gdalfrom osgeo import gdal_..原创 2021-07-23 09:40:10 · 1544 阅读 · 0 评论 -
基于gdal和osgeo库函数的矢量文件字段操作(不使用arcpy)
内容:shapefile文件字段添加、删除、赋值、面积计算(1)shapefile字段添加shp_file = r'' #shapfile绝对路径# 获取文件驱动driver = ogr.GetDriverByName('ESRI Shapefile')# 注意,如果这里第二个参数是0的话就是仅仅读取数据,不做写入或修改# 我们要新建字段,所以第二个参数要设置为1————1 is read/writedata_source = driver.Open(shp_file, 1)原创 2021-07-23 09:35:57 · 1060 阅读 · 0 评论