View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
geoimread

5.0
5.0 | 1 rating Rate this file 25 Downloads (last 30 days) File Size: 8.42 MB File ID: #46904 Version: 1.5
image thumbnail

geoimread

by

Aslak Grinsted (view profile)

 

10 Jun 2014 (Updated )

Reads a sub region of a geotiff or geojp2 image.

| Watch this File

File Information
Description

GEOIMREAD reads a sub region of a geotiff or geojp2 image.
  
  
 % Syntax
 
  A = geoimread(filename)
  A = geoimread(filename,xlim,ylim)
  A = geoimread(filename,latlim,lonlim)
  A = geoimread(...,buffer)
  [A,x,y,I] = geoimread(...)
  geoimread(...)
 
  
 % Description
 
  A = geoimread(filename) returns the full image given by a filename. This
  syntax is equivalent to A = geotiffread(filename).
  
  A = geoimread(filename,xlim,ylim) limits the region of the geotiff file to
  the limits given by xlim and ylim, which are map units (usually meters) relative
  to the data projection. For example, if the geotiff is projected in Texas Centric
  Mapping System/Lambert Conformal coordinates, xlim and ylim will have units of
  meters relative to the origin (100 W, 18 N). xlim and ylim can be multimensional,
  in which case the limits of the map will be taken as the limits of the limits of
  the distribution of all points in xlim, ylim.
  
  A = geoimread(filename,latlim,lonlim) if no values in xlim, ylim exceed
  normal values of latitudes and longitudes, geoimread assumes you've entered
  limits in geographic coordinates of degrees latitude and longitude. The first
  input is latitude, the second input is longitude.
 
  A = geoimread(...,buffer) adds a buffer in map units (usually meters or feet) to the
  limits of the region of interest. This may be useful if you want to load an image
  surrounding scattered lat/lon data. If you'd like an extra 2 kilometers of image around
  your data, enter 2000 as the buffer. If buffer is a two-element vector, the first
  element is applied to the left and right extents of the image, and the second element
  is applied to the top and bottom extents of the image.
 
  [A,x,y,I] = geoimread(...) also returns pixel center coordinates (x,y) of the
  output image and a geotiff info structure I. I is a useful input for projfwd and projinv.
 
  geoimread(...) without any outputs shows the output image A without loading
  any data into the workspace.
 
  
 % Examples:
  
  % Show a whole geotiff:
  geoimread('boston.tif');
  
  % Compare results from above to a subset geotiff:
  mapx = [765884 766035 766963]; % units are feet
  mapy = [2959218 2957723 2958972];
  geoimread('boston.tif',mapx,mapy)
  
  % Or if you have coordinates in lat/lon and you want a 500 foot buffer:
  lat = [42.3675288 42.3634246 42.3668397];
  lon = [-71.0940009 -71.0934685 -71.0900125];
  
  geoimread('boston.tif',lat,lon,500);

Acknowledgements

This file inspired Geotiffinterp, Ramp Radarsat Antarctic Mapping Project, Ibcso International Bathymetric Chart Of The Southern Ocean Toolbox, L8read Landsat 8 Level 1 Tiff Reader, Landsat, and How To Drape Landsat Images Over Bedmap2 Topography.

Required Products Mapping Toolbox
MATLAB release MATLAB 8.1 (R2013a)
MATLAB Search Path
/
/html
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (1)
28 Jan 2017 Constantin Silkin

Updates
08 Jul 2014 1.1

added ability to read geojp2 files

29 Jul 2014 1.2

changed description to better reflect the geo jpg2000 features.

11 Dec 2014 1.3

Renamed and added Misc improvements from Chad Greene:
* More input checking and error messages.
* Clarified syntax and description in header.
and more...

11 Dec 2014 1.4

Further improvements from Chad Greene: Added html example

05 Nov 2015 1.5

Added some workarounds to make some of the functionality work without mapping toolbox.

Contact us