weixin_42615630的博客

Python+GDAL | 读取矢量并写出txt 这篇文章主要描述了如何使用GDAL/OGR打开矢量文件、读取属性表，并将部分属性写出至txt。 代码 import ogr,sys,os import numpy as np os.chdir(r'E:\Python学习资料\GDAL\data\ospy_data1') #设置driver,并打开矢量文件 driver = ogr.GetDri...

这篇文章主要描述了如何使用GDAL/OGR创建矢量文件。 代码 import ogr,sys,os os.chdir(r'E:\') #设置driver driver = ogr.GetDriverByName('ESRI Shapefile') #打开输入的矢量 inDs = driver.Open('sites.shp', 0) if inDs is None: print("Co...

这篇文章主要描述了如何使用GDAL/OGR读取栅格数据，分为按波段、按行、按列读取。 代码 from osgeo import gdal from gdalconst import * import os, sys, time start_time = time.time() os.chdir(r'E:\') #获取所有注册类 gdal.AllRegister() #打开栅格 ds = gd...

这篇文章主要描述了如何使用GDAL/OGR打开矢量文件、读取点位的经纬度，并根据经纬度读取栅格数据。 代码 #读取矢量 from osgeo import ogr import os, sys #改变工作空间 os.chdir(r'D:\') #############获取矢量点位的经纬度 #设置driver driver = ogr.GetDriverByName('ESRI Shapefil...

这篇文章主要描述了如何使用GDAL设置掩膜，并计算平均值。 代码 from osgeo import gdal import os, sys import numpy os.chdir(r'D:\Python学习资料\GDAL\data\ospy_data4') #注册类 gdal.AllRegister() #打开栅格 ds = gdal.Open('aster.img', GA_ReadOn...

这篇文章主要描述了如何使用GDAL/OGR打开矢量文件、读取属性表头，并修改部分属性。 下列代码的实践目的是修改某矢量的唯一编码，将编码钱统一＋“BM-” 代码 import ogr import os import pandas as pd path = r'J:\' os.chdir(path) shapefile='Merge.shp' driver = ogr.GetDriverBy...