python 矢量数据转栅格数据代码实例

python 矢量数据转栅格数据代码实例,第1张

概述python 矢量数据栅格数据代码实例 这篇文章主要介绍了python 矢量数据转栅格数据代码实例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下 投影包osr与proj4的使用 osr投影转换示例 from osgeo import osr,ogr #定义投影 #wgs84 source=osr.SpatialReference() source.ImportFromEPSG(4326) #google target=osr.SpatialReference() target.Imp

这篇文章主要介绍了python 矢量数据转栅格数据代码实例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下

投影包osr与proj4的使用

osr投影转换示例

from osgeo import osr,ogr#定义投影#wgs84source=osr.SpatialReference()source.importFromepsg(4326)#Googletarget=osr.SpatialReference()target.importFromepsg(3857)#简单投影转换coordTrans=osr.Coordinatetransformation(source,target)#点转换coordTrans.transformPoint(117,40)#点数组转换coordTrans.transformPoints([(117,40),(117.5,39.5)])#SF几何对象转换g=ogr.CreateGeometryFromWkt("POINT(117 40)")#转换前wgs84print(g.ExportToWkt())print(g.GetX(),g.GetY())#转换后Googleg.transform(coordTrans)print(g.ExportToWkt())print(g.GetX(),g.GetY())2.投影转换示例from pyproj import Proj,Geod,transform# projection 1: UTM zone 15,grs80 ellipse,NAD83 datum# (defined by epsg code 26915)p1 = Proj(init='epsg:26915')# projection 2: UTM zone 15,clrk66 ellipse,NAD27 datump2 = Proj(init='epsg:26715')#点的转换(首先将地理坐标转换成p1投影坐标系下的平面直角坐标,再将x1,y1转换到p2投影坐标系下,最后将p2投影坐标系下的平面直角坐标转换成地理坐标)x1,y1=p1(-92.199881,38.56694)x2,y2 = transform(p1,p2,x1,y1)print('%9.3f %11.3f' % (x1,y1))print('%9.3f %11.3f' % (x2,y2))print('%8.3f %5.3f' % p2(x2,y2,inverse=True))#点数组的转换lats = (38.83,39.32,38.75)lons = (-92.22,-94.72,-90.37)x1,y1=p1(lons,lats)x2,y2=transform(p1,y1)xy=x1+y1print('%9.3f %9.3f %9.3f %11.3f %11.3f %11.3f' % xy)xy=x2+y2print('%9.3f %9.3f %9.3f %11.3f %11.3f %11.3f' % xy)lons,lats = p2(x2,inverse=True)xy=lons+latsprint('%8.3f %8.3f %8.3f %5.3f %5.3f %5.3f' % xy)p1 = Proj(proj='latlong',datum='wgs84')x1 = -111.5; y1 = 45.25919444444p2 = Proj(proj="utm",zone=10,datum='NAD27')x2,y1)print("%s %s" % (str(x2)[:9],str(y2)[:9]))

栅格数据投影转换

#栅格数据投影转换from osgeo import gdal,osrfrom osgeo.gdalconst import *#源图像投影source=osr.SpatialReference()source.importFromepsg(32650)#目标图像投影target=osr.SpatialReference()target.importFromepsg(3857)coordTrans=osr.Coordinatetransformation(source,target)#打开源图像文件ds=gdal.Open("fdem.tif")#仿射矩阵六参数mat=ds.GetGeotransform()#源图像的左上角与右下角像素,在目标图像中的坐标(ulx,uly,ulz)=coordTrans.transformPoint(mat[0],mat[3])(lrx,lry,lrz ) = coordTrans.transformPoint(mat[0] + mat[1]*ds.RasterXSize,mat[3] + mat[5]* ds.RasterYSize )#创建目标图像文件(空白图像),行列数、波段数以及数值类型仍等同原图像driver=gdal.GetDriverByname("GTiff")ts=driver.Create("fdem_lonlat.tif",ds.RasterXSize,ds.RasterYSize,1,GDT_UInt16)#转换后图像的分辨率resolution=(int)((lrx-ulx)/ds.RasterXSize)#转换后图像的六个放射变换参数mat2=[ulx,resolution,-resolution]ts.SetGeotransform(mat2)ts.SetProjection(target.ExportToWkt())#投影转换后需要做重采样gdal.Reprojectimage(ds,ts,source.ExportToWkt(),target.ExportToWkt(),gdal.GRA_Bilinear)#关闭ds = Nonets= None

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持我们。

总结

以上是内存溢出为你收集整理的python 矢量数据转栅格数据代码实例全部内容,希望文章能够帮你解决python 矢量数据转栅格数据代码实例所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/langs/1199722.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-06-04
下一篇 2022-06-04

发表评论

登录后才能评论

评论列表(0条)

保存