MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Question about SDTS format

  • To: mathgroup at smc.vnet.net
  • Subject: [mg102928] Re: [mg102882] Question about SDTS format
  • From: "Elton Kurt TeKolste" <tekolste at fastmail.us>
  • Date: Wed, 2 Sep 2009 04:04:48 -0400 (EDT)
  • References: <200908311033.GAA06953@smc.vnet.net>

I fiddled with the data provided by the government in SDTS format and
was able to extract a list of list of coordinates that produces a nice
map of the states of the US using Mathematica.ListLinePlot[].  Contact
me if you want a copy of either all or the continental US.  

There is a lot of other information available from the government
(climate, epidemiology, politics, ...) if you know how to use it.

drekt


On Wed, 02 Sep 2009 00:10 +0200, "Patrick Scheibe"
<pscheibe at trm.uni-leipzig.de> wrote:
> Hi,
> 
> 
> > Are there SDTS tools available for Mathematica that you know of?
> 
> Sorry, I don't know. I just checked your approach and recognized you
> unpack the files but I never used the stuff by myself.
> 
> Cheers
> Patrick
> 
> 
> > Kurt
> > 
> > 
> > On Tue, 01 Sep 2009 14:03 +0200, "Patrick Scheibe"
> > <pscheibe at trm.uni-leipzig.de> wrote:
> > > Hi,
> > > 
> > > don't unpack the tar.gz files. Just do something like
> > > 
> > > Import["blub.tar.gz", {"SDTS", "Graphics"}]
> > > 
> > > Cheers
> > > Patrick
> > > 
> > > On Tue, 2009-09-01 at 03:52 -0400, Elton Kurt TeKolste wrote:
> > > > The files on the US Government website nationalatlas.gov
> > > > (http://nationalatlas.gov/mld/statesp.html) claims to have files in SDTS
> > > > format.
> > > > 
> > > > In particular, I have downloaded and decompressed the files for the
> > > > state data, but Mathematica cannot import them.
> > > > 
> > > > For example, after downloading:
> > > > 
> > > > In[519]:= Import["Z:\BIGSIGMAA\Data\STAPAPPC.DDF"]
> > > > 
> > > > During evaluation of In[519]:= Import::fmterr: Cannot import data as
> > > > SDTSDEM format. >>
> > > > 
> > > > Out[519]= $Failed
> > > > 
> > > > It will read the "Element" field:
> > > > 
> > > > In[533]:= elements = 
> > > >  Import["Z:\BIGSIGMAA\Data\STAPAPPC.DDF", "Elements"]
> > > > 
> > > > Out[533]= {"Centering", "CentralScaleFactor", "CoordinateSystem", \
> > > > "CoordinateSystemInformation", "Data", "DataFormat", "Datum", \
> > > > "ElevationRange", "Graphics", "GridOrigin", "InverseFlattening", \
> > > > "LinearUnits", "Projection", "ProjectionName", "RasterSize", \
> > > > "ReferenceModel", "SemimajorAxis", "SemiminorAxis", "SpatialRange", \
> > > > "SpatialResolution", "StandardParallels"}
> > > > 
> > > > But not any of the elements
> > > > 
> > > > In[534]:= map = 
> > > >  Import["Z:\BIGSIGMAA\Data\STAPAPPC.DDF", #] & /@ elements
> > > > 
> > > > During evaluation of In[534]:= Import::fmterr: Cannot import data as
> > > > SDTSDEM format. >>
> > > > 
> > > > During evaluation of In[534]:= Import::fmterr: Cannot import data as
> > > > SDTSDEM format. >>
> > > > 
> > > > During evaluation of In[534]:= Import::fmterr: Cannot import data as
> > > > SDTSDEM format. >>
> > > > 
> > > > During evaluation of In[534]:= General::stop: Further output of
> > > > Import::fmterr will be suppressed during this calculation. >>
> > > > 
> > > > Out[534]= {$Failed, $Failed, $Failed, $Failed, $Failed, $Failed, \
> > > > $Failed, $Failed, $Failed, $Failed, $Failed, $Failed, $Failed, \
> > > > $Failed, $Failed, $Failed, $Failed, $Failed, $Failed, $Failed, \
> > > > $Failed}
> > > > 
> > > > I can think of three possiblities:
> > > > 1) I ammisusing Mathematica
> > > > 2) Mathematica does not actually import SDTS files as advertised
> > > > 3) The government did not publish this data in SDTS format despite the
> > > > data suffix.
> > > > 
> > > > Does anyone have any hints as to how to import SDTS data or information
> > > > about the formats used by the government?'
> > > > 
> > > > Regards,
> > > > Kurt Tekolste
> > > > 
> > > > 
> > > 
> > Regards,
> > Kurt Tekolste
> > 
> 
Regards,
Kurt Tekolste



  • Prev by Date: Re: Question about SDTS format
  • Next by Date: BoxWhiskerPlot - BoxLabels and FrameTicks won't work together
  • Previous by thread: Re: Question about SDTS format
  • Next by thread: new curve last night