2013年01月_箜_Kong

12月 11月 10月 09月 08月 07月 06月 05月 04月 03月 02月 01月

原创已知两点经纬度求球面最短距离的公式推导过程（几何法加向量法）

已知两点经纬度计算球面距离的公式，一搜一大堆，形式如下：可是至于这个公式为什么是这样的，今天推导了一下，详细推导过程如下。首先画个图（图1），要不然空间想象能力差的话容易犯糊涂。首先对图1做个大致的说明，红色的半圆表示赤道，蓝色的圆弧表示本初子午线（也就是经度为0的子午线）。球最上方是北极点，点A和点B分别为要计算的两个点，坐标分别为A（jA,wA）和B（jB，wB）。图1 示意图再开始推导之

2013-01-29 20:42:39 47705 18

翻译使用Mapnik生成地形图——thematicmapping.org译文（四）

原文地址：http://blog.thematicmapping.org/2012/07/terrain-mapping-with-mapnik.html作者：作者：Bj?rnSandvik（Google Plus：https://plus.google.com/118196887774002693676）在之前的三篇博文中，我们使用DEM数据分别创建了三个GroTiff的数据，分别是：山体阴影

2013-01-28 22:56:15 6578 3

翻译使用gdaldem创建彩色地形图和坡度阴影——thematicmapping.org译文（三）

原文地址：http://blog.thematicmapping.org/2012/06/creating-color-relief-and-slope-shading.html作者：Bj?rnSandvik（Google Plus：https://plus.google.com/118196887774002693676）在我之前的博客中，我们使用gdaldem传教的是灰度图像（山体阴影）。今

2013-01-28 21:54:08 7250 1

翻译使用gdaldem生成山体阴影——thematicmapping.org译文（二）

开始之前，先来几张DEM数据的处理效果图原文地址：http://blog.thematicmapping.org/2012/06/creating-hillshades-with-gdaldem.html 在上一篇中，我们创建了挪威Jotunheimen山区的DEM数据。接下来我们将使用这个DEM生成山体阴影和彩色地形图，这是一种很流行的地形图可视化技术，通过调节光线和阴影来进行地图渲染。 GD

2013-01-23 22:27:01 5426 1

翻译数字高程模型和地图——thematicmapping.org译文（一）

作者：BjørnSandvik（Google Plus：https://plus.google.com/118196887774002693676）原文地址：http://blog.thematicmapping.org/2012/06/digital-terrain-modeling-and-mapping.html 这是新博客的第一篇，该系列主要关注于数字地形数据的使用和处理。数字地形数据与

2013-01-22 23:30:58 5662 3

转载 Metamaterials perform image compression before light reaches the sensor（微波成像，拍摄与压缩同步进行）

转自果壳网上，地址：http://www.guokr.com/article/436242/ 英文原文地址：http://arstechnica.com/science/2013/01/metamaterials-perform-image-compression-before-light-reaches-the-sensor/ 异向介质材料（metamaterials）具有自然界中的材料所不具备

2013-01-22 20:50:08 1828

翻译使用阿富汗和巴基斯坦地区的SRTM数据生成山体阴影和彩色地形图

使用阿富汗和巴基斯坦地区的SRTM数据生成山体阴影和彩色地形图数据来自GIST，工具来自GDAL 原文地址：http://developmentseed.org/blog/2009/jul/01/generating-custom-map-tiles-rapidly-cloud/ 在阿富汗和巴基斯坦这样多山的地方，高程数据是非常重要的。GIST项目是使用航天飞机雷达地形测绘任务（STRM）获取各个

2013-01-21 22:10:23 7290 15

翻译使用GDAL对DEM渲染处理流程

下面是翻译的一篇国外的博客，没有严格按照原文的字面意思翻译，是按照我自己的理解来进行翻译的，如果有错误，还请大家指正，原文地址是：http://linfiniti.com/2010/12/a-workflow-for-creating-beautiful-relief-shaded-dems-using-gdal/。有时候我使用QGIS的Hillshade插件来生成山体阴影数据，然后和原来的DE

2013-01-20 20:05:14 12274

