GTOP030 Tutorial - Revision history

Dparks at 15:25, 4 June 2009

2009-06-04T15:25:19Z

Beiko at 15:20, 4 June 2009

2009-06-04T15:20:41Z

Dparks at 19:23, 23 May 2009

2009-05-23T19:23:44Z

Dparks at 19:18, 23 May 2009

2009-05-23T19:18:18Z

Dparks at 19:15, 23 May 2009

2009-05-23T19:15:46Z

Dparks at 19:15, 23 May 2009

2009-05-23T19:15:31Z

Dparks at 15:24, 22 May 2009

2009-05-22T15:24:50Z

Dparks at 15:19, 22 May 2009

2009-05-22T15:19:41Z

Dparks at 15:18, 22 May 2009

2009-05-22T15:18:13Z

Dparks at 02:55, 12 May 2009

2009-05-12T02:55:54Z

@@ Line 33: / Line 33: @@
   gdal_translate.exe -a_ullr 260 -10 300 -60 W100S10.DEM X260S10.DEM
-Will change the origin (upper left longitude) from -100 (as implied by the input file name) to +260. Note that I am using the prefix letter 'X' to indicate a remapped Western hemisphere file.
+Will change the origin (upper left longitude) from -100 (as implied by the input file name) to +260 (i.e., 360 + (-100)). Note that I am using the prefix letter 'X' to indicate a remapped Western hemisphere file.
 You also need to manually edit the .HDR file that is attached to the .DEM. The original header will contain consecutive lines like this:

← Older revision		Revision as of 15:20, 4 June 2009
Line 17:		Line 17:

	'''Note:''' GDAL does not overwrite files! If you wish to run either of these commands again with different parameters you either need to change the output file name or delete the previously made file.		'''Note:''' GDAL does not overwrite files! If you wish to run either of these commands again with different parameters you either need to change the output file name or delete the previously made file.
		+
		+	== Spanning the 180 degree line ==
		+
		+	A bit of extra work is required to build maps that span 180 degrees of longitude in the Pacific Ocean. In the usual coordinate system, this is a point of discontinuity, with -180 in the Western hemisphere and +180 in the Eastern.
		+
		+	The easiest solution to this is to keep the prime meridian as the zero coordinate, and monotonically increase the coordinate value as you go east, until you reach 360 degrees again at the prime meridian. Then there is a smooth transition in the Pacific, and North and South America will have longitudes in the 200-300 degree range.
		+
		+	GDAL can work in this coordinate system, and GenGIS has no problem with it either. But there are two things you will need to do before loading in your data:
		+
		+	# Use GDAL to convert any WESTERN HEMISPERE maps to the (180,360) coordinate system (note that EASTERN HEMISPHERE coordinates are unchanged)
		+	# Convert any Western hemisphere locations you may have from (-180,0) to (180,360)
		+
		+	Two things need to be done to address (1), if you are starting from a typical .DEM file. First, use 'gdal_translate' to redefine the origin of your map. For instance:
		+
		+	gdal_translate.exe -a_ullr 260 -10 300 -60 W100S10.DEM X260S10.DEM
		+
		+	Will change the origin (upper left longitude) from -100 (as implied by the input file name) to +260. Note that I am using the prefix letter 'X' to indicate a remapped Western hemisphere file.
		+
		+	You also need to manually edit the .HDR file that is attached to the .DEM. The original header will contain consecutive lines like this:
		+
		+	ULXMAP -99.99583333333334
		+	ULYMAP -10.00416666666667
		+
		+	You need to change ULXMAP to reflect the updated origin. Replace -99.99583333333334 with 260.00416666666667 to reflect the new point of origin. You do not need to change ULYMAP, but it's handy for copying over the numbers after the decimal point. Rename this file to make it consistent with the modified .DEM, in this case X260S10.HDR.
		+
		+	Secondly, in your location file any coordinates with a negative longitude x will need to be replaced with 360 + x. So a point at location (-120,30) will become (240,30).
		+
		+	That's it! Once the maps are individually modified, they can be manipulated as usual with GDAL, and loaded into GenGIS.

	[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]		[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]

