遥感影像读取和处理
我叫杨傲天
做一个谦卑的绅士
展开
-
python 矢量转栅格
import osimport gdalimport numpy as npfrom osgeo import ogrfrom geopandas import *#读取矢量数据元数据信息(要素包括 矩形坐标,字段及类型等)vdata_dir='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第七周\\Lesson7\\2-Second'filename=os.path.join(vdata_dir,'FirstLayer.shp')#print(fi.原创 2021-08-19 15:58:36 · 824 阅读 · 0 评论 -
python 矢量数据的读取
import osimport sysfrom osgeo import ogrfrom osgeo import gdaldriver=ogr.GetDriverByName('ESRI Shapefile')filename='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第六周\\Lesson6_osgeopy-data\\global\\ne_50m_populated_places.shp'dataSource=driver.Open(filen.原创 2021-08-19 11:25:54 · 905 阅读 · 0 评论 -
python 影像拼接
import gdal##拼接,不只是简单地把两张影像拼接,更重要的是地理位置配准#打开、读取影像filename1='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第五周\\Lesson5影像数据\\拼接影像\\test1.tif'filename2='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第五周\\Lesson5影像数据\\拼接影像\\test2.tif'ds1=gdal.Open(filename1)ds.原创 2021-08-14 09:51:10 · 1266 阅读 · 7 评论 -
python 裁剪影像
#栅格数据裁剪#添加包import gdalfilename='C:\\Users\\admin\\Desktop\\空间信息处理与系统开发资料\\第五周\\Lesson5影像数据\\裁剪影像\\test1.tif'dataset=gdal.Open(filename)#一共5个波段,读取每个波段band1=dataset.GetRasterBand(1)band2=dataset.GetRasterBand(2)band3=dataset.GetRasterBand(3)band4=d.原创 2021-08-11 19:48:30 · 809 阅读 · 0 评论 -
Python 循环读取影像 计算NDVI
from osgeo import gdalimport numpy as npfrom osgeo import gdalimport os#循环读取FilePath='D:\\影像练习'tif_name=os.listdir(FilePath) #读取FilePath下的目录名,包含扩展名SavePath='D:\\影像练习\\NDVI输出'print(tif_name)for i in range(len(tif_name)): #i从0开始到len(tif_name)-1 .原创 2021-08-09 08:33:57 · 587 阅读 · 1 评论 -
python 影像读取和计算NDVI
from osgeo import gdalfilename='D:\\影像练习\\001.tif'dataset=gdal.Open(filename)if dataset is None: print('could not open'+filename)cols=dataset.RasterXSize #列数rows=dataset.RasterYSize #行数 空间坐标 可以理解为行列号print(cols,rows)bands=dataset.RasterCount .原创 2021-08-08 16:21:35 · 4530 阅读 · 3 评论