原创 GDAL使用DEM数据计算山体阴影（Hillshade）

零、前言说起Hillshade，其实就是模拟太阳光照射地形所引起的明暗对比，然后来对地形图进行渲染，使之看起来具有立体效果的一种方式，常用于地图的渲染，如表1所示，具体的可以参考文献[1]，表1中的图均来自参考文献[1]。表1 DEM、山体阴影以及应用对比 DEM图像（使用颜色渲染）从左图的DEM图像中计算的山体阴影图 Paper Map Withou

2013-01-17 23:11:14 9947 4

原创 GDAL使用DEM数据计算地形指数

零、前言本文是接上文《GDAL使用DEM数据计算坡度坡向》，还是一篇关于DEM计算地形指数的一篇文章。这里所要计算的地形指数主要包括以下三个指数：地形耐用指数（Terrain Ruggedness Index，TPI）、地形位置指数（Topographic Position Index，TRI）和粗糙度（roughness）。上面三个地形指数都是在一个3×3的窗口中计算的，3×

2013-01-15 22:54:26 15324 4

原创 GDAL使用DEM数据计算坡度坡向

零、前言之前写过一个3×3的通用模板算子函数的博客《基于GDAL的一个通用的3×3模板函数》，网址：http://blog.csdn.net/liminlu0314/article/details/8316156。当时说是要基于这个函数写一个计算坡度坡向的函数。由于这段时间一直忙于别的事情，这件事情就拖着了，今天给大家补上。一、简介坡度（slope）是地表单元陡

2013-01-13 22:59:20 33984 42

原创 GDAL计算栅格图像统计值相关说明

一、简介在图像处理的过程中，会经常获取图像的统计值来进行一些计算。这里说的图像的统计值指：直方图，最大值，最小值，均值和方差这几个。下面分别对怎么使用GDAL来计算这些值做一个简单的说明，希望对大家有用。下面将分为三个部分进行说明，分别是：直方图统计、极值统计、均值标准方差统计。二、统计直方图关于直方图，这里不多介绍。这里主要对如何使用GDAL获取直方图进行说

2013-01-11 23:15:00 12655 21

原创使用GDAL生产jpg格式的快视图（拇指图）

使用GDAL生成快视图主要有两个步骤，第一重采样生产小图像；第二将生成的小图像CreateCopy生成jpg格式或者png格式。由于GDAL不支持直接创建Jpg或者png的压缩格式的图像，所以中间需要一个中转的过程。下面简单说明一下过程，代码就不贴了。首先通过获取采样比例计算输出图像的大小，这里可能需要根据需求自己定义。我遇到了一般有两种情况：第一种，将图像直接缩小到原来的10分之一（或者其他的

2013-01-08 22:15:11 7046 10

原创使用GDAL获取HDF等数据集中的图像

使用GDAL读取HDF，NetCDF等数据集中的数据时，一般需要两个步骤，第一，获取数据集中的子数据集；第二，通过第一步获取的子数据集读取图像数据。一般的HDF图像中会有很多个子数据集，比如常用的MODIS数据，使用Envi打开会弹出下面的对话框来让用户选择需要打开的子数据集（如图1）。图1 Envi打开Modis数据从图1中可以看出，Envi是把所有的子数据集的波段都进行了列举，不过这点和GD

2013-01-07 20:42:30 16145 27

FY2系列静止卫星圆盘数据校正经纬度查找表

2017-07-21

GDAL2.0.2 C# VS2013 X64 修改中文乱码

使用VS2013编译的GDAL2.0.2 C#的64位版本，修改了中文路径、中文图层名、中文字段名、中文属性值乱码问题

2017-01-06

GDAL1.11 X86 X64

GDAL1.11.0版本，包含X86和X64两个平台。具体修改信息如下： 1、增加OpenCL支持，可以使用GPU进行处理 2、所有的依赖第三方库全部编译为插件形式，需要的话，直接放进去，不需要就删掉 3、修改RPC模型支持像方仿射变化模型 4、增加CNSDTF的格网grd格式和CNSDTF的矢量vct格式（没经过严格测试）有问题请发邮件，或者CSDN私信，微博私信

