geotiffwrite

Write a 2D or 3D array to a single or multi-band GeoTIFF file

Description

This program can write a 2D or 3D array to a single or multi-band GeoTIFF file. For more information and examples, please visit:
http://research.baycrest.org/~jimmy/geotiffwrite
or type: help geotiffwrite

MATLAB's Mapping Toolbox only provides a "geotiffread" function, but it does not have a "geotiffwrite" function (Note). This is the MATLAB program to write a 2D or 3D array to a single or multi-band GeoTIFF file, where data can be either 1-bit monochrome data (i.e. binary or bilevel), 8-bit unsigned integer, 16-bit signed / unsigned integer, 32-bit signed integer, or 32-bit floating point.

Note: Starting from version R2011a, MATLAB also provides its own "geotiffwrite" function
(http://www.mathworks.com/help/toolbox/map/rn/bsq4us7-1.html#bsu1ro5-1).

This program is based on GeoTIFF Format Specification under:
http://www.remotesensing.org/geotiff/spec/geotiffhome.html
or http://download.osgeo.org/geotiff/spec

It does not need MATLAB's Mapping Toolbox, or any other library.

Usage:
[option bbox] = make_option([bbox]);
geotiffwrite(filename, bbox, image, [bit_depth, option]);

For more information on arguments and detail examples, please visit:
http://research.baycrest.org/~jimmy/geotiffwrite
or type: help geotiffwrite

Acknowledgements

This file inspired Utm2 Lat Long Geo Tiff Converter.

• MATLAB release: MATLAB 5.3.1 (R11.1)

[Sunke]

thanks for this export option! It helps sharing my MatLab data with colleagues that do not use MatLab, but ArcGIS, a lot !

[Kyle Mayer]

 

[Jesse Angle]

In order to make this code fully support the ModelTransformationTag the following changes are necessary:

line 873: comment out
error('ModelPixelScale is required...');

as this isn't necessarily a true statement and will cause problems when trying to use only the ModelTransformationTag.

line 1213 - 1214:
fwrite(fid, gifd.ModelPixelScaleTag,'double'); and
fwrite(fid, gifd.ModelTiePointTag, 'double');

should be surrounded by
if isfield(gifd,'Model...Tag')

end

statements for when these tags are not being used as is the case with ModelTransformationTag.

Finally, add in the code that actually writes the ModelTransformationTag portion (previously missing). Mine looks like this:

if isfield(gifd,'ModelPixelScaleTag')
fwrite(fid,gifd.ModelPixelScaleTag,'double');
end
if isfield(gifd,'ModelTiePointTag')
fwrite(fid,gifd.ModelTiePointTag,'double');
end
if isfield(gifd,'ModelTransformationTag')
fwrite(fid,gifd.ModelTransformationTag,'double')
end

That being said, great work Jimmy getting this put together and thank you for your efforts.

[Carlo M]

thank you very much Jimmy

[Jimmy Shen]

You need to use MATLAB's geotiffread to read L7 file into matrix, with all attributes. Then find out the index of your sub-region, and assign the sub-region to a new matrix. You can now use the make_option option in my package to select optional fields. However, you need to match the attributes that you read, and select the correct fields by yourself in order to create the optional fields to be used in my geotiffwrite. Since I don't have Mapping toolbox, I am afraid that I cannot provide you with more details.

[Carlo M]

thank you for your response... i'm trying to select a region of a Landsat 7 geotiff image and then to save it with all its atributes. Please, how can I select the rectangular window?

[Jimmy Shen]

My geotiffwrite can only write a rectangular dataset into the file. As long as your sub-region is a rectangular dataset, you should be able to write it into a separate file.

[Carlo M]

there is possible to select a sub-region of a geotiff image? i'm trying to do that but it seems not possible

[Jimmy Shen]

MATLAB finally has its own "geotiffwrite" function. So I suggest you use its geotiffread / geotiffwrite pair to load your Antarctica background, and superimpose your ship track data on top, and then re-write the entire data into GeoTIFF format using Matlab's geotiffwrite function.

I am currently no longer working in GeoScience field, so I am afraid that I cannot provide you with more useful advise.

[Kaycee]

Hi Jimmy, I found your post very informative, so I was hoping you could help me with a problem. I'm having issues attempting to plot a ship track on top of an already made geotiff image from Antarctica. The issue i believe is that it is stereographic, so i cannot use geotiff2mstruct. I am able to get the track on the ice image using profwd however the track is not exactly in the right spot. Also I cannot figure out how to get the latitude and longitude on the axes, the file has intrinsic or world x,y. Any advice?

[Evan]

Thanks for the fast response. I changed the colormap to include 256 elements and got it working. Seems to be a very useful script although it requires some troubleshooting with the 'option' inputs. Thank you!

[Jimmy Shen]

Please notice that the values in colormap range from 0 to 65535 (instead of from 0 to 1). So black is represented by [0 0 0], and white is represented by [65535 65535 65535]. If you put [1 1 1], it will be black. Please see more details from my web site: http://www.rotman-baycrest.on.ca/~jimmy/geotiffwrite

[Evan]

I've been playing around with this code for a bit, but I'm having trouble with colormaps. If I use a standard colormap (which I am currently using to write standard tifs of my data), I get a solid black image. Any advice of support on the specific data format/colormap format required to successfully created a colored geotiff?

[Jimmy Shen]

If Model Type is "projected", you cannot use "bbox" variable directly, and you have to use "make_option.m" program to generate "option" variable.

[Matthew]

great program. If the Model Type is not geographic, it clears the bbox variable ? My model type is "projected". So what should happen there? Thanks Matt

[E]

Excelent job! Works with Octave too. Just replace '&' with '&&' in lines 257 and 1165.

[Jimmy Shen]

Hi David:

SRTM (Shuttle Radar Topography Mission) data are nothing but Big-Endian, 16-bit Signed Integer binary numbers with either in 1201x1201 (SRTM3) or in 3601x3601 (SRTM1) dimension. You can simply use the following script to do the job:

function srtm=load_srtm(fn)
fid=fopen(fn,'rb','ieee-be');
srtm.dat=fread(fid,inf,'*int16');
if length(srtm.dat)==1442401
srtm.dat=reshape(srtm.dat,[1201 1201]);
elseif length(srtm.dat)==12967201
srtm.dat=reshape(srtm.dat,[3601 3601]);
end
srtm.dat=rot90(srtm.dat);

Please let me know if you encounter any problem.

Jimmy

[David Davenish]

Very helpful program, Thank you!

I have a question though, where can I find the "load_srtm" command that you examplhttp://www.mathworks.com/matlabcentral/fileexchange/27959-geotiffwrite#ed in link http://www.rotman-baycrest.on.ca/~jimmy/geotiffwrite in "Description"?

[Jimmy Shen]

Hi David:

Thanks for you comments and rating.

Yes, you can use some of the fields contents returned by geotiffinfo to provide information for my option struct when you run the "make_option.m" program. However, please be careful, since my struct fields are DIRECTLY correspond to GeoTIFF Format Specification (like IDL software does); while struct fields from geotiffinfo in MATLAB software are NOT DIRECTLY correspond to GeoTIFF Format Specification. Here's a simple example: you can use geotiffinfo.TiePoints.WorldPoints for TiepointXYZ; but you cannot use geotiffinfo.TiePoints.ImagePoints for TiepointIJK.

I have just submitted 1-bit GeoTIFF support, so please come back in 2~3 days to see if they make my program available for you.

Jimmy

[David]

Great capability that Matlab has needed for a long time. I found that I can use the struct returned by geotiffinfo as the option struct for the input. I'd like to write a 1-bit geotiff. Any pointers would be appreciated.

[Karl]

MATLAB has been needing this for years. Thanks Jimmy. I really like the flexibility of this program.

[Zina]