参考:将原GeoTIFF数据的投影坐标信息赋值到新创建的文件上 pro tiff_projection ;启动ENVI e = ENVI(/HEADLESS) ;打开文件 file = 'D:\01-Working\2017\20171204-IDL平均计算\GOCI_0510\GOCI_0510_01_GDPS_RC_GEO_ZJ_MASK_BIL_FLAASH_TSM.tif' raster = e.OpenRaster(file) ;读取影像的第一个波段 data = raster.Get

参考:Create Latitude/Longitude Arrays for GeoTIFF Image 用到的函数为 READ_TIFF,通过此函数可以获取 TIFF 数据的数组信息,同时可以获取地理坐标信息,包括最大纬度值.最小经度值.一个象元所代表的经纬度值,这些信息可以通过 ArcGIS 软件获取对比. 01   READ_TFF 读取单通道或多通道 TIFF 格式数据,返回图像数组及相关参数. 02   WRITE_TIFF 写入 TIFF 文件. 03   QUERY_TIFF 查

前言 在上两篇文章中我介绍了如何直接将Geotiff(一个或者多个)发布为TMS服务.这中间其实我遇到了一个问题,并且这个问题伴随Geotrellis的几乎所有使用案例,下面我进行详细讲述. 一.问题描述 无论在将TIff文件使用Geotrellis导入Accumulo中还是直接将其发布为TMS服务,其实这中间都存在一个问题:当多个Tiff文件存在重叠部分的时候如何接边以及在接边处的瓦片如何取出各Tiff文件中涉及到的数据. 这个问题可以说是一个问题也可以说是两个问题.当我们采用事先导入Accu

;+ ;:Description: ; Describe the procedure. ; ; Author: DYQ 2009-7-19; ; ;- PRO TEST_READHDF COMPILE_OPT idl2 file = FILE_DIRNAME(ROUTINE_FILEPATH('TEST_READHDF')) + $ '\data\MOD04_L2.A2004196.0140.005.2007023163311.hdf' ; IsHdf = HDF_ISHDF(file) IF

使用mrdfits函数 这是天文学标准库中的函数,下载地址:https://idlastro.gsfc.nasa.gov/homepage.html,下载后,将pro文件夹导入到IDL工程中. str = MRDFITS(file_or_unit, exten_no_or_name, [ header] )  file_or_unit:读取的文件名 exten_no_or_name:is either the extension number to be read (0 for the prim

