Select a category from the list below or use 'Edit', 'Find' on your
browser to search by key words. The Questions are in boldface and
the Answers are in lighter face.
3DEM
I have found your site to be an incredible resource
and introduction to
geospatial data acquisition and use.
I'm trying to follow along with your 3DEM tutorial but I'm having difficulty
getting usable data to load in the demo version (3dem70). I have acquired
through your suggested methods DEM, SDTS and DRG data for a practice site
(Greeley, CO which is in Weld County). I successfully (I think) ran the
SDTS
data through the SDTS2DEM converter application with the resulting file
named greeley.dem. Then I go to 3dem70 and click on "load terrain
-> digital
model" and select the greeley.dem file. In the window below a graphic
appears that doesn't seem to make much sense. It looks like this:
<<greeleydemscreenshot.jpg>>
Can you help me understand what is going on? Is this somehow corrupt data?
Furthermore, it seems that 3dem70 does not have the "operation ->
apply
overlay" function. Is this correct?
Your DEM file is apparently corrupt. I do not know how it got that
way. I
received an email about two weeks ago from a user who obtained good results
from the converter using his computer at work but somehow created a
corrupted DEM file using his computer at home. He sent me his DEM
file and
I looked at it in hexedit and could clearly see that all the records were
in
the right places but the elevation values were nonsensical (47,000 meters,
etc.)
SDTSDEM5 creates large temporary array data structures from heap memory
so
theoretically memory is limited only by the size of your hard disk.
However, 10m DEMs like that for Greely are very large. I am wondering
if
maybe there is a memory page size setting problem or some similar issue.
I
have an idea on how to fix this and may send you a revised program tomorrow
to test if you do not mind.
Anyway, I downloaded the Greely DEM and ran it through SDTSDEM5 on my
computer and then loaded it up in 3DEM70. Everything worked fine,
except
that it took forever to create the 3D view because of the size of the
file
(9MB). (30m dems are much more render-friendly.)
I have attached my output file called mygreely.dem. Try loading
this into
3DEM70 and make sure that it works. You could then try comparing
it to
yours using a text editor or alternatively send me your file and I could
diagnose the problem. (I would appreciate if you could do this so
I can
possibly figure out what is going on.)
You are technically correct about the overlay selection. In 3DEM70,
you can
do an overlay but you get at it a little differently. First generate
a 3D
view using 'View Scene'. Once the scene has rendered, choose 'Operation',
'Modify Textures'. A window will come up. The leftmost section
will be
labeled 'Terrain'. Under this section, choose 'Load'. You
will be prompted
for the name of the .bmp file to overlay. You will notice that the
registration features are not as good in 3DEM70 as in 3DEM.
That is the good news. The bad news is that I have never been able
to get
the overlay utility to work properly in 3DEM70. It always makes
funky
overlays that are sort of semi-opaque. The feature works perfectly
in 3DEM,
however. One used to be able to download a free 3DEM demo to try
but I
think the vendor has discontinued this offer. However, the new program
does
not cost much and works flawlessly.
Good luck and hope to hear from you soon.
Here is the modified converter to try. Just in
case there is some type of memory related problem I changed the sequencing
for opening and deleting the internal data structures. Instead of
opening them all at the beginning and closing them all at the end, I staggered
the opening and closing so that only two maximum are open at any one time.
This should reduce the memory burden on your computer, especially for
large files like the Greely cel0 file. See if this makes any difference.
Top
3DStudioMax
Hey there John, Nice site. I'm picking
up what I can from it, sifting through info to try and find what I'm looking
for. I'm hoping that you might expedite my quest. I need to
find/create some high resolution grayscale images of the grand canyon
to be used as displacement maps inside of 3dstudiomax. I downloaded
your DEM2TGA6
and
SDTSTGA2. To be honest, I'm not positive which I need.
I've used MicroDEM to spit out a bitmap of the DEM, but was left feeling
like there is a better way to get a higher resolution image. I could
use some advise on URLs to go to and what files I need. (some
of those files are so obscurely named that to the novice they might as
well say anything but something I'd recognize. That pretty much
sums up what I'm doing. Any help at all would be greatly appreciated.
I am not familiar with 3dsmax at all, so all I can do
is guess. I checked out the 3dsmax web page and there was a reference
to world xyz projections so maybe there is hope. What you want to
do is import xyz information into your application so that you can use
it for creating virtual terrains or whatever. Digital terrain data
consists of raster data but your application must map the data in order
to make an image. The raster itself is not the image. There
are dozens of formats that this data comes in, but most of it gridded
elevation data, i.e. an x coordinate, a y coordinate and an elevation
associated with the location on the xy plane. SDTS, USGS DEM, DTED0 and
GTOPO are a few of the formats that the raw data is presented to you in.
The trick then becomes to determine if your application contains a provision
to convert data in this form to something that you can use, for example
- a surface.
The only general-purpose rendering tool with which I
am familiar is POV-Ray. This application wants its gridded elevation
data presented to it in a special Truevision 24 byte TGA file format,
as explained on my web page. Each pixel represents a grid point
and the RGB color value is proportional to the elevation at that
point. Thus the utilities that I provide for converting the USGS
DEM (source) data available from the USGS website to this format.
These utilities will not do you any good unless 3ds max likewise wants
its data in that particular format.
So the first thing to do is to find out what your application
wants (if it can accept this type of data at all). Then you need
to find a source for elevation data ( several are listed on my website).
Then you need to find or create a utility to convert the data to the format
that your application will accept. If you take a look at POV-Ray
and see how they do it you will get a pretty good idea of what to
look for and what questions to ask for 3dsmax. An email to the 3ds
max website for instructions might be a good idea as well.
Top
3DViz
I am running 3D Viz and am attempting to support residential
and commercial land planning. AS far as local survey data no
problem but I also do presentation and "fly overs" and could
benefit from adding larger areas to each project. Is there any way
to convert the SDTS DEM files to any other AutoCAD file? We have
a good ray tracer in Viz but of course it doesn't recognize the .tga presented
by your file converter. Any help or advice would be appreciated.
I am embarrassed to admit that before your email I had
never heard of 3D Viz. I looked it up and could not find any indication
of a way to get DEM data, which is effectively raster information, into
3DViz. The 3DViz documentation indicated that Viz objects are solids,
(ACIS solids reference) which is not good news for importing raster data
although there was also a reference to 3D faces, which sounds like a surface
and would be good news. Gridded raster data like DEMs can be converted
to a faceted surface by "connecting the dots" of the grid and
calculating surface normals. DEM data could be converted to .tin,
.raw or .stl format this way if it would do any good.
I will keep checking and let you know if I find anything
later, Mark wrote...
John, thank you for the quick response. Perhaps
you work late as I. 3d Viz, being an AutoDesk product supports most
drawing files related to CAD. Viz will open or import the following
file types:
AutoCAD .DWG, DXF, Microstation .DGN,
... IGES, STL, .3ds, .max
Graphic, surface materials, and maps may be in the standard
formats including .tga, .jpg, .psd, .tiff, .rgb, .rla, .rpf, .mov, .png,
.bmp, .gif
Viz will create terrain similar to what you show by
loading the survey data as lines, set elevation, select the baseline,
create terrain, then select operands. As you select lines a mesh
is created and several options exist including color by level, multi-color
materials, and any map you can load thru the graphics programs.
I use Photoshop. Of course for a small area this works fine...larger
areas are difficult to make realistic. I have taken black and white
satellite photos, tiled them together, created bump maps, and applied
them to the terrain and I agree with your article that it is a bit tedious
to position and scale to fit. At this point, I don't know what
other questions to ask. If I have contours I can make it work.
There has to be similarity to all of this but right now we're talking
two different languages.
Top
ArcView
I'm interested in the work you've done in Wyoming
using digital terrain. I
work for bp out of Houston, Texas and my area of interest is Wyoming.
I
would like to utilize some of the resources you have listed but find them
mostly in a format I can't use. The mapping package I am currently using
can
accept shape files and this is why I'm contacting you. You list a number
of
converting programs on your web page but none of them convert to a shape
file. Are you aware of any way I can take the digital maps from the USGS
and
convert them to a shape file? Maybe there is a place where I can download
these files without converting them. Any help you can give would be
appreciated. I am mainly working around the Wind river Mountains and the
Uinta Mountains.
ArcView has a bundled extension that converts from SDTS
(which is what you
need) to shapfile, if you have that package. You just load the extension
in
order to import that type.
If not, here is the link to a
converter:http://www.engineering.usu.edu/cee/faculty/dtarb/SDTSext.html
Top
ASTER DEM
Thanks for your replies. I checked out your
updated
tutorial and it's as good as the others on your site.
Great step by step directions. Tried your software on
my downloaded ASTERs and it worked like a charm. It
was able to handle the conversion of the full files
without a problem on my machine (PIII500, 128MB),
producing wonderful results.
By the time I got your emails, I had already spent a
lot of time trying to find a free or inexpensive HDF
package that would allow me to open the original HDF
files and then transform them to georeferenced
DEMs...to no avail. I found several applications, but
they either wouldn't work on my machine or were way
too expensive. If I am not mistaken, the advantage of
having such a program is that it would open up access
to all the Aster HDF files available (there is limited
coverage for the ASTER DEM product), which I believe
can be converted to DEMs. I was hoping you could
provide your insight on the matter.
Sorry for the delay in replying to your email.
I got a flurry (actually a
ton) of emails following the article and somehow yours got missed.
Both HDF and Geotiff, carry embedded metadata. Both formats are
poorly documented but I was able to decipher geotiff first and already
had a
tiff reader and so I went with that format. I inspected the geotiff
proprietary tags and geokeys and know that there is very minimal embedded
metadata in the geotiff file. The HDF version may have more metadata
embedded.
If you want to get at the embedded data cleanly, I recommend buying
GlobalMapper from www.globalmapper.com.
Mike Childs there (no relation)
recently added ASTER Geotiff capability, and the program displays all
of the
embedded metadata. The program also writes the metadata to the USGS
ASCII
file automatically so that you not have to key it in.
HDF format becomes important for all of the other ASTER and other remotely
sensed satellite data, which is only available in HDF. This is my
next big
focus area after I finish with the geotiff DEMs and I will be publishing
my
findings then, hopefully around Easter time.
> 1) GEOTIF prompts for the UTM coordinates in decimal meters of
the south
> west corner. The metadata file lists four latitude and four
longitude
> values. Can you tell me definitively which ones correspond
to south west?
Here is object data from a .met file that I recently pulled:
OBJECT = GRingPointLatitude
CLASS = "1"
Value = (33.9744, 33.9561, 33.3185, 33.3364)
TYPE = "DOUBLE"
NUM_VAL = 4
END_OBJECT = GRingPointLatitude
OBJECT = GRingPointLongitude
CLASS = "1"
Value = (35.4294, 36.2276, 36.2039, 35.4116)
TYPE = "DOUBLE"
NUM_VAL = 4
END_OBJECT = GRingPointLongitude
I believe that the first
vertex of the gring (polygon) is the upper left (northwest) which would
be 33.9744N 35.4294E and they proceed in a clockwise direction
so the southwest coordinate would be the last pair, in this case 33.3364N
35.4116E. As far as I can tell, these coordinates do not coincide
with the
first valid data point but instead coincide with the first data
point of
any type, in this case a missing data point. That is, the coordinates
are
for the extent of the image, including the black border. I think.
The origin for a USGS ASCII DEM file is the lower
left corner with the rows increasing to the right and the columns increasing
up. This may give you a chance at getting oriented as far
as subsetting goes.
> 2) When we convert to UTM we use the NADCON program and functions.
Do we
> specify the input projection as Geographic coordinates with the WGS
1984
> ellipsoid and then enter the latitude and longitude of the southwest
corner
> in decimal degrees as reported in the Metadata file with the output
> projection being UTM, with the WGS 1984 ellipsoid and the appropriate
UTM
> zone? I just want to make sure we are georeferencing this correctly?
This is where I will really lose marks. Your question made
me realize that there is an error in GEOTIF .
ASTER DEMs are UTM WGS 84.
Unfortunately, I left
the horizontal and vertical datum flags set for NGVD29 and NAD27,
respectively when I wrote the program. The horizontal datum flag
should definitely be set to 3, (not 1like I had it in the program.)
I am not sure what
the vertical datum is for ASTER so I (now) set it to 1 for local
mean sea level. I
have corrected this and reposted the program. The two flags are
in bytes
889 and 891 of the DEM file. You can inspect the file and change
the state of these flags for existing files. The applicable information
from the USGS specification for ASCII DEM files is as follows:
Vertical datum INTEGER*1 I2 889 890
1=local mean sea level
2=National Geodetic Vertical
Datum 1929 (NGVD 29)
3=North American Vertical
Datum 1988 (NAVD 88)
Horizontal datum INTEGER*1 I2 891 892
1=North American Datum 1927 (NAD 27)
2=World Geodetic System 1972 (WGS 72)
3=WGS 84
4=NAD 83
5=Old Hawaii Datum
6=Puerto Rico Datum
7=NAD 83 Provisional (shifts in
horizontal coordinates are
computed, but old DEM nodes
are not resampled)
Now that this is corrected, inserting the proper coordinates and UTM zone
will allow the application that is reading the DEM file to display the
data accurately. I hope. To answer your original question,
since ASTER data is UTM WGS84 (unknown (to me) vertical datum) and now
the datum flags in the ASCII output file are set correctly (except for
possibly the vertical datum) the answer to your question is yes, I believe
that what you are doing is correct.
A couple of suggestions. For serious ASTER work you may want to
check out GlobalMapper. Although the
program is not free, it will read ASTER geotiff data and write USGS DEM
files just like GEOTIF, except probably better. It costs about $150
and is available at
www.globalmapper.com. Mike Childs (no relation) will probably
be happy to answer questions about that program. The second suggestion
is to contact the customer support people at EOS-EDC for information.
Brandy Adams has been very helpful and I am sure that she could answer
a lot of your questions better than I. Her email address is
bradams@usgs.gov and you can also talk to her by calling the telephone
number listed on the EOS Data Gateway page.
> 3) When we specify a subregion of the ASTER DEM to write out in UGS
DEM
> format will it be correctly oriented as long as we did 1 and 2 above
> correctly?
Well, mostly. However, when I wrote the ASCII data I assumed a rotational
angle of zero degrees, i.e. that the ASTER rows and columns were oriented
north and south. If what I stated above about the correspondence
of the Gring data to the corners of the image is true, that suggests that
there may be a slight rotation. Or maybe not if the polygon is not
a rectangle. I have not been able to find an indication of rotation
in the metadata (I have not looked too hard either) but there may or may
not be an error for this reason.
> 4) Can you think of problems using ASTER DEM's for hydrologic purposes?
They would seem to be ideal for this purpose.
>5) At what percent cloud cover should we be concerned?
You say not to use
>above
> 10, but is there a lot more error (for hydrologic practices) with
data at 7
> or 8 percent compared to data at 3 or 4?
Cloud cover is a dicey thing. Unfortunately, the ASTER bands are
not ideal
for either visible imaging or DEM production. My colleague Paul
Burkhardt
told me that he has had good luck with cloud cover values as high as 35%.
It all depends on what type of clouds and where they are located.
Heavy
dense clouds covering part of the image while leaving the rest clear will
yield good data for the clear zones while a thin uniform overcast may
cause
the entire DEM to be suspect. The best I can suggest is to look
at the
corresponding L1A image (or the browse image, but this is pretty small).
In
general arid areas will yield uniformly good DEMs while mountainous areas
in
the winter can be a problem.
> 6) Is there anything wrong with people from other countries logging
in and
> downloading these data?
No. If my website
traffic analysis is any indication, 25% to 50% of ASTER
usage comes from overseas. All of the emails that I have received
from
foreign users, express extreme gratitude to the
United States government for providing data that is unavailable in their
own
countries. My data indicates that Italy, France, UK, Spain and Canada
are
the largest foreign patrons.
I worked for Trojan Explosives in Spanish Fork from 1989 to 1993 in
environmental and process and ordnance engineering. I have lived
in various parts of the country working for the parent company, Ensign-Bickford.
Utah was my favorite place of all the states where I have lived.
To try to get some extra credit to make up for my
mistakes, I have posted a new article on ASTER L1A/L1B false color
composite images on the site. I must believe that this data would
be useful in your work as major and minor streams, rivers and lakes can
be seen easily. Since ASTER DEMs are made from L1A images, the ground
features should help in georeferencing as well.
Be sure to download the program MultiSpec if
you do not already have it. It is free and you definitely need some
type of data reader to wade through ASTER data. The program is very
easy to use and quite powerful. The link is in the article.
You don't by any chance know of a program that creates a DEM (relative
is
fine) from L1A data, do you? It sounds as though that must be ordered
through the EOS gateway. Several of the scenes that work for me
are
available in L1A format but there doesn't appear to be any DEMs available
or my area.
There are at least two capable programs, PCI Geomatics
at
http://www.pcigeomatics.com/
(this is the program the folks at EOS use to
produce the DEMs you order.) Virtuzo from Supresoft at is another
at
www.supresoft.com. Neither
program is particularly cheap.
I have not tried either but I have read good things about both.
Question: I was able to find about six raw ASTER
files, that will need to be
processed into DEMs. I'm a bit unclear about how to get the raw
data files
ground-truthed. Can you give me some hints?
The general method for georeferencing either ASTER DEMs or ASTER L1A images
requires the following steps:
1) Find the exact image locations of at least three
and preferably more ground control points distributed evenly across your
image. These points should ideally coincide with high and low points
of the DEM.
2) Derive an affine transformation matrix by running
a least-squares regression on the ground control points.
3) Apply the transformation matrix to each pixel
on the image to map the old location to the new location. In the
case of DEMs, each point has three dimensions while L1A data has two dimensions
and involves only planar transformations.
This is not the easiest thing to do. EOS-EDC
uses the same PCI Geomatics package used to produce the DEMs from the
stereo pairs to do the corrections to produce absolute DEMs. You
would need a similar tool to do the same, unless you request an absolute
DEM from EOS-EDC in which case they will do it for you.
I got the following information from the EOS-EDC
FAQ page regarding the latter procedure:
What is a Ground Control Point (GCP)?
Answer:
A precisely known UTM or Geographic coordinate of a point on the earth
used to geolocate the image. GCPs are needed to create absolute geolocated
DEMs. The requestor must supply these points, along with their precise
locations on the 3N and 3B images.
Examples:
Projection: Geographic Datum: WGS84
Feature Location in VNIR Band 3N: Line: 2434 Sample: 153
Feature Location in VNIR Band 3B: Line: 2104 Sample: 302
Coordinates: 44:32:36.42N, 116:26:35.27W
Elevation: 569 meters
Accuracy: X: 10 meters, Y: 10 meters, Z: 10 meters
Projection: UTM Datum: WGS84
Feature Location in VNIR Band 3N: Line: 2434 Sample: 153
Feature Location in VNIR Band 3B: Line: 2104 Sample: 302
UTM Zone: 14
Coordinates: Northing: 342567.34 Easting: 2349554.36
Elevation: 569 meters
Accuracy: X: 10 meters, Y: 10 meters, Z: 10 meters
Question:
How do I collect Ground Control Points for my image?
Answer
GCPs can come from several sources, including fieldwork where the location
the ground control point is determined using a Global Positioning System
(GPS), or from a map of known accuracy. These sources will provide the
exact
geographic coordinates of a certain location on the earth that is
identifiable on each of the stereo ASTER images (3N and 3B). Locations
that
make good GCPs include easily identifiable road intersections or other
ground level stationary objects. We require at least eight quality GCPs
and
they must be evenly distributed horizontally and vertically throughout
the
image to produce a quality DEM.
Question
What is the difference between a Relative DEM and an Absolute DEM?
Answer
An Absolute ASTER DEM is generated by using Ground Control Points (GCP)
supplied by the requestor. The software uses the GCPs to tie the DEM to
known points on the ground and yields a product with real ground elevations.
The Absolute DEM is geocoded using customer supplied Ground Control Points
(GCPs).
A Relative ASTER DEM is generated by using only the satellite ephemeris
data
and will yield a product with ground elevations that are quite near those
of
an Absolute DEM. The Relative DEM is not nearly as accurate as the Absolute
DEM horizontally since the Relative DEM is geocoded using only satellite
ephemeris data.
I read your Aster DEM! article on Spatialnews with great interest.
I am working on an archaeological project and Costa Rica and this seemed
perfect. I checked EOS and the Aster DEM product was not available
for this area, so I would like to be able to order request it. How
can I go about doing this? To order an on-demand DEM you have
to query the dataset: ASTER LEVEL 1A DATA SET - RECONSTRUCTED, UNPROCESSED
INSTRUMENT DATA V002 At the EOS Data Gateway. When your
L1A granules are presented there is a link on the search results page that
you select if you want a DEM produced from the ASTER L1A data.
There is an online help menu to aid this process. The official answer
from the EOS FAQ page is:
First, you must choose an appropriate ASTER L1A image. This is
done through the EOS Data Gateway. After the granule ID of the scene has
been obtained, you can order the production of the DEM on the ASTER On-Demand
Processing Requests page. To order a relative DEM, the granule ID is all
that is required. Absolute DEMs require the input of the granule ID along
with customer-derived Ground Control Points (GCPs). Note: The customer must
order the ASTER L1A image from the EOS Data Gateway to determine the Line
and Pixel values that must accompany the GCPs. You can probably
ignore the part about the absolute DEM unless you require the best possible
georeferencing. If you do not get any joy from this, I would suggest
contacting Brandy Adams at EROS South Dakota 605-594-6116 to see what other
options you might have, like perhaps ordering a fly-by if possible.
I have a question with regard to ASTER DEMs. First off, the
process that
is discussed on your website (using GEOTIF) to create a USGS quality DEM,
is this for a relative DEM or are there options in the program to add
GCPs
for the creation of an absolute DEM?
Second, through the EOS gateway, they will create a DEM for you, but as
you mentioned it takes a long time. The process they describe for
creating an absolute DEM requires at least 8 GCPs in addition to some
information from the ASTER 1A image (line of the sample location, etc.,
in
BOTH band 3N and 3B). When I downloaded the image and tried viewing
it in
ENVI Freelook, I was not able to find where you can visualize band 3.
Do
you have any experience with this or how I can "easily" visualize
band 3?
Thank you in advance for your help. Your website has been a great
help!
I have not done an absolute DEM myself but here is my opinion.
GEOTIF would not help you produce an absolute DEM as it only takes the
DEM
already produced and converts it to USGS ASCII format for subsequent
rendering. What you need to do is what you described in paragraph
2.
MultiSpec might be a useful tool for this job. I have been able
to view bands 3N and 3B for both L1A and L1B data easily. In order
to do this job right however, it would seem that you would
need to drop marker pixels in the image to check your referencing.
I do not
think that you can do this with MultiSpec, or any of the data readers
that I
have seen. I will keep on the lookout. Anyway the MultiSpec
link is given
on my web page (newest article) along with instructions on how to do exactly
what you are interested in.
The following is a copy of my letter to Brandy Adams of EOS-EDC regarding
the ASTER vertical elevation problem: Brandy:
I maintain a not-for-profit GIS webpage at www.terrainmap.com
where I have posted a few articles on ASTER. During the past week I have
received three emails from ASTER users reporting vertical errors on their
relative ASTER DEMs of up to 300m as compared to known standards like topo
map data. All the users are aware of the projection and datum used for ASTER
images so I do not think this is the problem. They have asked me for the
reason and I am at a loss for an explanation. Based
on the information on the FAQ page, I would have thought that the errors
due to uncertainties in satellite position would have been less than what
these users are observing. Do you think that this is within the normal range
for relative DEMs? I have written a program called
GEOTIF that translates ASTER geotiff files into USGS ASCII files. I would
like to add a translation/rotation algorithm to do something like what EOS-EDC
does for the production of absolute DEMs. Can you supply information on
how this is best done? Thanks for your help.
John Childs
Hi John, We believe that the errors in
elevation that your users are seeing are due to the DEMs being relative.
We have found our relative DEMs to have a decent vertical accuracy;
some do not have good horizontal accuracy, which is due to the accuracy
of the satellite ephemeris data (used to geocode the DEM). We have seen
DEMs up to 1000 meters off in the horizontal direction. If your users
are comparing our DEM to an existing DEM which is accurate in the horizontal
direction, the elevation differences could be great depending on the terrain.
We tested the vertical accuracy of our relative DEMs by first determining
the horizontal error. We did this by comparing the DEM to an existing map
source and then moving the DEM horizontally until it fit over the map. We
then did the vertical comparison. I forwarded this message to one
or your users as well. Your question about the translation/rotation
algorithm would best be answered by PCI software support. The software we
use to generate the DEMs is PCI/Orthoengine, and their contact information
can be found through the link below:
http://www.pcigeomatics.com/support/support.html If you have
other questions, please let me know. All the best,
Brandy
Land Processes DAAC User Services
Hello-I have a question about the DEM generated from ASTER data using
GEOTIFF4,I am using GEOTIFF4 to convert ASTER data of the Peruvian Andes,
then converting the DEM to an ArcInfo GRID, for some reason when I view
it in ArcInfo, the grid does not seem to be in the correct projection,
even though the correct projection found in the .met file is assigned
to the grid. In order to check the placement of where the images should
be, I generated a point coverage of the four corner points contained in
the .met file. These points project correctly, however the resulting DEM
from GEOTIFF4 is projecting 10,000km to the north of the corner points
where it should be projecting. Do you have any suggestions as to where
I am going wrong?
Sorry to bug you with such questions, but any help would be gratefully
appreciated!
Thank you
Todd
Todd:
Sorry for the delay in my reply. The UTM system is a little strange
in that within a given UTM zone the coordinates are not unique. The UTM
latitude coordinate in the northern hemisphere starts at 0 at the equator
and increases to 10,000,000m near the north pole. In the southern hemisphere,
the coordinates start at 0 near the south pole and increase to 10,00,000
at the equator. This leads to the strange situation where the same point
on the equator has two latitudes associated with it, 0 and 10,000,000.
ArcInfo is assuming that your southern hemisphere latitude coordinate
is in the northern hemisphere. When this error occurs, your latitude is
off by exactly 10,000,000m, regardless of the location of the point in
question. There must be a switch somewhere in ArcInfo that allows you
to specify that the coordinate transformation calculation shoud correspond
to a southern hemisphere location. There is no way to encode this in the
UTM coordinates themselves.
Thus, this is one of the rare instances where the problem is not my
fault (I think).
John
Hello, I'm trying to construct a DEM from an ASTER image (bands 3N
and 3B) using ERDAS 8.5.1 Otho-Base Pro. I've serious problems in finding
the values for side insidence (nadir and backwards) from the hdf or met
files.I found on the web that the keyword for side incidence is "POINTINGANGLE1"
but still no progress has been made.(In order to read hdf files I use
HDFView Version 1.2). Can you, please, help me out. Thank you in advance
for your time. Antoniou.
I do not have any experience with making DEMs from ASTER L1A/B but I
just opened a sample ASTER L1A file in my ASTER directory with my text
editor and found it right away in the metadata section at the very end
of the file: OBJECT = POINTINGANGLE CLASS = "3" NUM_VAL = 1 VALUE = -8.557000
END_OBJECT = POINTINGANGLE Use your text editor to search the string "=
POINTINGANGLE" and it should come up. Your reference to POINTINGANGLE1
seems to be in error, as the CLASS is indicated as 3. Let me know if this
is what you are looking for. John
Top
ASTER L1A/L1B
Thank you for your useful advice. I ordered
the unprocessed data and got the data within 12 days !
I got informed, that *.hdf.met and *.hdf were available in a certain
directory.
I was surprised to find out, that there were no corresponding *.tif files
on the other server (152.61.128.25/pub/asterdem/relative, as you mentioned
in your Aster Dem Download Procedure).
So I decided to download some of the huge *.hdf files (110 Megabytes each
) , hoping, I would be able to find (or program) a suitable converter
to *.tif format later. Although I have a highspeed DSL internet access
, there was not enough time to download all available *hdf files before
the deadline given by EOS. So I had to make a better selection.
Here I ran into a second obstacle: I was rather inconvenient to get the
"percent cloudcover" data for each granule. I had already tried
to make a good selection while ordering the granules, but the text-only
output (for Excel processing) of the selection pages did not contain any
data in the cloud-coverage attribute . I send a request to EOS about this
(see below) and got an answer, which I enclose. The cloudcoverage data
was available at the time, when browsing manually through the attribute
pages of each granule.
This data is vital for Germany, as I found out, that less than about 10%
of the granules returned have less than 10% cloud cover
As I am encouraged to dig for more high resolution elevation data of Germany
, here is my main question:
Do you know of a way to get to the *.tif data of my ordered granules,
or how I could convert the *.hdf files to *.tif format locally on my harddisk
?
Dear Sirs:
I found a flaw in your Search engine.
I wanted to get an overview of the available Aster Dem unprocessed LA1
data of central Germany and set up a query for lat 48 to 51 /longitude
7 to 11 and got a set of more than 100 granules returned. I transferred
this result to its own folder and modified the selected fields to id,
area coordinates and cloudcover in order to have a quality attribute for
the selection of the subsequent process orders . Although the cloudcover
is known in the individual metadata of each granule, in my above query
it remained blank. Is that a known problem and is there a workaraound
?
Second question :
Is there a way to request the view data (as another quality attribute)
of a set of granules, rather than requesting them individually ?
Thank you very much for your great website and your support !
Answer from edcdaac@usgs.gov
Subject Incorrect results returned for criteria entered
Thank you for your comments. There are certain data types, ASTER Level
1A and 1B included, that cannot display the cloud cover amount in the
Results List table. We are aware of this inconvenience, but unfortunately,
no changes are expected to be made to this feature at this time.
We have requested a multiple browse feature, but there is no estimated
time as to when it will be available. In the mean-time, it may be a little
more convenient if you FTP the browse images to your directory and preview
them in a HDF software. Select Request Sample from the Granule List and
fill in the appropriate information. You should only have to fill in your
information once, but you will have to click on Request Sample for each
granule. This may or may not be a faster method.
If you have any other questions or concerns, please feel free to contact
us.
Thank you and best regards,
Brandy
John added:
I think that the EOS reply sums it up. Do not forget
however that the ASTER DEM granules do offer both numerical cloud cover
data and browse images. One possibility might be to check the ASTER
DEM archive first, examine the browse there, and then locate the corresponding
L1A or L1B product.
Unfortunately, there is not a DEM for every L1A image
so this will not work all the time.
In answer to your second question, the program MultiSpec,
mentioned in the L1A article on my website can do exactly the conversion
that you require. That is how I prepared the images on the page.
The link is listed in the article.
Top
AutoCAD
I need to convert the stds dem to a format that can
be read in AutoCAD 2000i, any suggestions?
I do not know enough about AutoCAD's 3D capabilities
(particularly the way that it represents surfaces) to give you the answer
you are probably looking for. It may be possible to convert DEM
data to DXF format such that it can be rendered and manipulated directly
in AutoCAD but I doubt it.
However, it should be fairly straight forward to create
a raster image externally, import the image into a layer, and then use
the raster map as a base for further mapmaking by overlaying vector data
that you create in AutoCAD. (One problem is that you will not be
able to change the viewing perspective after you have imported the DEM.)
As you probably know, a DEM is a raster file, as opposed
to the vector data that AutoCAD and other CAD programs normally handle.
So the data structure is basically an array of elevation values arranged
on a rectangular pattern. The data object is the node. There
is one node per array location. The number in the array location
is proportional to or equal to the elevation of the terrain (z-value)
above the horizontal datum. The x and y coordinates are not included
in the array but can be inferred because the grid is regular and the grid
spacing and corner UTM coordinates are indicated in the DEM modules.
You must therefore import the raster data into AutoCad.
AutoCad 2000 apparently has the ability to import raster data in various
formats.
In order to create a useful map, you must first map
the numerical DEM data to a projection plane with the proper location,
lighting, etc. so as to create the image that you want. This image
must then be converted to a file format that AutoCad can handle.
The easiest way to do this is to use the program MicroDem
(see my web page for download information). Open your DEM as
a "New DEM" under the "File" menu. Choose the
"View" menu tab once your DEM has loaded and select "Open
GL" and "entire map". You should see a 3D rendering
of the DEM. Manipulate the viewpoint to get the proper perspective.
Check the "lighting" box if you want to get rid of the colors.
Adjust the sun position. Save the image as a .BMP file by right
clicking on the image. This file can now be imported into AutoCad
and serve as the basis of your map.
Alternatively, you could use DEMDRAPE and POV-Ray
or SDTSTER2 and Terragen to do the same thing with more flexibility but
it takes a little more expertise to work these programs correctly.
I have attached an image (Ticaboo Mesa, UT 7.5'
DEM) that I generated this way to show what I mean. A good link
that explains how to import raster data into AutoCAD if you do not already
know how to do this is:
http://www.cadenceweb.com/1998/1298/circles1298.html
Let me know if this does not make sense or if you need
help. Good luck.
We are trying to find a source for
mapping that we can use in AutoCAD 2002. What we are doing is taking the
mapping and putting things like city boundaries, sewer lines, water lines.
We need to be able to plot these out at various scales like 1” = 200 feet;
1” = 2000 feet; 1” = 50 feet. We will use these maps for design stages
of sewer and water projects. The contour lines need to adjust according
to the scale i.e. the line the same weight no matter what scale we plot
it. We need to plot it in color, similar to what a USGS map looks like.
These maps will be used from design and presentation to framing and putting
on the City hall wall. Can you help us? The key is it has to be editable
in AutoCAD and the line quality sharp.
I have check with ADCi and looked
through their stuff. It looks good but the printed product does not look
like the USGS maps.
We purchased DeLorne mapping software:
Xmap, Export module for XMap and the 3-D TopoQuads for West Virginia.
They almost work except that the print out is jagged and rough the lines
are not flexible with the scale used. They have big wide lines on smaller
scales such as 1” 50 feet.
Thank you for any help you can give.
I talked to the GIS specialist that comes into
our company once a week. I looked at his work and I am sure that
he is producing the type of publication-grade maps that you are interested
in. After talking to him, I am convinced that the only way you are
going to get the kind of results that you want is to use Arc View/Arc
Info from ESRI. I talked to him about doing what he did in AutoCad
and he said forget it, that AutoCad does not even come close to
the ESRI product in terms of capability.
As far as your idea of purchasing a map database and
then overlaying your vectors: I do not think this is possible at this
time. The reason is that there is a fundamental conflict between
the behavior of rasters and vectors upon scaling and other operations.
What our man does works approximately as follows:
First he obtains raster data that represents his mapping base. This
is typically a USGS topo map or an aerial photo. He imports the raster
into Arc Info and converts the raster image to a vector image using the
automated tools in Arc Info. By vectors, I actually mean Arc objects
such as points, lines and polygons which are mathematically speaking vectors
with spatial and other attributes. A road, for example will be contoured
and converted to an Arc polygon object. This object can then be
shaded and entered into the Arc Info database. The same for buildings,
vegetated areas, etc. The result is an array of Arc objects that
serve as a base for subsequent overlay of more vector data, such as fence
lines, utilities, etc.
He then imports the layers to Arc View which apparently
has a different set of tools that are used to create the finished map.
He said that it is possible to generate vector overlays in AutoCad and
then import them to Arc View or to work directly in Arc View.
The result is a vector map that is of flawless quality
and that behaves correctly upon scaling because every object in the map
is a vector object, even if it looks like a raster.
Anyway, I think I got this approximately correct.
I recommend visiting the ESRI website and talking to one of their sales
folks. The software is not cheap but it sounds like what you may
need to produce the kind of maps you described.
T hanks for your help. I called ADCi
and they assured me I could buy their maps and do the same in AutoCAD.
I think since they claim I have 30 days to try it that it may be cheaper
that buying new software. Have you had any experience with their Tiger
maps? I know that I have been a lot of trouble for you and again I want
to thank you for all your help.
Top
BLACKART
Hi John, I really enjoy your site. Keep it up! I was running through
your how to on using WinTopo and BLACKART to generate dems from topos. I
followed your instructions meticulously through WinTopo and successfully
saved out my vectorized topo with elevations as a arcView .shp file. When
I ran it through BLACKART , though, I kept getting the error "Invalid floating
point operation". I've tried various different combinations for the "LSQR
iterations" and "Laplacian Iterations", have regenerated my WinTopo data
sets several times (!), and have reread your how to several times in search
of my error, all to no avail. I was wondering if you had any insight on
what I might be doing wrong. Thanks in advance for any help. Cheers, Matt
Matt: I checked out your attached files as best I could.
Since you did not include the base .tif file I could not load it into
my demo copy of WinTopo. So I proceeded directly to testing it with BLACKART
and found the following: 1) The file produced the same error that you
reported when I attempted to run it. 2) I determined that the problem
was due to a mismatch between the number of records that were indicated
in the index file (this is the .shx file) and the number of records found
in the .dbf file. The .shx file indicted that there were 30 records. (The
program determines this by counting all the bytes in the .shx file, subtracting
100 bytes for the header section, and then dividing the difference by
8 as there are 8 bytes per index file record.) Inspection of the .dbf
file in a word processor will indicate that there are only 27 records
present (the record number is listed at the start of each record). As
a result, the .dbf file reader read past the end of the file and this
caused the immediate grief. (There also appear to be the correct 30 records
in the .shp file). 3) I added error handling for this case, instructing
the program to return zero if the dbf read attempt failed. This eliminated
the fp error that was occurring during the dbf read cycle. 4) Unfortunately,
the file did not get much further as another fp error occured later during
the lsqr computation for reasons that I do not understand right now. At
this point, I am attributing the problems to a faulty or corrupt input
file. As for the problems with the USGS ASCII output file, I can not offer
much help since I do not have the input file to test. I have attached
a simple ASCII input that you can run if this is any help. This file both
runs and produces a file readable by 3DEM. The ASCII file writer is one
that I have used in several programs previously includeing GEOTIFF4 so
it should work fine. It is not necessary that the file be any particular
size or even square in order to work properly. A couple of things to try:
1) I am assuming that you selected 'USGS ASCII' and not 'ASCII' for the
output. 2) I used WinTopo Pro v2.5 to generate my test files. If you have
an earlier version, you might try downloading the latest version. 3) I
uploaded a revised BLACKART to the web site in which I included the source
.tif file for the shape file. Try loading the test file and experimenting
with it. Perhaps this will help you get your own files running properly.
That is all that I have for now. If I think of anything else I will let
you know or post it. John
Hello,
The image is a simple example to try it. Sorry if my
english is not well, when I make a DEM from WinTopoPro and put in the
BlackArt throw and error. Do you know what I make bad? Thanks for all.
PD: And one another question, Do you know where I can download a Spain
DEM with 30m precision?. I search it but I only find GTOPO30, DTED 0 but
the precision of this files is badly. Regards.
Marc: Su ingles es mas mejor que mi espanol. Hablo solamente
unas pocas palabras de espanol y no puedo escribir tambien. The problem
is that your input file is corrupt. Specifically, the .shp file says that
there is only one row and one column in the input file, and sets up the
input array accordingly. When the program subsequently tries to write
all the data into the array, it over writes the array and causes the error.
The problem seems to be with WinTopo. This can be seen if you try to reload
your input .shp file back into WinTopo. A similar exception results and
WinTopo crashes. I do not know what it is about the geometry of your input
file that WinTopo does not like. Try another file and see if your luck
is better. Regarding your question about Spain 30m DEMs, the only source
is ASTER. Check the ASTER coverage map at http://edcdaac.usgs.gov/aster/dem_map_fullimage.html.
A few DEMs of Spain are available. See my website for full ASTER instructions.
Best regards, John.
Dear John
My name is David Maldonado and recently I have finished a
mesh terrain scenery of my country Venezuela, using data SMRT. I used your
program BLACKART to correct these data. I think that it's the best tool
that exists for this !!. In this respect I have some questions about your
program that I would thank a lot if you can respond them: What it´s the
meanings of LSQR? When I should use the LSQR Iterations? Is it good to use
LSQR Iterations and Laplacian Iterations at the same time? in which proportion?
How can I generate files DEM that can be processed directly by the Microsoft
SDK tools ? How can I visualize the coordinated X Y (lat, long) of a point,
while I draw an elevation?? Thank you, for your excellent program BLACKART
Greetings, David Maldonado
David: Thank you for your email. The point
you ask about is a little difficult to explain but is an excellent question
I will try to discuss in some detail. The following discussion relates mainly
to the problem of interpolating contour lines to DEMs. Filling holes in
DEMs is a different sort of problem and is discussed separately at the end
of my comments.
In general, BLACKART formulates the interpolation problem
according to the Franklin approximation algorithm by setting up a large
number of simultaneous linear equations, one per elevation node. The equations
are simply the Laplacian stencils (not to be confused with the Laplacian
solver discussed below), which simply says that each node is the average
of its four immediate neighbors. Then another set of equations is added
to the list, which says that each contour line node is equal to the elevation
of that contour line. This will always yield more equations than unknowns,
which of course has no exact solution. The Franklin algorithm then says
compute the "best" solution using a least squares approach. The LSQR algorithm,
which is a well-known iterative least squares solver developed by Dr. Saunders
and Dr. Paige at Stamford University is used to do this. LSQR is an iterative
solver so it starts with an initial estimate of the solution vector and
uses the initial estimate to compute the next iteration of the solution
vector, and then repeats until a solution of sufficient precision is achieved.
The closer the initial estimate is to the solution vector, the less iterations
it will take to compute a solution (actually an approximation rather than
an exact solution) if the process converges.
The solution computed by LSQR
yields a very good interpolation because of the way that the Franklin algorithm
formulates the equations. By this I mostly mean no terracing, although the
Franklin algorithm also avoids other unwanted artifacts as well. If no initial
estimate is supplied to the LSQR solver, it will simply use the zero vector
as its initial estimate and go from there. However, one of the problems
with this is that it can take a long time to converge to a high precision
solution. I determined during the course of my tests at RPI that the solution
time can be greatly reduced if a good initial estimate is supplied. I found
that the best way to estimate the initial solution is to use what is called
a Laplacian solver.
The Laplacian method formulates the problem similarly
to the Franklin formulation except that the nodes are not required to equal
the average of their four neighbors and to the contour elevations simultaneously.
This results in a system where the number of equations exactly equals the
number of unknowns, which can be solved (directly) much more quickly with
much less storage. The interpolation it computes is not very good for producing
a finished product when you are interpolating contour lines. However, it
is quite good for use as an initial estimate. (Many other interpolation
utilities use the Laplacian as their sole solution algorithm.) As a result,
the optimal solution time can be obtained by first running the Laplacian
solver to compute a good initial estimate, and then supplying this initial
estimate to the LSQR solver for the final interpolation. BLACKART will ask
you how many of each type of iteration you wish to run. There is no exact
formula for the optimum combination of each type. In general each Laplacian
iteration computes ten times faster than the corresponding LSQR iteration
so it usually pays to run a sufficiently high number of Laplacian iterations
and then run "enough" LSQR iterations, as judged by the result. You can
tell that you are running sufficient Laplacian iterations by running zero
LSQR iterations and then looking at the result. If it looks like a reasonable
DEM, you are probably running enough.
Because it is so fast, the Laplacian
solver is pretty good for jobs like patching holes in DEMs. In this case,
artifacts like terracing due to contour line remnants are not an issue and
the Laplacian solver often does very well with no LSQR iterations at all.
The LSQR may give a smoother interpolation but I would always try the Laplacian
first with zero LSQR iterations as this may be good enough. See additional
comments below. >
What it´s the meanings of LSQR? See above.
When I should use the LSQR Iterations?
Principally when interpolating
contour lines. When patching DEMs it may not be necessary to run any LSQR
iteratins. >
Is it good to use LSQR Iterations and Laplacian Iterations
at the same time? in which proportion?
When interpolating contour
lines it is always advisable to run both at the same time. BLACKART will
first compute a Laplacian solution using the specified number of iterations
and then automatically use the result as the initial estimate for the LSQR
solver.
How can I generate files DEM that can be processed directly
by the Microsoft SDK tools ?
I am not sure that I understand the question
but if you save in .bsq format your data will be stored in a flat binary
I16 format which is easy to machine process.
How can I visualize the
coordinated X Y (lat, long) of a point, while I > draw an elevation?
Right now BLACKART does not report this but you can infer it. An SRTM file
has 1201 rows and 1201 columns and spans one degree of latitude and longitude.
So each screen unit (which is reported) is exactly 1/1201 degrees. I may
add a feature to report this over the weekend. Editor's Note: This has since
been added.
bil Format
I am a student at the University of Delaware majoring
in computer science.
For my Object Oriented Software Engineering class project I put together
a
simple .bil file 3D viewer using OpenGL. THe viewer works in the sense
that it can dynamically display (rotate, zoom, translate accross the
file) the 3D data at a fairly reasonable resolution with smooth shading
and lighting. I am currently reworking the project a little for my own
enjoyment and I came across your site today looking for more info on file
formats. I first would like to say that your site has been more helpful
then the USGS as far as I am concerned. I wish I came across it earlier.
I was wondering if you would be interested in answering a few questions
I
have regarding the different file formats.
I used the .bil format for the project and I was successful at correctly
displaying the model but I'm unsure what elevation each point
represents. I understand from the metafile that (at least for the file
I
was using) each point is represented as a 16 bit integer in meters and
the
average elevation was around 1500. The file was one of the samples of
St.
Louis. From other sources I determined that the average elevation was
approx.
500 feet MSL so what point does the .bil file measured from? Have I
misunderstood the documentation? Also, the files seem to be also offered
in
a floating point format, where can I obtain exact representation of the
file.
I would like to make this program as robust as possible. Should I be using
a different file format other that the .bil? The program just needs simple
elevation in raster file format. Do you have any suggestions?
I had never heard of the .bil format before your email so I had to look
it
up to see what it was. I was unable to find anything that definitely
answered your question, unfortunately.
I do not think that your problem is the result of some odd vertical datum.
Elevations are almost always measured with respect to MSL or to some
geoid-based datum. These do not differ from each other by anything
near
1000 feet.
When I create image-type files from DEMs I always perform a linear
transformation of the elevation data so that the minimum elevation in
that
particular dataset maps to zero and the maximum maps to a value defined
by
the upper boundary of the integer space of the target file format.
For
example, if the elevations in the target file are to be represented by
two
byte unsigned integers, I map the minimum elevation to zero and
the maximum
to 65,535. (Actually, I might reserve 0x0000 and 0xffff for various
reasons
and span some subset of the available integer space.)
An unmapped image will be equivalent mathematically to the mapped image
but
the unmapped will look odd if it is a viewable format like GIF or TGA
if all
the elevations happened to be between 250 and 500, for example.
Also, any
smoothing function imposed by the target application will have more to
work
with if the elevations span the entire integer space. Often the
target
application is not particularly interested in the absolute elevations
anyway
and allows you to scale them arbitrarily for best appearance. Perhaps
this
is what happened to your data.
When I do this I usually report the slope and intercept so that the image
integers can be mapped back to the source elevations if this is necessary.
You might experiment if you have some known elevations and see if you
can
infer the slope and intercept of the transformation function. Hopefully
it
is not some more complicated offset. Or perhaps there are mapping
function
parameters hidden somewhere in the metadata of your file.
As for the file format question, USGS SDTS format is the most important
as
far as I am concerned because DEM data for all of the United States is
available for 30m and 10m resolution in this format. The problem
is that
SDTS is not a fixed-format file and is very unfriendly to programmers
as
compared to every other format. An SDTS profile consists of 18 .ddf
modules, many of which contain information that is required to display
an
image. You have to parse the Data Descriptive Record of each file
of
interest in order to get information about the Data Descriptive Area in
order to get information about the Data Record. The structure is
complicated and the USGS has done a horrific job of documenting the format
when you consider all the money that has been spent on creating it.
As for your question about floating point I assume you are referring to
the
.bil file. You need to find the file format specification that tells
you
how to decode the FPs because there are different ways to represent them.
Almost everybody uses IEEE 754 format but they can be 32 bit or 64 bit,
big
endian or little endian so you will have to do some digging on the
Internet
to get a definitive answer.
Top
CDED/GLOBE Format
My name is Joachim I am webmaster of German GPS Software
TOURATECH-QV. We integrated the possibility to add the Globe data to
our software, so one gets an elevation value as the cursor cruises
over a loaded topographical map.
I confess I know nothing about those various elevation file data
formats, and the amount of different formats available is more
confusing then enlightening me. SO maybe my question is obvious or
stupid, but I hope for an answer: We tried to integrate reading/using
those USGS/Canada DEM/CDED data too, but found out that this format has
significant disadvantages in online processing during use of a map.
Below you find a link to the sample files of Canadian CDED Data:
http://www.cits.rncan.gc.ca/siteCITS/servlet/CITS?site_id=01&sessionid=&page_id=1-006-002-001
So my question: Is there a tool (maybe commercially) available to
convert DEM/CDED data to Globe format including output of the
*.hdr(ESRI) header file ?
I hope this question is not too stupid and that you find the time to
send me a short answer.
Thank you very much for your help !
Dear Joachim:
I did some investigation into your request and discovered the following:
The CDED file format is the same as the USGS native DEM format so there
are
lots of programs that can read this type of file directly and output files
of various other formats. Unfortunately, the GLOBE format is not
one of
them. Building a GLOBE file is fairly easily and the directions
are given
on the GLOBE website. However I think that it is not practical
to
convert DEM files to the GLOBE format because of a basic compatibility
problem: USGS DEM files can have a variable number of rows and columns
of
elevation data, for example 1201 rows by 1201 columns is a common
arrangement. But according to the information at the GLOBE website
file
format page:
http://www.ngdc.noaa.gov/seg/topo/report/s11/s11C.html GLOBE
files cannot. There are only two types of GLOBE tiles allowed I
think. One
type has 10800 rows by 4800 columns and the other has 10800 rows by 6000
columns.
I did not see any information that would allow for varying the numbers
of
rows or columns. GLOBE readers probably decode the tile name (for
example
B10G) and determine which of the two possible tile types to expect from
a
lookup table. It is then a straight forward job to read the signed
16 bit
integer data stored in row major order (all the data for row one, followed
by all the data for row two, etc.) once you know how many columns there
are.
I can find no provision to store 1201 rows by 1201 columns in GLOBE format
unless the file is padded with 41,757,599 zeros, or unless lots of DEM
files
were packed into a single GLOBE tile which is probably not practical.
Perhaps you can convert the CDED format to some other format that is
friendly for data processing in your application. It is possible
to store
DEM elevation data very efficiently in TGA format, for example as one,
two
or three byte unsigned integers, either compressed or uncompressed.
I have
written a couple of programs to convert DEMS to this format because many
rendering applications can accept elevation data this way. The problem
with
this approach is that the metadata is lost. I believe GEOTiff format
can
overcome this problem by accommodating metadata, as can SDTS, MicroDEM
DEM,
.TER and some other formats.
If this does not confuse you more and you would like more information
on any
of these alternative file options, please write to me again and I will
send
information on a specific format.
I do have a problem. I downloaded a file of data
for Namibia from the GLOBE site. (not a tile). Unpacked (using microDEM
not WinZip) I got three files .hdr, .fmt .bin. Tried to import into microDem
but keep getting "Invalid floating point operation" error. I've
tried extracting various smaller chunks but get the same error. Data looks
OK in a binary editor.
I also downloaded the entire southern African tile.
(going to have a big phone bill) but when I try the compressed tile .gz
microDEM crashes altogether.
This file is OK as it opens fine in 3DEM.
Any Ideas?
I have no experience with GLOBE data as I prefer the
NIMA DTED0 dataset. I extracted the following instructions from
MicroDem v5.1:
GLOBE 30 sec DEM
GLOBE is a global digital elevation model (DEM) with
a horizontal grid spacing of 30 arc seconds (approximately 1 kilometer).
It used both DTED Level 0 and GTOPO30, so it might be superior to those
data sets.
-
You must download the data, and the GEO-VU headers,
to use the GLOBE data.
-
You must do an import of this data set before use.
-
You will have to subset the files for this data set
before use. The 100+ MB files GLOBE uses are too big to manipulate
and display efficiently. If you wanted to see the “big picture” you
should use a data set like ETOPO5. The import automatically does the
subset.
If the data is correctly subset, it will use a DIX index
file and support merge on the fly operations.
GLOBE is available on line at:
http://www.ngdc.noaa.gov/seg/topo/globeget.shtml
GTOPO30/GLOBE Data Set Import
GTOPO30 and GLOBE are global data set of land elevations
with 30” spacing.
-
Select File, Data Manipulation.
-
Select Import, GTOPO30.
-
Select the file you want. They are named based on
the NW corner; you will pick the HDR file and the program will open
both the header and the accompanying DEM file. The window will show
the corners of this data set, and the maximum size DEM that you can
import. If you need a DEM covering a larger area, you should consider
a smaller scale DEM like ETOPO5.
-
Select the NW corner you want.
-
Select the SE corner you want.
-
Give the file a name.
-
The import will be done, and the program will show
you the directory where it was written.
If you need a larger area than will fit in one DEM, you
can extract pieces using the steps above, thin them, and then merge the
thinned files.
If you need an area that spans several files from GTOPO30
or GLOBE, you can extract pieces using the steps above, and then merge
the resulting files.
ETOPO, unlike GTOPO30 and GLOBE, has elevation data over
both land and water. GTOPO30 and GLOBE cover only the continents.
Have you tried these instructions already? If so
let me know and I will try to think of something else.
What an amazing web site you have. I am interested
in producing maps in
western Canada that show terrain relief and shading. I have a software
program called Surfer that allows me to import DEM's and then create DTMs.
These can be shaded to build relief maps etc.
Have you been successful in building the maps you show in Canada?. Are
DEM's
available here? What I would like to be able to is build the kind of map
(in
Canada) that you have shown in your section on 3DEM Overlays.
My objective is to build an map of an area with UTM co-ordinates and then
provide a trail through the area for backcountry use with GPS.
Thanks for considering my questions and for the great site!
You can build digital terrain maps for Canada. Unfortunately the
Canadian
government, like most of the world does not yet make its spatial data
available for free. You can purchase DEMs in CDED format, which
is exactly
the same as USGS ASCII native DEM format at :
http://www.cits.rncan.gc.ca/siteCITS/servlet/CITS/site_id=01&sessionid=&page_id=1-006-002.html
_id=1-006-002-001
Any application that accepts USGS ASCII DEM should accept CDED format,
including some of the converters on my page.
Here is another link for free DEM data for the Yukon Territory:
http://renres.gov.yk.ca/pubs/rrgis/data/data_desc/90m_dem.html
You can always get GTOPO30 or NIMA DTED0 data for free, but the resolution
is poor at 1000m postings. Look for information on my page.
You can probably find more sources by exploring the web.
Top
Converters
Hello, I'm having a problem with all of the converters
I have downloaded from your site, and they all respond with cannot open
file. I'm running win98 with a 500mmx processor,128 meg ram. Can you help
me out.
None of the converters will accept a file name with
greater than 8 characters. Rename the input file using Explorer
if necessary.
Top
DEMs from Imagery
I was wondering if you could steer me to information
on how to create
DEMs from imagery. Are there any software products out there, for a
relatively low cost, that can perform that kind of conversion? (I'm
guessing that the big GIS packages might have an interactive module for
creating elevation products, or DEMs).
My friend Paul recommends the software from PCI Geomantic. You can
take a
look at http://www.pcigeomatics.com/.
He told me the software runs $800,
which is a lot but not out of line with first rate GIS packages.
I think
this package is pretty much the standard.
I tried to look at Virtuozo, an application offered by a Chinese company.
Their website is at http://www.supresoft.com/.
They appear capable and
offer a free demo but I was never able to get set up with their license
key
because of a network card id problem so I do not know how it works.
Top
DEM to Contours
I frequent your interesting and informative
website from time to time on the search for useful software and new developments.
I wonder if you might be able to point me in the right direction in my
search for a freeware / shareware software pack which could do the following:
derive elevation contour sets from DEM data
output these contours to .SHP files compatible with Arc View
Will look forward to a reply and please do keep up the great work
on your webpage
It's not exactly free (~$130.00) but GlobalMapper is the only program I
know that will (probably) do the job. You just load up the dem and
click "generate contours", then export shapefile, selecting "export
lines". You might send Mike Childs (no relation) at GlobalMapper
an email and check this out with him. I have attached a sample contour
and shape file generated by the program.
Top
DTED Data
I need to obtain DTED level 1 and DFAD level 1 data
for Japan. Specifically,
N39 E139 to
N43 E144.
I work for a company contracted by the US government to produce terrains
for
flight simulators. Getting the required government approval is not a
problem, paying for the data is not a problem, FINDING the data in the
proper format (DTED, DFAD level 1) is proving to be a problem.
Any assistance you provide, (a point in the right direction?), would be
greatly appreciated.
Thank you for your time, in the mean time I'll be digging through your
web
pages and references.
I'm still a newbie here, but I successfully retrieved some SDTS data
of
interest and want to convert it to DTED, because I have some sample OpenGL
programs that will take DTED as input.
How do I get from SDTS to DTED? What's the relationship between DTED and
DEM?
Oh, on another topic, is there elevation data for Mt. Everest?
I do not know of any SDTS to DTED converter. DTED and
USGS DEM are similar in that they both contain elevation data and metadata,
but in different formats.
Both DEM and DTED are somewhat messy formats to read/write to. While
I already have a DEM reader that I use in other converters, I do not have
a DTED writer that I could combine it with to make a quick converter.
The file format specifications for both are available. Let me know
if you are interested in them and I will send you the links
There are many free and low cost viewers out there
that read SDTS directly. Why don't you see if one of
these will do the job for you. MicroDEM and 3DEM are
two examples. Their links are on my web page. A
little searching on the web will produce others. Let
me know if you need help in finding them.
As far as Everest data, the best available at this
moment is NIMA DTED0. I usually import the DTED into
MicroDEM and then go from there. Detailed
instructions for getting at this data are also shown
on the web page, in the International DEMs and Sinai
sections.
I am trying to find dted 1 data for the manchester
area of the uk. to use it for planning a private wireless network set
up between friends.
Do you know of any sites that I can download
this information for free? I have been trying to use radio mobile but
the data is not accurate
can you help?
NIMA DTED1 (90m resolution) is currently classified.
The best that is available is DTED0 (1km resolution) from NGA Raster Roam.
You might also check out the UK Ordnance Survey.
The data is not free but may not be very expensive. The advantage
is that it will be of excellent quality and accessible. I have included
the text of an earlier request regarding Ordnance Survey data for reference:
Do you know of a converter to convert DTED into
DEM files? I can't seem to find one.
The only application that I know that is able to read
DTED data is MicroDem,
available for free at
http://www.nadn.navy.mil/Users/oceano/pguth/website/microdemdown.htm.
Unfortunately, this program only exports in a native MicroDem format.
I
wrote a converter called MDEM2DEM that is available from my website that
can
convert MicroDem to USGS ASCII DEM.
A big problem is that MicroDem "loses" most
of the metadata upon export.
The converter inserts dummy data just to get the image to render.
You must carefully
inspect the type A record in the USGS DEM file and insert the metadata
manually.
You might also look at NIMAMUSE at
http://www.nima.mil/geospatial/SW_TOOLS/NIMAMUSE/doc/apps/intro/intro.htm
I
have not looked at it myself but it might have an export format you can
use.
Top
ETOPO5
Would you know if there is a freely available DEM
of the entire US, excluding Hawaii and Alaska? I am aware of the GTOPO30
site, but those are tiled, and several are needed to get the entire country.
I don't need high resolution.
If you want a DEM, the GTOPO may be optimal. I
believe that the ETOPO dataset is coarser (5') so you may want to check
this out. Here is a starting link:
http://evlweb.eecs.uic.edu/pape/vrml/etopo/mapurl.html.
If you just want a rendering, check out the National Elevation Dataset.
http://edcnts14.cr.usgs.gov:81/Website/store/viewer.
Top
Flight Simulation
Like your site! Just wondering if you knew where
one could obtain some terrains in following file extensions, .DEM, . MAP,
and .PTT files of the same geographic area? I am part of flight
sim "development team" and we were curious to see if other terrains
could be brought into the game (Jane's USAF - Electronic Arts).
I noticed that some of the files that are available through the USGS were
brought into Microsoft Flight Sim. I also noted there were several
translator programs available to convert from one file type to the other.
If you are interested, the following thread is currently being weaved
at the following web-address:
http://www.simhq.com/simhq3/sims/boards/bbs/Forum49/HTML/001864.html
We are basically looking for new terrains to fly over.
The game comes supplied with parts of Nevada, Vietnam, Germany and Iraq.
For each terrain type there is a) terrain tiles, b) the elevation component,
and c) an aviation grade map complete with airport locations, etc. for
larger scale navigation. We are looking for all three. Good
luck??!!
This is just a hobby for me, I am not a programmer
by any stretch of the imagination. Anyway, if you are at all familiar
about the terrains used in flight sims (or at least with the file extensions
I have given above), I would love to hear from you!
I took a quick look at the link you enclosed. From
what I gathered there and from your email the terrain data is in at
least two files. The DEM of course is the digital elevation model
which is just a regular grid in the xy plane mapped to a z coordinate.
It is typical to "drape" a (2D) raster image over the DEM in
order to add features to the otherwise stark landscape. This information
may be in the MAP and PTT files. If they are TGAs this would make
sense because this is a 2D raster format commonly used for this purpose.
There are several examples of this type of draped image on my webpage.
DEM data for the United States is readily available from
gisdatadepot.com. (again, see my site for the link). Getting good
raster data for the drape is a more difficult challenge. USGS land
use DLGs are commonly used, as are aerial and satellite images.
Artificially generated texture maps are also used.
You might try contacting Mike Heir at heir@home.net.
He is a Flight Simulator type and he and I worked together to write the
DEM2BSQ converter to support the Flight Simulator community. Perhaps
he can offer some insight to your problems.
Thank you for your message and understanding.
Not knowing how vast and complex the problems with DEMs, I have been using
a
simple method of making dems in grey scale with PSP7. What
a wonderful new
world you have opened up.
I enclose one of my test .raw files which I use with Grises50 to turn
into
the DEM.
This one is 257x257 pixels and has 50 different grays. Darker at
the lower
altitudes and getting lighter as the hill rises. I test mine
with the
dropper in PSP.
As I only do my own local scenery, different .raw files seem to work well
but I am up against a big problem in that FS2002 does not accept LODs
above
11. To overcome this I make a .raw 257x257 with a very wide
darkest grey
and the hill rises near the centre. This I turn into the DEM
and make the
lower level just under the Microsoft ground level where I am locating
the
hill. This seems to work. In experimenting with
different ground heights
one can lower the scenery in FS2002 with this type of procedure.
Thanks for your message. I do not know a lot about
FS but I worked with a
fellow named Mike Heir to develop MDEM2BSQ. The idea is to import
DEM data
(from the USGS for example) into MicroDem, merge files if necessary, and
then save in flat binary format that he called .BSQ and may be the same
as
your .RAW format. He then uses these files as input to FS for his
terrains.
Perhaps this would work for you. The advantage is that MicroDem
fairly
easily ingests many DEM formats.
Top
Hexedit
I have been searching the Internet, and download
sites upon it, for a type
of bitmap editor program, with little success, when I arrived at your
website, and it showed programs that are closer to what I have been looking
for than anything else I have found, and I feel that you may be able to
understand the type of program that I am looking for and may be able to
help
me or suggest something.
I have been searching for a very basic type of bitmap editor which can
perform the sort of functions as these examples:
* Convert numerical data (either with direct input or text files) into
simple digital bitmap (.BMP) files: An input such as (255,255,255)
or
"FFFFFF" would result in a white pixel; "FFOOOOO",
a red pixel; a text data
could be loaded, in which text consisting of continuous "000000"
field with
a "FFFFFF" field located somewhere within would ultimately result
in the
program creating a bitmap consisting of a single white pixel on a black
background
* Convert bitmap files into numerical data; a blank white bitmap would
ultimately return data in a series of "FFFFFF" fields or usable
(255,255,255) data
* Allow editing of individual RGB values of pixels and changing the color
of
pixels by changing red, green, and blue values independently; selecting
a
red pixel would bring up a dialogue box with information like this: "red=255
green=0 blue=0" and provide the ability to increase the green value
to 255,
resulting in that individual pixel changing color from red to yellow
* Allow merging of bitmap files so that individual color values in pixels
would produce a pixel equal to the sum of the RGB values of the
corresponding pixels of each bitmap; if two bitmaps of like size (100x100)
were superimposed upon one another, and the RGB value at location 30,30
on
each bitmap was a dim red "770000", the result on the final
bitmap would
result in a bright red "FFOOOO" at location 30,30
The type of program that I am looking for is certainly much more basic
than
the ones shown on your web page, but it is because it is so basic that
I
have become as frustrated trying to find such a program. If you
know of any
simple program that offers any of the features I explained, please reply
with any suggestions. If you cannot help but know of someone who
may be
able to, please feel free to forward this e-mail intact to anyone who
may
have suggestions or recommendations.
I hear what you are saying. This kind of integrated
editing environment
would be a handy thing to have. I have never seen anything like
what you
are describing. I wonder if the higher end Graphics programs like
3DStudio
Max, etc. have this functionality.
I have been able to do pretty much everything you described using some
simple tools. The first is the book called "Encyclopedia of
Graphics File
Formats by James D. Murray and William van Ryper". I use this
book to help
decrypt various graphics file formats. I can then use a hex editor
like
UltraEdit, hexedit or XVI32 to display each byte of the file in as hex
number. I can
then edit the file by hand. I often have to do this when converting
an
under-specified GIS format to a graphics format where I do not know what
the
orientation of the decoded file is, whether it is mirror-imaged, etc.
I
usually tag the file by setting a few dozen pixels by hand in my hex editor
and then inspect it using PaintshopPro.
PaintshopPro does not have all of the tools you are probably looking for
either, but it does have the ability to display the color palette, along
with the hex value of each color. You can also perform arithmetical
operations on files by combining color values of pixels and so on.
Sorry I could not be of much help but that is how I handle these tasks.
Let
me know if you find anything better.
Actually, that program you told me about, XVI32, has been very useful
to me,
and I thank you for letting me know about it. I have been able to
manipulate numerical data using QBasic, output this data to a text file
and
cut-and-paste to a blank bitmap file, and I have this process working
fine.
It's allowed me to do some things I didn't think I'd be able to do (I've
found editing a bitmap using a text editor is impossible, for example).
It
also got me back to programming again, and I almost forgot how much I
liked
to program.
I'm thinking that I will need a viewer to check the data against my topo,
and photo data. Then I will need a 3D tool of some sort to manipulate the
terrain elevation data, preferably some kind of program where I can see
the changes I am making in real time, (a GUI) rather than entering huge
volumes of tabular data into a database is preferable.
My questions to you are as follows: What kind of
viewers would you recommend for this project? What
3D rendering/ imaging/ manipulation tool would you recommend for making
slight adjustments to the USGS data? I'm already going to have to
buy one for making 3D objects for the game, (hangers dams, roads, runways,
bridges) so recommend several if you would, and I'll try to find one that
will work for both issues for me. Is there a single
tool that will extract, view, and manipulate the USGS data that I want to
import into the game ? The short answer to your question is that
I do not know of an application that does exactly what you are looking for.
I use hex editors like Ultraedit and Hexedit to do the kind of manipulations
you are talking about (but I do not do this too much because it is so cumbersome).
To do this I simply navigate around the USGS ASCII DEM file or the SDTS
xxxxcel0.ddf file making adjustments using the hex editor and then view
my results using an application like GlobalMapper, MicroDEM, etc.
There are some applications that try to do what you are talking about, applications
such as WebWinds and Noesys. These are aimed at processing satellite images.
I am not sure they will handle SDTS format. I have experimented with
them using geotiff and hdf files and they and the data are extremely
difficult to work with. The data is problematic because DEM files
are just so huge. Wish I could be of more help. I will think about
this a little more and if I come up with anything I will let you know.
Top
IKONOS
Thanks for the beautiful and comprehensive article
about Digital Terrain Models that is posted on geocities.com.
In the part about IKONOS, you mentioned that you
have not been able to determine if IKONOS is capable of generating stereo
orthophoto pairs suitable for the production of DEMs. The answer
is that yes it can and companies and governments are indeed using it to
do so.
I have to correct myself here, IKONOS is capable of
collecting stereo pairs
of satellite images in two ways. The first one is the stereo coverage
resulting from the overlapping of two consecutive strips of imagery from
two
adjacent passes of the satellite over the same area. This kind of images
is
available to the public. The other kind of stereo coverage that
IKONOS
sensors offer is a stereo that is within the same strip. This is
produced
by collecting the image captured by the first line sensor, then separately
collecting the image captured by the third line sensor. This has
been
possible since the IKONOS camera has 3 bush-broom line sensors; the first
one is placed at a slight angle to look forward, the middle sensor
is
placed in a nadir position so it is always orthogonal to the terrain,
and
the last sensor is a back-ward looking sensor. The second type of
stereo
coverage is available only to Space Imaging and its government customers
and
can not be sold to the private industry.
Now we come to the correction part. I the geocities article you
mentioned
"stereo orthophoto pairs". In fact stereo pairs of any
photography can be
utilized to produce digital terrain models (DTMs) that can be used to
ortho
rectify the imagery and produce ortho photos. Two orthos of the
same area,
even if they had the same scale, will not give you stereo coverage because
all the distortions due to relief displacement have been removed from
them.
Top
Kitterman's Map
I'm working with ESRI software and Illustrator on
a map of a sub range of the Andes (Cordillera Huayhuash, Peru).
I am using mapublisher to import all my layers into illustrator. No
problem there. But I also want to create a shaded relief layer. Any suggestions
how to get that layer back into a drawing program and register it.
I've tried exporting an EPS greyscale DEM into
Photoshop and rendering it from there with OK results, I can get
a smooth looking terrain model that I can then use with hypsometric tinting
and shading underneath all my contours and line work, it looks superb
but... The problem is registration and the resultant size of this layer,
even as an 8bit grey scale it is extremely huge(50+ MEGs) unless
I drop the resolution and then it becomes very pixelated and not so good
for print output. Any idea how they got it to look so good on the Rainier
map. I think they are working in the same software environment as
me (MAPublisher). But maybe Op.Sys makes a difference PC vs. MAC? Or maybe
I am way of track, any advice would be appreciated and acknowledged in
final product.
I am not familiar with the applications that you are
using and your description of the DEM data type did not ring a bell either.
I am sure it is not the OS or computer that is the key to making a great
image. This is the way I see it:
1) A DEM is raster data of course and each elevation
posting is represented by a number that takes from one to four bytes to
represent. The size of the DEM as far as your application is concerned
is equal to the number of elevation postings times the number of bytes
per posting, plus a small amount of file overhead for any header information.
2) The overlay image is also raster data. Each pixel
is represented by from one to four bytes (usually) depending on the color
depth per pixel. The input file may be in a compressed state (RLE, etc.)
but it needs to be uncompressed by the application in order to process
it. The pixel topology of the overlay will not match the elevation posting
topology of the DEM, unless you manipulate things to cause this to be
so.
3) To make the image, the DEM and overlay need to be
mapped one to the other. The application will maintain the overlay and
DEM data separately, and then should create a third object which is the
DEM/overlay rendering. This should be represented by a data structure
with size of, at most the number of bytes of the overlay, depending on
whether the DEM is mapped to the overlay or the overlay is mapped to the
DEM.. Presumably, some interpolation or sampling algorithm is used to
accomplish this. If the DEM is mapped to the overlay, this is good. If
the overlay is mapped to the DEM, this is bad because sampling or especially
interpolation routines cause distortion of your image.
4) If this image is exported, the file size will depend
on the characteristics of the image and the type of file chosen to represent
it (jpeg, gif, etc.).
5) As layers inside an application, a large amount of
data may have to be maintained in order to allow arbitrary perspectives
to be displayed. For example, a 1201X1201 DEM will require about 2.88
Mbytes at two bytes/posting. A 32 bit uncompressed TGA image of 1024X768
would require 3.14 Mbytes for a total of about 6 Mbytes of data.
Mr. Kitterman's image was probably 4MB.
6) There are other factors that come into play as well.
For example, what is the quality of the source overlay image? An inferior,
muddy image rendered using a large pixel count will merely faithfully
reproduce the inherent infidelity of the source image.
This is what I suggest trying in order to make the
best quality image:
1) Plan out the size of the target image carefully
and choose the DEM and overlay appropriately. Don't use such
large source datasets that manipulating them is too cumbersome.
Remember however that the final image will be somewhat to much smaller
if output in compressed format.
2) Scan your own topo map using the highest quality
source possible.
3) Scan the image in at exactly the size that you
intend to use for the overlay. Any resizing of the image after scanning
will degrade it.
4) Match the pixel topology of the DEM exactly to that
of the scanned image. There are several ways to do this. One way would
be to write a program to manipulate and rewrite the DEM file but this
would take programming skills. Another way would be to convert the DEM
to a TGA, resize the TGA , import the TGA to POV-Ray as a height field
and then perform the overlay by using the POV-Ray utility for this purpose.
Your application may have a similar utility.
5) You are going to be working with large files for these
steps. Notice that Mr. Kitterman's image is only 448,422 bytes in its
final form but is 4.13 Mbytes when rewritten as a 24 bit TGA file..
6) Save the resulting image as a jpeg for maximum compression.
As far as registration of the overlay to the DEM, this
can be done by imposing small cues onto the DEM and matching them with
data on the overlay. (I give an example of this on my web
page under the Extended DEM). Hopefully things match well enough
because of your file prep so that these tweaks are minimal. At all
costs avoid "rubber sheeting" the overlay because you will lose
image quality. section.
Making an image like this is a lot of work but is probably
necessary to get the kind of results that Mr. Kitterman got.
Hello, my name is Ovidiu I.
and I visited your webpage. The area which
commanded my attention most was Mr. Kitterman's overlay. Could you tell
me
what software and data set he used? Any information would be greatly
appreciated. Thank you, Ovidiu.
I don't know much about the image other than the info
at the link
http://www.avenza.com/MPcomp/
. The software is apparently MAPublisher 4.0.
Charles Kitterman works for Kulshan Cartographic Services. The topo
map
used was apparently one generated by that company. The DEM dataset
was most
probably USGS.
Thanks for your very informative and useful site on
terrain modeling. I
have been experimenting with 3dem and have a question: In the second
to
last paragraph of your Mt. Rainier Challenge section you mentioned looking
at the image in a straight overhead view rather than in perspective
view. How do you view it from directly overhead and also generate
hill
shading? I can only seem to generate the hill shade showing terrain
when I
select >operation >3d view but then I have to view it from an angle.
The overhead view requires a little trick to work easily. First,
select 3D
view. Under the "terrain position" box select the "far"
radio button. Then
generate your image normally.
When the image has generated, select "Operation" and then "change
position".
Use the "tilt forward" button to rotate the image to the overhead
view
position.
Top
National Elevation Dataset
I am working with 3DEM and would like to use NED
data that I can get from USGS (mentioned in your DEM starter kit).
The
problem is I can't get the NED data into 3DEM, what format do you recommend
I get from USGS? the options are Grid-Float, ArcGrid, or Bil.
It seems much simpler to use NED data where you can get a large area rather
than piecing together multiple dems based on quads.
.bil format would definitely be the one to go with. It is a very
simple two
byte integer format that would be a breeze to convert to something useful,
like USGS DEM for example. The problem is that there is no .bil
to USGS
converters in existence that I know of. It is on my list to write
but I have
slowed down writing converters since school started for me.
Kim Kisner of the USGS told me via email that the USGS will be making
NED
data available in SDTS format "soon". In that case, 3DEM
will be able to
read NED data directly. I guess you will have to wait a bit for
now but
since you asked I might get motivated to whip up the converter in a hurry.
Keep checking my website and the USGS website for a solution one way or
the
other. If you want, you could contact Kim directly at Kisner@usgs.gov
for
more information regarding SDTS availability.
Top
Mars DTM
Any idea where one could obtain a specification for
the DTM format (Mars Digital
Terrain Model) or GTDR (Venus Global Terrain Data Record) format?
The document at
http://www-pdsimage.jpl.nasa.gov/jbcache/viking/vo_2007/document/volinfo.txt
seems to be a good starting place.
If I was deciphering the format I would first read through this document,
then load an appropriate data file from the corresponding ftp at
http://www-pdsimage.jpl.nasa.gov/jbcache/viking/vo_2007/
into a hex editor
and see if I could figure it out. If I could not figure out the
format from
these I would contact someone associated with the ftp and ask for help.
The
scientists that maintain these sites are usually pretty cooperative.
I got these links from
http://www.ccm.net/~jrsmith/dtm.html so contacting
Mr. Smith, who seems to have a lot of experience would probably also yield
valuable information.
Another possible link is the Mars
Polar Lander website.
Let me know if this leads anywhere in case someone else asks.
Top
Moon DTM
Do you know where I can find a DEM of the moon ??!!
I don't know of any source offhand but you could contact Dr. Margot at
this
link http://www.gps.caltech.edu/~margot/
and ask about how he obtained the
radar interferometery data from
http://dsnra.jpl.nasa.gov/
Top
MicroDEM
Thank you for your excellent Digital Terrain Modeling
site. It's resources are excellent.
I'm hoping you can help me solve a *somewhat* unique
issue I'm having with converting DEM files into a
useable Photoshop format.
I promise I'll try and keep this as short as possible
but I must provide a little history for you.
Over the past two weeks I've been trying to solve the
issue of how to convert DEM files into a useable
Photoshop format so I can a create shaded relief
overlay for the trail maps I produce here in Colorado.
I'm on a Windows 98 machine and I use Macromedia
Freehand to produce my maps.
So far I've collected nearly every DEM utility out
there and so far no luck. I obtain my DEM data from
GeoComm.com in the form of the current USGS
SDTS-format 7.5' DEM files. Using SDTSDEM2 I convert
those pesky files into the "native" USGS DEM format.
From that point I've used the following apps with the
corresponding results.
1) MicroDEM - Opens the DEM files perfect! Only they
show up in this little window that's 361x464 pixels
and I haven't been able to figure out how to export a
BMP, Tiff or TGA file that will match the physical
dimensions of a DRG file (GeoTiff) when opened in
Photoshop.
2) DEM2TGA7 - Seems to run fine but when I open the
resulting *.TGA file in Photoshop it looks like orange
snowstorm.
3) Wilbur - Same issue as MicroDEM
Am I going about this the wrong way?! Most of my maps
are output at a minimum size of 11"x17"! I must not be
doing something right. Can you give me a pointer or
two?
I am going to be shooting in the dark a little because I have no experience
with Macromedia Freehand or Photoshop but here goes.
1) The typical SDTS 7.5' DEM has about 480 rows by 380 columns (plus or
minus) of elevation postings. Each elevation posting gets mapped
to one
graphics pixel in a TGA or other graphics file. So the image size
that
MicroDem produces is exactly what you would expect. The larger scale
1:250,000 USGS DEMs have typically 1201 by 1201 elevation postings so
the
TGA produced by converting one of these would be about that size.
You can
of course increase the size of these by image processing (resizing) but
then
an interpolation occurs between pixels to create the new pixels.
2) You can export graphics files from MicroDem by selecting 'Save Image'
from the main menu. You get to select from BMP, JPG or GIF format.
3) This is the third report I have received of DEM2TGA7 output files
not
being read by Photoshop, so I believe it now. I think this is a
problem
with the Photoshop reader but have no easy way of resolving this because
I
do not have a copy of this program. (Anyone who has written a file
reader
knows how difficult it is to anticipate all possible flavors of an input
file allowable under the file specification.) I do know that PaintshopPro
will read the output files from DEM2TGA7. You can download a demo
copy of
PaintShop for free (see link on my page), import the file, export it to
any
format (including TGA if that is what you need), and then read it into
Photoshop. I have recommended this a couple of times before when
this
problem came up and have had no reports of problems in doing it this way.
If none of these help, please contact me again.
I have been trying for the last 2 days to get 4 24K
DDF into Bryce 4 (Bryce is able to import DEM's). I finally figured
out how to merge the 4 24K DDF's in MicroDEM and it looks great.
Also, I exported the merged file into a DEM. Verified this by reopening
the merged file into MicroDEM. But, Bryce complains with "Bad
I/O" error when I try to import the new DEM. Bryce will import
any of the original DDF's . . .
Any idea on how to rectify this? Any help would
be appreciated.
I would rather use a DEM as opposed to some
other file format . . .
While typing this, I realized I should just try downloading
some DEM modelers from ZDNET and see if they have a hard time importing
the merged DEM . . .
This is project for the chamber. I want to
start out with this small scene. Then, if I really have the ambition
and time, I want to do an animation/fly through of this area to put on
the website . . .
Question: You probably know this already but you
cannot import MicroDEM output files into Bryce directly, because MicroDEM
DEM format is different from USGS DEM format.
The second question is, if you converted the MicroDEM
output to USGS DEM format, how did you do it? Did you use my program?
If you could clarify these points I will supply more
information
Thanks for putting your web site together; it is a
wonderful resource!
You mention that MicroDEM will Merge two 7.5 min DEMs and indeed it behaves
as if it does. However, I am using a program recommended by the USGS (DEM
Works) to open the resulting file. DEM Works just ignores the file with
no
errors.
3D Studio Max is really where I would like to end up but when it tries
to
open the resulting MicroDEM 4.0 file it replies with "Improper file
format."
Finally there are another set of programs that written by Sol Katz that
list
the information in a DEM called D2X1.exe. For the MicroDEM merged file
the
parameters all come back as zero and the program hangs on the MircoDEM
generated file.
MicroDEM will open its own Merged file and display fine.
Do you know of a way to get a merged MicroDEM file in a valid format for
other programs?
The problem is that a MicroDem .dem file is not the same as a USGS .dem
file. MicroDem outputs files in its own native file format.
(I am always
amazed by the number of different GIS file formats that exist. It
is hard
to get serious in this field without becoming somewhat knowledgeable about
file format conversions, or as one of my former programming teachers called
it "shoveling dirt", i.e. you take it from here and shovel it
over there.)
According to my research, only two other programs are able to read a
MicroDem .dem directly: Daylon Leveller and DEMScape.
Regardless, I am hoping that the simple answer to your question is to
use
the "Save Image" selection in Microdem rather than the "Save
DEM". When you
use "Save Image", you get the choice of saving as a .BMP, .JPEG
or .GIF
format.
This is what I have done in the past to create POV-Ray input files.
I start
in MicroDem, import the DEM, save as a .BMP, import to Paintshop Pro,
and
then save from Painshop to .TGA format which POV-Ray can convert into
a 3D
image. You should probably save the image in grayscale to get the
best
results. You may have to fiddle with the source image depending
on your
initial MicroDem settings. Try to perform a "color merge"
under the
"Modify" and "Display Parameter" settings and then
pick "Grayscale Map" from
the "File" menu.
If your target programs can handle one of these formats, you are home
free.
If this does not solve your problem, let me know. I looked at the
MicroDem
.dem file with a binary editor and it looks pretty straight forward.
I can
probably write a quick conversion from this format to any target format
if
necessary.
Here is a link that has more MicroDem information from someone more expert
than I:
http://www.3dartist.com/3dao/r/allenbil/md/md4art.htm
Good luck.
Top
NDEM (South Africa) File Format
It was interesting reading about your Digital Terrain
Mapping and your
enviable knowledge on your deep understanding of the formats involved
as
highlighted on your Website.
I use the Metacust Ultra Version 2.05 3D environment. I am about to purchase
a particular terrain model from a supplier here in South Africa. Because
I
do not have a clue how these different formats will impact on my system
I
decided to ask the expert. Tell me, what are the specific format details
to
look out for in order for me to make a wise investment. I have a Silicon
Graphics SGI 02 Hardware running the IRIX version 6.5.7M Software and
I need
to create animated graphics for weather presentations. I need to purchase
a
base map of South Africa onto which I should be able to use available
terrain data and weather data to achieve an excellent graphical
presentation. Now, please excuse me for my lack of understanding of the
various applicable formats, I need to know, what must be my prescribed
format? In other words, which format and or applicable rules must I watch
out for during this purchase? All I know is that my manual stipulates
that
"The DEM files usable for Metacust Ultra, need to be either the USGS
1:250
000 scale DEM standard format or they should be in run-time performance
optimized proprietary format , nDem." Please explain this.
Awaiting your most urgent response.
Please accept my apologies for the delayed reply.
I was on
travel this week and had a small backlog of emails related to the site
to
work through when I returned. Also, please call me John. Also, please
do
not call me an expert, which I am not. "In the land of the
blind, the
one-eyed man is king" applies in this field.
From what I can tell from the information you related, you have few worries.
The application that you are considering appears to like two formats.
I
know the first one well so let me start with that.
USGS 1:250,000 DEM format (also known as USGS native format) is an old,
well-established format for digital topological data in the United States.
Although it is being superceded in the U.S. by the newer SDTS format,
many,
many GIS applications read this format. In the U.S. at least it
serves as
something of a neutral translation format, i.e. you may not speak my
language and I may not speak yours but we both speak USGS native DEM.
(This
is ironic considering that the USGS wanted the hugely complicated SDTS
format to serve this purpose.)
The format is ASCII meaning that it is human-readable. This makes
conversion from this format to another, if this should be necessary, a
lot
easier. It will also be easy to transfer data across platforms (like
from
Windows to IRIX) since you do not have to worry about issues like
machine-dependant byte order preferences, etc. With all this
going for it,
there must be a down side, right? One is file size. ASCII
USGS DEM files
are going to be a lot larger than binary files. This used to be
an issue
before the era of cheap computing but is pretty much a non-issue today.
The
other is that it is becoming extinct as eventually the USGS will convert
their data sets completely to SDTS. (However, it will live for a
long time
as a translation medium.)
If your application can read data in USGS native DEM format, you are off
to
a good start. It will not be able by any means to accommodate all
the GIS
data out there but you will be able translate fairly easily anything that
you cannot read directly.
The link to the file specification can be found at:
ftp://mapping.usgs.gov/pub/ti/DEM/demguide/dugdem.txt.
I have no experience with the second so the following is conjecture.
NDEM
appears to stand for National Digital Elevation Model. This is a
raster
file format used by the Chief Directorate to Surveys and Mapping (South
Africa) for their DEM data set. You apparently purchase this DEM
data
directly from CDSM. When you do, you get a CD with your map data
in NDEM
format.
While I was not able to locate the file specification on the CDSM website,
it is very likely that it is there or available from them. It is
described
as a raster file format, which essentially means that conversion of this
format, if ever necessary, should be straight forward. It is probably
a
binary format.
The only thing that bothers me a little about your letter is the qualifier
"run-time optimized proprietary format" that precedes the reference
to NDEM.
I hope that this is redundancy but you would do well to check with the
software supplier to make absolutely sure that the application can read
the
file format of your source data set directly and that this does not denote
some "improved" version of the format that your application
will choke on.
This is important because data for South Africa is going to be in NDEM
format for the most part. Conversions are pretty easy but you should
know
up front about any possible problems in reading your source data.
The link to CDSM where you can get information is:
http://w3sli.wcape.gov.za/Default.htm.
I hope this helps. Write with further questions if it does not.
Your
software supplier should be your most important source of information.
I
think you are in as good of shape as you can be in this field. Your
application sounds very interesting. I would like to learn of your
future
experiences so drop me a line sometime if you are able.
Top
Ordnance Survey Data
My brother has just sent me details of your site,
which I was very impressed with. I am a landscape architect working out
of Yorkshire in the UK.
I am interested in doing some modeling with terrain
data available in the UK from Ordnance Survey. The data comes in 2 forms
- DXF and NTF. Which (if at all) of these two forms can be handled by
the software you use for modeling?
I am particularly interested in being able to use
the data to demonstrate change in the landscape - either land raising
or land extraction
You win the award for the most interesting question of
the week. I am not very familiar with the Ordnance Survey database
other than knowing of its fabled history.
I looked at the website, however and found a helpful
document at
http://www.ordsvy.gov.uk/downloads/height/profile/profil_w.pdf.
It looks like the NTF format is analogous to the SDTS format used by the
USGS here in the states. It has a versatile, free-format file
structure that is correspondingly complex. The DXF format is well-known.
The version used by the OS is the 3D variety, meaning it carries its DEM
data in xyz format. This is of very wasteful since the xy data constitutes
a fixed-spacing grid. It is of course only necessary to specify
the number of rows, number of columns and the spacing (along with at least
two geocoordinates). This is why data in the DXF format is so much larger
than the NTF version.
None of the GIS applications that I use, including MicroDEM
and Global Mapper (which are my most versatile) will accept DEMs in either
NTF or DXF format. I did some research and found two applications
that should get you started. Leveller at
http://www.daylongraphics.com/products/leveller/ and Maprender
at http://www.maprender3d.com/ claim
to accept NTF format so I would check these out and go with NTF
format if they worked. Leveller only costs around US$50 and there
is a free demo that you can try out to make sure it works. Maprender
looks like a nice program but it is more expensive. Hopefully one
or both will be able to do some of what you have in mind. (There
may be others more widely known in the UK. Check with the OS and see what
they recommend.)
I hope that this will get you started. Write again
if it does not and we will think of something else.
Top
Ortho Imagery
I am currently working on a research project using
a DEM (Winfield, CO), and
aerial photos of the area, that have been scanned. I would like to
georeference the photos in order to drape them over the DEM. Do you have
any
suggestions on how to georeference the photos? I have not been able to
find
any existing georeferenced photos for the area.
Your website is very helpful! Thanks for the insight and direction.
You get the prize for asking the most difficult question
of the week. Or
maybe you win for picking the most out of the way location. Winfield?
That's a ghost town! I hope you get extra credit for your selection.
Only kidding. The location that you are working with certainly is
an
interesting one from a topographical standpoint. Close to Buena
Vista home
of the annual FibArk whitewater kayak race, the famous Numbers rapids,
Mt.
Elbert, etc.
I checked the USGS DOQ directory to see if DOQs were available for your
area
of interest. Unfortunately, there is coverage for almost all of
the state
except for a big hole over the Collegiate range - i.e. your area of
interest. You can take my word for it or check the coverage map
yourself at
http://mcmcweb.er.usgs.gov/status/rmmc/co/co_doq.html.
It would have been
nice if they were available because they are orthorectified, georeferenced,
and cheap.
I next checked Mapquest. They offer free aerial photos of some areas of
the
country, but not for Chaffe county, CO.
I then did a web search for commercial sources. There is lots available
for
the front range (probably recycled USGS DOQs), but not for where you are
looking.
I then tried satellite imagery. NIMA offers 10m DOIs but coverage
stops
just south of I70, apparently outside your area of interest. You
can check
this for yourself at the NIMA geospatial engine at www.nima.mil
(same link
as on my website.)
Then I looked for ASTER L1A satellite imagery from the EOS Data Gateway
http://edcimswww.cr.usgs.gov/pub/imswelcome/.
They do not have any images
of your area of interest in their archive. You can order a satellite
overflight and get an image of Winfield, but it will take a couple of
weeks.
Let me know if you are interested and I can tell you how to do it.
(However, ASTER has visible spectrum sensors so there is no predicting
the
cloud coverage you will get. ASTER may not even produce a useful
image if
taken at this time of year.)
Then, I looked at the Space Imaging Quicklook site at
http://www.spaceimaging.com/gallery/quicklook/collection3_geo.asp?vFIPS_CODE
=08015 and found some images very close to your target. There are
13
available. I did not check them all but there does not appear to
be an
exact match. However, depending on the extent of coverage of the
images you
have, there may be some overlap and you may be able to infer some
coordinates by looking at the Space Imaging photos. They are all
georeferenced and orthorectified.
This brings us back to using the images that you have and trying to
georeference them using ground features. This is generally not as
easy as
it sounds because one mountain can appear quite similar to another.
I have used Paintshop Pro (Photoshop would probably do just as well.)
to do
ground referenced overlays in the past. In this method you
convert your
DEM to .tga format and open it in Paintshop. Then you paste
your overlay
as a new layer over the DEM tga. You decrease the opacity of the
overlay so
that you can kind of see through it. Then by sliding the overlay around,
referencing ground features and resizing you can usually get an
OK match.
(It is important to (re)load the original overlay each time you resize
it
because the interpolation algorithms will lose information each time if
you
do not.) When you are satisfied with the match you crop both
the DEM and
the overlay and save them separately. You then use an application
like
POV-Ray that can map the now perfectly matched raster onto the tga-DEM
height field to complete the job.
Alternatively, you might try using 3DEM and guess the corner coordinates.
If you are wrong, try keep trying until you get it right. This should
work
OK if you know the corner coordinates approximately. (There is a
nudge
command in 3DEM that is helpful in doing this job.) This should
be possible
by carefully comparing your aerial photos to the USGS topo map and guessing
at some approximate corner coordinates.
Anyway, I have not been too helpful but the short answer is that you can
definitely use your images if you work hard enough at it. There
is no one
correct way in this case but you should succeed with some patience.
Let me
know if anything is interesting to you and I can send more details.
Top
Overlay Problem
I went through your web page.It is very informative.
I have a problem in
registration of satellite image with DEM images. The problem is my DEM
image
size is small and the satellite image size is big and my reference
image is the
DEM image. So if I do registration of DEM with satellite image then
I am
getting a satellite image with reduced spatial resolution. How can
I register
the DEM with the satellite image without reducing the spatial resolution
of satellite
image?
I am guessing that the problem is that your overlay software works by
mapping one pixel
of the overlay image to exactly one elevation value (and effectively one
pixel) of your DEM. When there are more pixels in the overlay than
there
are DEM pixels to map to, your software solves the problem by thinning
the
overlay image. This in turn results in a degraded overlay image.
The only solution is to increase the number of pixels in the DEM.
This
means either using a higher resolution DEM, or, if such a thing is not
available then artificially increasing the number of pixels by some type
of
resizing algorithm.
I looked around for some combination of standard applications that would
do
the job and did not find anything. I tried saving a DEM as a BMP,
resizing
the DEM image upward in Paintshop Pro, then converting it to a tiff, then
importing it into GlobalMapper and then exporting it to an ASCII DEM.
This
did not work because GlobalMapper refused to take the tiff because it
was
not a Geotiff with embedded georeferencing information.
I think the only solution would be to write a program to resize the DEM.
This is a big enough job so that it would have to be a pretty important
overlay to justify the effort.
Sorry I could not think of anything better. If I do, I will let
you know.
Top
R2V
We are in a process of evaluating the stability of
rock slopes is a part
of Saudi Arabia. We are interested in constructing a DEM map, a 3D DEM
and a slope category map from a 1:50,000 contour map. What software can
we use? Will ArcView or MapInfo help?
Thanks.
I looked into this problem as a way of getting DEM data
for hard-to-obtain
international locations. I am in the process of exploring this area
myself
at the moment, and here is what I know.
In order to convert a contour map into DEM data you need to first convert
the raster data (bitmapped contours) into vector data (tagged vector
contours.) In order to do this, you must first scan your topo map
to
convert it to digital data, and then use a raster to vector conversion
program to create the vector data. There are several software packages
available to automate the task. I have used R2V from Able software.
I do
not know what ArcView or MapInfo can do, but I would strongly recommend
that
you look at software that is designed to do this specific task, rather
than
a general purpose package.
After you do the conversion, you have to clean the image, repair the contour
lines and then tag them with the appropriate elevation values manually
or in
a semi-automated fashion.
When this is done your program will automatically generate a DEM file
by
interpolating to grid points.
Let me tell you the bad news. In order to get good results you must
do a
formidable job of cleaning the image. A typical contour map has
all sorts
of information superimposed on it that will corrupt your results.
You must
remove all conflicting lines that might be interpreted as contour lines,
such as roads, political boundaries, text, shaded regions,etc.
The contour lines, no matter how good your source image is will be broken
in
many places and will bleed together in many others. These must all
be
repaired (manually) or otherwise rationalized before you can use the data.
I am in the process of writing some preprocessing software that will help
clean and repair images automatically but am nowhere near done.
I have concluded that this task is virtually impossible and it is much
more
practical to locate or produce DEM data by other means. I will post
an
article in Spatial News in the near future to explain how this can be
done.
Unfortunately, It will be a couple of weeks before I have all the details
worked out so that the process I have in mind is feasible. But check
SN or
my website for more on this problem in the near future.
Thanks for your prompt response on my email. I reached the stage you
stated by screen digitizing the raster map as I work on large scale
maps. I have the contours in a tagged vector form. Will the R2V create
the DEM then?
Thanks again.
Well, R2V will generate the DEM. I just have the
free demo version of R2V
so my test bitmap was limited to 250X250 pixels. When I ran the
program, it
generated what looked like a proper USGS native DEM file. Since
I never got
as far as you are at the moment in terms of a high-fidelity source file,
I
had to go back and work on that problem. However, R2V is set up
to take a
raster image, convert it to vector form automatically, and then it has
utilities for tagging the contours. (these utilities are great if
your
vectorized contours are continuous and do not bleed.) If you already
have
this work done by other means, I am not sure that you can easily import
the
vector data into R2V.
Your best way to answer the question is to download the R2V demo, subset
your input file, run it through the program and see what you get.
The link
is http://www.ablesw.com/.
I emailed some questions to Dr. Wu already so
if you need some pointers in using the program I will tell you what I
know.
I would appreciate any tips you could give me in the area of preparing
the
input file. I am trying to write a section for my web page to help
amateur
cartographers prepare DEMS from topos. The literature is sparse
in this
area. I recently bought a book called "Machine Interpretation
of Line
Drawing Images" by Ablameyko but it does not seem to have too much
more than
I know already.
I got the R2V and could not open the digitized map.
So I started from
scratch, scanned the image then use the R2V to create a vector file then
a 3D DEM. I found it easier to clean the unwanted lines in the scanned
image using the Photoshop software. Further cleaning was also necessary
using the RTV before the creation of the DEM. My problem now is that I
have a 3D DEM, a fancy one that I can animate, but could not have a 2D
DEM file. Any suggestions?
I am very well acquainted (frustrated) with certain applications
reading
legal USGS DEM files and others choking on them. If you want a few
tips,
here they are:
Fortunately USGS DEMs are ASCII files so they are easy to read.
Try
inspecting the file manually and comparing it to one that does work for
you
and try to find small but important differences. I find it better
to use a
hex editor like HexEdit (free download) rather than a word processor because
it shows you the exact position of each character in the file. It
is best
to use small example files for this test.
Sometimes the application will reject the file because of missing or
inconsistent data in the type A record. This is the front part of
the file
that has all the metadata. This is especially a problem when you
are
building a file from scratch as opposed to converting from another format
where you may not have all the metadata. Try to insert all the metadata
that R2V allows you to put in. You may want to try inserting other
data
manually (see HexEdit above).
The DEM file specification is available at
http://rockyweb.cr.usgs.gov/nmpstds/acrodocs/dem/2DEM0198.PDF.
You can use
this to as an aid in examining your file.
How about opening the DEM in another application and then exporting from
there to your target. MicroDEM seems to be pretty forgiving in terms
of
reading DEMs, although its export capabilities are limited. (MicroDEM
will
open the R2V file, as will DEMworks). You could then save it as
a MicroDEM
DEM and then use my converter MDEM2DEM to change it back to a DEM and
see if
your target will accept that version.
I told you this would be a pain. But if you really enjoy abuse,
try
extracting DEM data from satellite imagery.
Regarding your last email: Is the 'join' operation automated in R2V or
do
you have to manually draw all the little lines between each break?
In order to tag the contours, you need first to join the segments of
each contour in one line. You need to turn the line editor ON, select
all the line segments forming on contour and join them. Click View,
Overlays, Line IDs and start assign value for each contour. Sometimes
you still see few spots on each contour having a zero value. You need
to
clean these. After that from the file 3D data, you can create the 3D
DEM. I suggest you check the R2V tutorial file.
I tried to export the 3D DEM to IDRISI, but it wouldn't accept it. I'm
still trying with other software.
Before I start with the editing, I'll try other software.
I'll also ask
one of our GIS Department staff to help in the editing process. The join
operation is an automated one for the selected segments and does draw
line between each break. But you should be careful because it sometimes
creates un-needed lines. Each time you select a segment to join, check
and make sure that no other un-needed lines are added to your contour.
After you finish joining and set value for the contour line, you will
still see some points on the contour line that have zero values. You
need to delete them.
Hi I just read your web page and found it interesting. I have some DTED
level 1 files that have been created by digitizing bathymetry contour
lines
but without inputting any depth data in between the depth contour lines.
So
what I am left with are DTED format files with depth steps let say -20,
-30,
-50, -100 and -200 meters like instead real gridded depth data. Many of
them
only have few depth contour lines (let say 4 or 5) over the whole 1x1
degree
cell. I want to generate real gridded DTED level 1 without loosing the
location of the depth contour lines and shorelines. Meaning replacing
the
flat depth steps between contour lines by interpolated depths. What would
you think could be the best tool or approach to do so. I look at some
stuff
and I have installed MicroDem 5.1 beta version but didn't the good one.
It sounds like R2V from Able Software might help.
It will allow you to
import your data in raster or bitmap format, add information for example
tagging contour lines, and then output a DEM file.
The program will show up in your search engine. I am in the process
of
writing a web page section on how this can be done.
I found a lot of interesting information in your
site.
I am just going to make some 3D models of several
bottom areas of the Baltic sea (as a part of my research project). I am
ecologist so GIS is more like a hobby, but it is very powerful thing to
visualize my ecological data. Unfortunately getting DEM or whatever data
on relief of the Baltic sea is a BIG problem.
So I decided to make my own. Just from simple paper
maps. But, I don't know what software I should use (keeping in mind the
fact it should be not too expensive). Could you please give me an advice.
Thank you in advance.
I have to do a DEM map in Idrisi. I have paper
maps with level curves.
May you tell me what is the best way to digitalize
the curve levels and transform them to DEM?
The best way that I know of is to use the tools in ArcInfo
to convert the raster contours to tagged vectors. This program is
very expensive (about USD 3000) and not easy to use but is apparently
the standard for this type of work.
As I mentioned on my web page, RV2 is less expensive
(USD 900) but requires more work to make the conversion.
I am not an expert in this area and the article on my
web page www.terrainmap.com tells
most of what I know about this subject.
Top
Seafloor Bathymetry
Please could you help direct me to topo`s of the ocean
area around Long Island.
I do not know much about sea floor topography. The NOAA is probably
the
best source of data. Here is a link for you to explore:
http://www.ngdc.noaa.gov/mgg/bathymetry/relief.html
GlobalMapper will read ETOPO2 data files which are available from the
NOAA.
These might get you started. You have to pay for both the ETOPO
data and
GlobalMapper.
MicroDEM will read NOS EEZ Bathymetry. A source for this data and
more
coastal bathymetry can be found at:
http://www.ngdc.noaa.gov/mgg/bathymetry/multibeam/multibeam_products.html
It seems like none of this data is free but it is not too expensive.
Top
SDTS
Hi, my name is Jody. I am trying to use some
DEM data I have downloaded off the web and am having a few issues.
Your web site has been most hopeful with information regarding DEM's as
well as file conversions, but I am still having issues with files I am
getting.
The highlighted file is the original, I have
unzipped the files and as you can see they all have the file ext. .DDF?
Do you know of a utility that will allow me to convert these to .DEM?
Any input on this issue will be appreciated. Thanks
for your help.
Download the utility SDTSDEM2 listed under the "News
and Information" section of my web page at
www.terrainmap.com Run it and type in the path and then the first
four numbers in the file name of the profile for example c:\DEM\9769
(unless the utility is in the same directory, in which case just type
in 9769) . The utility will query you for the output file name,
to which you should respond <path><filename>.dem, i.e. c:\dem\outfile.dem.
A file in USGS native format file will be produced.
Several of the .ddf files in the SDTS profile that
you downloaded contain information required by SDTSDEM2 to (re)construct
the USGS native DEM. After converting you can read the .dem
file with a text editor because it will be in ASCII format.
Can you tell me the difference between your two versions
of the sdts to
dem converter and the SDTS2DEM that comes with source code at version
.016 that was updated on 9/5/01 by Gregg Townsend (based on the original
sol katz code) and is available at www.gisdatadepot.com/dem/?
As far as the dem data goes, they are equivalent in that they both convert
the gridded elevation data in the xxxxcel0 file to a series of USTS DEM
type B records.
Differences would be in the metadata contained in the type A record.
As you
may know, there are 18 SDTS files in the profile, each of which contains
information, some more relevant than others. My program dips into the
.cel0
.ldef, lspdm, .iref, and ddsh files to extract such data and
insert it
into the proper fields in the type A record of the output .dem file.
Gregg's revision of the original Sol Katz program may be better than mine,
in that he may extract more of this metadata from the SDTS profiles.
The
problem is that even though his source code is available, it uses an obscure
(to me) standard library for handling the SDTS files. I was unable
to
figure out what was going on in that program so I decided to write
my own. I wrote my program is an extension of an SDTS reader that
I needed
for other purposes anyway.
Like I say in the disclaimer on my page, if your need is for
mission-critical output you had better write your own conversion program
so
that you know exactly what is going on. The weakness of my program
is that
it does not extract all the available metadata from the SDTS profile.
In
some cases I insert default information in type A fields that I consider
of
secondary importance and that I know to be generally true for the USGS
data.
However, there is a possibility that one of these more obscure fields
could
be in error of some input file throws it a curve. These errors are
obviously not important if all you are doing is creating images like I
do.
Top
STL Format
I'm interested in combining GIS with 3D prototyping.
Are you aware of any
converters that are freeware that will produce STL (StereoLithography)
files
from DEM files?
Here's a high end version of what I'd like to do: http://www.stm-usa.com/
you may have seen these guys already. Pretty cool synthesis of a
couple of
technologies. I really just want to mill hardwood topo models.
Thanks for any pointers you might offer. And thanks for the freeware
converters. I'm a big POV fan.
I had communicated with Mike G. about STL format a while ago and
recently received a follow up email from him. I had been planning
to write
an STL converter but I keep getting interested in other things and never
get
around to starting it. You might try contacting Mike and see where
he ended
up. Here is a copy of the email that he sent me.
If you don't have any luck let me know and maybe I can write a converter
quickly over the semester break.
His email:
Sorry about the delay in responding to your letter.
I found a shareware (fairly cheap) converter called "AccuTrans3D"
which can convert to/from many different 3D file formats, including
from DEM to STL. I have been using that successfully in making
some 3D topographic maps, starting out with Mt. Saint Helens.
Unfortunately, now it seems that the AccuTrans3D web site
is unavailable. In any case, I'm going to get Rhino3D soon.
Now my big problem is dust collection! Yech.
I just released a product which you might be interested in.
It allows extremely large images to be viewed on the web.
It's now available in a beta release at http://www.xippix.com/cooker
Top
Terragen
I am doing a final year project on the Sao Miguel
island of Azores.
I have been using arc/view and recently came across your website and thought
it
would be a good idea to do an animation of the area that I am studying.
however I am not to sure which extensions I am currently using. I did
not
digitize the contour lines, I had to scan in the image and proceed to
'trace' the contours and other features using ArcFiew.
therefore I would like to use the package Terragen, to create an animation
flyby, or just 3d visualizations of the area.
The problem being I don't know what file extensions are compatible with
Terragen, and what extension my files are currently in.
I can't find any data about this island on the Internet to carry out my
idea.
I hope that you can enlighten me on this topic and help me with a suitable
solution.
Your digitized image is probably in some type of vector format that is
incompatible with Terragen. In fact, Terragen requires its own unique
file
format for input data. This is called a TER file.
Your best method is probably to export your file from ArcView in some
type
of DEM format. Usually USGS ASCII DEM is offered as an option (although
I
do not know for sure whether Arc View can do this.)
Then import the file into MicroDem (see my web page for the link if you
do
not have it). Save the file as a MicroDem DEM.
Then use my utility MDEM2TER to convert the MicroDem DEM to TER format
using
MDEM2TER. I realize that this is cumbersome but it is the only way
that I
know.
Firstly, I apologies for the nature of this question
since it no doubt has
been answered before, but being a 3D newbie I'm still grappling with some
of
the techniques.
Like a lot of people, I have started my 3D terrain experience with Terragen.
I've played around and am beginning to produce some nice images
and
animations. However, what I really want to do is build a model of
Grenada,
West Indies.
Now I am aware that I need to download the relevant data. From your
website
I had a look at the NIMA website and seemed to find the co-ordinates I
would
need to create my model. The problem is, which format should I download
and
then how do I get the file into Terragen?
From viewing the pages on your website, I understood that the SDTSTER4
utility you created would allow for the export of .ter files. What
I'm
having problems with is using the information I get from NIMA so that
I can
use the utility to convert it to a .ter file. Can this be done?
As far as I understand it I need a USGS 7.5' DEM file for the utility
to
work, but is that what the NIMA file is?
If you could provide me with any pointers to put me back on the right
track
I'd really appreciate it.
Terragen only accepts its own TER file format.
So you have to translate whatever your input is into that. NIMA
data must be first imported into MicroDem, merged, saved in MicroDem native
DEM format, and then converted to TER using my utility MDEM2TER.
Look at the Sinai Challenge section of my web page for further instructions.
You can go directly from SDTS to TER using SDTS2TER4
for US DEM data as you noted. Merging is not usually required in
this case.
Thanks so much for those pointers. I had completely
missed the Sinai
Challenge' section for some reason!!
Anyway, it all worked well, but I was really disappointed with the quality
of the data in that I could hardly even get the land forms to appear in
the
render. Is this just simply because the detail just isn't there
from the
source? As a result I went to have a look at GTOPO30 and downloaded
a 10MB
that might be better.
Is this a fairly typical problem with international DEMs? Or do
you know of
any other sites that I might try? The problem is, the island of
Grenada is
only 17 by 23 KM!!!
Top
TopoUSA
Hi!
I have been looking at your page, it's wonderful!
I purchased TopoUSA regional per your suggestion but
it is entirely too "rough" for my liking. I would like
to produce much more refined maps. My question is: I
can save alot of time if I can use the data cd in
TopoUSA since it has done the downloads and has the
data on it. Is there any way I could access the
information (it seems very cryptic)? If not I am going
to return the product since it is useless to me. I am
evaluating a product called Surfer 7 made by Golden
Software to use for my map creations. Do you know
anything about it? Thanks for your time!
I am sorry that I steered you wrong on TopoUSA.
I agree about the
limitations and do not use the program at all for my work, but I definitely
would use the maps thus created for hiking and camping. Perhaps
you could
sell/give it a hiker. If it is any consolation I will insert your
feedback
into the page so that others will be forewarned. In my defense do
get a lot
of emails from people who are not exactly power users and really thought
it
offered a good solution for them.
I had the same idea about the data when I bought it. Unfortunately,
it is
not stored in any type of conventional file structure that I have been
able
to get at. When I tried to open a file using my hex editor it just
locked
up my computer. Even if I could open a file, the format is unknown
and the
data is almost certainly compressed using an unknown compression scheme.
I
did a web search to see if anyone had decrypted it yet but came up empty.
I got an email once asking for a converter from DEM to Golden Soft Surfer
Grid format, but have no experience with using it at all. My experience
is
pretty much limited to the packages discussed on the web page.
Top
VistaPRO
I'm trying to import the SDTS files from www.gisdatadepot.com
into Vistapro. The program has a utility to convert the "old style"
ASCII USGS DEMs. When using your program mdem2dem the result looks different
than the ASCII DEMs I know from a while ago, and I can't convert it any
longer (3DEM doesn't like them either). The biggest difference I see so
far, is that the converted file has lots of columns with -9999 in it.
I hope you have an idea what I'm doing wrong....
Thank you for sending the files. Unfortunately,
I was not able to solve your problem right away. I unzipped your
file named "Bonita.dem" and looked at it using the text editor.
It looked normal. I tried to open the file using MicroDem, DEMWorks
and 3DEM. The file opened correctly in all three, as normal.
I have attached a rather large BMP screen shot of the 3DEM image in
case you do not believe me. I do not have a copy of VistaPro so
I cannot run tests with this program.
I am using the 3DEM 8.4.1 demo version that I downloaded
from http://www.davecentral.com/5457.html very
recently. Are we using the same version of 3DEM? I think I
can get 3DEM70 for free. I will download this version and try the
import to see if this is the problem. I notice that the ny-e
file that you sent me is a 1:250K DEM. DLGV-32 from the USGS has
the same problem of not reading converted 1:24,000 SDTS DEMS (both from
MDEM2DEM and the Sol Katz program SDTS2DEM) but reading the 1:250,000
unconverted DEMS successfully. DEM readers are finicky and sometimes
they will not read the file if only one flag in the ASCII DEM is set to
a wrong integer. Since MicroDem loses a lot of the Metadata I had
to set a lot of the internal flags for a 1:24000 UTM file. Perhaps
this is the problem.
For the moment, I am stuck until I can run some additional
tests. Try downloading the 3DEM demo from the link above and loading
Bonita.dem. I will write back to you when I have more information.
later...
I downloaded 3DEM70 from
http://www.simtel.net/pub/pd/6700.shtml. This program also rendered
your file Bonita.dem successfully. I just entered <File><Load
Terrain><Digital Model>, made sure that the file type was set
to <USGS DEM or GTOPO30 DEM> and then clicked on the file.
It loaded perfectly. The rendered image is attached.
I also downloaded the demo version of VistaPro.
This program is looking for a VistaPro native DEM (which is different
from a USGS DEM) or a V4s file. The demo version would not load
any of the USGS DEMS, including ny-e.
later, Klaus wrote...
First of all, you're right, I can read the dem file your
program creates with 3dem. Sorry for the mistake.
For my (unsuccessful) test reading the SDTS dems I used
3dem Version 8.4.4. The program starts reading the sdts and asks for a
dem file to save the converted data to. After the conversion 3dem tries
to read in the dem file and gives a "Incompatible File Format"
error message. When looking into the file, it has some values like 4.512D+006
in it. I thought, that this might be the problem I have. I downloaded
version 10.0.5 like you suggested and this one could read the sdts,
but doesn't write an ASCII dem....
So there still seems to be a difference between the
"old time" ASCII dems and the new ones. Those created by your
program can be used with 3dem, but not with Vistapro. The way of reading
a sdts with 3dem to get an ASCII dem, Vistapro likes, results in this
"Incompatible File Format" message.
So I still need to find out, what has changed. I attached
the conversion program, which reads the ASCII dems and writes the Vistapro
dems to this mail. Maybe you can figure out, what the tool doesn't like.
Hello, could you say me how can I convert a VISTAPRO
4.0 DEM
to a USGS DEM ?.
Thank you.
Arthur (Spain)
Arturo:
Muchas gracias para sus dos cartas.
Busque en el Internet por el dato tecnico del archivo
de datos de VistaPro.
Noticias buenas: descubri el dato tecnico basico
facilmente y rapidamente. Noticias malas: el dato
tecnico no inluyo la informacion del algoritmo para la
compression del archivo de datos, se llama "run length
compression" en ingles. Busque muchas horas en el
Internet por este informacion pero no tuve exito. Es
necessario a obtener este informacion para hacer la
programma a leer el archivo. Se aparece que el
algoritmo es un secreto commercial.
Escribi a el autor de la programa y lo pregunte por el
algoritmo que desiero. Espero que el lo respondera
pronto a mi peticion.
Muchas gracias a me permite a practica me espanol con
usted. Usted sabe, hablo solamente unas pocas
palabras de espanol.
John
y Arturo escribe...
Hola de nuevo, antes de nada felicitarle por su buen
nivel
de Castellano (Español).
Cuando le pregunte por un conversor de VistaPro a USGS DEM,
era porque necesitaba introducir un DEM que cree de una
parte de España en el programa World Construction Set, ya he
solucionado ese problema, ya que he encontrado en la WEB
http://amber.rc.arizona.edu/lw/demscape.html
un conversor de
VistaPro DEM a STM, y World Construction Set acepta este
formato sin problemas.
Ahora mi problema es saber como hacer funcionar el World
Construction Set ya que no tengo ningun manual.
Le envio un archivo que me dejo un amigo de las Islas
Canarias (España)en formato XYZ, tengo problemas con el ya
que World Construction Set no me lo acepta, podria darme
alguna idea de como poder trabajar con el.
Gracias.
Arturo.
I just started digesting you Digital Terrain Modeling and Mapping using
DEM,SDTS, DRG, DLG and DTED Data and think its great.
The first terrain modeler I used was Vistapro and I still like it. The
problem with it is finding data. I would like to take advantage of the
current SDTS data and can do it indirectly by producing a shaded contour
map which vistapro can read. This is really pain and I would like to go
directly from SDTS to the VistaPro DEM file or the old USGS DEM file for
which there is a Vistapro DEM converter. Can any of the SDTS converters
do
this
Top
XYZ To ArcView
Are you aware of any xyz to dem converters. I have
three columns, east, north and elevation and would like to generate contours.
Thanks in advance,
I am not aware of any. The problem in general is that standard DEM
file formats require a bunch of metadata. A raw xyz grid file will
not have this data and so a converter will not know what to put in the DEM
file fields. This question came up from another person just last
week. As a test I tried converting an xyz grid into a tiff format
and then asked GlobalMapper to read it in as if it was a Geotiff.
The read went OK but I was not permitted to write the file because there
was no metadata keys.
Top
|