The software described in this document is furnished under a license agreement. The software may be used
or copied only under the terms of the license agreement. No part of this manual may be photocopied or
reproduced in any form without prior written consent from The MathW orks, Inc.
FEDERAL ACQUISITION: This provision applies to all acquisitions of the Program and Documentation
by, for, or through the federal government of the United States. By accepting delivery of the Program
or Documentation, the government hereby agrees that this software or documentation qualifies as
commercial computer software or commercial computer software documentation as such terms are used
or defined in FAR 12.212, DFARS Part 227.72, and DFARS 252.227-7014. Accordingly, the terms and
conditions of this Agreement and only those rights specified in this Agreement, shall pertain to and govern
theuse,modification,reproduction,release,performance,display,anddisclosureoftheProgramand
Documentation by the federal government (or other entity acquiring for or through the federal government)
and shall supersede any conflicting contractual terms or conditions. If this License fails to meet the
government’s needs or is inconsistent in any respect with federal procurement law, the government agrees
to return the Program and Docu mentation, unused, to The MathWorks, Inc.
Trademarks
MATLAB and Simulink are registered trademarks of The MathWorks, Inc. See
www.mathworks.com/trademarks for a list of additional trademarks. Other product or brand
names may be trademarks or registered trademarks of their respective holders.
Patents
The MathWorks products are protected by one or more U.S. patents. Please see
www.mathworks.com/patents for more information.
Summary by Version ...............................1
Contents
Version 3.1 (R2010a) Mapping Toolbox Software
Version 3.0 (R2009b) Mapping Toolbox Software
Version 2.7.2 (R2009a) Mapping Toolbox Software
Version 2.7.1 (R2008b) Mapping Toolbox Software
Version 2.7 (R2008a) Mapping Toolbox Software
Version 2.6 (R2007b) Mapping Toolbox Software
Version 2.5 (R2007a) Mapping Toolbox Software
Version 2.4 (R2006b) Mapping Toolbox Software
Version 2.3 (R2006a) Mapping Toolbox Software
Version 2.2 (R14SP3) Mapping Toolbox Software
......4
......10
....16
....17
......23
......34
......49
......55
......61
.....65
Version 2.1 (R14SP2) Mapping Toolbox Software
Version 2.0.3 (R14SP1) Mapping Toolbox Software
Version 2.0.2 (R14) Mapping Toolbox Software
Version 2.0.1 (R13SP1+) Mapping Toolbox Software
Version 2.0 (R13SP1+) Mapping Toolbox Software
.......86
.....68
....82
...87
....89
iii
Version 1.3.1 (R13SP1) Mapping Toolbox Software ....105
Version 1.3 (R13) Mapping Toolbox Software
Version 1.2 (R12) Mapping Toolbox Software
Compatibility Sum mary for Mapping Toolbox
.........107
.........109
........112
ivContents
SummarybyVersion
This table provides quick access to what’s new in each version. For
clarification, see “Using Release Notes” on page 2 .
Mapping Toolbox™ Release Notes
Version
(Release)
Latest Versi
V3.1 (R2010a
V3.0 (R2009b)
V2.7.2 (R2009a)
V2.7.1 (
V2.7 (R2008a)
V2.6 (R2007b)
R2007a)
V2.5 (
V2.4 (R2006b)
on
)
R2008b)
New Features
and Changes
Yes
Details
Yes
Details
NoYes
Yes
Details
Yes
Details
Yes
ls
Detai
Yes
Details
Yes
Details
Version
Compatibilit
Consideratio
Yes
Summary
Yes
Summary
Summary
Yes
Summary
Yes
Summary
Yes
Summa
Yes
Summary
Yes
Summary
ry
y
ns
Fixed Bugs
and Known
Problems
Bug Reports
includes fix
Bug Reports
includes fixes
Bug Repor
includes
Bug Reports
includes fixes
Bug Reports
includes fixes
Bug Rep
des fixes
inclu
Bug Reports
includes fixes
Bug Reports
includes fixes
es
ts
fixes
orts
Related
Documentation
at Web Site
Printable R elease
Notes: PDF
Current product
documentation
No
No
No
No
No
No
No
V2.3 (R2006a)
.2 (R14SP3)
V2
V2.1 (R14SP2)
Yes
ails
Det
Yes
Details
Yes
Details
Yes
mary
Sum
Yes
Summary
Yes
Summary
Reports
Bug
ludes fixes
inc
Bug Reports
includes fixes
Bug Reports
includes fixes
No
No
No
1
Mapping Toolbox™ Release Notes
Version
(Release)
V2.0.3 (R14SP1)
V2.0.2 (R14)
V2.0.1 (R13SP1+)
V2.0 (R13SP1+)
V1.3.1 (R13SP1)
V1.3 (R13)
V1.2 (R12)
New Features
and Changes
Yes
Details
Yes
Version
Compatibility
Considerations
Yes
Summary
No
Fixed Bugs
and Known
Problems
Fixed bugs
No bug fixes
Related
Documentation
at Web Site
No
No
Details
Yes
Details
No
“Fixed Bugs” on
page 87
Version
2.0.1 product
documentation
Yes
Details
Yes
Details
Yes
Yes
Summary
No
No
“Fixed Bugs” on
page 103
“Fixed Bugs” on
page 105
Fixed bugs
No
No
No
Details
Yes
No
No bug fixes
No
Details
Using Release Notes
Use release notes when upgrading to a newer version to learn about:
• New features
• Changes
• Potential impact on your existing files and practices
Review the release notes for other MathWorks™ products required for this
product (for example, MATLAB
®
or Simulink®). Determine if enhancements,
bugs, or compatibility considerations in other products impact you.
If you are upgrading from a software version other than the m ost recent one,
review the current release notes and all interim versions. For example, when
you upg rade from V1.0 to V1.2, review the release notes for V1.1 and V1.2.
2
SummarybyVersion
What Is in the Rel
New Features and
• New functional
• Changes to exi
Version Compa
When a new fea
versions, th
impact.
Compatibil
Reports at
in incompa
compatibi
Fixed Bug
The Math
view Bug
time and
provisi
availa
is not a
ity issues reported after the product release appear under Bug
The MathWorks™ W eb site. Bug fixes can sometimes result
tibilities, so review the fixed bugs in Bug Reports for any
lity impact.
sandKnownProblems
Works offers a user-searchable Bug Reports database so you can
Reports. The development team updates this database at release
as more information becomes available. Bug Reports include
ons for any known workarounds or file replacements. Information is
ble for bugs existing in or fixed in Release 14SP2 or later. Information
vailable for all bugs in earlier releases.
ity
sting functionality
tibility Considerations
ture or change introduces a reported incompatibility between
e Compatibility Considerations subsection explains the
ease Notes
Changes
s Bug Reports using your MathWorks Account.
Acces
3
Mapping Toolbox™ Release Notes
Version 3.1 (R2010a) Mapping Toolbox Software
This table summarizes what’s new in Version 3.1 (R2010a):
New Features
and Changes
Yes
Details below
• “WMS Database Modified” on page 4
• “The etopo Function Now Supports the ETOPO1 and ETOPO2v2 Data
Sets” on page 5
• “Now Possible to Retrieve Legend for WMS Map” o n page 6
• “The shaperead and shapewrite Functions Now Suppo rt Non-ASCII
Characters” on page 6
• “Display Range Increased for eqdazim and eqaazim Projections ” on page 6
• “The GUIs mlayers and mobjects No Longer Support EraseMode” on page 7
• “scatterm Now Returns a Handle to an hggroup Object” on page 7
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of changes,
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Bug Reports at
Web site
Related
Documentation
at Web Site
Printable
Release Notes:
PDF
Current product
documentation
• “mdistort Now Returns a Handle to a contourgroup Object” on page 8
• “polybool No Longer Errors when Given Empty Input Vertex Arrays” on
page 8
• “Functions Being Removed” on page 9
WMS Database Modified
The WMS database changes on a release-to-release basis, as some new
servers are added and some unavailable servers are removed. A total of 199
servers (10.2% of the original number of servers) and 66,270 layers have
4
Version 3.1 (R2010a) Mapping Toolbox™ Software
been removed from the database. A total of 602 servers, with 207,269 layers,
have been added.
Some notable new servers in the database are:
• The OnMars (
(
http://onmoon.jpl.nasa.gov/wms.cgi) servers from the Je t Propulsion
Some servers are no longer accessible and have been rem ov ed from the
database. If you have code based on theseservers,itwillnolongerrun.Tofix
this problem, search the WMS Database for another comparable s erver.
Some examples from the R2009b documentation have been modified due to
server inaccessibility. Some notable servers that have been removed are:
• Two of the servers from the European Space Agency (ESA)
(
http://mapdev.eo.esa.int/mapServer/mapServer and
http://mapdev.esrin.esa.int/mapServer/mapServer)
®
The etopo Function Now Supports the ETOPO1 and
ETOPO2v2 Data Sets
Before R2010a, the etopo function supported ETOPO2–2001 (2-minute)
and ETOPO5 (5-minute) data. Support has been added for ETOPO2v2c
(2-minute) and ETOPO1c (1-minute) data. The ETOPO1 model, released in
March 2009, is the most recent and contains the highest resolution data. For
information on downloading the ETOPO data sets, see Technical Note 2101:
5
Mapping Toolbox™ Release Notes
Accessing Geospatial Data on the Internet for the Mapping Toolbox on the
MathWorks Web site.
Compatibility Considerations
The etopo function still works with ETOPO2 and ETOP O5 data.
Now Possible to Retrieve Legend for WMS Map
The Details property of the WMSLayer classcontainsaStyle field. A
LegendURL structure has been added to this field. The information in the
LegendURL structure, if provided by the server, enables you to retrieve a
legend image for a specific WMS map.
The shaperead an d shapewrite Functions Now
Support Non-ASCII Characters
You can now use the shaperead and shapewrite functions to import and
export attributes with non-ASCII characters. To use this feature, you must
set your character encoding scheme to match that used by the shapefile. For
example, if your session is configured to support
and you want to import a shapefile with Japanese characters, you m ust first
change your configuration to support
US_ASCII character encoding
Shift_JIS.
Display Range Increased for eqdazim and eqaazim
Projections
In previous releases , the Equidistant Azimuthal (eqdazim)andLambert
Azimuthal Equal-Area (
display of data points farther than 160 degrees from the projection origin.
Now you can set the projection radius for either of these projections to up
to 179.5 degrees.
Use the
the largest possible value, pass the following p roperty name-value pair to
axesm or setm:
FLatLimit property to control this setting. For example, to choose
'FLatLimit',[-Inf 179.5]
6
eqaazim) projections did not allow projection or
Version 3.1 (R2010a) Mapping Toolbox™ Software
Compatibility Considerations
The d efault value of the projection radius remains 160 degrees, so you do not
need to update any code that relies on the default value.
BeforeR2010a,youcouldusethefollowinglineofcodetoinitializethe
equatorial aspect of an azimuthal projection:
The EraseMode property, represented by the Emode button, has been
removed from the
mlayers and mobjects GUIs.
Compatibility Considerations
Before R2010a, the mlayers and mobjects GUIs had an EraseMode option,
whichmadeitpossibletosettheerasemodeonaparticularmaplayeror
object. This property controlled the technique MATLAB used to draw and
erase hggroup child objects.
scatterm Now Returns a Handle to an hggroup
Object
The function scatterm in the syntax
h = scatterm(...)
7
Mapping Toolbox™ Release Notes
now returns a handle to an hggroup.
Compatibility Considerations
In previous releases, sc atterm returned a vector of patch handles.
mdistort Now Returns a Handle to a contourgroup
Object
The output of the syntax
is now a scalar handle to a contourgroup object containing the contours and
text.
Compatibility Considerations
In previous releases, mdistort returned handles to the line and text objects
andusedthesyntax
h = mdistort(...)
[h,ht] = mdistort(...)
with two output arguments. The second output of mdistort is now redundant
because these handles will be available as children of
h.
polybool No Longer Errors when Given Empty Input
Vertex Arrays
When one or both pairs of input vertex arrays is empty, the polybool function
now returns either empty values or the values of the non-empty input pair,
depending on the requested operation.
In previous releases, if one or both pairs of input vertex arrays were empty,
the
polybool function would issue an error.
Functions Being Removed
Function Name
etopo5
What Happens
When You Use
the Function?
Errors
Use This
Instead
etopo
Compatibility
Considerations
Replace all
existing
instances of
etopo5 with
etopo.
9
Mapping Toolbox™ Release Notes
Version 3.0 (R2009b) Mapping Toolbox Software
This table summarizes what’s new in Version 3.0 (R2009b):
New Features
and Changes
Yes
Details below
• “New Features for Creating Web Map Service Maps” on page 10
• “New makerefmat Syntax for Constructing Referencing Matrices” on page
11
• “Some Functions Now Accept Referencing Matrices as Input” on page 11
• “New Angle Conversion Functions” on page 11
• “Expanded Support for GSHHS Global Coastline Data” on page 12
• “New Behavior for polymerge when Three or More Line Segments Have
Common End Point” on page 12
• “Automatic Conversion of Latitude Limits to Ascending Order” on page 12
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of changes,
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Bug Reports at
Web site
Related
Documentation
at Web Site
Printable
Release Notes:
PDF
Current product
documentation
10
• “Second Input Argument of roundn No Longer Supports Complex Numbers,
Non-integers, or Default Values” on page 13
• “Functions Removed” on page 13
New Features for Creating Web Map Service Maps
New functions and classes now make it possible to interact with Web Map
Service (WMS) servers and render WMS maps. Use the new features to
search a built-in database of pre-qualified WMS servers and layers. Retrieve
customized geographic data sets and related metadata from WMS servers.
The new classes encapsulate WMS servers, data layers, metadata, and map
Version 3.0 (R2009b) Mapping Toolbox™ Software
requests. See the Creating Web Map Service Maps chapter in the User’s
Guide and related “Class Reference” for more information.
New makerefmat Syntax for Constructing
Referencing Matrices
A new parameter name-value pair syntax makes it easier to construct
referencing matrices with
image or raster grid that is referencedtoandalignedwithageographic
coordinate sy stem but not for one that is referenced to a 2-D map coordinate
system. Use parameters to set the number of rows (M) and columns (N) of
the raster or image to be used with the referencing matrix; the latitude and
longitude limits of the geographic quadrangle bounding the georeferenced
raster; and the edges from which row and column indexing start, designating,
for example, columns that run either south-to-north or north-to-south.
makerefmat. You can use the new syntax for an
Some Functions Now Accept Referencing Matrices
as Input
The functions below now accept referencing matrices as input, and some of
them (
output. The functions that generate referencing matrices as output do so only
in cas es where referencing matrices are used as input. If referencing vectors
are used as input, referencing vectors are also generated as output. Note that
the functions in this table work exclusively with data grids or images that are
referenced to geographic (latitude or longitude) coordinates.
maptrims, resizem,andvec2mtx) also generate referencing matrices as
areamatimbedmmaptrimssetltlnusamap
filtermlimitmmeshgratsetpostnworldmap
findmlos2meshlsrmresizem
getseedsltln2valmeshmvec2mtx
gradientmmapprofileneworigviewshed
New Angle Conversion Functions
Two new functions, radtodeg and degto rad, replace older functions rad2deg
and deg2rad. They are functionally identical.
11
Mapping Toolbox™ Release Notes
Compatibility Considerations
The older functions, rad2deg and deg2rad,willcontinuetowork,butwhen
writing new code the newer versions are recommended.
Expanded Support for GSHHS Global Coastline Data
The gshhs function now suppo rts Versions 1.4 and later of the Global
Self-Consistent, Hierarchical, High-Resolution Shoreline Database (GSHHS).
New B ehavior for polymer ge when Three or More
Line Segments Have Common End Point
The behavior of the polymerge function has changed in cases of three or more
distinct parts with a common end point. In such cases, the choice of which
parts to merge is ambiguous; therefore, none of the corres ponding parts are
connectedatthatcommonpoint.
Compatibility Considerations
In previous releases, if three or more parts shared a common end point, the
polymerge function attempted to merge them. T he result, however, was
unspecified and sometimes obviously wrong.
12
Automatic Conversion of Latitude Limits to Ascending
Order
The functions axesm and setm require that the latitude limits in the
'MapLatLimit' property be provided in ascending orde r . If you enter the
limits in descending order, these functions will now automatically convert the
limits to ascending order, and return a warning message notifying you of
this change.
Compatibility Considerations
In prev io us releases, if you entered the latitude limits of the 'MapLatLimit'
property in descending order when using axesm or setm, you could end up
with a map axes that was in t ern ally inconsistent, possibly resulting in
unexpected errors during subsequent operations. This is no longer the case.
Version 3.0 (R2009b) Mapping Toolbox™ Software
Second Input Argument of roundn No Longer
Supports Complex Numbers, Non-integers, or
Default Values
roundn no longer accepts certain types of input for the second input argument,
N, which is supposed to be a real, integer-valued exponent. Now, if you use
a complex number or non-integer as the second input to
receive a n error; and if you omit
receive an error if you call
error message strings.
N, you will receive a warning. You will also
roundn with a second output argument to capture
Compatibility Considerations
In previous releases, if you used a complex number or non-integer as the
second input to
you called
strings, you would receive an obsolete syntax warning.
roundn, th is number would be converted into a real integer. If
roundn with a second output argument to capture error message
roundn,youwill
In R2009b,
default is being p h ase d out, however, and in the future you will receive
an error if you fail to supply the second input argument. If there are any
instances in your code with the usage
with
The two-output syntax option was pre v iously deprecated and has resulted in
a warning in the past several releases.
-2 is still the default value for the second input argument. This
roundn(x), you should replace them
roundn(x,-2).
Functions Removed
Functions Being Removed in a Future Release
Name
colorui
eastof
Stage
Still
runs
Warns
Compatibility Considerations
Replace all existing instances of colorui with
uisetcolor.
If you are using degrees, replace
eastof(lon, meridian, 'degrees')
with
13
Mapping Toolbox™ Release Notes
Functions Being Removed in a F uture Release (Continued)
Name
imagem
smoothl
tgrline
str
unit
ong
Stage
Errors
Warns
Still
runs
Warns
Compatibility Considerations
meridian + mod(lon - me ridi an, 360)
and if you are
eastof(lon, meridian, 'radians')
using radians, replace
with
meridian + mod(lon - me ridi an, 2*pi)
Replace all existing instances of
grid2image.
Use
require
in degre
More r
in sha
shape
apMultipart
unwr
s its input to be in radians. When working
es, use
radtodeg(unwrapMultipart(degtorad(lon)))
ecent Tiger/Line
pefile format and can be imported using
read
.
instead. This function
®
data sets are available
imagem with
The syntax
14
westof
Wa
str =
unitstr(str,'times')
has already been removed.
rns
If you are using degrees, replace
stof(lon, meridian, 'degrees')
we
with
eridian - mod(meridian - l on, 360)
m
and if you are using radians, replace
Version 3.0 (R2009b) Mapping Toolbox™ Software
Functions Being Removed in a F uture Release (Continued)
Name
Stage
Compatibility Considerations
westof(lon, meridian, 'radians')
with
meridian - mod(meridian - lon, 2*p i)
Functions Removed in R2009b
deg2dmhms2hmhr2secsec2hmscmapui
deg2dmshms2hrmat2dmssec2hrtigermif
dms2deghms2matmat2hmstime2strtigerp
dms2dmhms2secrad2dmtimedim
dms2mathr2hmrad2dmscontorm
dms2radhr2hmssec2hmcontor3m
The functions above have been completely removed from the toolbox and
error if used.
15
Mapping Toolbox™ Release Notes
Version 2.7.2 (R2009a) Mapping Toolbox Software
This table summarizes what’s new in Version 2.7.2 (R2009a):
New Features
and Changes
No
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of changes,
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Bug Reports at
Web site
Related
Documentation
at Web Site
Printable
Release Notes:
PDF
Current product
documentation
geoshow and mapshow Now Construct Ordinary
Patch Objects
When displaying polygons, instead of constructing graphics objects whose
classes derive from patch,
objects.
Compatibility Considerations
This change has no effect on the display, and it is unlikely that you will need
tochangeanycodeinresponsetoit. Infact,itischallengingtodetectthis
change at all.
geoshow and mapshow now construct ordinary patch
16
The older (derived) class included several extra properties used only for
internal bookkeeping. You might notice them if you call
In R2009a, these properties have b een removed, and the output of
different in terms of both layout and property order. (It now looks the same
as for any ordinary patch.)
get on the handle.
get looks
Version 2.7.1 (R2008b) Mapping Toolbox™ Software
Version 2.7.1 (R2008b) Mapping Toolbox Software
This table summarizes what’s new in Version 2.7.1 (R2008b):
New Features
and Changes
Yes
Details below
• “Using the Map Axes Map Limit Properties with axesm, setm, and
defaultm” on page 17
• “Changing Projection Type of an Existing Map Axes with setm” on page 18
• “Other Bug Fixes with Compatibility Considerations” on page 19
• “coast.MAT Data File Revised” on page 19
• “Map Limit Syntaxes Removed” on page 20
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of changes,
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Bug Reports at
Web site
Related
Documentation
at Web Site
Printable
Release
Notes: PDF
Current product
documentation
UsingtheMapAxesMapLimitPropertieswith
axesm, setm, and defaultm
Changes and enhancements have been made to axesm, setm,anddefaultm
with respect to map axes properties that affect the fundamental display
geometry:
•
MapProjection
• Zone
• Origin
• FLatLimit
• FLonLimit
17
Mapping Toolbox™ Release Notes
• MapLatLimit
• MapLonLimit
The changes result in the following improvements:
• The use of the map limit properties to set up a map axes is mo re intuitive.
• The way in which
properties is now consistent with the behavior of
• The map limit properties (
properties (
and
Zone propertie s interact in a more clear and predictable fashion.
For more information, see the section “Using the Map Limit Properties” in the
Mapping Toolbox™ User’s Guide and bug report 319891 on the MathWorks
Web site.
FLatLimit and FLonLimit), and the Origin, MapProjection,
defaultm resolves possible inconsistencies betw een these
axesm and setm.
MapLatLimit and MapLonLimit), the frame limit
Changing Projection Type of an Existing Map Axes
with setm
In previous releases, calling the setm function to change the MapProjection
property of a map axes, especially when switching between an azimuthal and
non-azimuthal projection (e.g., a conic or cylindrical projection), often resulted
in the following types of problems:
• The modified map axes mig ht cover a different part of the Earth.
• The map frame and graticule might fail to update properly.
• Map limit properties changed at the same time as the projection might
not have the proper effect.
18
The
setm function now more effectively resets the projection, clearing out
settings that were specific to the earlier projection, updating the map frame
and graticule, and staying in the same general part of the world (even when
switching between azimuthal and non-azimuthal projections).
Version 2.7.1 (R2008b) Mapping Toolbox™ Software
Compatibility Considerations
You may need to change the way in wh ich you reset various map axes
properties, such as
Origin, FLatLimit,andFLonLimit after changing
projections, as discussed in the section “Switching Between Projections” in the
Mapping Toolbox User’s Guide. In many cases it will no longer be necessary
to reset as many properties.
Other Bug Fixes with Compatibility Considerations
• The default FLatLimit for lambert and lambertstd has been changed to
[-45 45]. In previous releases, axesm produced huge map frames, due to
the
FLatLimit default of [-90 90].
• The function
the parallels and meridians. In previous releases, a call to
linespec or property name/property value syntaxes returned empty.
• The function
maplonlimit fields.
• A reference ellipsoid set to a non-default value (via the
gridm now returns handles to the line objects used to display
gridm using
geotiff2mstruct no longer sets the maplatlimit and
geoid property)
no longer reverts to the default when the UTM zone is reset. For more
information, see bug report 459353 on the MathWorks Web site.
• The
daspectm function now works for azimuthal projections and units of
radians.
coast.MAT Data File Revised
Portions of the global coastline l ati tud e-longitude vectors in the coast.MAT
data file have been revised to ensure proper polygon topology. The data edits
comprise the following:
• Replacingorremovingvarious"bow-tie" and degenerate linear
(non-polygonal) island features.
• Opening a "pinched" section in the middle of Lake Balkhash in Central
Asia.
• Merging the eastern and western sections of Wrangel Island ne ar the
Bering Strait (cut by the 180-degree m e ridian) into a single polygon with
longitudes ranging from slightly less than 180 to slightly greater than 180.
19
Mapping Toolbox™ Release Notes
• Eight additional edits to pull apart landmasses w ith points of contact and
remove coastal "spikes."
Map Limit Syntaxes Removed
The f ollow ing syntaxes are obsolete. An error occurs if you u se them.
•
pcolorm(Z)
• pcolorm(Z,gratsize)
• surfacem(Z)
• surfacem(Z,gratsize)
• surflm(Z)
• surflm(Z,s)
• surfm(Z)
• surfm(Z,gratsize)
20
These syntaxes displ ay ed a data grid with geographic limits that matched
the map latitude and longitude limits in the current map axes. Using the
old syntaxes correctly involved knowing the latitude and longitude limits of
your data and matching them to the values listed under
maplonlimit in the map axes properties. We have replaced these syntaxes
maplatlimit and
with a more dire ct approach that requires you to enter the latitude and
longitude limits for the data grid.
Compatibility Considerations
The table below suggests alternative code to replace the obsolete syntaxes.
In the following table,
and
gratsize is a two-element vector specifying the size of the graticule
isahandletothesurfacethatisdisplayed.Andlatlim and lonlim are the
h
Z displays:
geographic limits of the data grid (in degrees):
latlim = [southern_limit northern_limit]
Z is a regular data grid (a 2-D array of cl a ss double)
Version 2.7.1 (R2008b) Mapping Toolbox™ Software
lonlim = [western_limit eastern_limit]
Original SyntaxReplacement Syntax
h = pcolorm(Z)
[lat,lon] = meshgrat(latlim,lonlim,size(Z));
constructs a surface using the regular
data grid
(using
size(Z) and with geographic limits
Z and a graticule mesh
meshgrat) with size equal to
that match the map latitude and
longitude limits in the current map
axes.
h = pcolorm(Z,gratsize)
uses a graticule mesh with size equal
to
gratsize.
h = surfacem(Z)
constructs a surface using the regular
data grid
(using
Z and a graticule mesh
meshgrat) of size 50-by-100.
The geographic limits match the map
latitude and longitude limits in the
current map axes.
h = surfacem(Z,gratsize)
uses a graticule mesh with size equal
to
gratsize.
h = surflm(Z)
h = pcolorm(lat,lon,Z)
[lat,lon] = meshgrat(latlim,lonlim,gratsize);
h = pcolorm(lat,lon,Z)
h = surfacem(latlim,lonlim,Z)
[lat,lon] = meshgrat(latlim,lonlim,gratsize);
h = surfacem(lat,lon,Z)
h = surflm(latlim,lonlim,Z)
constructs a surface using the regular
data grid
(using
size(Z) and with geographic limits
Z and a graticule mesh
meshgrat) with size equal to
that match the map latitude and
longitude limits in the current map
axes. It is displayed with a default
light source.
21
Mapping Toolbox™ Release Notes
Original SyntaxReplacement Syntax
h = surflm(Z,s)
h = surflm(latlim,lonlim,Z,s)
specifies the direction of the light
source.
s is a two- or three-element
vector that specifies the direction
fromthesurfacemaptothe
light source as defined in the
documentation for
h = surfm(Z)
surfl.
h = surfm(latlim,lonlim,Z)
constructs a surface using the regular
data grid
(using
size(Z) and with geographic limits
Z and a graticule mesh
meshgrat) with size equal to
that match the map latitude and
longitude limits in the current map
axes.
h = surfm(Z,gratsize)
[lat,lon] = meshgrat(latlim,lonlim,gratsize);
uses a graticule mesh with size equal
to
gratsize.
22
h = surfm(lat,lon,Z)
Version 2.7 (R2008a) Mapping Toolbox™ Software
Version 2.7 (R2008a) Mapping Toolbox Software
This table summarizes what’s new in Version 2.7 (R2008a):
New
Features
and
Changes
Yes
Details b elo w
• “Functions for Working with Geographic Quadrangles” on page 24
• “Fixes and Improvements to Function avhrrgoode” on page 25
• “Improved Accuracy for the limitm and setpostn Functions” on page 26
• “New Point Location Demo Data for Tsunami Events” on page 27
• “Better Trimming Benefits fillm Function” on page 28
• “Restored units Options for Function angl2str” on page 28
• “New Longitude-Wrapping Option in the closePolygonParts Utility ” on
page 29
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of new features
and changes,
below. See also
Summary.
Fixed Bugs
and Known
Problems
Bug reports at
Web site
Related
Documentation
at Web Site
Printable Release
Notes: PDF
Current product
documentation
• “Changes to Terminology for Geographic Data Structures” on page 29
• “Identifiers Provided for all Warnings” on page 30
• “Documentation for Functions tigermif and tigerp Removed” on page 31
• “Removed Syntaxes that Returned Error Messages in Optional A rgument”
on page 31
23
Mapping Toolbox™ Release Notes
Functions for Working with Geographic Quadrangles
A geographic quadrangle is an area on the surface of a sphere or ellipsoid
bounded on the east and west by a pair o f meridians and on the north and
south by a pair of parallels. In many ways, such an object is similar to
a bounding rectangle in the plane, but they can be difficult to work with
because of the way longitudes w rap around and the way meridians converge
at the p oles. For example,
• The western longitude limit can have a larger numerical value than the
eastern longitude limit.
• If one of the bounding latitudes is +90 or -90 degrees, the quadrangle has
three sides rather than four.
• As noted below, the intersection o f two geographic quadrangles might
possibly comprise two separate parts—with the eastern e nd of the first
quadrangle intersecting the western end of the second quadrangle, and
vice versa.
Mapping Toolbox software typically represents a geographic quadrangle in
terms of its latitude and longitude limits, stored i n 1-by-2 vectors having
the forms
Vectors like these have been used in various Mapping Toolbox functions
since its inception, and can appear in the input or output argument lists of
over dozen functions.
In R2008a, three new functions let you query, intersect, and display
geographic quadrangles, and account for subtleties such as those described
above:
•
ingeoquad — Returns true for points inside or on latitude-longitude
quadrangle
•
intersectgeoquad — Returns intersection(s) of two latitude-longitude
quadrangles
•
outlinegeoquad — Returns sampled p olyg on vertices for a
latitude-longitude geographic quadrangle
Version 2.7 (R2008a) Mapping Toolbox™ Software
Use inge oqua d, for example, to check whether a geographic point is located
within the area covered by a regular data grid, given the latitude and
longitude limits computed by
Use
intersectgeoquad to compute overlap, if any, between two quadrangles.
limitm.
Interestingly, three general results are possible : no intersection, an
intersection that is itself a geographic quadrangle, and an intersection the
comprises two distinct geographic quadrangles. (The intersection can have
two parts if the input quadrangles wraparoundinlongitudetooverlapon
both their eastern and western sides. This case, of course, is not possible for
bounding boxes in the plane.)
Use
outlinegeoquad to generate a pair of latitude and longitude coordinate
vectors that define a polygon that traces the outline of a geographic
quadrangle. This can be useful for displaying the quadrangle graphically
using
geoshow, for example, especially on a projection where the meridians
and/or parallels do not project to straight lines, because in addition to
connecting the four corners
outlinegeoquad lets you interpolate additional
vertices along parallels, meridians, or both.
Fixes and Improvements to Function avhrrgoode
Function avhrrgo ode has been rewritten to improve its efficiency and to
remove a number of problems and limitations:
• Fixed a spatial referencing problem when a nonglobal region has been
specified which caused locationstobeoffsetbyhalfapixel.
• The function no longer returns incorrect NaN coordinate values at the
equator when given certain latitudelimitsthatcrosstheequator.
• The function no longer errors when attempting to read a file name with
certain legal latitude and longitude limits.
• The new version executes at least five times faster.
Compatability Considerations
• The nonfunctional syntaxes avhrrgoode and avhrrgoode(r egio n) have
been removed from the documentation.
25
Mapping Toolbox™ Release Notes
• The function now returns empty when the user-supplied limits are outside
data limits.
• The function no longer permits longitude limits to be specified outside
the interval [-180 180].
• Parameters other than
to use their default values.
• In versions prior to R2008a, when reading from the global data set and
a smaller region data set, the size of the outputs differed by one column
from each other when given identical latitude and longitude limits. Now
thesizesarethesame.
region and filename can be specified as empty
Improved Accuracy for the limitm and setpostn
Functions
In previous releases, after calculating the latitude and longitude limits of
the geographic quadrangle bound in g a regular data grid , function
arbitrarily rounded those limits to the nearest one millionth o f a degree
(equivalent to about 10 cm in latitude or equatorial longitude). Although it is
small, this rounding operation in effect applied an arbitrary shift to points on
or very near the edge of the grid. The direction of the shift and its m agnitude
were arbitrary because rounding can either increase or decrease a value. In
any given case, the shift depended on the specific referencing vector and the
number of columns and rows in the data grid. This behavior unnecessarily
degraded the numerical accuracy of
on it, and it has now been rem oved. For more information, see bug report
420038 on the MathWorks Web site.
limitm and those functions which depend
limitm
26
In the
Additional changes eliminate a problem for certain input points near
boundaries between grid cells that caused row and column subscripts returned
by
of the data grid—but still within the grid—returned subscript values could
exceed the corresponding grid size. For more information, see bug report
173338 on the MathWorks Web site.
setpostn function , an identical rounding step has been removed.
setpostn to be off by 1. For points near the northern and eastern edges
Version 2.7 (R2008a) Mapping Toolbox™ Software
Compatibility Considerations
These corrections can cause subtle changes in the behavior of other functions
that work with regular data grids referenced to latitude-longitude, for
example,
imbedm.
If your referencing vector contains approximations to rational numbers that
do not have an exact a 64-bit floating point representation (e.g., for cells that
are 1.5 degrees wide,
refvec(1) is 0.666666...), you may still find that certain
points that are extremely close to a grid cell boundary cross into a neighboring
cell just across the boundary. Such numerical ambiguity is inevitable given
how the information in a referencing vector is encoded. Although it cannot
be eliminated within
within a few factors of
setpostn, the inexactness only affects points that fall
eps (very much less than a millionth of a degree) away
from a given cell boundary.
New Point Location Demo Data for Tsunami Events
The Mapping Toolbox demo data in the $MATLABROOT/toolbox/map/mapdemos
directory now includes a global tsunami data set in shapefile format with
'Point' geometry. The data set comprises four files:
is not part of the shapefile set. It is a text file docum enting
the data set.
The data includes tidal wave events for which the maximum water height
was at least one meter, ranging for the years 1950 to 2006, inclusiv e. The
Global Tsunami Databa se, U.S. National Geospatial Data Center (NGDC),
National Oceanic and Atmospheric Administration (NOAA), available at
http://www.ngdc.noaa.gov/seg/hazard/tsu.shtml, is the source of the
data. All the files consist of U.S. Government information that is in the public
domain and is not subject to copyright protection.
The approximate location of each event is a single point in geodetic coordinates
(latitude-longitude) with an unspecified datum. The
.dbf file contains 18
27
Mapping Toolbox™ Release Notes
separate text or numeric attributes for most events, including wave height,
causes and seismic magnitudes, and location and country names.
TheshapefileswerecreatedatTheMathWorksfromqueryingtheonline
source data, importing the results into the MATLA B workspace, and
exporting them using the Mapping Toolbox
information, type
edit tsunamis.txt
at the M ATLAB prompt.
Better Trimming Benefits fillm Function
The changes described in the Versin 2.6 (R2007b) release note “Improvements
to Data Trimming in patchm and patchesm” on page 40 resulting from
improved polygon trimming also apply to the
Restored units Options for Function angl2str
The angl2str function once again can format strings for angles in
degrees-minutes (DM) and degrees-minutes-second (DMS) notations. These
options were removed in Version 2.6 (R2007b), and have now been restored.
In addition to the
obtain DM- and DMS-formatted strings by specifying
shapewrite function. For more
fillm function.
'degrees' and 'radians' units optio ns, you can now
28
•
'degrees2dm' — for degrees-decimal minutes formatting
'degrees2dms' — for degrees-minutes-decimal seconds formatting
•
To use these options, input angles must be in degrees. That is,
the string
units to indicate both the units in which the angle argument is
angl2str uses
provided and to control the output format.
This change restores the behavior of
angl2str prior to V ers ion 2.6 in a
slightly different form. Before V. 2.6, the DM and DMS options were specified
by a
units strings of 'dm' and 'dms', respectively. The new strings that
replace them signify that the functions
degrees2dm and degrees2dms,
introduced in Version 2.5 (R2007a), perform the conversions of inputs given in
degrees to DM and DMS notation.
Version 2.7 (R2008a) Mapping Toolbox™ Software
New Longitude-W
closePolygonPa
The close Pol ygo
angleunits,th
'radians'. If
first two
for longitude
latitude with
alongitudeo
is not added.
at must be string-valued and can be either
you include this argument with a value appropriate for the
, lon)
(lat
wrapping. For example, a polygon that begins at a given
a longitude of -180 degrees, and ends at the same latitude with
f 180 degrees is regarded as closed and an a dditional vertex
Changes to T
rts Utility
nParts
arguments, cl oseP olygonParts can correctly account
erminology for Geographic Data
Structures
From Versi
to “versio
structure
shorthan
the toolb
defined t
• Geograp
feature
on 2.0 onward, the Mapping Toolbox docume ntation has referred
n 1 geographic data structures” and “version 2 geographic data
s,” using the terms “geostruct1” and “geostruct2” respectively as
d for them. T o reflect current usage, starting with this version of
ox, these terms are obsolete; new terms and distinctions hav e been
o help clarify what these structures are and can be used for:
hic data structure arrays, introduced in Version 2.0, contain vector
s and are called either
rapping Option in the
function now accepts an optional third argument,
'degrees' or
- Geostr
- Mapstr
• Displ
geogr
geoda
Due t
form
age raste r geodata is with regular or geolocated data grids (2-D numeric
pack
ays accompanied by referencing matrices or vectors). There are only a
arr
Mapping Toolbox functions that can still generate display structures (by
few
orting data from external file formats):
imp
wdata
•
dc
ucts, if they contain geographic coordinates (latitudes and
udes)
longit
ucts, if they contain projected map/planar coordinates (x and y)
ay structure arrays, dating from Version 1, also used to be called
aphic data structures, and can contain either vector features or raster
ta.
o their greater generality, geostructs and mapstructs are the preferred
in which to represent vector features in the toolbox. The preferred way to
— Returns line/patch display structures
29
Mapping Toolbox™ Release Notes
• dcwgaz — Returns line/patch display structures
demdataui — Returns “regular”—as in regular data grid, that is—display
•
structures
•
mlayers — GUI to control plotting of display structure elements
tgrline — Returns line/patch display structures
•
vmap0data — Returns line/patch display structures
•
vmap0ui — GUI for selecting data from Vector Map Level 0
•
Even fewer functions accept display structures as inputs:
•
displaym — Displays elements of a display structure
extractm — Extracts lat-lon coordinates from line/patch display structure
•
In addition to
displaym and extractm,theupdategeostruct function
converts a line or patch display structure to a geostruct.
For m ore information, see “Mapping Toolbox Geographic Data Structures”.
Identifiers Provided for all Warnings
All warnings issued from within Mapping Toolbox functions now include
identifiers, enabling you to suppress them at your own discretion. Previously,
this was possible for only certain warnings, but with the addition of new
identifiers in over two dozen functions in R2008a, all warnings are now
covered. For example, you canturnoffthewarningthat
when given a latitude-longitude position outside the limits of the specified
data grid. In this case, the warning identifier is
'map:setpostn:pointOutsideLimits'
You can suppress it using the following statement:
SeetheR2007breleasenote“FunctionstigermifandtigerpAreObsoleteand
Error if Used” on page 48 for alternatives to
tigermif and tigerp.
Removed Syntaxes that Returned Error Messages in
Optional Argument
In earlier versions, the following Mapping Toolbox functions supported
syntaxes that included an optional output argume nt called
argument was included in a call to one of these functions, and certain error
conditions were encountered while executing the function, then instead of
issuing a n error, the function would return the corresponding error message
in
msg. The following functions are affected:
msg. If this output
•
axesm
• defaultm
• displaym
• gcm
• handlem
• lightm
• linem
• maps
31
Mapping Toolbox™ Release Notes
• meshm
• namem
• patchesm
• roundn
• surfacem
• surflsrm
• textm
• unitstr
• utmzone
• utmzoneui
For example, even with no map axe s present, the command
[mstruct, msg] = gcm
32
returned without error in R2007b and earlier, but placed an error message
string in
msg.
These syntaxes have been disabled in R2008a. If you try to use them, a
warning is issued. The warning may be followed by an error, depending on
whether or not an error condition is encountered within the function. For
example, if a map axes is present, the command above results in
Warning: Function GCM no longer returns error message strings in
output argument MSG. Inst ead any errors are thrown where they occur.
You should remove the la st output argument (MSG) from your call to
GCM in order to avoid this warning. If you want to handle errors
yourself, call GCM in a try-catch block.
> In mapdisp/private/warnObsoleteMSGSyntax at 6
In gcm at 20
Ifthereisnomapaxes,itresultsin
Warning: Function GCM no longer returns error message strings in
output argument MSG. Inst ead any errors are thrown where they occur.
You should remove the la st output argument (MSG) from your call to
Version 2.7 (R2008a) Mapping Toolbox™ Software
GCM in order to avoid this warning. If you want to handle errors
yourself, call GCM in a try-catch block.
> In mapdisp/private/warnObsoleteMSGSyntax at 6
In gcm at 20
??? Error using ==> gcm> chec kaxes at 41
No axes in current figure.
Select a figure with map axes or use AXESM to define o ne.
Error in ==> gcm at 24
h = checkaxes(varargin{:});
Compatibility Considerations
As suggested by this warning, if you have any scripts or functions of your own
that depend on the old syntax, you should remove the
the function call in a
try-catch block instead.
msg argument and place
33
Mapping Toolbox™ Release Notes
Version 2.6 (R2007b) Mapping Toolbox Software
This table summarizes what’s new in Version 2.6 (R2007b):
New
Features
and
Changes
Yes
Details b elo w
• “Exporting Vector Geodata to Earth Browsers” on page 35
• “Improved Conversion Between Angle Units” on page 36
• “ImprovementsinHandlingLengthUnits”onpage38
• “New Angle Wrapping Functions” on page 39
• “New Function to Unwrap Sequences of Angles” on page 39
• “Improvements to Data Trimming in patchm and patchesm” on page 40
• “Higher Quality boston.tif GeoTIFF Satellite Image” on page 41
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of new features
and changes,
below. See also
Summary.
Fixed Bugs
and Known
Problems
Bug reports at
Web site
Related
Documentation
at Web Site
Printable Release
Notes: PDF
Current product
documentation
34
• “Map Axes Now Display Transparent Objects More Easily” on pag e 42
• “The arcgridread Function Now Imports Noninteger Data Grids” on page 43
• “Change to avhrrlam bert Function Behavior When No Data Is Available”
on page 43
• “Enhancements to Mapping Toolbox User’s G uid e” on page 43
• “Functions deg2rad and rad2deg No Longer Convert Complex to Real”
on page 43
Version 2.6 (R2007b) Mapping Toolbox™ Software
• “Degrees-Minutes-Seconds Conversion Functions Are Obsolete and Error
if Used” on page 44
• “Time Conversion Functions Are Obsolete and Error if Used” on page 47
• “cmapui GUI is now Obsolete” on page 47
• “Functions tigermif and tigerp Are Obsolete and Erro r if Used” on page 48
Exporting Vector Geodata to Ear th Browsers
kmlwrite is a new function for exporting vector point data to a file in KML
format. KML stands for Keyhole Markup Language; it is an XML dialect
used to structure geographic data for display in an Earth browser, such as
Google™ Earth, Google Maps, and Maps for Google Mobile. KML has a
hierarchical structure of nested elements and attributes.
simpleAPIthatletsMappingToolboxuserswritevectordatatoaKMLfilein
order to subsequently display the data onto an Earth browser.
kmlwrite has a
When used with Google Earth, files output from
kmlwrite can be seen
immediately in Google Earth, if that application is available to the user.
If the files are uploaded to a publicly accessible Web server, they can be
viewed by anyone on the Internet via Google Maps or other Web sites and
browser utilities that can read and display KML files. Google Maps and
GoogleMapsformobiledonotsupporttherangeofKMLmarkupthatGoogle
Earth supports (for example, placemark locations m ust be specified to them
as coordinates, not as addresses). See the Google KML documentation at
http://code.google.com/apis/km l /d o cu mentation/mapsSupport.html for more
information.
kmlwrite accepts latitude and longitude point vectors, passed either in
geostructs or as column arrays. It also accepts addresses, w hich can be
as general as a country’s name or as specific as a street address. When
geostructs are the input, the attribute data in the geostruct can be formatted
as HTML tables and included in the KML output. When latitude-longitude
arrays are input, you can pass attributes to
kmlwrite with strings. When
addresses are the input, geostructs are not used.
To customize p lacemarks, you can control formatting of geostruct attributes in
the KML file with an attribute specification, a struct used to format them (for
example, to add units to length attributes or to control the number of decimal
places for numeric values). A new support function,
makeattribspec lets you
35
Mapping Toolbox™ Release Notes
change the names used as labels in placemarks (geostruct field names are
used by default), omit fields from placemarks, and add HTML markup to the
attributes displayed in placemark tables.
See “Exporting Vector Geodata” in the Mapping Toolbox User’s Guide and
the mapexkmlexport demo, “Exporting Vector Point Data to KML” for more
information.
Improved Conversion Between Angle Units
The ang ledi m function has been replaced by four, more specific, functions:
in “Four New Angle-Unit Conversion Functions” on page 37). However,
angledim has been retained in Version 2.6 for backward compatibility.
The functions
alternatives.
degtorad, radtodeg,andunitsratio provide additional
Because it must resolve both the input and output units,
excessive for most applications. It works only for class
angledim is
double and it quietly
discards the imaginary part of any complex input. You can use any of several
more efficient alternatives:
If you are working from the command line, you can often replace
angledim
with degtorad or radtodeg. If you are converting angle units within a script
or function and you know both the
coding, then you can also replace
If you know either
Also note that the functions in the fromRadians family can convert multiple
variables in a single function call. For example, you can replace this code
angle1 = angledim(angle1InRadians,'radians',to);
36
Version 2.6 (R2007b) Mapping Toolbox™ Software
angle2 = angledim(angle2InRadians,'radians',to);
with
[angle1,angle2] =
fromRadians(
to,angle1InRadians,angle2InRadians);
If you do not know either f rom or to atthetimeofcoding,thenyoucancall
unitsratio to obtain the correct conversion factor, then multiply the values
of one or more variables. For example, you can replace:
r = unitsratio(to, from);
angle1Out = r * angle1In;
angle2Out = r * angle2In;
Four New Angle-Unit Conversion Functions
The following functions have been added for efficient conversion of angle
units (degrees or radians) when either the target or destination units (but not
both) are unknown before run time.
•
toDegrees — Convert angles to degrees
toRadians — Convert angles to radians
•
fromDegrees — Co nvert angles from deg rees
•
fromRadians — Convert angles from radians
•
If the output units match the inputs units, as in
angle2,...)
,whereunits turns out to equal ' degr ees', then the input
angles are simply copied to the output angles.
Use these functions in place of
angledim. The new functions are more
efficient, especially when the value of either the
angledim is known in advance and the value of the other angle-unit argument
is not.
toDegrees(units, angle1,
from or to argument of
37
Mapping Toolbox™ Release Notes
Improvements in
Handling Length Units
Alternatives to the distdim Function
There are now mo
than the
distd
compute multi
units of dista
For other alte
page for deta
re efficient ways to convert length and distance units
im
function. In place of distdim,youcanuseunitsratio to
plicative factors to apply when converting between different
nces and angles, which you can use in subsequent calculations.
rnatives, see “Replacing distdim” in the
ils.
The unitstr function Is Obsolete
The unitstr
distance, a
The syntax
Instead, s
Compatib
unitsrat
function, which validates names and abbreviations for units of
ngle, and time, is obsolete and will be removed in a future release.
str = unitstr(str,'times') has already been removed.
ee the documentation for
unitsratio for a list of val id unit strings.
ility Considerations. Thereisnoreplacementfor
io
recognizes all the unit strings known to the toolbox.
Interpretation of “Miles” Units has Changed
As of R20
'miles'
07b, the following functions interpret distance units specified as
as International Miles,notStatute Miles:
distdim reference
unitstr,but
38
•
almana
• daspec
• eleva
• mappr
• pape
• scal
rscale
eruler
c
tm
tion
ofile
Version 2.6 (R2007b) Mapping Toolbox™ Software
Compatibility Conside rations. T his will not materially affect the
accuracy of results in most cases; the lengths of the two types of miles only
differ by about two parts per million (three millimeters). The
distdim
function’s interpretation of miles has not changed. However, there are better
alternatives to it; see the release note “Alternatives to the distdim Function”
on page 38.
New Angle Wrapping Functions
Four new low-level functions have been added that force longitudes, azimuths,
orphaseanglestospanintervalsof
2*pi]
•
•
•
•
or [-pi pi] radians.
wrapTo180 — Wrap angle in degrees to [-180 180]
wrapTo360 — Wra p angle in degrees to [0 360]
wrapToPi — Wrap angle in radians to [-pi pi]
wrapTo2Pi — Wra p angle in radians to [0 2*pi]
[0 360] or [-180 180] degrees or [0
The first two functions work in degrees, the next two in radians. None of
them perform argument checking.
You can use the new wrapping and functions in place of
npi2pi and zero22pi
for greater efficiency. The older functions will eventually be removed from
the toolbox.
New Function to Unwrap Sequences of Angles
The new unwrapMultipart function unwraps vectors of angles similarly to
the MATLAB function
NaN separators, unwrapping each section separately. Use it to remove
discontinuities from vectors of longitudes, azimuths, or phase angles that
contain NaN-delimited sequences and as a replacement for the obsolete
function
smoothlong.
unwrap,exceptthatithandlesvectorsthatinclude
39
Mapping Toolbox™ Release Notes
Improvements to
Data Trimming in patchm and
patchesm
The patchm and pa
parts of polygo
also affects
f
shifted coord
where they wou
map frame. Thi
constructed
use of
by
AlphaD
Now Display T
tchesm
ns th a t fall ou tside your current map limits. This improvement
illm
, which calls patchm. Previously the patch functions simply
inates inward so that vertices collected at the edge of the limits,
ld appear as lines along map borders, unless obscured by the
s change allows OpenGL to better render the patch objects
patchm and patchesm, making them more compatible w ith the
ata
to achieve transparency. Seethereleasenote“MapAxes
ransparent Objects More Easily” on page 42 for more details.
Compatibility Considerations
Themoreco
circumsta
the data pr
objects c
display f
whenever
map limi
will bec
are trim
separat
will be
map lim
by
patc
map ax
mplete trimming in
nces under which automatic reprojection can no longer display all
ovided to these functions. Automatic reprojection causes map
reated with
unctions (but not
acallto
setm changes certain map axes properties, including the
ts and projection type. In the case of
ome unavailable for automatic reprojection if all of the polygons
med away completely. In the case of
e object for each polygon, any polygon that is trimmed away completely
unavailable for reprojection, even if it would lie within newly defined
its. In either of these cases, you should delete the handle(s) returned
hm
or patchesm, then repeat the original calls after changing your
es properties.
functions now completely trim away polygons and
patchm and patchesm means that there are
plotm, linem, patchm, patchesm, and certain other
geoshow) to be remov ed, projected, and redisplayed
patchm, a set of polygons
patchesm, which constructs a
40
Other
•
• You m
See
Eas
potential compatibility issues:
hm
patc
poly
and patches m exhibit greater sensitivity to incomplete or incorrect
gon topology.
ight need to manually set the renderer f or proper display of some
ch data
pat
the release note “Map Axes Now Display Transparent Objects More
ily” on page 42 for information about rendering and the Mapping Toolbox
Version 2.6 (R2007b) Mapping Toolbox™ Software
demo mapexgshhs for an example of a situation where polygon topology
necessitates manual setting of the renderer.
Higher Quality boston.tif GeoTIFF Satellite Image
The original boston.tif GeoTIFF satellite image has been replaced by a
higher resolution image, created by and provided courtesy of GeoEye™.
Thenewimagehasthesamenameastheoldone,
boston.tif file, and an ov erview image in JPEG format, boston_ovr.jpg,
latitude-longitude, with a ground pixel size of approximately 16 meters. For
further information, refer to the text files
boston_metadata.txt in toolbox/map/mapdemos.
boston.txt, boston_ovr.txt,and
boston.tif. The new
Compatibility Considerations
Older satellite images of Boston and a demo have been removed from
Mapping Toolbox directories. The new
images replace the images having the same names previously included
in
toolbox/map/mapdemos. In addition, several older images related to
boston.tif have been removed:
•
boston_red.tif
• boston_green.tif
• boston_blue.tif
• boston_pan.tif
• boston_enhanced_pan.tif
boston.tif and boston_ovr.jpg
41
Mapping Toolbox™ Release Notes
The mapexenhance demo (“Enhancing Multispectral GeoTIFF Images”), which
used several of these images, has also been removed.
Map Axes Now Display Transparent Objects More
Easily
It is now much easier to achieve transparency effects from the toolbox by
setting the
lightm, contourm,andcontour3m se t the figure’s Renderer property: axesm
and lightm set it to 'zbuffer',whilecontourm and contour3m set it to
'painters'. You then had to manually reset Renderer to 'opengl' in order
for transparency to take effect.
AlphaData property of an object. Previously, functions axesm,
Now the
of
'auto', causing M ATLAB to select the most appropriate renderer for
you; it will use OpenGL when appropriate, given your
RendererMode of the figure retains the default MATLAB value
AlphaData settings.
Using OpenGL not only enables transparency effects, it also can make use of
hardware graphics acceleration capabilities should they be available.
Compatibility Considerations
If you need a particular map display to look the same as it did in Mapping
Toolbox Version 2.5 (R2007a), in most cases you can just issue the command
set(gcf,'Renderer','zbuffer')
after you construct your map axes. If you are calling contourm or contour3m,
issue the command
set(gcf,'Renderer','painters')
after you call the contouring function.
Theconsequenceofdoingthisisthatyouwillnotbeabletousetransparency
with that map figure until you reset its renderer to
'RendererMode' back to 'auto', which is its default state.
'opengl' or set its
42
Version 2.6 (R2007b) Mapping Toolbox™ Software
The arcgridread
Function Now Imports Noninteger
Data Grids
In pr evious rele
that had intege
been removed,
notation.
Change to avh
Data Is Avail
In previous
quadrangle
in the appro
intersect t
release,
matrices.
a
Compatibility Considerations
If you dep
quadran
gle, you will need to change your code.
Enhance
Severa
been re
topics
and th
new se
use th
l sections of the chapter “Understanding G eo spatial Geometry” have
written and new material has been added to better explain critical
such as ellipsoid models, units of and notations for angles and length,
e conversions that are possible between various units. There is also a
ction, “Exporting Vector Geodata”, explaining and illustrating how to
enew
ases of the toolbox,
r values (often of meters or feet). This limitation has now
such that input grids can contain arbitrary values in decimal
rrlambert Function Behavior When No
able
releases of the toolbox,
defined by
priate Lambert Equal Area Azimuthal projection) failed to
he bounding box of the data in the projected coordinates. In this
vhrrlambert
end on
latlim and lonlim (when projected to form a polygon
does not error when this occurs but returns empty
avhrrlambert to error when there is no data in your
ments to Mapping Toolbox User’s Guide
kmlwrite and makeattribspec functions.
arcgridread could only import data grids
avhrrlambert would error if the
Func
Comp
In p
iss
lon
fun
Fo
tions deg2rad and rad2deg No Longer Convert
lex to Real
rior versions, when given complex inputs, functions
ued a warning and then converted their inputs to real. Now they no
ger do either of these things; in the unlikely event of complex input, these
ctions simply scale the imaginary part by the same factor as the real part.
r example, in R2007a and earlier releases, they behave as follows:
deg2rad and rad2deg
43
Mapping Toolbox™ Release Notes
>> deg2rad(180i)
Warning: Imaginary parts o f complex ANGLE argument ignored
> In deg2rad at 16
ans =
Going forward from this release, the result is
>> deg2rad(180i)
ans =
Degrees-Minutes-Seconds Conversion Functions Are
Obsolete and Error if Used
The following functions, which accepted or produced double scalars to
represent degrees, minutes, and seconds no w error when used, and will b e
removed completely from the toolbox in a future release:
0
0 + 3.1416i
44
•
deg2dm
• deg2dms
• dms2deg
• dms2dm
• dms2mat
• dms2rad
• mat2dms
• rad2dm
• rad2dms
The scalar DM and DMS encodings are being eliminated from the toolbox
because they were never used for internal computations, and alw ay s had
the potential to generate serious numerical errors if passed accidentally to
functions that expected normal latitude-longitude tuples. They also made the
functions that accepted them less efficient due to the need to convert from DM
or DMS to fractional latitudes and longitudes before processing the input data.
Version 2.6 (R2007b) Mapping Toolbox™ Software
In every case, an alternative that does not use the old degrees-minutes-seconds
scalar encoding exists. See the following section on compatibility for
replacements and “New Functions for Degrees-Minutes-Seconds Conversions”
on page 51 in the V2.5 Release Notes for descriptions of replacement
functions, and the compatibility considerations below for descriptions of
alternative syntaxes and expressions you can use for degrees-minutes-seconds
conversions.
Compatibility Considerations
DM and DMS representations are widely used in published reports and
can occur in geodata that you want to read into the MATLAB workspace.
You can still import and export DM and DMS data, but Mapping Toolbox
functions no longer accepts the old encodings as alternatives to floating-point
representations of latitude and longitude for internal manipulations.
The following functions (which all use scalar DMS encoding) are being re tired.
They remain in the product for R2007b, but now generate errors when used.
They will be removed completely in the next version. Use the alternative
suggested in lieu of these functions.
•
deg2dm — Instead use degrees2dm to convert degrees to degrees-minutes
vector.
•
deg2dms —Insteadusedegrees2dms to convert d egre es to
degrees-minutes-seconds vector.
•
dms2deg — Instead use dms2degrees to convert degrees-minutes-seconds
vector to degrees.
•
dms2dm —Insteadcombinedms2degrees and degre es2dm,asin
degrees2dm(dms2degrees([-29 42 18.7])) to remove the seconds
component from a degree-minutes-second v ector.
•
dms2mat —Insteadusedegrees2dms to convert d egre es to
degrees-minutes-seconds vector.
•
dms2rad — Instead use dms2degrees to convert degrees-minutes-seconds
vector to degrees and call
mat2dms — Instead use dms2degrees to convert degrees-minutes-seconds
•
degtorad or multiply by pi/180.
vector to degrees.
45
Mapping Toolbox™ Release Notes
• rad2dm — Instead, call radtodeg or multiply input arguments by 180/pi,
and then call
rad2dms —Instead,callradtodeg or multiply input arguments by 180/pi,
•
and then call
degrees2dm.
degrees2dms.
In addition, the
axesm and setm functions no longer accept the strings 'dms'
and 'dm' for setting either the AngleUnits or LabelUn its properties of
amapaxes.
Many other Mapping Toolbox functions optionally accept angle strings for
their
units parameter; the following 57 functions now only accept 'degrees'
and 'radians', whereas in prior versions they would also accept 'dm' and
'dms' as values for units:
angl2strdistancehistrputpolestdist
angledimeastofimbedmrcurvestdm
antipodeelevationinterpmreckontimezone
areaintellipse1intrplatrhxrhtrack
areamatepsmintrplonrotatemtrack1
areaquadeqa2grnmapprofilerspheretrack2
axesmgc2scmeanmscalerulerunitstr
azimuthgcxgcmeshgratscircle1westof
convertlatgcxscneworigscircle2zero22pi
crossfixgradientmnewpolescxsc
daspectmgrn2eqanpi2pisetm
departurehistaorg2polsmoothlong
46
These functions now error when provided 'dm' or 'dms' for their units
argument.
Version 2.6 (R2007b) Mapping Toolbox™ Software
Time Conversion Functions Are Obsolete and Error if
Used
The following functions, which converted time representations, now error
when used and will be rem ove d completely from the toolbox in a future release:
•
hms2hm
• hms2hr
• hms2mat
• hms2sec
• hr2hm
• hr2hms
• hr2sec
• mat2hms
• sec2hm
• sec2hms
• sec2hr
• time2str
• timedim
Compatibility Considerations
These functions now raise errors when they are invoked. They will be
completely removed in a future versionofthetoolbox. Nosubstitutes have
been provided, as no operations of the toolbox have ever depended on them.
cmapui GUI is now Obsolete
The cmapui GUI now errors if you attempt to use it. It will be completely
removed from the next Mapping Toolbox version. Use the MATLAB
colormapeditor GUI instead, which provides better functionality. You can
also use the Colormap drop-down menu in the Property Editor (part of
the MATLAB plotting tools and available via the
select a built-in colormap; the
colormapeditor. To set up a colormap for terrain displays, you can use the
custom option on that drop-down menu ope ns
propedit command) to
47
Mapping Toolbox™ Release Notes
demcmap function. To generate an appropriate (but random) colormap for
political maps, use the
Functions tigermif and tigerp Are Obsolete and Error
if Used
The following functions error and issue an error message when you attempt
to use them:
•
tigerp — Read TIGER p and pa thinned boundary files (ArcInfo format)
• “New Functions for Degrees-Minutes-Seconds Conversions” on page 51
• “Time Conversion Functions to be Removed” on page 54
Version
Compatibility
Considerations
Yes — Details
labeled as
Compatibility
Considerations
in descriptions
of new features
and changes,
below. See also
Summary.
Fixed Bugs
and Known
Problems
Bug reports at
Web site
Related
Documentation
at Web Site
Printable Release
Notes: PDF
Current product
documentation
Performance Improvements for los2 and viewshed
This release includes a faster los2 function (which com p utes intervisibility
between locations on or above a terrain grid). The
computes the portions of a terrain gridthatcanbeseenfromagiven
viewpoint) has also been accelerated as a result.
viewshed function (which
49
Mapping Toolbox™ Release Notes
Utility Functions for Computing Distance and Position
Along Meridians
Two functions that reckon position and distance along a meridian on the
ellipsoid are now available:
•
meridianarc — Computes distance along a meridian between two latitudes
meridianfwd — R eckons position along meridian given a starting point
•
and distance
Some GUIs Are No Longer Available from the
Command Line
In prior releases, when you typed certain Mapping Toolbox function names
with no argument list, a specialized GUI appeared that enabled you to
interactively set parameters related to the function. This feature was seldom
used and sometimes raised errors when users attempted to operate the GUIs.
Starting in this release, a GUI will no longer appear when you issue the
following commands:
50
•
comet3m
• cometm
• contourfm
• contour3m
• contourm
• demcmap
• fill3m
• fillm
• lightm
• limitm
• linem
• meshlsrm
• meshm
• patchesm
• patchm
• pcolorm
• plot3m
• plotm
• quiver3m
• quiverm
• scatterm
• stem3m
• surfacem
• surflm
• surfm
• surflsrm
• symbolm
Version 2.5 (R2007a) Mapping Toolbox™ Software
• textm
Compatibility Considerations
Use the above functions with arguments to avoid raising errors. Their GUIs
will continue to be available via
containing map axes), but they are not being actively supported and will be
eliminated in a future release.
maptool (which places menus on a figure
New Functions for Degrees-Minutes-Seconds
Conversions
Four new functions have been added to convert to and from decimal degrees
and degrees-minutes-seconds (DMS):
•
dms2degrees — Convert degrees-minutes-seconds to degrees
dm2degrees — Convert degrees-minutes to degrees
•
51
Mapping Toolbox™ Release Notes
• degrees2dms — Convert degrees to degrees-minutes-seconds
degrees2dm — Convert degrees to degrees-minutes
•
The DMS inputs and outputs of these functions are vectors of one row and
three columns for each row in the decimal degrees input or output. The first
column contains the “degrees” element and is integer-valued. The second
column contains the “minutes” element and is integer-valued. The third
column contains the “seconds” element, and may have a nonzero fractional
part. Similarly, DM inputs and outputs are two-column vectors with integer
degrees and fractional minutes parts.
The new conversion functions dispense with the DMS encoding used in prior
versions of the toolbox. These represented DMS angles by a single real
number, the format of which is
internally by M apping Toolbox functions, as it is not self-documenting and
can lead to erroneous computations. For example, two DMS-encoded real
numbers cannot be added to obtain a meaningful result.
dddmm.ss. Such an encoding is no longer used
52
Compatibility Considerations
DM and DMS representations are widely used in published reports and
can occur in geodata that you want to read into the MATLAB workspace.
You can still import and export DM and DMS data, but Mapping Toolbox
functions no longer accepts the old encodings as alternatives to floating point
representations of latitude and longitude for internal manipulations.
The scalar DM and DMS encodings are being eliminated from the toolbox
because they were never used for internal computations, and always had the
potential to generate serious numerical errors if passed accidently to functions
that expected normal latitude-longitude tuples. They also made the functions
that accepted them le ss efficient d ue to the need to convert from DM or DMS
to fractional latitudes and longitudes before processing the input data.
The following existing functions (which all use scalar DMS encoding) are
being retired. They remain available but now issue warnings that they are
obsolete when used:
•
deg2dm — Instead use degrees2dm to convert degrees to degrees-minutes
vector
Version 2.5 (R2007a) Mapping Toolbox™ Software
• deg2dms —Insteadusedegrees2dms to conv ert degrees to
degrees-minutes-seconds vector
•
dms2deg — Instead use dms2degrees to convert degrees-minutes-seconds
vector to degrees
•
dms2mat —Insteadusedegrees2dms to convert d egre es to
degrees-minutes-seconds vector
•
dms2rad — Instead use dms2degrees to convert degrees-minutes-seconds
vector to degrees and call
degtorad or multiply by pi/180
• mat2dms — Instead use dms2degrees to convert deg rees -minute s-seconds
vector to degrees
•
rad2dm — Instead, call radtodeg or multiply input arguments by 180/pi,
and then call
degrees2dm
• rad2dms —Instead,callradtodeg or multiply input arguments by 180/pi,
and then call
degrees2dms
In addition, the axesm and setm functions no longer accept the strings 'dms'
and 'dm' for setting either the AngleUnits or LabelUn its properties of
amapaxes.
Many other Mapping Toolbox functions optionally accept angle strings for
their
units parameter; the following 57 functions now only accept 'degrees'
and 'radians', whereas in prior versions they would also accept 'dm' and
'dms' as values for units:
angl2strdistancehistrputpolestdist
angledimeastofimbedmrcurvestdm
antipodeelevationinterpmreckontimezone
areaintellipse1intrplatrhxrhtrack
areamatepsmintrplonrotatemtrack1
areaquadeqa2grnmapprofilerspheretrack2
axesmgc2scmeanmscalerulerunitstr
azimuthgcxgcmeshgratscircle1westof
53
Mapping Toolbox™ Release Notes
convertlatgcxscneworigscircle2zero22pi
crossfixgradientmnewpolescxsc
daspectmgrn2eqanpi2pisetm
departurehistaorg2polsmoothlong
These functions now issue warnings when provided 'dm' or 'dms' for their
units argument.
Time Conversion Functions to be Re m oved
The following functions to convert between time units and encodings will be
removed from a future release of the toolbox:
•
hms2hr
• hms2hm
• hms2mat
54
• hms2sec
• hr2hm
• hr2hms
• hr2sec
• mat2hms
• sec2hm
• sec2hms
• sec2hr
• time2str
• timedim
Compatibility Considerations
These functions remain available, but when they are invoked now issue
warnings that they are obsolete.
Version 2.4 (R2006b) Mapping Toolbox™ Software
Version 2.4 (R2006b) Mapping Toolbox Software
This table summarizes what’s new in Version 2.4 (R2006b):
New
Features
and
Changes
Yes
Details b elo w
New features and changes introduced in this version are
• “Standard Formulations of Five Major Map Projections” on page 55
• “Two New Geodetic/Geocentric Latitude Conversion Functions” on page 56
• “Accelerated Performance for geoshow, mapshow, and bufferm” on page 56
• “Changes in Behavior of mapshow and geoshow” on page 56
• “dted Automatically Fixes Incorrectly Specified Longitude Directions in
DTED Data” on page 60
Version
Compatibility
Considerations
Yes — Details
in “Changes
in Behavi or
of mapshow
and geoshow”
on page 56
below. See also
Summary.
Fixed Bugs
and Known
Problems
Bug reports at
Web site
Related
Documentation
at Web Site
Printable Release
Notes: PDF
Current product
documentation
Standard Formulations of Five Major Map Projections
New formulations of five conic map projections are provided. The existing
implementations remain available under their old names. The new versions
use the same names as the ones they supplement, appended with “std”:
• Cassini Transverse Cylindrical (
• Albers Equal-Area Conic (
• Equidistant Conic (
• Lambert Conformal Conic (
eqdconicstd)
cassinistd)
eqaconicstd)
lambertstd)
55
Mapping Toolbox™ Release Notes
• Polyconic (polyconstd)
Computations used for the new versions differ from the old ones only when the
latitude origin (the first element of the origin vector)isnonzero.Inthiscase,
the old versions shift the origin off the equator through a solid body rotation
of the sphere (or, for an ellipsoidal earth model, a suitable auxiliary sphere).
This is technically correct, but differs from accepted industry standards
for these projections. The new versions use the standard formulations and
give results that are consistent with projection results from other software
packages, regardless of the latitude origin. The old versions are retained in
the toolbox, with no change in behavior, to ensure backward compatibility.
See the Projections Reference documentation for more information.
Two New Geodetic/Geocentric Latitude Conversion
Functions
Two new functions provide a more direct route to functionality already
available via the
convertlat function:
56
•
geocentric2geodeticlat converts an array of geocentric latitude in
radianstogeodeticlatitudeinradianson a reference ellipsoid given a first
eccentricity
•
geodetic2geocentriclat converts an array of geodetic latitude in radians
to geocentric latitude in radians on a reference ellipsoid given a first
eccentricity.
Accelerated Performance for geoshow, mapshow,
and bufferm
Functions geoshow, mapshow,andbufferm run substantially faster in many
cases, especially when vector display is being controlled via symbol specs in
mapshow and geoshow.
Changes in Behavior of mapshow and geoshow
In addition to operating faster, the mapshow and geoshow functions now
behave slightly differently regarding their defaults, handles returned,
warnings issued, and several other aspects:
Version 2.4 (R2006b) Mapping Toolbox™ Software
Default Symbols and Colors
• Point marker type changes from 'X' to '+'
• Point marker color changes from 'black' to 'red'
• Line color changes from 'black' to 'blue'
• Polygon facecolor changes from 'black' to pale yellow
Polygon edgecolor remains
'black'
Contour DisplayType Behavior Changes
• The DisplayType option 'contour' now returns an hggroup handle. The
children of the hggroup are patches. In prior versions, an array of line
handles was returned.
• You can specify any contourgroup property as a parameter value pair. In
previous versions,
no other contour properties.
• Both
• In R2006b, when plotting contours on a regular axes (not a map axes),
mapshow and geoshow might return a different number of contour
levels by default than in previous versions, in which you could not specify
contour intervals; in R2006b, you can control contour intervals and levels
via the
geoshow projects the contour lines using a Plate Carree projection; in
previous versions it simpl y displayed longitudes as x and latitudes as y
without doing any trimming or longitude wrapping.
LevelStep or LevelList contourgroup properties, among others.
mapshow allowed you to set the LineStyle property, but
Graphic Objects and Return Values for Vector Inputs
• Vector coordinate array input (x-y or latitude-longitude pairs) with a
DisplayType of 'Line' or 'Point' now generates an ordinary line object
insteadofamapgraphicsline.
• For geostruct input, an hggroup object is constructed; its handle is returned
instead of an array of handles to map graphic objects:
- For polygon geostructs, map graphics polygon objects are still
constructed, but become children of the hggroup.
57
Mapping Toolbox™ Release Notes
- For point, multipoint, and line geostructs, the children of the hggroup
In both cases each child of the hggroup, rather than each element in an
array of handles, corresponds to a distinct feature in the geostruct.
Handles Returned for G raphic Objects
• Geostruct inputs result in an hggroup handle containing either line objects
(for point, multipoint, and line inputs) or modified patch objects (for
polygon inputs) as their children.
• Coordinate arrays (x,y pairs) displayed as lines now result in ordinary
line objects.
• Geostructs containing lines result as ordinary line objects within hggroups.
New Warnings Issued
• mapshow and geoshow now warn when given a geostruct within which the
Geometry field differs from a specified 'Dis playType' parameter.
are ordinary line objects; map graphics objects are no longer constructed.
58
mapshow will warn if it is given a geostruct containing Lat and Lon fields
•
instead of
geoshow will warn if it is given a geostruct containing X and Y fields instead
•
of
Lat and Lon fields.
X and Y fields.
geoshow Supplies Default Projection
geoshow now p rojects vector and raster in puts using a default projection
(Plate Carree) if the parent axes is not a map axes. The axes itself is
unchanged (it is not modified to become a map axes), but the scale factor of
the projection is set such that latitudes and longitudes in degrees can be read
directly from the axes ticks and grid lines.
Duplicate Parameter/Value Pair Inputs Allowed
mapshow no longer errors or warns if giv en duplicate Parameter/Value pair
inputs; in such circumstances,
SymbolSpecs)
mapshow now u ses the last value (even with
Version 2.4 (R2006b) Mapping Toolbox™ Software
geoshow Supports True Surface Display
geoshow now creates a true 3-D surface if given a 'surface' DisplayType
rather than setting the ZData values to 0.
Texturemap DisplayType Behavior Changes
The 'texturemap' DisplayType now uses the pixel edges to create XData and
YData grids rather than using the pixel centers, which correctly registers the
display to map coordinates. The
same dimensions as the
XData and YData arrays, which exceed the input
grid in size by one in both the x and y dimensions.
ZData contains an array of zeros having the
You should use
coded by color (i.e. 2-D displays); use
show data w ith relief (nonzero
'texturemap' displays when the a ttribute being displayed is
'surface' displays when you n eed to
ZData).
More General Support for Graphics Properties
All Handle Graphics®patch properties are now supported for polygon inputs.
All Handle Graphics line properties are now supported for point and line
inputs, except that
'linestyle' is ignored for point inputs.
Limitations on Referencing Matrices for Geoshow Removed
geoshow is now capable of accepting any referencing matrix. Previously
it could only accept those referencing matrices that were convertible to
referencing vectors.
mapshow and geoshow Ignore Empty Inputs Rather than
Erroring
In previous versions, mapshow and geoshow would throw errors when provided
with empty (
these functions via scripts. The new behavior is also more consistent with
that of MATLAB plotting functions such as
[]) arrays. This behavior could be inconvenient when running
plot, surf, mesh,andcontour.
59
Mapping Toolbox™ Release Notes
dted Automatically Fixes Incorrectly Specified
Longitude Directions in DTED Data
Some DTED level 0 files available via the Nation al Geospatial-Intelligence
Agency’s (NGA) web interface may have minor errors. Specifically, Level
0 data for cells just to the east of the prime meridian may have longitude
coordinate strings with ’W’ substituted for ’E’. The
and automatically corrects this data error.
dted function now detects
60
Version 2.3 (R2006a) Mapping Toolbox™ Software
Version 2.3 (R2006a) Mapping Toolbox Software
This table summarizes what’s new in Version 2.3 (R2006a):
New Features
and Changes
Yes
Details below
New features and changes introduced in this version are
• “Full Support for 64-Bit Windows” on page 61
• “Third-Party Library and Code Upgrades” on page 61
• “Support for 32-Bit Floating-Point GeoTIFF Images” on page 62
• “Utility Functions for NaN-Separated Polygons and Lines” on page 63
• “Standardized Vector Topology in coast.mat” on page 63
• “Three New Demos” on page 63
Version
Compatibility
Considerations
Yes — Details
in “GeoTIFF
Compatibility
Considerations”
on page 62
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Bug reports at
Web site
Related
Documentation
at Web Site
Printable
Release
Notes: PDF
Current product
documentation
Full Support for 64-Bit Windows
Version 2.3 adds support for the mex- and library-based functions
geotiffinfo, geotiffread, sdtsinfo,andsdtsdemread on this new
MATLAB platform via library upgrades (described below) and a custom port
of STDS++.
Third-Party Library and Code Upgrades
Third-party libraries and software packages have been upgraded to their
current versions to ens ure best performance and compatibility with external
geospatial data s ou rces and applications software:
61
Mapping Toolbox™ Release Notes
• General Polygon Clipper (GPC) upgraded to Version 2.32
• PROJ.4 library upgraded to Version 4.4.9
• SDTS++ library u pgraded to Version 1.5.1
• GeoTIFF library upgraded to Version 1.2.2
Support for 32-Bit Floating-Point GeoTIFF Images
The MATLAB function imread can now import TIFF images containing 32-bit
floating-point data. As a result,
variety of GeoTIFF.
GeoTIFF Compatibility Considerations
The structure returned by geotiffinfo in V. 2.3 has changed. The following
table describes the differences between the current and previous versions:
Version 2.3Previous Versions
geotiffread now reads the corresponding
The TiePoints structure contains
two substructures,
The CornerCoords structure
contains six
respectively,
Lon.
[1-by-4] row vectors,
X, Y, Col, Row, Lat,and
The Zone field contains [] if the
UTM zone is not applicable or was
missing from the metadata.
TiePoints structure contained
The
two [
3-by-1] arrays, ImageP oint s
and WorldPoints.
The
CornerCoords structure
contained six
vectors:
LON.
[4-by-1] column
PCSX, PCSY, X, Y, LAT,and
The Zone fie ld contained 32767 if the
UTM zone was not applicable or was
missing from the metadata.
62
Version 2.3 (R2006a) Mapping Toolbox™ Software
Utility Functions for NaN-Separated Polygons and
Lines
closePolygonParts
Closes all rings in a multipart polygon to ensure proper analysis and
rendering.
isShapeMultipart
Boolean-valued function that returns true if a polygon or line has multiple
parts.
removeExtraNanSeparators
Eliminates redundant NaN separators that might exist in polygons and lines.
Standardized Vector Topology in coast.mat
Polygons in the low-resolution coastline sample data file coast.mat now
follow the convention used by
polygons with “holes” (inner rings representing lakes, inland seas, and islands
within them). Outer contours now always run clockwise and inner contours
run counterclockwise. These edits, which reversed the order of vertices in
some rings, enable th e display functions to f ill outer rings properly while
leaving inner rings blank.
geoshow, mapshow,andmapview to display
Three New Demos
If you are view ing these release notes using the Help browser, clicking any of
thedemolinksbelowwillopenthedemoin a browser window. Click the links
at the top of that window to view o r run the code for the demo.
Converting Coastline Data (GSHHS) to Shapefile Format
Shows how to extract coastlines from the Global Self-consistent Hierarchical
High-resolution Shorelines (GSHHS) data set, manipulate the polygon
features, and save the result to a polygon shapefile.
63
Mapping Toolbox™ Release Notes
Plotting a 3-D Dome as a Mesh Over a Globe
Illustrates how to construct a 3-D feature in a system of local vertical
coordinates, then transform and combine it with a globe display in
Earth-Centered, Earth-Fixed (ECEF ) coordinates.
Unprojecting a Digital Elevation Model (DEM)
Shows how to unproject a georeferenced terrain elevation grid from Universal
Transverse Mercator (UTM) into a regular latitude-longitude grid having
comparable spatial resolution.
64
Version 2.2 (R14SP3) Mapping Toolbox™ Software
Version 2.2 (R14SP3) Mapping Toolbox Software
This table summarizes what’s new in Version 2.2 (R14SP3):
New Features
and Changes
Yes
Details below
New features and changes introduced in this version are
• “Geodetic-Geocentric Coordinate Conversion Functions” on page 65
• “Additional User Control Over Shapefile Content” on page 65
• “Shapefile Read/Write Efficiency E nhanced” on page 66
• “Improved Rendering of Polygons with Inner Rings” on page 66
• “Map Viewer Now Georeferences Images It Saves” on page 66
New three-dimensional coordinate conversion functions (geodetic2ecef,
ecef2geodetic, ecef2lv, lv2ecef) transform 3-D point locations between
geodetic (latitude, longitude, height), geocentric Cartesian (Earth Centered,
Earth Fixed), and local vertical Cartesian coordinate systems.
Additional User Control Over Shapefile Content
Function shapewrite now allows user control over field names, lengths, and
decimal precision w hen writing feature attributes to the DBF file, via a
“DBF specification.” The new function
specification from a geographic data structure (geostruct2) array. Users can
customize the output and pass it to
makedbfspec constructs a default DBF
shapewrite.
65
Mapping Toolbox™ Release Notes
Shapefile Read/Write Efficiency Enhanced
Improved implementations of functions shaperead and s hapewrite process
data substantially faster (about four times faster for a 10-MB shapefile of
major roads in Massachusetts).
Improved Rendering of Polygons with Inner Rings
The M ap Viewer (function mapview) and functions mapshow and geoshow now
properly render polygons containing inner rings (e.g., lakes and inland seas
within a continent, islands within a pond). Features in underlying layers
“show through” inner rings because they are not obscured by the patch faces
used to render the polygons.
Compatibility Considerations
Polygon Vertex Ordering Is Now Significant for Properly Rendering
Filled Polygons. The map display functions
now require that coordinate vectors representing polygons have consistent
directionality, such that
geoshow, m apsh ow,andmapview
66
• Vertices defining outer rings (to be filled) be encoded in a clockwise
direction.
• Vertices defining inner rings (often termed “lakes” or “islands,” to be
rendered as transparent holes) be encoded in a counterclockwise direction.
If you have vector map data sets that violate these conditions, the map display
functions
them as filled polygons. To determine the directionality of polygon vertices,
use the logical function
CDs or printed documentation should substantially mitigate recurrent
problems with stale links. Please report any stale links that you
might find in the technical note to MathWorks Technical Support
(
http://www.mathworks.com/contact_TS.html), so that it can be updated
promptly.
Compatibility Considerations for Atlas Data and
Associated Functions
Through Version 2.0.3, the toolbox include d a set of atlas data with global
geopolitical data embedded as MATLAB arrays in four MAT-files:
worldhi, worldmtx,andworldmtxmed. How ever, g eopolitical data is difficult
to keep current, and is subject to inaccuracies and interpretations that can
cause contention. Therefore, starting with Version 2.1, Mapping Toolbox demo
data now excludes geopolitical data that would specify national sovereignty
over specific regions of the Earth. The only exceptions are the boundaries of
the 50 U.S. states and the District of Columbia.
worldlo,
This change means that the
worldlo, worldhi, worldmtx,andworldmtxmed
MAT-files are no longer part of the toolbox. However, the nonpolitical data
on global coas tlines, major lak es and inland seas, major rivers, and major
cities and populated places that was in
worldlo.mat has been retained
in the toolbox and transformed into shapefile format. This includes the
addition of name attributes for many previously unnamed features. There
are four new shapefiles in this category:
worldrivers.shp,andworldc itie s.shp.
landareas.shp, worldlakes.shp,
For consistency, the atlas data for the United States that was originally
stored in the
usalo and usahi MAT-files has also been transformed, although
none has been removed. Thes e data sets now reside in the following
shapefiles and MAT-files:
and
greatlakes.mat.
usastatelo.shp, usastatehi.shp, conus.mat,
The toolbox originally included four functions dedicated to extracting data
from the atlas data MAT-files:
worldlo, world hi, usalo,andusahi.Withthe
data removal/transformation described above, these functions are no longer
needed and h av e been removed from the toolbox in Version 2.1. You can
easily access the new shapefiles using the
shaperead function, which includes
powerful and flexible options for selecting features and even controlling which
attributes are read. In addition, function
country2mtx, whose sole purpose
70
Version 2.1 (R14SP2) Mapping Toolbox™ Software
was to rasterize the country boundary polygons in worldlo.mat, has been
removed.
Related changes extend to the
worldmap function, which formerly combined
two purposes:
• Select an appropriate map projection and parameters with which to display
a given latitude-longitude area.
• Automatically display atlas data for that area.
In Version 2.1,
worldmap constructs a map axes object and can easily be followed with a
worldmap supports only the first of these actions. A call to
variety of Mapping Toolbox commands todisplaythemapdataofyourchoice.
Because the
usamap function is so similar to worldmap, corresponding changes
have been made there as well.
To help those who have relied heavily on
worldmap and usamap to plot base
maps w ith automatically selected vector map data, examples thro u ghout the
User’s Guide, reference pages, and MATLAB function help entries have been
updated to illustrate the new behavior of
worldmap and usamap,andtoshow
how to create maps including vector shapefile data layers. These examples
cover a wide variety of ways to read and subset data with
display data with
A good p lace to start is with the examples for the
geoshow and other Mapping Toolbox display functions.
worldmap and usamap
shaperead and
functions. Also see example code in “Changes to worldmap and usamap
Display Types” on page 72.
To help you update commands, scripts, and data for constructing and
maintaining base maps, a recently published technical note on the MathWorks
Web site provides links to data and documentation for many sources of vector
and raster digital map data that you can access over the Internet:
Compatibility Considerations for worldmap and
usamap
worldmap and usamap have been simplified to construct appropriate map axes
for a given area without displaying any map data.
71
Mapping Toolbox™ Release Notes
In all cases, map frames, latitude-longitude grid lines, meridian labels, and
parallel labels are turned on. You can use the following command sequence
to remove them:
framem off; gridm off; m labe l off; plabel off
Other changes include the following:
•
usamap now accepts two-letter U.S. Postal S ervice abbreviations for state
names (e.g., AL, AK, AR, etc.).
• The following input opti o n s are n ow obsolete (if used, a warning is issued):
- A first argument equal to 'lo' or 'hi'
- The regiononly and stateonly syntax: a state or country name with
- All type options: 'none', 'line', 'lineonly', 'patch', 'patchonly',
the string
'mesh', 'meshonly', 'dem', 'demonly', 'dem3d', 'dem3donly',
'lmesh3d', 'lmesh3donly', 'ldem3d',and'ldem3donly' (the new
behavior matches the
'only' appended
'none' option)
72
Changes to worldmap and usamap Display Types
As of this release, the worldmap and usam ap functions no longer supports
the
type input argument. This argument provided an easy way to control
display behavior.
The
type option in worldmap was a single argument that could be one of
the following strings:
'mesh', 'meshonly', 'dem', 'demonly', 'dem3d', 'dem3donly ', 'lmesh3d',
'lmesh3donly', 'ldem3d',and'ldem3donly'.Inusamap, type was a subset
of the above names (the 3-D options were not supported).
In the current release, the various
following a call to
Toolbox commands. The following table specifies how you can achieve the
effects of the old
'none', 'line', 'lineonly', 'patch', 'patchonly',
type display options can be simulated by
worldmap or usamap with an appropriate set of Ma p ping
worldmap type argument using such auxiliary meth od s:
Version 2.1 (R14SP2) Mapping Toolbox™ Software
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load topo
load topo
worldmap(topo, ...
topolegend, 'dem')
load topo
worlmdap(topo, ...
topolegend,
'demonly')
load topo
worldmap(topo, ...
topolegend,
'dem3d')
worldmap(topo, topolegend)
meshm(topo, topolegend)
demcmap(topo)
land = shaperead('landareas.shp',...
'UseGeoCoords', true);
geoshow([land.Lat], [land.Lon])
load topo
worldmap(topo, topolegend)
meshm(topo, topolegend)
demcmap(topo)
load topo
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 7.5*pba(3)/da(3);
daspect(da);
demcmap(topo)
land = shaperead('landareas.shp',...
'UseGeoCoords', true);
geoshow([land.Lat], [land.Lon])
73
Mapping Toolbox™ Release Notes
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load topo
load topo
worldmap(topo, ...
topolegend,
'dem3donly')
load korea
worldmap (map,
refvec,...
'mesh')
load korea
worldmap(map,
refvec,...
'meshonly')
load topo
worldmap(topo,...
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),topo)
demcmap(topo)
load korea
worldmap(map, refvec)
meshm(map, refvec)
land = shaperead('landareas.shp',...
'UseGeoCoords', true);
geoshow([land.Lat], [land.Lon])
(Text North Korea and South Korea will be
missing)
load korea
worldmap(map, refvec)
meshm(map,refvec)
load topo
worldmap(topo, topolegend)
74
topolegend,
'mesh3d')
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 7.5*pba(3)/da(3);
daspect(da);
Version 2.1 (R14SP2) Mapping Toolbox™ Software
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load topo
load topo
worldmap...
(topo,topolegend,...
'mesh3donly')
load topo
worldmap(topo,...
topolegend,
'ldem3d')
load topo
worldmap(topo, ...
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),top
load topo
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 7.5*pba(3)/da(3);
daspect(da);
demcmap(topo)
camlight(90,5);
camlight(0,5);
lighting phong
material([0.25 0.8 0])
load topo
worldmap(topo, topolegend)
topolegend,'ldem3donly')
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 7.5*pba(3)/da(3);
daspect(da);
demcmap(topo)
75
Mapping Toolbox™ Release Notes
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load topo
load topo
worldmap(topo,...
topolegend,
'lmesh3d'
load topo
worldmap(topo,...
topolegend,...
'lmesh3donly')
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 2*pba(3)/da(3);
daspect(da);
camlight(90,5);
camlight(0,5);
lighting phong
material([0.25 0.8 0])
load topo
worldmap(topo, topolegend)
meshm(topo, topolegend,size(topo),topo)
da = daspect;
pba = pbaspect;
da(3) = 2*pba(3)/da(3);
daspect(da);
76
Version 2.1 (R14SP2) Mapping Toolbox™ Software
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load korea
load korea
worldmap(map,
refvec,...
load korea
worldmap(map,
refvec,...
'lineonly')
load korea
'line')
worldmap(map,
refvec,...
worldmap(map, refvec)
land =...
shaperead('landareas.shp',...
'UseGeoCoords', true);
geoshow([land.Lat], [land.Lon])
(Text North Korea and South Korea will be
missing. Land area boundaries resolution is
lower.)
load korea
worldmap(map, refvec)
land = ...
shaperead('landareas.shp',...
'UseGeoCoords', true);
geoshow([land.Lat], [land.Lon])
load korea
worldmap(map, refvec)
'none')
77
Mapping Toolbox™ Release Notes
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
load korea
load korea
worldmap(map,
refvec,...
'patch')
load korea
worldmap(map,
refvec,...
'patchonly')
worldmap(map, refvec)
land = ...
shaperead('landareas.shp',...
'UseGeoCoords', true);
faceColors = ...
makesymbolspec('Polygon',...
{'INDEX', [1 numel(land)],...
'FaceColor', ...
polcmap(numel(land))});
geoshow(land,'SymbolSpec',...
makesymbolspec('Polygon',...
faceColors)
(Text North Korea and South Korea will be
missing. Country coloring will be missing.)
load korea
worldmap(map, refvec)
land = ...
shaperead('landareas.shp',...
78
'UseGeoCoords', true);
faceColors =...
{'INDEX', [1 numel(land)],...
'FaceColor', ...
polcmap(numel(land))});
geoshow(land,'SymbolSpec',...
Version 2.1 (R14SP2) Mapping Toolbox™ Software
Mapping 1.x to
2.0.x UsageMapping 2.1 Usage
faceColors)
(Country coloring will be missing.)
DataFilesAddedinThisRelease
Thefollowingfileswereaddedtothemapdemos directory, for use in toolbox
demos and examples:
•
landareas — Polygon shapefile: global coastlines, both exterior and
interior, including names for larger land masses
•
worldlakes — Polygon shapefile: coastlines and names of major lakes
and inland seas worldwide
•
worldrivers — PolyLine shapefile: major world rivers and their names
worldcities — Point shapefile: locations and names of major cities and
•
populated places worldwide
•
usastatelo — Polygon shapefile: low-resolution outlines and names of
the 50 U .S. states plus D.C.
•
usastatehi — Polygon shapefile: moderate-resolution outlines and names
of the 50 U.S. states plus D .C.
•
conus — MAT -f il e: Low-resolution latitudes and longit ud es, in degrees,
for the perimeter of the conterminous United States (CONUS), the Great
Lakes, and interstate borders
•
greatlakes — MAT-file: A Version 1 geographic data structure
(geostruct1) with outlines and names for the Great Lakes of North America
Atlas Data MAT-Files Removed in T his Release
MAT-files containing Atlas Data have been removed in Version 2.1. Some of
the data has been retained in a different form. The disposition of these data
sets and variables is described below.
79
Mapping Toolbox™ Release Notes
World MAT-File Data
• worldlo.mat, which contained the following variables:
- DNline — Data moved to worldrivers.shp
- DNpatch —Datamovedtoworldlakes.shp
- POline — Data removed from toolbox
- POtext — Data removed from toolbox
- PPpoint — Data moved to worldcities.shp
- PPtext — Data moved to worldcities.shp
- gazette — Data removed from toolbox
worldhi.mat — D ata removed from toolbox
•
worldmtx.mat — D ata removed from toolbox
•
worldmtxmed.mat — D ata removed from too lbox
•
80
United States MAT-File Data
• usalo.mat, which contained the following variables (all retained):
- conus —Datamovedtoconus.mat
- greatlakes — Data moved to great lakes.mat
- gtlakelat —Datamovedtoconus. mat
- gtlakelon —Datamovedtoconus. mat
- state — Data moved to usastatelo.shp
- stateborder — Data moved to conus.mat
- statelat — Data moved to conus.mat
- statelon — Data moved to conus.mat
- uslat —Datamovedtoconus.mat
- uslon —Datamovedtoconus.mat
• usahi.mat — Data moved to usastat ehi.shp
Version 2.1 (R14SP2) Mapping Toolbox™ Software
Functions Delet
The following fu
sets, have been r
•
country2mtx —C
data
usahi —Return
•
•
usalo —Retur
•
worldhi —Ret
•
worldlo —Re
Functions M
• etopo5 —Re
•
tigerp —Do
format an
•
tigermif
format a
nctions, which performed specific operations on Atlas Data
emoved in Version 2.1:
high-resolution vector atlas data for the United States
n vector atlas data for the United States
urn high-re solution vector atlas data for the world
turn v ector atlas data for the world or oceans
ade Obsolete in This Release
placed by
wnload U.S. Census cartographic boundary files in shapefile
duse
shaperead instead.
— Download U.S. Census cartographic boundary files in shapefile
nd use
shaperead instead.
ed in This Release
reate a raster map grid of a country from
etopo.Useetopo instead.
worldlo atlas
81
Mapping Toolbox™ Release Notes
Version 2.0.3 (R14SP1) Mapping Toolbox Software
This table summarizes what’s new in Version 2.0.3 (R14SP1):
New Features
and Changes
Yes Details
below
The new features and changes introduced in this version are
Version
Compatibility
Considerations
Yes Details in
“Compatibility
Considerations”
on page 83
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
Fixed bugs
Related
Documentation
at Web Site
Current product
documentation
Polygon Vertex Ordering Functions
Four new functions have been added to Mapping 2.0.3. These are called by or
can be used in conjunction with the
Has New Helper Functions” on page 83.) Three of the functions check or
change the ordering of vertices that define a polygon, and the fourth one splits
polygons with holes in a consistent fashion.
These functions are
•
ispolycw — True if vertices of polygonal contour are clockwise ord ere d
polybool function (see “polybool Revised,
82
poly2cw — Convert polygonal contour to clockwise ordering
•
poly2ccw — Convert polygonal contour to counterclockwise ordering
•
poly2fv —Convertpolygonalregiontoface-vertex form for use with patch
•
in order to properly render polygons containing holes
Version 2.0.3 (R14SP1) Mapping Toolbox™ Software
Compatibility C
onsiderations
axesm Does Not Reorder Longitude Limits
Function axes m
[west east],b
The form
now requires that 'MapLonLimit' values have the form
ut east can be less than west numerically (e.g.,
[east
west]
is no longer supported.
contourm Returns C ontourgroup Handle
Function con
an array of ha
tourm
now returns a handle to a contourgroup object instead of
ndlestolineobjects.
elevation Has New Azimuth Output
Function e
correspon
This elimi
with func
levation
ding to the elevation angle (output 1) and slant range (output 2).
nates the need to use a geodesic or rhumb line azimuth, computed
tions
can now return a third output,theline-of-sightazimuth
azimuth or distance, as an approximate substitute.
gtopo30 Has New Options
The func
• Data can
tionality of
be read from a more general directory structure.
gtopo30 has been extended so that:
[120 -60]).
• Empty m
• Latitu
Instea
NaNs.
ay be supplied in place of any input argument.
de and longitude limits are no longer restricted to lie within a tile.
d, areas not covered by available data are automatically filled with
inputm Has New Button Output
tm
inpu
func
now has an optional third output argument, 'button',which
tions like the third output of
ginput.
polybool Revised, Has New Helper Functions
ping 2.0.3 fixes major problems with the
Map
putes set operations on polygonal regions. The fixes include incompatible
com
ta structure and syntax changes.
da
polybool function, which
83
Mapping Toolbox™ Release Notes
• The 'cutvector' option has been removed from the syntax. This option
facilitated plotting of polygonal regions containing holes, but the cutvector
computation was not robust, nor was it an effective way to represent
polygonal regions containing multiply nested holes and/or discontiguous
regions. To replace the cutvector functionality, a new function,
has been added.
holes and discontiguous regions, into a faces and a vertices matrix that
can be used with the patch function to display the region. The MATLAB
function help for
how to use
poly2fv,
poly2fv converts a polygonal region, possibly including
polybool and poly2fv contain examples illustrating
poly2fv.
• The
'cell' and 'vector' options have been removed from the syntax.
polybool now returns the output polygonal region using the same format
as the input polygonal region. Use
polysplit and polyjoin to convert
from one format to the other.
•
polybool now uses the “clockwise” rule to distinguish between external
polygonal contours and internal polygonal contours. This is the same rule
used by ESRI shapefiles.
polybool assumes that a polygonal contour
whose vertices are arranged in clockwise order is an external contour, and
that a polygonal contour whose vertices are arranged in counterclockwise
order is an internal contour. Use
polygonal contour is clockwise ordered, and use
ispolycw to determine whether a
poly2cw and poly2ccw to
convert polygonal contours to be clockwise or counterclockwise ordered,
respectively. If a polygonal region input has no external contours according
to this rule,
An input polygonal region can either take the form of a pair of
polybool issues a warning message.
NaN-separated
vectors, or it can take the form of a pair of cell arrays. In the cell array format,
each cell must contain the vertices of a single polygonal contour. Representing
a polygonal region as a cell array whose cells can contain
NaN-separated
vectorsisnolongersupported,andan error message will be issued.
Summary of new polybool-related helper functions. You can use the
helper functions that support
polybool foryourownpurposes. Theyare
84
•
ispolycw — True if vertices of polygonal contour are clockwise ordered
poly2cw — Convert polygonal contour to clockwise ordering
•
poly2ccw — Convert polygonal contour to counterclockwise ordering
•
Version 2.0.3 (R14SP1) Mapping Toolbox™ Software
• poly2fv —Convertpolygonalregiontoface-vertex form for use with patch
Summary of known problems addressed. The issues solved by the
new implementation of
•
polybool no longer errors out or produces incorrect results for proper
polybool are
inputs.
•
polybool now produces consistent output when computing intersections,
with results no longer depending on the order of its inputs.
•
polybool now operates correctly on polygons with collinear edges.
polymerge Output Vertex Order Changes
Improvements in the algorithm that polymerge uses to chain together
adjacent polygon or line segments, required to address a bug and improve
efficiency, in some cas es also chang e the ordering of the output vertices.
85
Mapping Toolbox™ Release Notes
Version 2.0.2 (R14) Mapping Toolbox Software
This table summarizes what’s new in Version 2.0.2 (R14):
New Features
and Changes
Yes
Details below
New features and changes introduced in this version are
Version
Compatibility
Considerations
No
Fixed Bugs
and Known
Problems
No bug fixesCurrent product
Related
Documentation
at Web Site
documentation
Mapping Toolbox Compilability Enhanced
Almost all Mapping Toolbox command line functions can be compiled using
the new MATLAB Compiler. However, large Atlas Data MAT-files need to be
added explicitly to use them in applications.
SDTS Functions Now Available on HP-UX
The input functions sdtsinfo and sdtsdemread are now available on the
HP-UX platform.
86
Version 2.0.1 (R13SP1+) Mapping Toolbox™ Software
Version 2.0.1 (R13SP1+) Mapping Toolbox Software
This table summarizes what’s new in Version 2.0.1 (R13SP1+):
New Features
and Changes
Yes
Details below
Thenewfeatureintroducedinthisversionis
Version
Compatibility
Considerations
No
Fixed Bugs
and Known
Problems
Fixed bugs,
below
Related
Documentation
at Web Site
V2.0.1 product
documentation
mapexreg Demo (Georeferencing an Image to an
Orthotile Base Layer) Generalized and Clarified
Calls to imtransform and makerefmat arenolongerhard-codedforanoutput
pixel size of one meter. The choice of output pixel size now is explained
in detail.
Fixed Bugs
Version 2.0.1 contains the following bug fixes:
almanac
Function almanac now gives a precise inverse flattening (1/f) value for Bessel
1841 ellipsoid. Previously,
of the full precision value (299.1528128).
almanac used a truncated value (299.1528) instead
axesm
Function axesm now accepts the 'meridianlabel' parameter. Previously,
axesm would throw an error if the 'meridianlabel' parameter was specified.
distance
Function distance now computes correct results for rhumb line distances
along a parallel (including the e quator) on a (nonspherical) ellipsoid.
87
Mapping Toolbox™ Release Notes
distance
Function distance now uses the robust “aversine formula” to compute great
circle distances on a sphere. For certain inputs roundoff sensitivities in the
previous implementation resulted in a small, but nonzero, distance between
identical points.
distance
Function distance now correctly obtains correct results for rhumb lines that
follow the e quator, or any other parallel, and cross the 180-degree meridian; it
gives the distance along the short arc of the parallel. Previously the length of
the long arc could be returned instead.
distance, azimuth
Functions distance and azimuth now produce accurate results even for long
geodesics on an ellipsoid, up to near-antipodal distances.
projfwd, projinv
Functions projfwd and projinv now use NaN separators on all platforms.
Previously, under certain circumstances,
coordinate arrays with
projfwd and projinv returned
Inf separators rather than NaNs.
88
reckon
Function reckon now supports geodesics on an ellipsoid, in addition to rhumb
linesonanellipsoidandbothgeodesicsandrhumblinesonasphere.
Version 2.0 (R13SP1+) Mapping Toolbox™ Software
Version 2.0 (R13SP1+) Mapping Toolbox Software
This table summarizes what’s new in Version 2.0 (R13SP1+):
New Features
and Changes
Yes
Details below
New features and changes introduced in this version are
• “Release Summary” on page 90
• “New Demos and Sample Data” on page 90
• “New Map Viewer” on page 92
• “Spatially Referenced Image Formats” on page 93
• “Working with Images Referenced to Map Coordinates” on page 94
• “SDTS Terrain Data Access” on page 95
• “Shapefiles and Vector Features” on page 96
• “Geographic Data Structure Enhancements” on page 96
Version
Compatibility
Considerations
Yes—Detailsin
“Compatibility
Considerations”
on page 101
below. S e e als o
Summary.
Fixed Bugs
and Known
Problems
“Fixed Bugs” on
page 103, below
Related
Documentation
at Web Site
V2.0.1 product
documentation
• “Attribute-Driven Vector D ata Symbolization” on page 97
• “Access to PROJ.4 Map Projections” on page 97
• “Minor Enhancements” on page 98
• “Summary of Functions Added in Version 2.0” on page 99
• “Compatibility Considerations” on page 101
• “Fixed Bugs” on page 103
89
Mapping Toolbox™ Release Notes
Release Summary
Mapping Toolbox
graphical user i
data analysis i
toolbox, new fu
Version 2.0 provides a compre hensiv e set of functions and
nterfaces for building map displays and performing geospatial
n the MATLAB environment. In this major release of the
nctions a nd capabilities include:
• Support for st
- ESRI shapefil
only)
- “Worldfiles
• New display
grids in x-y
• New functi
• Transvers
• Anewinte
- Support f
- Data imp
- Zooming
- Cursor c
- Data ti
- Graphi
- Map ex
• Amajo
conce
and e
rial to mak e the documentation easier to access
mate
map coordinates (as well as unprojected latitude-longitude)
ons for vector symbolization based on feature attributes
e Mercator projection a nd PROJ.4 projection library support
ractive M ap Viewer with the following features:
or multiple raster and vector layers
ort from file or workspace
/panning/map navigation tools
oordinateandmap-scalereadout
pandinfotools
cal overlays
porting/printing
r update of the Mapping Toolbox User’s Guide, focused on clarifying
pts and terminology, explaining features and functions, adding new
nhancing existing examples, and organizing and cross-referencing
• mapexrefmat — Creating and Using Referencing Matrices
• mapexreg — Georeferencing an Image to an Or th o tile Base Layer
• viewmaps — GUI Demonstrating Map Projections (not new in 2.0)
Note that the above commands run the demo scripts to produce figures,
whereas
mapdemos describes and illustrates the demos in the Help Browser.
In addition, a number of new sample data sets containing vector features,
digital elevation models, and georeferenced images have been added for
use in demos and examples. The new data includes satellite images o f the
Boston area, topographic grids of the White Mountains in New Hampshire,
and vecto r data for roads and hydrographic features in the Boston area. The
data sets are provided in standard geospatial and GIS formats: GeoTIFF,
TIFF images with worldfiles, SDTS Raster Profile, Arc ASCII Grid, and ESRI
shapefiles. They are accompanied by ASCII text files (with suffix
.txt),
containing descriptive metadata.
The new data files are in addition to existing sample data sets stored
in MAT-files, such as
geoid, russia,andkorea. In addition, metadata
indicating the source and describing the contents of most sample data and
atlas data has been added to their respective MAT-files.
To read the metadata for a Mapping Toolbox MAT-file data set,
inspect the
source and description workspace v aria bles.
load it and
To read brief descriptions of the demos and sample data sets, type
help mapdemos
91
Mapping Toolbox™ Release Notes
New Map Viewer
The toolbox includes a new interactive tool for displaying and examining map
data, called the Map Viewer. The Map Viewer helps you work with data
that is already in a projected map coordinate system, which is the case for
many high-resolution satellite and aerial imagery products, as well as many
vector map data sets that cover small areas of the Earth in substantial detail.
The Map Viewer gives you a view (or views) of an x-y map coordinate plane
where, for example, x and y may correspond to the easting and northing
coordinates of a given UTM or State Plane zone. If some of your d ata is
referenced to geographic (latitude-longitude) coordinates, you first need to use
an appropriate projection to transform it to map coordinates.
Key features of the Map Viewer include:
• A two-dimensional view of the x-y map plane with interactive tools for
navigating in that plane: magnification, de-magnification, panning, setting
the map scale, and zooming the view to the extent of a giv en data set, or
all data sets
• Vector and image data import from files or from the MATLAB workspace
92
• Management of each data set in a separate layer. You can control layer
ordering and visibility to bring what you need into view.
• Multiple, simultaneous views (differing in scale or view extent) of the same
data layers
• A data tip tool to identify vector shapes or query a feature attribute that
you specify
• An info tool to display, in a separate window, all the attributes of a selected
feature
• Tools to annotate the map with lines, arrows, and text
• Ability to expo rt a raster snapshot of the map view
To start the Map Viewer, use the
mapview
will open a Map Viewer window. You can then import standard format files
or load workspace data into the Map Viewer. The following figure shows the
mapview function. Typing
Version 2.0 (R13SP1+) Mapping Toolbox™ Software
Map Viewer’s main window after loading image, line, and point data sets, also
illustrating user-specified symbology for streets and landmarks, and Data
Tip labels for selected point features:
You can print maps as they appear in the Map Viewer w indow, to either a
printerorafile.EitherclickthePrint button on the left end of the toolbar or
select Print in the File menu. A standard figure print dialog appears. You
should consider selecting a driver appropriate to your printer; this may be
necessary if you attempt to print in color.
Spatially Referenced Image Formats
The toolbox adds support for several industry-standard formats for spatially
referenced images and data grids: Arc ASCII G rid files, GeoTIFF files, and
worldfiles. You do this by using the following new functions:
arcgridread,
93
Mapping Toolbox™ Release Notes
geotiffinfo, geotiffead, ge tworldfilename, worldfileread,and
worldfilewrite.
The grid and image file-reading functions construct referencing matrices,
enabling a general and consistent approach to managing the relationship
between map coordinates and pixel columns and rows (see below).
Working with Images Referenced to Map Coordinates
The toolbox includes a family of function s supporting images and gridded data
sets that are referenced to two-dimensional projected map coordinates. These
functions represent the relationship between a point at position (
an image and point in map coordinates (x,y) with a referencing matrix.A
referencing matrix is a 3-by-2 matrix,
transformation such that
[x y] = [row col 1] * R
Most often the image (or data grid) consistsofsquarepixelsandhasedgesthat
align with the map coordinate axes, but you are not limited to this situation.
row,col)in
R, that represents a general affine
94
The new function
makerefmat provides several options for constructing
referencing matrices from scratch. See its MATLAB function help for further
information. However, you can most often obtain a referencing matrix from a
world file, via
worldfileread, or from a GeoTIFF file via geoti ffre ad.If
you have created a new referencing matrix for an image in a generic image
format such as TIFF or JPEG, you can save this information in a world file
using
worldfilewrite.
New functions
defined above, along with its inverse. With functions
pix2map and map2pix implement the basic transformation
mapbbox and mapoutline
you can calculate the image bounds or perimeter in map coordinates. You
can u se the function
pixcenters to calculate the coordinates of the center of
each image pixel. This is helpful if you have a data grid that you want to
display using
The new function
surface.
mapshow accepts a referencing matrix in order to correctly
position an image in a standard Handle Graphics axes (i.e., not a map axes).
Version 2.0 (R13SP1+) Mapping Toolbox™ Software
Because the affine transformation is extremely general, it can equally well
represent the registration of an imagetogeographiccoordinates,with
longitude taking the place of map
[lon lat] = [row col 1] * R
x and latitude replacing map y:
A referencing matrix is a more general version of the 1-by-3 map legend
vector, perhaps better characterized as a referencing vector, that many toolbox
functions already use. To convert a referencing vector to a referencing
matrix, you can use function
refvec2mat. The inverse transformation exists
only for referencing matrices having a special form, but if it does you can
compute it with
refmat2vec. Functions latlon2pix and pix2latlon support
transformations for referencing m atrices in geographic coordinates, allowing
for longitude-wrapping differences.
Use the new function
geoshow along with map axes to display data grids
referenced to latitude and longitude via either a referencing vector or a
referencing m atrix. Or, you can use the new function
grid2image to quickly
display an unprojected grid with latitude as the ordinate (x-axis) and
longitude as the abscissa (y-coordinate).
SDTS Terrain Data Access
Mapping Toolbox functions can now import terrain grids stored in United
States Spatial Data Transfer Standard (SDTS) raster profile format, giving you
access to a wide selection of public domain terrain data sets, principally those
prepared by the U.S. Geo logical Survey (USGS). USGS DEM data is no longer
distributed from government Internet sites. Instead, these and other USGS
files are now available from a commercial server at no charge. For further
information, see
and http://edc.usgs.gov/geodata/.
Two new toolbox functions support SDTS data handling:
•
sdtsinfo returns a structure containing metadata for an SDTS data set.
sdtsdemread returns elevation data and a referencing matrix from an
•
SDTS DEM data set.
For details, type
prompt.
http://edc.usgs.gov/products/elevation/dem.html
help sdtsinfo or help sdtsdemread at the MATLAB
95
Mapping Toolbox™ Release Notes
Shapefiles and Vector Features
The toolbox now reads the shapefile format defined by Environmental
Systems Research Institute (ESRI) and widely used as a data exchange format
for nontopological GIS data. You can query a shapefile for its shapetype,
attribute names and types, number of records, etc., using
read 2-D shapes usin g
geostruct2 geographic data structure array (see “Geographic Data Structure
Enhancements” on page 96 ), tailored to the contents read in from shapefiles.
The s hape data is represented in the MATLAB workspace as a 1-D structure
array, with one e lem ent per feature. Attribute values are stored in structure
fields,orinaseparate,parallelstructure array. Because attributes are
defined according to the needs of the shapefile author and might be unknown
until the time a given shapefile is read, the new Mapping Toolbox geographic
data structures (des cribed below) are necessarily more general than the
existing toolbox geographic data structures.
Geographic Data Structure Enhancements
Certain Mapping Toolbox functions introduced in Version 2 use an enhanced
geographic data structure layout (called a geostruct2) to store and manipulate
vector geodata. This layout has the flexibility to store any kind and number of
attributes, and handles either geographic (latitude and longitude) or plane
(x and y ) coordinates. In contrast, a V ersion 1 geographic data structure
(geostruct1) — which is still supported — is limited to a fixed set of fields and
can contain geographic coordinates only.
shapeinfo,and
shaperead.Theshaperead function constructs a
96
One way to create a Version 2 geographic data structure is to input vector
geodata to the workspace from a shapefile. The function
shaperead returns
a geostruct2 that encapsulates s ome or all of the data stored in a shapefile
and its supporting index and dBASE files. To determine what kinds of data a
shapefiles contains, you can use the
The new functions
in geostruct2 form. Use
the geostruct, and use
mapshow, geoshow,andmapview each display vector data
geoshow when Lat and Lon coordinate fields exist in
mapshow or mapview when X and Y fields are present.
shapeinfo function to query it.
You can transform a geostruct1 into a geostruct2 (but not the reverse). Use
the function
updategeostruct for this purpose. See “Mapping Toolbox
Loading...
+ hidden pages
You need points to download manuals.
1 point = 1 manual.
You can buy points or you can get point for every manual you upload.