2014-05-05

GDAL110 全套 C\C++ C# Java Python

1：修改C#版本的dll，可以支持所有的中文路径 2：增加其他第三方库支持（GEOS，PROJ4，FileGDB、HDF4，HDF5，NetCDf，CURL，EXPAT，ICONV，ECW，OCI，MrSID，JPEG2000，SQLite3，xerces-c等）使用VS2008SP1编译，如果不能执行请安装VS2008SP1的分发库。 C++文件夹下为常用C和C++库以及dll文件 C#文件夹下为C#使用的八个dll和GDAL自带的C#程序 Python文件夹下为Python接口 Java文件夹下为Java库所有语言的使用都需要把C++文件夹下所有的dll拷贝至自己的程序所在目录，否则会出现不能运行的情况对于C#程序和C++程序，为了能够正常访问中文路径，请将GDAL_FILENAME_IS_UTF8设置为NO。代码分别为： C++： CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO"); C# : OSGeo.GDAL.Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO"); 如果程序不能执行，或者提示0xc1000052错误，请到这里下载VS2008 SP1的运行库，地址为： http://www.microsoft.com/zh-CN/download/details.aspx?id=5582

2013-07-24

GDAL110版本对应的C#完美支持中文路径

GDAL110版本修改支持中文路径后的C#版本，需要与原来的GDAL110的C++库支持，C++库下载地址为： http://download.csdn.net/detail/liminlu0314/5748203

2013-07-24

GDAL 1.10全套增加FileGDB支持

使用GDAL1.10正式版本编译，编译环境VS2008 SP1 内部将FileGDB、HDF4、HDF5、GEOS、PROJ4、ECW、Jpg2000、libexpt、libcurl、libiconv、netcdf、xerces_c、sqlite3等常用的库全部编译进去语言 C\C++、C#、Java、Python 含有帮助文档如果程序不能执行，或者提示0xc1000052错误，请到这里下载VS2008 SP1的运行库，地址为： http://www.microsoft.com/zh-CN/download/details.aspx?id=5582

2013-07-13

GDAL 1.10开发库（VS2008SP1编译）

1：修改默认文件路径为非UTF8编码，可以直接打开中文路径数据 2：修改增加OPENCL支持，对于warp操作可以使用GPU运算 3：增加其他第三方库支持（HDF4，HDF5，NetCDf，GEOS，PROJ4，CURL，EXPAT，ICONV，ECW，OCI，MrSID，JPEG2000，SQLite3，xerces-c等）使用VS2008SP1编译，如果不能执行请安装VS2008SP1的分发库。 C++文件夹下为常用C和C++库以及dll文件 C#文件夹下为C#使用的八个dll和GDAL自带的C#程序 Python文件夹下为Python接口 Java文件夹下为Java库所有语言的使用都需要把C++文件夹下所有的dll拷贝至自己的程序所在目录，否则会出现不能运行的情况

2013-05-10

GDAL1.92 修改S57中文乱码

6：修改读取S57海图数据中文属性值乱码 e:\3rdPart\gdal-1.9.2\ogr\ogrsf_frmts\s57\s57reader.cpp line 883 // Add By liml 2013-04-25 Convert UCS-2 to Utf-8 const char* pachBuffer = poRecord->GetStringSubfield("NATF",0,"ATVL",iAttr); if( EQUAL(pszAcronym, "NOBJNM")) { size_t nLength = strlen(pachBuffer); char* pszValue = new char[nLength]; memcpy(pszValue, pachBuffer, nLength); pszValue[nLength-1] = '\0'; char *pszGetUTF8 = CPLRecodeFromWChar( (const wchar_t*)pszValue, CPL_ENC_UCS2, CPL_ENC_UTF8); pachBuffer = CPLRecode( pszGetUTF8, CPL_ENC_UTF8, CPL_ENC_LOCALE); delete []pszValue; } // Add By liml 2013-04-25 poFeature->SetField( pszAcronym, pachBuffer);

2013-04-25

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

GDAL专栏