目录
序言：
? ? ? ? 本文学习内容是来自Utah State University的openGIS课程：
相关资料，包括讲义、源码、数据样例
https://www.gis.usu.edu/~chrisg/python/
- OGR矢量库：简单的矢量数据读写，是GDAL的一部分
- GDAL地理空间数据抽象库： a、读写栅格数据 b、ArcGIS也是基于GDAL开发的 c、C++库，但是可以用python调用
需要的模块：
1、进入https://www.lfd.uci.edu/~gohlke/pythonlibs/? 搜索osgeo，根据自己python版本下载需要的gdal版本
?2、numpy（高速的数组处理，对栅格数据尤其重要）
cmd：pip install numpy
读矢量数据：
导入库：
import ogr或者：
from osgeo import ogr万能方法：
try:
import ogr
except:
from osgeo import ogr读取数据之前必须先载入数据驱动，初始化一个对象，让它‘知道’某种数据结构：
try:
import ogr
except:
from osgeo import ogr
driver = ogr.GetDriverByName('ESRI Shapefile')数据驱动driver的open方法中，返回一个数据源对象, update为0是只读，为1是可写，file_name（绝对路径)
driver.opem(<file_name>, <update>)读取数据层:
ESRI的shapefile都是填0的，如果不填的话默认也是0.
layer = dataSoure.GetLayer(0)获取数据层有多少个点：
n = layer.GetFeatureCount()
print(n)
读取上下边界：
extent = layer.GetExtent()
print('extent:', extent)
print('ul:', extent[0], extent[3])#上下
print('lr:', extent[1], extent[2])#左右读取某一要素，此处读取的是一个点
feat = layer.GetFeature(41)
fid = feat.GetField('id')
print(fid)循环读取要素，遍历所有：
feat = layer.GetNextFeature() #读取下一个
while feat:
feat = layer.GetNextFeature()
later.ResetReading()提取feature的几何形状：
geom = feat.GetGeometryRef()
geom.GetX()
geom.GetY()
print(geom)释放内存：
feature.Destroy()
dataSource.Destroy()写矢量数据：
创建一个新的文件：
这个文件如果已经存在就会报错：
driver.CreateDataSource(<filename>)创建一个新的layer：
dataSource.CreateLayer(<name>,CreateLayer(<name>, geom_type=<OGRwkbGeometryType>, [srs])删除一个shp文件：
driver.DeleteDataSource(‘test.shp’)要添加一个新字段，只能在layer里面加，而且还不能有数据添加的字段如果是字符串，还要设定宽度：
fieldDefn = ogr.FieldDefn('id', ogr.OFTString)
fieldDefn.SetWidth(4)
layer.CreateField(fieldDefn)然后从layer中读取相应的feature类型，并创建feature：
featureDefn = layer.GetLayerDefn()
feature = ogr.Feature(featureDefn)设定几何形状：
feature.SetGeometry(point)设定某字段的值：
feature.SetField('id', 23)将feature写入layer：
layer.CreateFeature(feature)