影像拼接工具使用说明 一.影像像素位深检查 1.采用开源库GDAL的gdalinfo.exe读取GeoTiff文件的信息,如类型.投影,范围等. 2.采用DOS脚本遍历所有GeoTiff文件,输出各个GeoTiff文件的信息到日志. 3.统计GeoTiff文件个数,Type=Byte的文件的个数.数目一致表明全部GeoTiff像素类型一致. 运行示例: for /f %i in ('dir /b /s d:\tifdata\*.tif ') do @( e:\gdal1.4.2\bin\gdal

1.打开文件 IDL从磁盘上的文件读写数据,必须首先把一盒逻辑设备号连接到一个指定的文件,然后进行文件操作,如打开,关闭和读取等.IDL中的逻辑设备号的范围是-2——128,其中1-99是用户可以任意指定的,100-128是根据get_lun和free_lun来管理的,其他事特殊函数专用的. 0是标准输入流,一般是键盘,也就是说下面命令是一只的. read,x read,,x -1是标准的输出流,一般显示屏幕,下面代码是一样的 printf,-,x print,x -2是错误流,一般显示屏幕 文

前言 传统上我们需要先将Tiff中存储的影像等数据先切割成瓦片,而后再对外提供服务.这样的好处是服务器响应快,典型的用空间来换时间的操作.然而这样造成的问题是空间的巨大浪费,一般情况下均需要存储1-18级左右的瓦片数据.我一直在思考有没有办法不存储瓦片而直接发布TMS服务,当然这样响应速度肯定是要受一点影响,但是基于Geotrellis的分布式计算对这一点提供了巨大帮助,大大缩短了瓦片临时切割(存储于内存中)所用的时间.而且这样不仅仅是节省了存储空间的问题,何况我们有时可能只是为了查看数据情况(

1.下标方式 IDL> array=indgen(8) IDL> print,array 0 1 2 3 4 5 6 7 IDL> print,array[3] IDL> array=indgen(8) IDL> print,array[-1] 7 IDL> print,array[-5:-1] 3 4 5 6 7 2.向量方式 IDL> ;读取数组中第一,二,四和第六个元素的代码 IDL> array=indgen(8) IDL> indices=[

前言 上一篇中简单介绍了 COG 的概念和 Geotrellis 中引入 COG 的原因及简单的原理,本文为大家介绍如何在 Geotrellis 中使用 COG 来写入和读取 GeoTIFF数据. 一.写入数据--ETL 1.1 实现方案 其实这与之前的普通 ETL 操作在概念上是相似的,都是将原始数据转换成系统能用的数据的过程,这是宽泛的 ETL 的定义.在 Geotrellis 中实现很简单,与之前代码基本一致,只要切换一下 writer 类型以及最后建立金字塔额时候略有不同.实现方案如下:

参考:I/O - NetCDF Routines —— NetCDF 操作函数 01   NCDF_OPEN 打开一个 NetCDF 文件. 02   NCDF_CLOSE 关闭一个 NetCDF 文件. 03   NCDF_VARID 返回一个变量 ID,varid. 04   NCDF_VARGET 读取一个变量. 05   NCDF_ATTGET 读取一个属性. 06   NCDF_ATTINQ 返回一个关于 NetCDF 属性信息的结构体. 07   NCDF_ATTRENAME 属性重

参考:I/O - HDF Routines —— HDF 操作函数 01   HDF_SD_START 打开一个 SDS 模式的 HDF 文件. 02   HDF_SD_END 关闭一个 SDS 模式的 HDF 文件. 03   HDF_SD_NAMETOINDEX 返回变量索引. 04   HDF_SD_SELECT 返回变量标识符. 05   HDF_SD_GETDATA 读取变量数据. 06   HDF_SD_ENDACCESS 结束到一个变量的通道. 07   HDF_SD_ATTRFI

HDF(Hierarchical Data Formats)数据格式由 NCSA 开发.HDF 提供了大量的数据模式,包括多维数组.表格.图像.注解和调色板.在下面的章节中,将描述 HDF 科学数据系列(SDS)的数据模式,因为它是 HDF 中最具灵活性的,并且它和 NetCDF 具有相似性.也就是说,HDF SDS 的基本组成也是变量.属性和维数. 注意:IDL读取数据与在其他软件上显示的数据位置相反,左上对右下! 注意:HDF记录数据与实际数据是通过一个数量关系获取的,如下图所示: 实际结果

NetCDF(network Common Data Form)由位于科罗拉多州波尔市的 Unidata 程序中心开发,主要应用于大气科学的研究.NetCDF 的数据模式具有简单性和灵活性的特点.NetCDF 文件的基本组成为变量.属性和维数: 变量为标量或多维数组.NetCDF 所支持的 IDL 数据类型有 string.byte.int.long.float 和 double 属性包含一个变量或这个文件的附加性质.包含变量信息(如单位.有效范围.尺度因子等)的属性称作变量属性:包含文件信息的

  源码: 1: PRO TIFF_NDVI,F1,F2,FOUT 2: F1 = DIALOG_PICKFILE(TITLE = 'B4 TIFF',FILTER='*.TIF',/READ) 3: F2 = DIALOG_PICKFILE(TITLE = 'B3 TIFF',FILTER='*.TIF',/READ) 4: FOUT =DIALOG_PICKFILE(TITLE = 'RESULT PATH',FILTER='*.TIF',/WRITE) 5: IF(FILE_TEST(F1

IDL Reference 操作符号 数组 字符及字符串 结构体 指针 链表 & 哈希表 程序控制(循环.条件.跳转语句) 过程 & 函数 输入与输出 系统变量 文件系统操作 直接图形系统 参考:IDL部分系统函数(全) 1. 操作符号 Relational Operators:[比较运算符]EQ(=), GT(>), LT(<), GE(≥), LE(≤), NE(!=). Mathematical Operators:[数学运算符]+.-.*./.++.--.MOD Mat

Geotrellis系列文章链接地址http://www.cnblogs.com/shoufengwei/p/5619419.html 目录 前言 需求说明 实现方案 总结 一.前言        最近一段时间比较忙,没能继续推进Geotrellis项目开发,周末和这两天抽空又实现了一个功能--导出自定义的Tiff文件.又恰巧碰上今天这么重要的日子,当然要写点东西来纪念一下,所以就有了这篇文章,闲话莫说,进入正题. 二.需求说明        很多时候我们需要从一块(或者很多块)大的Tiff中根

参考: 01   IMAGE 将图像数据以图形窗体的形式显示. 02   COLORBAR 在已经存在的IDL图形中增加一个colorbar或创建. 03   MAPGRID 在已经存在的IDL地图图形上面绘制地图网格. 序号 名称   功能说明   语法 & 举例 01 IMAGE   将图像数据以图形窗体的形式显示.-----------------------------------------------graphic = IMAGE( Data [, X, Y ] [, Keyword

C# (C Sharp)是微软公司在2000年6月发布的一种新的编程语言.C#与Java有很多的相似之处,包括了诸如单一继承.界面.与Java几乎同样的语法,和编译成中间代码再运行的过程.它又借鉴了Delphi的一个特点,与COM(组件对象模型)是直接集成的,而且它是微软公司.NET windows网络框架的主角. IDL则一直是应用程序开发和科学家进行可视化与分析的首选语言.因为它功能强大,简单易学,很少的几行代码就能实现其他语言很难实现的功能,所以它是进行科学数据分析.可视化表达和跨平台应用

先列后行 arr=indgen(3,4) SIZE(arr,/dimensions) print ,where(arr gt7) print,arr[where(arr gt 7)] print,where(arr gt 7,cout,COMPLEMENT=com,ncomplement=ncom) array[WHERE(array GT 5, /NULL)] = 5 arr=indgen(2,8); arrReform=reform(arr,1,8);维数一定,改变个数 arrRebin=R