当前位置:首页 > 代码相关 > 正文内容

GDAL使用(二)

admin4年前 (2021-02-04)代码相关12565

本篇说一下CreateCopy()和Create()

CreateCopy()是从一个源数据拷贝到一个新的数据中。

有关说明写到了代码的注释中。C++代码和Python代码实现的功能是一致的。

C++代码:

#include <iostream>
#include "gdal_priv.h"
#include "cpl_conv.h" 
#include <cpl_string.h>

int main()
{
	clock_t start = clock();
	GDALAllRegister();
	const char* srcDs = "E:\testdata\CZIOrigin.tiff";
	const char* dstDs = "O:\landsatdata.tif";
	const char* pszFormat = "GTiff";
	//打开数据
	GDALDataset* poSrcDs = (GDALDataset*)GDALOpen(srcDs, GA_ReadOnly);
	//创建驱动
	GDALDriver* poDriver = GetGDALDriverManager()->GetDriverByName(pszFormat);
	if (poDriver == NULL)
		exit(1);

	//相关创建选项可以到栅格驱动程序中查看
	char** papszOptions = NULL;
	papszOptions = CSLSetNameValue(papszOptions, "TILED", "YES");
	papszOptions = CSLSetNameValue(papszOptions, "COMPRESS", "LZW");
	papszOptions = CSLSetNameValue(papszOptions, "NUM_THREADS", "ALL_CPUS");

	//创建拷贝数据集
	GDALDataset* poDstDs = poDriver->CreateCopy(dstDs, poSrcDs, FALSE, papszOptions, GDALTermProgress, NULL);

	//释放资源
	if (poDstDs != NULL)
	{
		GDALClose((GDALDatasetH)poDstDs);
	}
	GDALClose((GDALDatasetH)poSrcDs);
	CSLDestroy(papszOptions);

	clock_t end = clock();
	double endtime = (double)(end - start) / CLOCKS_PER_SEC;
	//以秒为单位输出时间
	std::cout << "Total time:" << endtime << std::endl;	

	printf("over!");
	system("pause");
}

Python代码:

# -*- coding: utf-8 -*-
"""
Created on Thu Feb  4 11:47:06 2021

@author: cyhu
"""


from osgeo import gdal
import time

time_start=time.time()

# 文件路径
src_file = "E:\testdata\landsatdata.tiff"
dst_file = "O:\landsatdata.tif"
# 打开数据
src_ds = gdal.Open(src_file,gdal.GA_ReadOnly)
if not src_ds:
    print("open failed!")

driver = gdal.GetDriverByName("GTiff")
# 相关的创建选项可以到栅格驱动程序中查看
dst_ds = driver.CreateCopy(dst_file,src_ds,strict = 0,
                           options=["TILED=YES", "COMPRESS=LZW","NUM_THREADS=ALL_CPUS"])

src_ds = None
dst_ds = None

time_end = time.time()

print("time ",time_end - time_start)


Create()本身并不复杂,但是需要进行一些元数据的写入,所以代码量稍多一些。

C++代码:

#include <iostream>
#include "gdal_priv.h"
#include "cpl_conv.h" 
#include <cpl_string.h>

int main()
{
	//clock_t start = clock();
	GDALAllRegister();
	//const char* srcDs = "E:\testdata\CZIOrigin.tiff";
	const char* dstDs = "O:\landsatdata.tif";
	const char* pszFormat = "GTiff";
	//创建栅格驱动程序
	GDALDriver* poDriver = (GDALDriver*)GDALGetDriverByName(pszFormat);

	//创建栅格
	char** papszOptions = NULL;
	GDALDataset* poDstDS = poDriver->Create(dstDs, 512, 512, 1, GDT_Byte, papszOptions);

	/* 创建上面的栅格数据后写入元数据和栅格数据 */
	//要写入栅格的数据块
	GByte abyRaster[512 * 512];
	//创建仿射变换参数
	double adfGeoTransform[6] = { 444720, 30, 0, 3751320, 0, -30 };
	//创建OGR空间参考
	OGRSpatialReference oSRS;
	char* pszSRS_WKT = NULL;
	//设置数据集的变换参数
	poDstDS->SetGeoTransform(adfGeoTransform);
	//设置投影分带
	oSRS.SetUTM(11, TRUE);
	//设置基准面
	oSRS.SetWellKnownGeogCS("NAD27");
	//输出到WKT格式
	oSRS.exportToWkt(&pszSRS_WKT);
	//设置投影
	poDstDS->SetProjection(pszSRS_WKT);
	//获取波段
	GDALRasterBand* poBand = poDstDS->GetRasterBand(1);
	//写入数据
	poBand->RasterIO(GF_Write, 0, 0, 512, 512, abyRaster, 512, 512, GDT_Byte, 0, 0);
	//释放资源
	CPLFree(pszSRS_WKT);
	GDALClose((GDALDatasetH)poDstDS);

	std::cout << "over!" << std::endl;
	system("pause");
	return 0;
}

Python代码:

# -*- coding: utf-8 -*-
"""
Created on Thu Feb  4 17:21:33 2021

@author: cyhu
"""


from osgeo import gdal
from osgeo import osr
import numpy


# 获取栅格驱动
driver = gdal.GetDriverByName("GTiff")
dst_ds = "O:\test.tif"
# 创建栅格
dataset = driver.Create(dst_ds, xsize = 512, ysize = 512, bands = 1, eType = gdal.GDT_Byte)
# 设置仿射变换参数
dataset.SetGeoTransform([444720, 30, 0, 3751320, 0, -30])
# 设置空间参考
spatialref = osr.SpatialReference()
spatialref.SetUTM(11,1)
spatialref.SetWellKnownGeogCS("NAD27")
dataset.SetProjection(spatialref.ExportToWkt())
# 创建要写入的数据
raster = numpy.zeros((512, 512), dtype=numpy.uint8)
# 写入数据
dataset.GetRasterBand(1).WriteArray(raster)
# 释放资源
dataset = None


扫描二维码推送至手机访问。

版权声明:本文由lovedm.club发布,如需转载请注明出处。

本文链接:https://lovedm.club/?id=116

标签: GDAL
分享给朋友:

“GDAL使用(二)” 的相关文章

C#的类型系统

C#的类型系统

C#的五大数据类型:    类(Class)    结构体(Structures)    枚举(Enumerations)    接口(In...

C# Lambda表达式

简单用法,一句一句来,便于理解 Func<int, int, int> func = new Func<int, int, int>((int a, int b) => { return a * b; });(int a, int b) => { ret...

C# 正则表达式(2)

// pattan = @"[^ahou]"; 表示匹配除ahou之外的字符,^在表示反义 string res4 = Regex.Replace(s, @"[^ahou]",&...

C# 与文件相关的几个类(2)

Directory类:静态类,主要处理文件目录。方法:CreateDirectory(String)在指定路径中创建所有目录和子目录,除非它们已经存在。返回值是一个DirectoryInfo对象Delete(String)从指定路径删除空目录。无返回值。Exists(String)确定给定路径是否引...

C、C++、C#交换变量

C、C++、C#交换变量

最近重新看了看C和C++,觉得有些地方挺有意思。作为一开始不管什么资料都会用来做例子的一个程序,交换变量。不管在哪,常用的int,float,double类型的变量都是值类型的,作为参数传到函数(方法)中时,其实是复制了一个值进去,也就是说通过函数是无法直接更改这些值的,只能通过一些间接的方法来更改...

C语言 异或运算符 ^

C语言中异或运算符^表示参见运算的二进制运算符相同为0,不同为1,如下:1 ^ 1 = 01 ^ 0 = 10 ^ 1 = 10 ^ 0 = 0下面举例说明∧运算符的应用:  (1)使特定位翻转  假设有01111010,想使其低4位翻转,即1变为0,0变为1。可以将它与00001111进行∧运算,...