@@ Line 3: / Line 3: @@
 * Download and install [http://fwtools.maptools.org/ FWTools]. This will place a ''FWTools Shell'' icon on your desktop. From within this command line environment, a number of useful GDAL functions can be accessed.
 * Goto [http://edc.usgs.gov/products/elevation/gtopo30/gtopo30.html GTOPO30 site] and download the four DEM maps which span that area we are interested in: w100n90.tar.gz, w100n40.tar.gz, w060n40.tar.gz, w060n90.tar.gz
-* Start the FWTools Shell by double clicking the ''FWTools Shell'' icon on your desktop.
+* Decompress the archives into a folder.
 * To extract the area spanning 100W50N to 50W5N from these four DEMs use <tt>gdal_merge</tt>:

@@ Line 1: / Line 1: @@
 These instructions indicate how to construct a custom digital elevation map (DEM) using the GTOPO30 data. In this tutorial, we will create a map of the North Atlantic Seaboard which runs from 100W50N to 50W5N (see image below).
-* Download and install [http://www.gdal.org/ GDAL]. This will place a ''FWTools Shell'' icon on your desktop. From within this command line environment, a number of useful GDAL functions can be accessed.
+* Download and install [http://fwtools.maptools.org/ FWTools]. This will place a ''FWTools Shell'' icon on your desktop. From within this command line environment, a number of useful GDAL functions can be accessed.
 * Goto [http://edc.usgs.gov/products/elevation/gtopo30/gtopo30.html GTOPO30 site] and download the four DEM maps which span that area we are interested in: w100n90.tar.gz, w100n40.tar.gz, w060n40.tar.gz, w060n90.tar.gz
 * Start the FWTools Shell by double clicking the ''FWTools Shell'' icon on your desktop.

@@ Line 1: / Line 1: @@
 These instructions indicate how to construct a custom digital elevation map (DEM) using the GTOPO30 data. In this tutorial, we will create a map of the North Atlantic Seaboard which runs from 100W50N to 50W5N (see image below).
-* Download and install [http://www.gdal.org/|GDAL]. This will place a ''FWTools Shell'' icon on your desktop. From within this command line environment, a number of useful GDAL functions can be accessed.
+* Download and install [http://www.gdal.org/ GDAL]. This will place a ''FWTools Shell'' icon on your desktop. From within this command line environment, a number of useful GDAL functions can be accessed.
 * Goto [http://edc.usgs.gov/products/elevation/gtopo30/gtopo30.html GTOPO30 site] and download the four DEM maps which span that area we are interested in: w100n90.tar.gz, w100n40.tar.gz, w060n40.tar.gz, w060n90.tar.gz
 * Start the FWTools Shell by double clicking the ''FWTools Shell'' icon on your desktop.

← Older revision		Revision as of 15:24, 22 May 2009
Line 17:		Line 17:

	[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]		[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]
−
−	~~<html>~~
−	~~<head>~~
−	~~<!-- saved from url=(0025)http://www.techsmith.com/ -->~~
−	~~<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />~~
−	~~<meta name="DC.title" content="H1N1_SwineFlu_Spread" />~~
−	~~<meta name="DC.date" content="2009-05-21 04:19:41 PM" />~~
−	~~<meta name="DC.language" content="ENU" />~~
−	~~<meta name="DC.subject" content="Influenza" />~~
−	~~<meta name="itunes.keywords" content="Influenza, H1N1, Swine Flu" />~~
−	<meta name="DC.description" content="Geographic visualizations built using GenGIS (http://kiwi.cs.dal.ca/GenGIS) Automated data retrieval and tree construction: Norm MacDonald GenGIS visualizations and data handling: Donovan Parks Mike Porter Suwen Wang The development of GenGIS wa" />
−
−	~~<title>Created by Camtasia Studio 6</title>~~
−	~~<script type="text/javascript" src="swfobject.js"></script>~~
−	~~<script type="text/javascript">~~
−	~~swfobject.registerObject("csSWF", "8", "expressInstall.swf");~~
−	~~</script>~~
−	~~<style type="text/css">~~
−	~~body~~
−	{
−	~~background-color: #1a1a1a;~~
−	~~font: .8em/1.3em verdana,arial,helvetica,sans-serif;~~
−	~~text-align: center;~~
−	}
−	~~#media~~
−	{
−	~~margin-top: 40px;~~
−	}
−	~~#noUpdate~~
−	{
−	~~margin: 0 auto;~~
−	~~font-family:Arial, Helvetica, sans-serif;~~
−	~~font-size: x-small;~~
−	~~color: #cccccc;~~
−	~~text-align: left;~~
−	~~width: 210px;~~
−	~~height: 200px;~~
−	~~padding: 40px;~~
−	}
−	~~</style>~~
−	~~</head>~~
−	~~<body>~~
−	~~<div id="media">~~
−	~~<object classid="clsid:D27CDB6E-AE6D-11cf-96B8-444553540000" width="768" height="514" id="csSWF">~~
−	~~<param name="movie" value="H1N1_SwineFlu_Spread.swf" />~~
−	~~<param name="quality" value="best" />~~
−	~~<param name="bgcolor" value="#1a1a1a" />~~
−	~~<param name="allowfullscreen" value="false" />~~
−	~~<param name="scale" value="noscale" />~~
−	~~<param name="allowscriptaccess" value="always" />~~
−	~~<param name="flashvars" value="autostart=false&color=0x1A1A1A,0x1A1A1A" />~~
−	~~<!--[if !IE]>-->~~
−	~~<object type="application/x-shockwave-flash" data="H1N1_SwineFlu_Spread.swf" width="768" height="514">~~
−	~~<param name="quality" value="best" />~~
−	~~<param name="bgcolor" value="#1a1a1a" />~~
−	~~<param name="allowfullscreen" value="false" />~~
−	~~<param name="scale" value="noscale" />~~
−	~~<param name="allowscriptaccess" value="always" />~~
−	~~<param name="flashvars" value="autostart=false&color=0x1A1A1A,0x1A1A1A" />~~
−	~~<!--<![endif]-->~~
−	~~<div id="noUpdate">~~
−	<p>The Camtasia Studio video content presented here requires a more recent version of the Adobe Flash Player. If you are you using a browser with JavaScript disabled please enable it now. Otherwise, please update your version of the free Flash Player by <a href="http://www.adobe.com/go/getflashplayer">downloading here</a>.</p>
−	~~</div>~~
−	~~<!--[if !IE]>-->~~
−	~~</object>~~
−	~~<!--<![endif]-->~~
−	~~</object>~~
−	~~</div>~~
−	~~<!-- Users looking for simple object / embed tags can copy and paste the needed tags below.~~
−	~~<div id="media">~~
−	~~<object id="csSWF" classid="clsid:d27cdb6e-ae6d-11cf-96b8-444553540000" width="768" height="514" codebase="http://download.macromedia.com/pub/shockwave/cabs/flash/swflash.cab#version=9,0,115,0">~~
−	~~<param name="src" value="H1N1_SwineFlu_Spread.swf"/>~~
−	~~<param name="bgcolor" value="#1a1a1a"/>~~
−	~~<param name="quality" value="best"/>~~
−	~~<param name="allowScriptAccess" value="always"/>~~
−	~~<param name="allowFullScreen" value="false"/>~~
−	~~<param name="scale" value="noscale"/>~~
−	~~<param name="flashVars" value="autostart=false"/>~~
−	<embed name="csSWF" src="H1N1_SwineFlu_Spread.swf" width="768" height="514" bgcolor="#1a1a1a" quality="best" allowScriptAccess="always" allowFullScreen="false" scale="noscale" flashVars="autostart=false&color=0x1A1A1A,0x1A1A1A" pluginspage="http://www.macromedia.com/shockwave/download/index.cgi?P1_Prod_Version=ShockwaveFlash"></embed>
−	~~</object>~~
−	~~</div>~~
−	~~-->~~
−	~~</body>~~
−	~~</html>~~

@@ Line 18: / Line 18: @@
 [[Image:NorthAtlantic.jpg|frame|Map of North Atlantic Seaboard from 100W50N to 50W5N.|center]]
-<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
+<html>
      <head>
          <!-- saved from url=(0025)http://www.techsmith.com/ -->

← Older revision		Revision as of 15:18, 22 May 2009
Line 17:		Line 17:

	[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]		[[Image:NorthAtlantic.jpg\|frame\|Map of North Atlantic Seaboard from 100W50N to 50W5N.\|center]]
		+
		+	<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
		+	<!-- saved from url=(0014)about:internet -->
		+	<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en">
		+	<head>
		+	<!-- saved from url=(0025)http://www.techsmith.com/ -->
		+	<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />
		+	<meta name="DC.title" content="H1N1_SwineFlu_Spread" />
		+	<meta name="DC.date" content="2009-05-21 04:19:41 PM" />
		+	<meta name="DC.language" content="ENU" />
		+	<meta name="DC.subject" content="Influenza" />
		+	<meta name="itunes.keywords" content="Influenza, H1N1, Swine Flu" />
		+	<meta name="DC.description" content="Geographic visualizations built using GenGIS (http://kiwi.cs.dal.ca/GenGIS) Automated data retrieval and tree construction: Norm MacDonald GenGIS visualizations and data handling: Donovan Parks Mike Porter Suwen Wang The development of GenGIS wa" />
		+
		+	<title>Created by Camtasia Studio 6</title>
		+	<script type="text/javascript" src="swfobject.js"></script>
		+	<script type="text/javascript">
		+	swfobject.registerObject("csSWF", "8", "expressInstall.swf");
		+	</script>
		+	<style type="text/css">
		+	body
		+	{
		+	background-color: #1a1a1a;
		+	font: .8em/1.3em verdana,arial,helvetica,sans-serif;
		+	text-align: center;
		+	}
		+	#media
		+	{
		+	margin-top: 40px;
		+	}
		+	#noUpdate
		+	{
		+	margin: 0 auto;
		+	font-family:Arial, Helvetica, sans-serif;
		+	font-size: x-small;
		+	color: #cccccc;
		+	text-align: left;
		+	width: 210px;
		+	height: 200px;
		+	padding: 40px;
		+	}
		+	</style>
		+	</head>
		+	<body>
		+	<div id="media">
		+	<object classid="clsid:D27CDB6E-AE6D-11cf-96B8-444553540000" width="768" height="514" id="csSWF">
		+	<param name="movie" value="H1N1_SwineFlu_Spread.swf" />
		+	<param name="quality" value="best" />
		+	<param name="bgcolor" value="#1a1a1a" />
		+	<param name="allowfullscreen" value="false" />
		+	<param name="scale" value="noscale" />
		+	<param name="allowscriptaccess" value="always" />
		+	<param name="flashvars" value="autostart=false&color=0x1A1A1A,0x1A1A1A" />
		+	<!--[if !IE]>-->
		+	<object type="application/x-shockwave-flash" data="H1N1_SwineFlu_Spread.swf" width="768" height="514">
		+	<param name="quality" value="best" />
		+	<param name="bgcolor" value="#1a1a1a" />
		+	<param name="allowfullscreen" value="false" />
		+	<param name="scale" value="noscale" />
		+	<param name="allowscriptaccess" value="always" />
		+	<param name="flashvars" value="autostart=false&color=0x1A1A1A,0x1A1A1A" />
		+	<!--<![endif]-->
		+	<div id="noUpdate">
		+	<p>The Camtasia Studio video content presented here requires a more recent version of the Adobe Flash Player. If you are you using a browser with JavaScript disabled please enable it now. Otherwise, please update your version of the free Flash Player by <a href="http://www.adobe.com/go/getflashplayer">downloading here</a>.</p>
		+	</div>
		+	<!--[if !IE]>-->
		+	</object>
		+	<!--<![endif]-->
		+	</object>
		+	</div>
		+	<!-- Users looking for simple object / embed tags can copy and paste the needed tags below.
		+	<div id="media">
		+	<object id="csSWF" classid="clsid:d27cdb6e-ae6d-11cf-96b8-444553540000" width="768" height="514" codebase="http://download.macromedia.com/pub/shockwave/cabs/flash/swflash.cab#version=9,0,115,0">
		+	<param name="src" value="H1N1_SwineFlu_Spread.swf"/>
		+	<param name="bgcolor" value="#1a1a1a"/>
		+	<param name="quality" value="best"/>
		+	<param name="allowScriptAccess" value="always"/>
		+	<param name="allowFullScreen" value="false"/>
		+	<param name="scale" value="noscale"/>
		+	<param name="flashVars" value="autostart=false"/>
		+	<embed name="csSWF" src="H1N1_SwineFlu_Spread.swf" width="768" height="514" bgcolor="#1a1a1a" quality="best" allowScriptAccess="always" allowFullScreen="false" scale="noscale" flashVars="autostart=false&color=0x1A1A1A,0x1A1A1A" pluginspage="http://www.macromedia.com/shockwave/download/index.cgi?P1_Prod_Version=ShockwaveFlash"></embed>
		+	</object>
		+	</div>
		+	-->
		+	</body>
		+	</html>

@@ Line 6: / Line 6: @@
 * To extract the area spanning 100W50N to 50W5N from these four DEMs use <tt>gdal_merge</tt>:
-gdal_merge -o north_atlantic.tif -of GTiff -ul_lr -100 50 -50 5 W060N40.DEM W060N90.DEM W100N40.DEM W100N90.DEM
+ gdal_merge -o north_atlantic.tif -of GTiff -ul_lr -100 50 -50 5 W060N40.DEM W060N90.DEM W100N40.DEM W100N90.DEM
 * This will result in a large (~62MB) GeoTiff file. Generally, it is best to work with smaller files within GenGIS. We can subsample this GeoTiff file using <tt>gdal_translate</tt>: