write geotiff

From DavinciWiki
Jump to: navigation, search

Description

This function writes a GeoTIFF file. The input data must be a structure from read_geo()

This function relies on GDAL > 1.8.

Arguments and Return Values

Arguments: A Geo-Struct and file name

Return Value: none

Usage

Syntax: write_geotiff(Geo-Struct, filename [,force=BOOL][,subdata=INT][,a_srs=STRING])

Geo-Struct - the davinci Geo-Struct

filename - full path to data

force - force the creation/overwrite of the specified file

a_srs - source projection argument (either a proj string or wkt), assigned by default values in Geo-Struct

subdata - write out only the subdataset specified. Currently we only support the writing of 1 subdataset.

Examples

dv> a=read_geo($DV_EX+"/I01001001.cub")
Downloading...........................................done.
   anc: struct, 11 elements
       wkt: Text Buffer with 12 lines of text
           1: PROJCS["SIMPLE_CYLINDRICAL MARS",
           2:     GEOGCS["GCS_MARS",
           3:         DATUM["D_MARS",
           4:             SPHEROID["MARS",3396190,0]],
           5:         PRIMEM["Reference_Meridian",0],
           6:         UNIT["degree",0.0174532925199433]],
           7:     PROJECTION["Equirectangular"],
           8:     PARAMETER["latitude_of_origin",0],
           9:     PARAMETER["central_meridian",0],
           10:     PARAMETER["standard_parallel_1",0],
       proj4: "'+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=3396190 +b=3396190 +units=m +no_defs '"
       origin: 1x2x1 array of double, bsq format [16 bytes]
       pixel_size: 1x2x1 array of double, bsq format [16 bytes]
       x: 186699.0000
       y: -3541.000000
       lat: 1x5x1 array of float, bsq format [20 bytes]
       lon: 1x5x1 array of float, bsq format [20 bytes]
       lpo: 1x5x1 array of float, bsq format [20 bytes]
       spo: 1x5x1 array of float, bsq format [20 bytes]
       center: 2x1x1 array of float, bsq format [8 bytes]
   ignore: -3.402822655e+38
   projection: "Equirectangular"
   spheroid: struct, 3 elements
       name: "MARS"
       radius: "3396190"
       iflattening: "0"
   resolution: 100.0000000
   x: 186699.0000
   y: -3541.000000
   ul: 2x1x1 array of float, bsq format [8 bytes]
   ur: 2x1x1 array of float, bsq format [8 bytes]
   ll: 2x1x1 array of float, bsq format [8 bytes]
   lr: 2x1x1 array of float, bsq format [8 bytes]
   w: 661
   h: 2391
   data: 661x2391x10 array of float, bsq format [63,218,040 bytes]

dv> write_geotiff(a,"test.tiff",force=1)
Writing /tmp/dv_3044/test.tiff.tiff: 1503x1561x1 TIFF file.

Attaching Ancellary Data
Writing test.tiff: 1503x1561x1 GeoTIFF File.
Done.

DavinciWiki Mini-Nav Bar

Contents


Contact Developers

  • davinci-dev [AT] mars.asu.edu

All other topics

  • See navigation on the left

Major Sub-Functions

Related Functions

Recent Library Changes

Created On: 04-11-2013
Modified On: 05-14-2013

Personal tools