Python读取表格数据,python 栅格转矢量
安装方法:
手动下载whl安装包https://www.lfd.uci.edu/~gohlke/pythonlibs/pip安装安装包#-*-编码:utf-8-*-from osgeo导入gdal dataset=gdal . open(r f:\ mod 11a 1 . a 2018 152 . lst _ day _ 1km . TIF )#栅格矩阵的列数im_width=dataset。栅格矩阵im_height=dataset的RasterXSize#行数。RasterYSize#波段数im_bands=数据集。RasterCount#仿射矩阵,大地坐标和左上角像素的像素#有六个参数,子表表示左上角X坐标;图像在东西方向的分辨率;如果北面朝上,地图旋转角度0,表示图像的线条平行于X轴;左上角的y坐标;#如果北面朝上,地图的旋转角度,0表示图像的列平行于Y轴;地图在南北方向的分辨率。im _ geo trans=dataset . getgeotransform()#地图投影信息im _ proj=dataset . get projection()#读取一个像素的值#(1)读取一个波段,其参数是波段的索引号,波段索引号从1开始(我打开的这张图片只有一个波段)band=dataset . getrasterband(1)#(2)用READA SARAY (xoff,yoff)从(xoff,yoff)开始读出size (xsize,ysize)的矩阵,为了读取整个图像im _ data=band.readasarray (0,0,im _ width,im _ height) # (3),获取一个或几个像素的值data=im _ data [0,1046] print (data) #释放内存。如果不释放它,当你在arcgis或envi中打开图像时,显示的文件已经被del数据集占用。
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。