python使用gdal将shp文件转为TIF

方法一

# 缺少获取shp文件坐标系的步骤
def vector2raster(inputfilePath, outputfile, resp):sf = shapefile.Reader(inputfilePath)# 读取shp四至min_x, min_y, max_x, max_y = sf.bboxresp = 100  # 设置每个cell的大小tifrow = int((max_x - min_x) / resp)tifcol = int((max_y - min_y) / resp)vector = ogr.Open(inputfilePath)layer = vector.GetLayer()targetDataset = gdal.GetDriverByName('GTiff').Create(outputfile, tifrow, tifcol, 1, gdal.GDT_Byte)transform = (min_x, resp, 0, min_y, 0, resp)targetDataset.SetGeoTransform(transform)# targetDataset.SetProjection(data.GetProjection())band = targetDataset.GetRasterBand(1)NoData_value = -999band.SetNoDataValue(NoData_value)band.FlushCache()gdal.RasterizeLayer(targetDataset, [1], layer, )inputfilePath = "input_shp.shp"  # 输入的待转换为tif的shp文件
outputfile = 'ouput_tif.tif'  # 转为TIF后的存储地址
resp = 100  # 栅格中每个cell的大小vector2raster(inputfilePath, outputfile, 100)

方法二

# 003 通过模板设置坐标系
import gdal
from osgeo.gdal import gdalconst
from osgeo import gdal, ogr
import numpy as np
import os
import glob
import shapefiledef vector2raster(inputfilePath, outputfile, tempaltefile):inputfilePath = inputfilePathoutputfile = outputfiletempaltefile = tempaltefiledata = gdal.Open(tempaltefile, gdalconst.GA_ReadOnly)x_res = data.RasterXSizey_res = data.RasterYSizevector = ogr.Open(inputfilePath)layer = vector.GetLayer()targetDataset = gdal.GetDriverByName('GTiff').Create(outputfile, x_res, y_res, 1, gdal.GDT_Byte)targetDataset.SetGeoTransform(data.GetGeoTransform())# targetDataset.SetProjection(data.GetProjection())band = targetDataset.GetRasterBand(1)NoData_value = -999band.SetNoDataValue(NoData_value)band.FlushCache()gdal.RasterizeLayer(targetDataset, [1], layer, )inputfilePath = "input_shp.shp"
outputfile = 'output_tif.tif'
tempaltefile = 'tempaltefile .tif'
vector2raster(inputfilePath, outputfile, tempaltefile)

python使用gdal将shp文件转为TIF相关推荐

java使用GDAL将shp文件转为geojson文件
java使用GDAL将shp文件转为geojson文件代码如下(示例): //将shp文件转为geojson文件//shp文件所在的位置String strVectorFile = "E: ...
python使用gdal读取shp文件
from __future__ import division import pandas as pd from osgeo import ogr, osr, gdal import osdef re ...
使用NetTopologySuite将shp文件转为GeoJson格式
使用NetTopologySuite将shp文件转为GeoJson格式需要引入的库(NetTopologySuite可从Nuget中下载): using NetTopologySuite.Featu ...
python shp文件_对python 读取线的shp文件实例详解
如下所示: import shapefile sf = shapefile.reader("e:\\1.2\\cs\\dx_csl.shp") shapes = sf.shapes ...
批量将NC文件转为tif格式
批量将多年的NC文件转为tif格式进行处理(来源:https://www.geek-share.com/detail/2763962738.html) 所参考和借鉴的文章的链接如下: https:// ...
Python 画降雨图，mash 数据掩膜裁剪出想要的区域，根据shp文件裁剪tif数据
1. 背景介绍 2. 导库 import numpy as np import xarray as xr import matplotlib.pyplot as plt from matplotlib ...
MATLAB：读取nc文件并将nc文件转为tif文件输出
NC(NetCDF network Common Data Form,译为"网络通用数据格式") 文件是将多个时段的数据存储在一个文件当中. 内容格式如下: •从数学上来说,net ...
python转csv_python脚本如何将Excel文件转为csv文件（代码）
本篇文章给大家带来的内容是关于python脚本如何将Excel文件转为csv文件(代码),有一定的参考价值,有需要的朋友可以参考一下,希望对你有所帮助.#!/usr/bin/env python __ ...
【数据转换】ascii文件转换为tif文件的函数
目的:脱离arcpy环境(ArcGIS平台),在python环境下实现ascii文件转换为tif文件格式原因:在利用numpy库进行数值运算的过程中,多次涉及tif文件与ascii文件的转换.为了提 ...

python使用gdal将shp文件转为TIF

python使用gdal将shp文件转为TIF

方法一

方法二

python使用gdal将shp文件转为TIF相关推荐

最新文章

热门文章