[GRASSweb-list]scott: web/faq importing_usgs_7_5min.html,1.2,1.3
grass at intevation.de
grass at intevation.de
Wed Dec 17 12:47:24 EST 2003
Author: scott
Update of /grassrepository/web/faq
In directory doto:/tmp/cvs-serv31191
Modified Files:
importing_usgs_7_5min.html
Log Message:
commenting page according to request tracker #1871
Index: importing_usgs_7_5min.html
===================================================================
RCS file: /grassrepository/web/faq/importing_usgs_7_5min.html,v
retrieving revision 1.2
retrieving revision 1.3
diff -u -d -r1.2 -r1.3
--- importing_usgs_7_5min.html 30 Oct 2002 11:14:16 -0000 1.2
+++ importing_usgs_7_5min.html 17 Dec 2003 17:47:21 -0000 1.3
@@ -1,302 +1,151 @@
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
-<title>Importing USGS 7.5min SDTS DEM into GRASS GIS</title>
-<body bgcolor=#FFFFE0>
- <strong>
-Importing USGS 7.5min SDTS DEM into GRASS GIS
-</strong><p>
-
-This shows how to import USGS 7.5min SDTS DEM data into grass. It includes
-where to find those pesky SDTS programs and FIPS libraries for working
-with the USGS data sets. Another method of importing the USGS 7.5min SDTS DEM files is by using the <a href="../gdp/html_grass5/html/r.in.gdal.html">r.in.gdal</a> command followed by the <a href="../gdp/html_grass5/html/r.in.ascii.html">r.in.ascii</a> command.
-<p>
-
-<b>1) Acquire the fips123 library for the sdts utilities</b>
-<p>
-The fips123 library (~619k) is available from:
-<p>
-ftp://sdts.er.usgs.gov/pub/sdts/software/fips123/fips123.ver104.tar.gz
-<p>
-Extract and compile the fips123 library:
-<p>
-<pre>
- tar xzvf fips123.ver104.tar.gz
- cd fips123
- make
-</pre>
-Once the make is complete, you should have a fips123.a library in the
-fips123 directory. You should move the library and include files to a
-more suitable location:
-<p>
-<pre>
- cp f123inc/*.h /usr/local/include/
- cp fips123.a /usr/local/lib/libfips123.a
-</pre>
-<p>
-
-<b>2) Acquire the sdts utilities</b>
-<p>
-A multitude of utilities for converting the SDTS format data into
-formats more readily importable into grass are available (~38k) from:
-<p>
-ftp://ftp.blm.gov/pub/gis/sdts/dem/sdtscode.zip
-<p>
-Extract the utilities:
-<p>
-<pre>
- unzip -L sdtscode.zip
-</pre>
-<p>
-Applying the patch listed in Appendix A to the sdtscode source will
-help with case sensitivity issues on filenames. The datasets from
-USGS use all uppercase filenames. These patches force the use of
-all uppercase filenames.
-<p>
-<pre>
- patch < sdtscode.patch
-</pre>
-<p>
-Compile the utilities. Appendix B shows the Makefile I used.
-<pre>
- make
-</pre>
-Once the make is complete, you should have several utility programs in
-this directory. You should copy the utilites to a directory in your path:
-<p>
-<pre>
- cp demout sdts2dem sdts2sur sdtsgrid sdts2arc sdts2rcz sdts2xyz /usr/local/bin/
-</pre>
-<p>
-
-
-<b>3) Acquire data from the USGS website</b>
-<p>
-Browse available data at:
-<p>
-http://edcwww.cr.usgs.gov/doc/edchome/ndcdb/ndcdb.html
-<p>
-In this case, our quest for Washington, DC data leads us to the
-"Washington West, DC" quad:<br>
-7.5 Minute Digital Elevation Model Data for N38.94 W77.08
-and the actual data set:
-<p>
-http://edcftp.cr.usgs.gov/pub/data/DEM/7.5min/W/washington_west_DC/30.2.1.1191443.tar.gz
-<p>
-Download the data to a temporary directory, and extract the data files:
-<p>
-<pre>
-tar zxvf 30.2.1.1191443.tar.gz
-</pre>
-<p>
-NOTE: the files are extracted into the current directory.
-<p>
-
-<b>4) Preprocess the USGS data</b>
-<p>
-Convert the data to an ascii format (nearly) suitable for importing:
-<p>
-<pre>
-sdts2arc 9304 washdc L0 | tee washdc.txt
-</pre>
-<p>
-9304 - the numeric prefix on all files in the data set<br>
-washdc - the base of our output filename (the filename will be washdc.grd)<br>
-L0 - (ell zero) is the "cell id", usually L0, follows the CE in nnnnCExx.DDF<br>
-| tee washdc.txt - save a copy of the program output to washdc.txt
-<p>
-The output of sdts2arc will show the latitude and longitude of the quad:
-<p>
-<pre>
-LAT/LONG of the 7.5 USGS Quad
- SW 315655.156250 4304843.500000
- NW 315979.093750 4318716.500000
- NE 326804.687500 4318471.500000
- SE 326500.250000 4304598.500000
-</pre>
-<p>
-These are really the UTM coordinates for the corners of the quad. This
-will be needed when setting up location information for the data. Full
-output from the sdts2arc run for Washington West, DC in Appendix C.
-<p>
-Filter the void values from the arc file:
-<p>
-<pre>
-sed -e 's/ 32770/ 0/g' washdc.grd > washdc2.grd
-</pre>
-<p>
-[NOTE: grass5.0 may handle void values differently]
-<p>
-The argument after the -e is:
-<pre>
-quote-slash-space-three-two-seven-seven-zero-slash-space-zero-quote
-</pre>
-<p>
-Edit the header information from the file to include the grass ascii
-raster file header information. Remove the original header data. Use
-your favourite plain text editor. The header should contain the
-followinng data (explained below):
-<p>
-<pre>
-north: 4318716.500
-south: 4304598.500
-west: 315655.1562
-east: 326804.6875
-rows: 471
-cols: 372
-</pre>
-<p>
-The values are obtained from the LAT/LONG of USGS Quad of sdts2arc
-output. NOTE: the coordinates are printed in long/lat order.
-<p>
-north is chosen as the larger lat of NW/NE corners<br>
-south is chosen as the smaller lat of SW/SE corners<br>
-west is chosen as the smaller long of SW/NW corners<br>
-east is chosen as the larger long of SE/NE corners<br>
-rows is the nrow parameter from the sdts2arc output<br>
-cols is the ncol parameter from the sdts2arc output
-<p>
-Save this file as washdc.ascii
-<p>
-
-<b>5) Set up a new location in grass for the data</b>
-<p>
-Start up grass and give it a new location name. When it asks for the
-parameters for your new location, use the following:
-<p>
-Coordinate System: 1 (UTM)
-<p>
-UTM Zone for Location: 18 [get this value from sdts2arc output]
-<p>
-One Line Description: WASHINGTON WEST, DC-MD-VA [anything appropriate]
-<p>
-Default Region: [pick values from LAT/LONG of USGS Quad of sdts2arc
-output, note that the coordinates are printed in long/lat order]
-<pre>
- North Edge: 4318716.500 [larger lat of NW/NE corners]
- South Edge: 4304598.500 [smaller lat of SW/SE corners]
- West Edge: 315655.1562 [smaller long of SW/NW corners]
- East Edge: 326804.6875 [larger long of SE/NE corners]
-
-Grid Resolution:
- East-West: 30 [the xhrs value from the sdts2arc output]
- North-South: 30 [the xhrs value from the sdts2arc output]
-
-Ellispoid Name: wgs84 [not sure where this came from, but seems to work]
-</pre>
-<p>
-Tweak the DEFAULT_WIND created by grass:
-<p>
-change the e-w resol and n-s resol to 30 by editing the file
-<p>
-Propagate those changes into the WIND file:
-<p>
-<pre>
- cp DEFAULT_WIND WIND
-</pre>
-<p>
-
-<b>6) Import the elevation data as an ascii file and take a look</b>
-<p>
-Importing as an ascii raster file is done as follows:
-<p>
-<pre>
-r.in.ascii input=washdc.ascii output=elev title=Elevation mult=1.0
-</pre>
-<p>
-To view the data:
-<p>
-<pre>
-d.mon start=x0 select=x0
-d.rast map=elev -o
-</pre>
-You should see a lovely rendering of Washington, DC, prominently
-featuring the Potomac River (heading to the NW), and Rock Creek
-(heading N)
-<p>
-<hr>
-<b>Appendix A:</b> sdtscode.patch - case sensitivity patches for sdtscode.zip utilities
-<p>
-<pre>
-diff -ubBp ./temp/sdts2sur.c sdtsutils/sdts2sur.c
---- ./temp/sdts2sur.c Sat Aug 27 17:07:56 1994
-+++ sdtsutils/sdts2sur.c Wed Sep 22 18:19:12 1999
-@@ -266,7 +266,7 @@ void dem_rc(int status)
- {
+<head>
+ <title>Importing USGS 7.5min SDTS DEM into GRASS GIS</title>
+</head>
+ <body bgcolor="#ffffe0">
+ <strong> Importing USGS 7.5min SDTS DEM into GRASS GIS<br>
+<br>
+The instructions below are probably out of date. They are retained
+here for now pending confirmation of this assertion. New users are
+directed towards Moritz Lennert's <a
+ href="http://grass.itc.it/gdp/grass5tutor/HTML_en/book1.html">Quickstart
+to GRASS 5</a> guide, instead.<br>
+<br>
+---------------------------<br>
+ </strong>
+<p> This shows how to import USGS 7.5min SDTS DEM data into grass. It includes
+where to find those pesky SDTS programs and FIPS libraries for working with
+the USGS data sets. Another method of importing the USGS 7.5min SDTS DEM
+files is by using the <a href="../gdp/html_grass5/html/r.in.gdal.html">r.in.gdal</a>
+command followed by the <a
+ href="../gdp/html_grass5/html/r.in.ascii.html">r.in.ascii</a> command. </p>
+<p> <b>1) Acquire the fips123 library for the sdts utilities</b> </p>
+<p> The fips123 library (~619k) is available from: </p>
+<p> ftp://sdts.er.usgs.gov/pub/sdts/software/fips123/fips123.ver104.tar.gz
+</p>
+<p> Extract and compile the fips123 library: </p>
+<p> </p>
+<pre> tar xzvf fips123.ver104.tar.gz<br> cd fips123<br> make<br></pre>
+ Once the make is complete, you should have a fips123.a library in the fips123
+directory. You should move the library and include files to a more suitable
+location:
+<p> </p>
+<pre> cp f123inc/*.h /usr/local/include/<br> cp fips123.a /usr/local/lib/libfips123.a<br></pre>
- strcpy (file_name,base_name);
-- strcat (file_name,"ldef.ddf");
-+ strcat (file_name,"LDEF.DDF");
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
- {
- printf ("\nCAN'T OPEN 'dem' LPEF FILE %s",file_name);
-@@ -349,9 +349,9 @@ void degenlines(int status)
- {
- /* degenerate lines sored as points */
- strcpy (file_name,base_name);
-- strcat (file_name,"ne");
-+ strcat (file_name,"NE");
- strcat (file_name,cellid);
-- strcat (file_name,".ddf");
-+ strcat (file_name,".DDF");
- p=0;
- /* strcat (file_name,module);*/
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
-@@ -843,9 +843,9 @@ void demnodes(int status)
- }
- */
- strcpy (file_name,base_name);
-- strcat (file_name,"no");
-+ strcat (file_name,"NO");
- strcat (file_name,cellid);
-- strcat (file_name,".ddf");
-+ strcat (file_name,".DDF");
- p = 0;
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
- {
-@@ -1029,9 +1029,9 @@ void demareas(int status)
- {
+<p> <b>2) Acquire the sdts utilities</b> </p>
+<p> A multitude of utilities for converting the SDTS format data into formats
+more readily importable into grass are available (~38k) from: </p>
+<p> ftp://ftp.blm.gov/pub/gis/sdts/dem/sdtscode.zip </p>
+<p> Extract the utilities: </p>
+<p> </p>
+<pre> unzip -L sdtscode.zip<br></pre>
- strcpy (file_name,base_name);
-- strcat (file_name,"na");
-+ strcat (file_name,"NA");
- strcat (file_name,cellid);
-- strcat (file_name,".ddf");
-+ strcat (file_name,".DDF");
- /* strcat (file_name,module);*/
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
- {
-@@ -1211,7 +1211,7 @@ void demmbr(int status)
- double fl;
- /* int label ;*/
- strcpy (file_name,base_name);
-- strcat (file_name,"spdm.ddf");
-+ strcat (file_name,"SPDM.DDF");
- /* strcat (file_name,module);*/
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
- {
-@@ -1565,9 +1565,9 @@ void loadpoints(int status)
- {
- /* these are the degenerate lines */
- strcpy (file_name,base_name);
-- strcat (file_name,"ne");
-+ strcat (file_name,"NE");
- strcat (file_name,cellid);
-- strcat (file_name,".ddf");
-+ strcat (file_name,".DDF");
- nxy = 0;
- /* strcat (file_name,module);*/
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
-@@ -1766,7 +1766,7 @@ void cells_out(int status)
- strcpy (file_name,base_name);
- strcat (file_name,"CE");
- strcat (file_name,cellid);
-- strcat (file_name,".ddf");
-+ strcat (file_name,".DDF");
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
- {
- printf ("\nCAN'T OPEN CELL VALUES FILE %s...",file_name);
-<p>
-<p>
+<p> Applying the patch listed in Appendix A to the sdtscode source will help
+with case sensitivity issues on filenames. The datasets from USGS use all
+uppercase filenames. These patches force the use of all uppercase filenames.
+</p>
+<p> </p>
+<pre> patch < sdtscode.patch<br></pre>
+
+<p> Compile the utilities. Appendix B shows the Makefile I used. </p>
+<pre> make<br></pre>
+ Once the make is complete, you should have several utility programs in this
+directory. You should copy the utilites to a directory in your path:
+<p> </p>
+<pre> cp demout sdts2dem sdts2sur sdtsgrid sdts2arc sdts2rcz sdts2xyz /usr/local/bin/<br></pre>
+
+<p> <b>3) Acquire data from the USGS website</b> </p>
+<p> Browse available data at: </p>
+<p> http://edcwww.cr.usgs.gov/doc/edchome/ndcdb/ndcdb.html </p>
+<p> In this case, our quest for Washington, DC data leads us to the "Washington
+West, DC" quad:<br>
+ 7.5 Minute Digital Elevation Model Data for N38.94 W77.08 and the actual
+data set: </p>
+<p> http://edcftp.cr.usgs.gov/pub/data/DEM/7.5min/W/washington_west_DC/30.2.1.1191443.tar.gz
+</p>
+<p> Download the data to a temporary directory, and extract the data files:
+</p>
+<p> </p>
+<pre>tar zxvf 30.2.1.1191443.tar.gz<br></pre>
+
+<p> NOTE: the files are extracted into the current directory. </p>
+<p> <b>4) Preprocess the USGS data</b> </p>
+<p> Convert the data to an ascii format (nearly) suitable for importing:
+ </p>
+<p> </p>
+<pre>sdts2arc 9304 washdc L0 | tee washdc.txt<br></pre>
+
+<p> 9304 - the numeric prefix on all files in the data set<br>
+ washdc - the base of our output filename (the filename will be washdc.grd)<br>
+ L0 - (ell zero) is the "cell id", usually L0, follows the CE in nnnnCExx.DDF<br>
+ | tee washdc.txt - save a copy of the program output to washdc.txt </p>
+<p> The output of sdts2arc will show the latitude and longitude of the quad:
+</p>
+<p> </p>
+<pre>LAT/LONG of the 7.5 USGS Quad<br> SW 315655.156250 4304843.500000<br> NW 315979.093750 4318716.500000<br> NE 326804.687500 4318471.500000<br> SE 326500.250000 4304598.500000<br></pre>
+
+<p> These are really the UTM coordinates for the corners of the quad. This
+will be needed when setting up location information for the data. Full output
+from the sdts2arc run for Washington West, DC in Appendix C. </p>
+<p> Filter the void values from the arc file: </p>
+<p> </p>
+<pre>sed -e 's/ 32770/ 0/g' washdc.grd > washdc2.grd<br></pre>
+
+<p> [NOTE: grass5.0 may handle void values differently] </p>
+<p> The argument after the -e is: </p>
+<pre>quote-slash-space-three-two-seven-seven-zero-slash-space-zero-quote<br></pre>
+
+<p> Edit the header information from the file to include the grass ascii raster
+file header information. Remove the original header data. Use your favourite
+plain text editor. The header should contain the followinng data (explained
+below): </p>
+<p> </p>
+<pre>north: 4318716.500<br>south: 4304598.500<br>west: 315655.1562<br>east: 326804.6875<br>rows: 471<br>cols: 372<br></pre>
+
+<p> The values are obtained from the LAT/LONG of USGS Quad of sdts2arc output.
+ NOTE: the coordinates are printed in long/lat order. </p>
+<p> north is chosen as the larger lat of NW/NE corners<br>
+ south is chosen as the smaller lat of SW/SE corners<br>
+ west is chosen as the smaller long of SW/NW corners<br>
+ east is chosen as the larger long of SE/NE corners<br>
+ rows is the nrow parameter from the sdts2arc output<br>
+ cols is the ncol parameter from the sdts2arc output </p>
+<p> Save this file as washdc.ascii </p>
+<p> <b>5) Set up a new location in grass for the data</b> </p>
+<p> Start up grass and give it a new location name. When it asks for the
+parameters for your new location, use the following: </p>
+<p> Coordinate System: 1 (UTM) </p>
+<p> UTM Zone for Location: 18 [get this value from sdts2arc output] </p>
+<p> One Line Description: WASHINGTON WEST, DC-MD-VA [anything appropriate]
+</p>
+<p> Default Region: [pick values from LAT/LONG of USGS Quad of sdts2arc output,
+note that the coordinates are printed in long/lat order] </p>
+<pre> North Edge: 4318716.500 [larger lat of NW/NE corners]<br> South Edge: 4304598.500 [smaller lat of SW/SE corners]<br> West Edge: 315655.1562 [smaller long of SW/NW corners]<br> East Edge: 326804.6875 [larger long of SE/NE corners]<br><br>Grid Resolution:<br> East-West: 30 [the xhrs value from the sdts2arc output]<br> North-South: 30 [the xhrs value from the sdts2arc output]<br><br>Ellispoid Name: wgs84 [not sure where this came from, but seems to work]<br></pre>
+
+<p> Tweak the DEFAULT_WIND created by grass: </p>
+<p> change the e-w resol and n-s resol to 30 by editing the file </p>
+<p> Propagate those changes into the WIND file: </p>
+<p> </p>
+<pre> cp DEFAULT_WIND WIND<br></pre>
+
+<p> <b>6) Import the elevation data as an ascii file and take a look</b>
+</p>
+<p> Importing as an ascii raster file is done as follows: </p>
+<p> </p>
+<pre>r.in.ascii input=washdc.ascii output=elev title=Elevation mult=1.0<br></pre>
+
+<p> To view the data: </p>
+<p> </p>
+<pre>d.mon start=x0 select=x0<br>d.rast map=elev -o<br></pre>
+ You should see a lovely rendering of Washington, DC, prominently featuring
+the Potomac River (heading to the NW), and Rock Creek (heading N)
+<p> </p>
+<hr> <b>Appendix A:</b> sdtscode.patch - case sensitivity patches for sdtscode.zip
+utilities
+<p> </p>
+<pre>diff -ubBp ./temp/sdts2sur.c sdtsutils/sdts2sur.c<br>--- ./temp/sdts2sur.c Sat Aug 27 17:07:56 1994<br>+++ sdtsutils/sdts2sur.c Wed Sep 22 18:19:12 1999<br>@@ -266,7 +266,7 @@ void dem_rc(int status)<br> {<br> <br> strcpy (file_name,base_name);<br>- strcat (file_name,"ldef.ddf");<br>+ strcat (file_name,"LDEF.DDF");<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br> {<br> printf ("\nCAN'T OPEN 'dem' LPEF FILE %s",file_name);<br>@@ -349,9 +349,9 @@ void degenlines(int status)<br> {<br> /* degenerate lines sored as points */<br> strcpy (file_name,base_name);<br>- strcat (file_name,"ne");<br>+ strcat (file_name,"NE");<br> strcat (file_name,cellid);<br>- strcat (file_name,".ddf");<br>+ strcat (file_name,".DDF");<br> p=0;<br> /* strcat (file_name,module);*/<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br>@@ -843,9 +843,9 @@ void demnodes(int status)<br> }<br> */<br>
strcpy (file_name,base_name);<br>- strcat (file_name,"no");<br>+ strcat (file_name,"NO");<br> strcat (file_name,cellid);<br>- strcat (file_name,".ddf");<br>+ strcat (file_name,".DDF");<br> p = 0;<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br> {<br>@@ -1029,9 +1029,9 @@ void demareas(int status)<br> {<br> <br> strcpy (file_name,base_name);<br>- strcat (file_name,"na");<br>+ strcat (file_name,"NA");<br> strcat (file_name,cellid);<br>- strcat (file_name,".ddf");<br>+ strcat (file_name,".DDF");<br> /* strcat (file_name,module);*/<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br> {<br>@@ -1211,7 +1211,7 @@ void demmbr(int status)<br> double fl;<br> /* int label ;*/<br> strcpy (file_name,base_name);<br>- strcat (file_name,"spdm.ddf");<br>+ strcat (file_name,"SPDM.DDF");<br> /* strcat (file_name,module);*/<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,
&fpin))<br> {<br>@@ -1565,9 +1565,9 @@ void loadpoints(int status)<br> {<br> /* these are the degenerate lines */<br> strcpy (file_name,base_name);<br>- strcat (file_name,"ne");<br>+ strcat (file_name,"NE");<br> strcat (file_name,cellid);<br>- strcat (file_name,".ddf");<br>+ strcat (file_name,".DDF");<br> nxy = 0;<br> /* strcat (file_name,module);*/<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br>@@ -1766,7 +1766,7 @@ void cells_out(int status)<br> strcpy (file_name,base_name);<br> strcat (file_name,"CE");<br> strcat (file_name,cellid);<br>- strcat (file_name,".ddf");<br>+ strcat (file_name,".DDF");<br> if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))<br> {<br> printf ("\nCAN'T OPEN CELL VALUES FILE %s...",file_name);<br><p>
+</p><p>
diff -ubBp ./temp/sdts2xyz.c sdtsutils/sdts2xyz.c
--- ./temp/sdts2xyz.c Thu Apr 30 14:59:56 1998
+++ sdtsutils/sdts2xyz.c Wed Sep 22 10:34:18 1999
@@ -306,7 +155,7 @@
strcpy (file_name,base_name);
- strcat (file_name,"ldef.ddf");
+ strcat (file_name,"LDEF.DDF");
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
+ if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
{
printf ("\nCAN'T OPEN 'dem' LPEF FILE %s",file_name);
@@ -640,7 +640,7 @@ void dem_mbr(int status)
@@ -316,7 +165,7 @@
- strcat (file_name,"spdm.ddf");
+ strcat (file_name,"SPDM.DDF");
/* strcat (file_name,module);*/
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
+ if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
{
@@ -965,7 +965,7 @@ void cells_out(int status)
strcpy (file_name,base_name);
@@ -324,59 +173,27 @@
strcat (file_name,cellid);
- strcat (file_name,".ddf");
+ strcat (file_name,".DDF");
- if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
+ if (! beg123file (file_name,'R',&int_level,&ice,ccs,&fpin))
{
printf ("\nCAN'T OPEN CELL VALUES FILE %s...",file_name);
+</p><p>
+</p><hr>
+<b>Appendix B:</b> Makefile for sdtscode.zip utilities<br><p>
+</p><pre>all: demout sdts2dem sdts2sur sdtsgrid sdts2arc sdts2rcz sdts2xyz<br><br>%: %.c<br> gcc -o $@ $@.c -Wl,/usr/local/lib/libfips123.a<br></pre>
<p>
-<hr>
-<b>Appendix B:</b> Makefile for sdtscode.zip utilities
-<p>
-<pre>
-all: demout sdts2dem sdts2sur sdtsgrid sdts2arc sdts2rcz sdts2xyz
-
-%: %.c
- gcc -o $@ $@.c -Wl,/usr/local/lib/libfips123.a
-</pre>
+</p><hr>
+<b>Appendix C:</b> sdts2arc for Washington West, DC DEM elevation data<br><p>
+</p><pre> Spatial Data Transfer Standard (SDTS)<br> SDTS to ARCINFO ASCII GRID Utility<br> SDTS2ARC BETA ver .003<br> Another gis translation tool from<br> Sol Katz(skatz at blm.gov), Mar 1998.<br><br>Title: WASHINGTON WEST, DC-MD-VA - 24000 <br>Data ID: LAT:: 38.875 LONG:: -77 SCALE:: 24000<br>Data set creation date: 19980708<br><br>cell width: 30.000000<br>cell height: 30.000000<br>rsnm: UTM<br>zone: 18<br>nrow= 471, ncol= 372 <br>Range: max= 423, min= 0, void= -32767, fill= -32766<br>LAT/LONG of the 7.5 USGS Quad <br> SW 315655.156250 4304843.500000<br> NW 315979.093750 4318716.500000<br> NE 326804.687500 4318471.500000<br> SE 326500.250000 4304598.500000<br>processing cells <br><br> END OF FILE <br><br>End of Program<br></pre>
<p>
+</p><hr>
+<a href="index.html">Go back to GRASS GIS FAQ...</a>
<hr>
-<b>Appendix C:</b> sdts2arc for Washington West, DC DEM elevation data
-<p>
-<pre>
- Spatial Data Transfer Standard (SDTS)
- SDTS to ARCINFO ASCII GRID Utility
- SDTS2ARC BETA ver .003
- Another gis translation tool from
- Sol Katz(skatz at blm.gov), Mar 1998.
-
-Title: WASHINGTON WEST, DC-MD-VA - 24000
-Data ID: LAT:: 38.875 LONG:: -77 SCALE:: 24000
-Data set creation date: 19980708
-
-cell width: 30.000000
-cell height: 30.000000
-rsnm: UTM
-zone: 18
-nrow= 471, ncol= 372
-Range: max= 423, min= 0, void= -32767, fill= -32766
-LAT/LONG of the 7.5 USGS Quad
- SW 315655.156250 4304843.500000
- NW 315979.093750 4318716.500000
- NE 326804.687500 4318471.500000
- SE 326500.250000 4304598.500000
-processing cells
-
- END OF FILE
-
-End of Program
-</pre>
-<p>
-<HR>
-<A HREF="index.html">Go back to GRASS GIS FAQ...</A>
-<HR>
<em>Last Modified: 30th Dec. 1999</em>
<address>
-<A HREF="mailto:emitchell at altaira.com">Eric Mitchell</A>
+<a href="mailto:emitchell at altaira.com">Eric Mitchell</a>
</address>
+</pre>
+<br>
</body>
</html>
More information about the grass-web
mailing list