缘起 对于一般的遥感影像文件,金字塔文件默认都是与影像文件放在同一个目录下,金字塔文件名一般与源影像文件名相同,但后缀名不同.或者金字塔内建于影像内部,但这不是这里所涉及的. 在使用ArcGIS桌面版或者Erdas遥感影像处理软件打开遥感影像文件的时候,如果影像不含(带有)金字塔,则会提示是否创建金字塔. 我们碰到这么一种情况,就是原始影像是不带金字塔的,并且所在的目录不允许修改,也就是不能将金字塔创建在其中,原始影像文件更是不允许修改,所以也不能更新内建金字塔. 解决 对于这个问题,为了加速影

#!/usr/bin/env python # -*- coding: utf-8 -*- import sys, gdal, os from gdalconst import GA_ReadOnly from os.path import realpath from shapely.geometry import LineString #根据坐标点计算该点所在栅格上的值 def get_elevation(x_coord, y_coord, raster, bands, geo_trans):

gdalwarp输出tif图像的时候,默认如果没有使用BIGTIFF=YES选项,则会根据输出影像的大小进行判断,低于4G则不适用bigtiff格式. 对于非bigtiff图像,如果这时候使用gdaladdo创建金字塔(内建模式),如果出现了文件总大小超过4G的情况,则会导致图像文件出错. gdaladdo会输出相关信息,但不会终止金字塔的创建,于是导致金字塔是黑色的. C:\Program Files\GDAL>gdaladdo D:\data\fs\ID_Guangdong.shp.tif

主要转自http://blog.csdn.net/liminlu0314?viewmode=contents 部分GDAL工具功能简介 gdalinfo.exe 显示GDAL支持的各种栅格文件的信息. gdal_translate.exe 在不同的格式间进行转换.同时,潜在的执行了一些切割.重采样和使像素比例变化的任务. gdalwarp.exe 投影转换和投影绑定.同时也可以进行图像镶嵌.这个程序可以重新投影所支持的投影,而且如果图像("raw" with)控制信息也可以把GCPs

开工之前要先了解一下瓦片地图,瓦片地图金字塔模型是一种多分辨率层次模型,从瓦片金字塔的底层到顶层,分辨率越来越低,但表示的地理范围不变.实现原理就是,首先确定地图服务平台所要提供的缩放级别的数量N,把缩放级别最低.地图比例尺最大的地图图片作为金字塔的底层,即第0层,并对其进行分块,从地图图片的左上角开始,从左至右.从上到下进行切割,分割成相同大小(比如256x256像素)的正方形地图瓦片,形成第0层瓦片矩阵;在第O层地图图片的基础上,按每2x2像素合成为一个像素的方法生成第1层地图图片,并对其进

GDAL原生支持超过100种栅格数据类型,涵盖所有主流GIS与RS数据格式,包括•  ArcInfo grids, ArcSDE raster, Imagine, Idrisi, ENVI, GRASS, GeoTIFF •  HDF4, HDF5•  USGS DOQ, USGS DEM •  ECW, MrSID •  TIFF, JPEG, JPEG2000, PNG, GIF, BMP 完整的支持列表可以参考http://www.gdal.org/formats_list.html 导入

读写影像可以说是图像处理最基础的一步.关于使用GDAL读写影像,平时也在网上查了很多资料,就想结合自己的使用心得,做做简单的总结. 在这里写一个例子:裁剪lena图像的某部分内容,将其放入到新创建的.tif文.以此来说明GDAL读写影像的具体实现. 1.打开图像 用GDAL打开lena.bmp,实现如下.注意这里打开图像,指的是获取图像的头文件,以此得到图像的一些信息,没有涉及到读取像素操作. GDALAllRegister(); //GDAL所有操作都需要先注册格式 CPLSetConfigO

初次使用visual studio 对工具不熟悉,有些步骤可以配置的 (1)在VS2005中新建win32控制台程序 testGDALconsole,(向导中附加选项不能选为空项目) 将C:\GDAL\bin\gdal14.dll拷贝到testGDALconsole目录下的debug目录中.(否则运行时会提示找不到gdal14.dll) (2)在工程的Library files(库目录)中和Include files(包含目录)中分别添加GDAL的LIB文件目录和头文件目录, (3)配置附加依赖

在遥感图像处理中,GDAL库不仅能读取和处理大部分的遥感图像数据,而且还能够实现图像处理后将数据保存为图像的功能. 本文就详细介绍如何将内存中的图像数据保存为.tif格式. 首先,遥感数据处理完,保存在一维数组data中,图像的宽和高为width和height,图像有三个波段.在保存时要使图像的每一行对其,保证位数为32的倍数 先上实现的代码,源程序如下: +)/;//字节对齐 unsigned char *data=new unsigned char[bytesPerLine*height];

以下是设定一个矩形框,用IPixelBlock将256*256瓦片tile拼接成一个整块影像的代码,row1, col1, row2, col2是一个矩形框行列号范围.level是瓦片的金字塔等级.这里的瓦片已经下载完毕,位于domSavePath文件夹下. //选择的Google瓦块的行列号范围 int row1, col1, row2, col2; ; row1 = topLeft.Row; col1 = topLeft.Col; row2 = bottomRight.Row; col2 =

部分源码选自GDAL库的官方网址:www.gdal.org,其余的代码为笔者自己编写. // readfile.cpp : 定义控制台应用程序的入口点. // /* part of the codes were cite from: http://www.gdal.org/gdal_tutorial.html and remaining of code were created :by www.cnblogs.com/AmatVictorialCuram/‎ and please mark th

GDAL is a useful command line tool to process spatial data, if you haven’t heard of the tool before some examples of what it can do are: Create contours from a DEM Create a TMS tile structure Rasterize vector into a raster file Build a quick mosaic f

         为了能快速的显示大影像,最近一直在学习GDAL,GDAL确实是一个功能强大的开源库,其核心部分数据集和波段,下面这个图很详细的描述了它们之间的关系,还有其中的细节:     GDAL的RASTERIO功能非常强大,简短的一句话就能实现图像的显示,但也是这个简单函数,把我折腾的半死.在学习使用GDAL的过程中,非常感谢李林大哥和貟建明大哥,他们不厌其烦地解答我的疑问,提示关键性要点,使我在解决问题的过程中事半功倍.      我现在要把最近学习过程中的心得写下来,留给自己以后看看

作者:朱金灿 来源:http://blog.csdn.net/clever101 最近使用gdal库比较多,就谈谈gdal库的一些使用心得. 第一个是GDALOpen的访问权限参数会影响图像的创建金字塔方式.比如你是这样打开图像和创建金字塔: std::string strImgPath = _T(“C:\\1.tif”); GDALDataset* mGdalDataset=(GDALDataset*)(GDALOpen(strImgPath.c_str(),GA_Update)); mGda

文章版权由作者李晓晖和博客园共有,若转载请于明显处标明出处:http://www.cnblogs.com/naaoveGIS/ 1.背景 为利于项目实施,团队用AE写过一个插件式的工具集,主要包括了数据入库.数据统计.数据修改.数据发布等功能. 由于arcgis环境的不同,数据库的环境不同,导致工具的版本出现了arcgis10.0,10.1,10.2,以及相对的Oracle版本和mysql版本,从而使维护量直线上升. 痛定思痛,我们决定采用GDAL来将各工具逐一开源化,并且将不同数据库操作封装成

Landsat系列卫星提供的数据,一般都是每个波段一个tif文件,然后外加一个MTL.txt的元数据文件,使用gdal可以直接打开每个波段的tif文件,但是有时候想在打开tif数据的同时能够自动读取MTL文件里面的元数据信息,这个时候就只能自己再解析这个文件了.比较麻烦,下面就是针对这种情况,直接在gdal库里面扩展一种支持landsat的mtl的格式,直接打开mtl文件的同时自动打开tif数据以及元数据. 关于扩展GDAL库可以参考我之前写的CNSDTF格式扩展,以及gdal的官网说明,网址为

0.目的 GDAL库中提供了四五种读取Jpeg2000的驱动,但是各个驱动读取数据的效率各不相同,下面就针对三种读取jpeg2000的效率进行测试. GDAL库中提供的读取Jpeg2000的驱动有下面六个,具体的说明请移步GDAL官网查看.地址为:http://www.gdal.org/formats_list.html. JPEG2000 - Implementation of the JPEG-2000 part 1 JP2ECW – ERDAS JPEG2000 (.jp2) JP2KAK

在上一篇博客中,讲了一下使用GDAL从文件流中构造一个GDAL可以识别的数据来进行处理.原以为这个接口在C#中没有,仔细看了下GDAL库中源码,发现C#版本也有类似的函数,下面是GDAL库中的一个C#使用示例代码. using System; using System.IO; using System.Runtime.InteropServices; using OSGeo.GDAL; /// <summary> /// 一个基于C#语言的简单示例,使用内存文件 /// </summar

扩展GDAL,支持CNSDTF格式(一) 一.        简介 本文主要根据<中华人民共和国国家标准GB/T17798-2007--地理空间数据交换格式(Geospatialdata transfer format)>(http://www.bzxz.net/bzxz/129804.html)中定义的格网数据和矢量数据格式内容,对GDAL库进行扩充,使之可以直接打开这两种格式. 此外GDAL格式扩展就参考GDAL官方文档中的扩展方式.网址为:www.gdal.org/gdal_driver

GDAL——命令使用专题——gdalinfo命令  前言 GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库.它利用抽象数据模型来表达所支持的各种文件格式.它还有一系列命令行工具来进行数据转换和处理. 主要针对:ogrinfo命令.gdalinfo命令.gdalsrsinfo命令.gdallocationinfo命令.gdal_translate命令.gdaltransform命令 等.这一节讨论gdalinfo