Re: [mapserver-users] Using e00 files

2010-02-04 Thread Jon Miner


On 02/02/2010 11:32 AM, Jon Miner wrote:
> Now, I'm trying to layer on information from the GRANIT website:
>
> http://www.granit.unh.edu/
>
> ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00
>
> These files are "e00" files. Based on some other mailing list entries,  I 
> used the ogr2ogr command to convert the file into Shapefiles:
>
> ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00
>

I got it working!  The solution was to re-run the ogr2ogr command like this:

ogr2ogr -t_srs EPSG:26719 -f "ESRI Shapefile" temp rrnh.e00

My question is... if tools like ogrinfo can read all the projection data out of 
the file, why can't mapserver do this on the fly?  Or is it just the case where 
we don't want to have to do all the work every time for performance reasons?

Thanks for everyone's help!

- Jon



  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-03 Thread Peter N. Schweitzer

On 02/03/2010 09:29 AM, Jon Miner wrote:

I had to specify the raster layer as:

LAYER
   NAME "usgs"
   TYPE RASTER
   STATUS on
   DATA o43071a4.tif
   PROJECTION
 "init=epsg:26719"
   END
END

So, I think the answer to your question is, no they're not in the same projection.  What 
do I need to do to get mapserver to translate between the two? Seems to me that mapserver 
should be able to read the information out of the files (like ogrinfo and gdalinfo do) 
and do the "right thing". I admit this might be a total misunderstanding on my 
part... still learning about mapserver's capabilities.


Yes, epsg:26719 refers to UTM zone 19; the NH state plane is, I think,
somewhat different.

The ARC.prj file generated by ogr2ogr contains this (expanded for clarity):

PROJCS[
  "NAD83_New_Hampshire",
  GEOGCS[
"GCS_North_American_1983",
DATUM[
  "D_North_American_1983",
  SPHEROID["GRS_1980",6378137,298.257222101]
  ],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]
],
  PROJECTION["Transverse_Mercator"],
  PARAMETER["latitude_of_origin",42.5],
  PARAMETER["central_meridian",-71.67],
  PARAMETER["scale_factor",0.7],
  PARAMETER["false_easting",984250],
  PARAMETER["false_northing",0],
  UNIT["Foot_US",0.30480060960121924]
  ]

epsg:3437 would use

PROJECTION
  "proj=tmerc"
  "lat_0=42.5"
  "lon_0=-71.67"
  "k=0.7"
  "x_0=30.0"
  "y_0=0"
  "ellps=GRS80"
  "datum=NAD83"
  "units=us-ft"
  end

but because the false_easting in the prj file is not 30,
it looks to me like you should try this:

PROJECTION
  "proj=tmerc"
  "lat_0=42.5"
  "lon_0=-71.67"
  "k=0.7"
  "x_0=984250.0"
  "y_0=0"
  "ellps=GRS80"
  "datum=NAD83"
  "units=us-ft"
  end

I wouldn't claim to be expert in any of this, so I welcome a review
by someone more familiar with these projection parameters.

Peter
--
Peter N. Schweitzer (MS 954, U.S. Geological Survey, Reston, VA 20192)
(703) 648-6533  FAX: (703) 648-6252  email: pschweit...@usgs.gov

___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-03 Thread Jon Miner


- Original Message 
From: Peter N. Schweitzer 
To: mapserver-users@lists.osgeo.org
Sent: Tue, February 2, 2010 12:24:44 PM
Subject: Re: [mapserver-users] Using e00 files

On 02/02/2010 11:32 AM, Jon Miner wrote:
> Now, I'm trying to layer on information from the GRANIT website:
>
> http://www.granit.unh.edu/
>
> ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00
>
> These files are "e00" files. Based on some other mailing list entries,  I 
> used the ogr2ogr command to convert the file into Shapefiles:
>
> ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00
>
> then copied the files into my data directory, and added a new layer to the 
> map file:
>
> LAYER
>NAME "rail"
>TYPE LINE
>STATUS on
>DATA nhrail
>CLASS
>  NAME "Railroads"
>  SYMBOL 0
>  COLOR 168 0 184
>END
> END
>
> Something is wrong since the line data does not get rendered. Maybe there is 
> a better approach instead of converting to Shapefiles?

Jon,

The rr data are in NH state plane coordinates--are the other data you
are using in the same projection?


--

Thanks Peter,

I had to specify the raster layer as:

LAYER
  NAME "usgs"
  TYPE RASTER
  STATUS on
  DATA o43071a4.tif
  PROJECTION
"init=epsg:26719"
  END
END

So, I think the answer to your question is, no they're not in the same 
projection.  What do I need to do to get mapserver to translate between the 
two? Seems to me that mapserver should be able to read the information out of 
the files (like ogrinfo and gdalinfo do) and do the "right thing". I admit this 
might be a total misunderstanding on my part... still learning about 
mapserver's capabilities.



  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-02 Thread Jon Miner
I've added the output from `ogrinfo -summary` at the bottom off this thread...


-Original Message-
From: mapserver-users-boun...@lists.osgeo.org 
[mailto:mapserver-users-boun...@lists.osgeo.org] On Behalf Of Jon Miner
Sent: Tuesday, February 02, 2010 10:32 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] Using e00 files

Making progress... Thanks to everyone's help, I have a simple app displaying a 
USGS Topo sheet using GeoTIFF.

Now, I'm trying to layer on information from the GRANIT website:

http://www.granit.unh.edu/

ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00

These files are "e00" files. Based on some other mailing list entries,  I used 
the ogr2ogr command to convert the file into Shapefiles:

ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00

then copied the files into my data directory, and added a new layer to the map 
file:

LAYER
  NAME "rail"
  TYPE LINE
  STATUS on
  DATA nhrail
  CLASS
NAME "Railroads"
SYMBOL 0
COLOR 168 0 184
  END
END

Something is wrong since the line data does not get rendered. Maybe there is a 
better approach instead of converting to Shapefiles?
-
$ ogrinfo -summary nhrail.shp nhrail

INFO: Open of `nhrail.shp'
  using driver `ESRI Shapefile' successful.

Layer name: nhrail
Geometry: Line String
Feature Count: 2207
Extent: (745250.25, 70273.328000) - (1231517.80, 915630.00)
Layer SRS WKT:
PROJCS["NAD83 / New Hampshire",
GEOGCS["GCS_North_American_1983",
DATUM["North_American_Datum_1983",
SPHEROID["GRS_1980",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",42.5],
PARAMETER["central_meridian",-71.67],
PARAMETER["scale_factor",0.7],
PARAMETER["false_easting",984250.031003876],
PARAMETER["false_northing",0],
UNIT["Foot_US",0.3048006]]
UserId: Real (11.0)
FNODE_: Real (11.0)
TNODE_: Real (11.0)
LPOLY_: Real (11.0)
RPOLY_: Real (11.0)
LENGTH: Real (12.3)
RRNH#: Integer (5.0)
RRNH-ID: Integer (5.0)
GRANITID: Integer (8.0)
RRA: Integer (2.0)


  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-02 Thread Jon Miner
Hi David,

The layer name is explicitly called out in the URL; changing it to DEFAULT has 
no effect.

- Jon



- Original Message 
From: "Fawcett, David (MPCA)" 
To: Jon Miner ; "mapserver-users@lists.osgeo.org" 

Sent: Tue, February 2, 2010 12:07:29 PM
Subject: RE: [mapserver-users] Using e00 files

Jon, 

If you are not calling the layer name explicitly in the URL, try setting the 
layer status to DEFAULT.

David.

-Original Message-
From: mapserver-users-boun...@lists.osgeo.org 
[mailto:mapserver-users-boun...@lists.osgeo.org] On Behalf Of Jon Miner
Sent: Tuesday, February 02, 2010 10:32 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] Using e00 files

Making progress... Thanks to everyone's help, I have a simple app displaying a 
USGS Topo sheet using GeoTIFF.

Now, I'm trying to layer on information from the GRANIT website:

http://www.granit.unh.edu/

ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00

These files are "e00" files. Based on some other mailing list entries,  I used 
the ogr2ogr command to convert the file into Shapefiles:

ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00

then copied the files into my data directory, and added a new layer to the map 
file:

LAYER
  NAME "rail"
  TYPE LINE
  STATUS on
  DATA nhrail
  CLASS
NAME "Railroads"
SYMBOL 0
COLOR 168 0 184
  END
END

Something is wrong since the line data does not get rendered. Maybe there is a 
better approach instead of converting to Shapefiles?



  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users



  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-02 Thread Brent Fraser
You haven't defined the LAYER's PROJECTION (we really should enhance mapserver 
to read .prj files!)




Jon Miner wrote:

Making progress... Thanks to everyone's help, I have a simple app displaying a 
USGS Topo sheet using GeoTIFF.

Now, I'm trying to layer on information from the GRANIT website:

http://www.granit.unh.edu/

ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00

These files are "e00" files. Based on some other mailing list entries,  I used 
the ogr2ogr command to convert the file into Shapefiles:

ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00

then copied the files into my data directory, and added a new layer to the map 
file:

LAYER
  NAME "rail"
  TYPE LINE
  STATUS on
  DATA nhrail
  CLASS
NAME "Railroads"
SYMBOL 0
COLOR 168 0 184
  END
END

Something is wrong since the line data does not get rendered. Maybe there is a 
better approach instead of converting to Shapefiles?



  
___

mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


Re: [mapserver-users] Using e00 files

2010-02-02 Thread Peter N. Schweitzer

On 02/02/2010 11:32 AM, Jon Miner wrote:

Now, I'm trying to layer on information from the GRANIT website:

http://www.granit.unh.edu/

ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00

These files are "e00" files. Based on some other mailing list entries,  I used 
the ogr2ogr command to convert the file into Shapefiles:

ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00

then copied the files into my data directory, and added a new layer to the map 
file:

LAYER
   NAME "rail"
   TYPE LINE
   STATUS on
   DATA nhrail
   CLASS
 NAME "Railroads"
 SYMBOL 0
 COLOR 168 0 184
   END
END

Something is wrong since the line data does not get rendered. Maybe there is a 
better approach instead of converting to Shapefiles?


Jon,

The rr data are in NH state plane coordinates--are the other data you
are using in the same projection?

Peter
--
Peter N. Schweitzer (MS 954, U.S. Geological Survey, Reston, VA 20192)
(703) 648-6533  FAX: (703) 648-6252  email: pschweit...@usgs.gov

___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


RE: [mapserver-users] Using e00 files

2010-02-02 Thread Fawcett, David (MPCA)
Jon, 

If you are not calling the layer name explicitly in the URL, try setting the 
layer status to DEFAULT.

David.

-Original Message-
From: mapserver-users-boun...@lists.osgeo.org 
[mailto:mapserver-users-boun...@lists.osgeo.org] On Behalf Of Jon Miner
Sent: Tuesday, February 02, 2010 10:32 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] Using e00 files

Making progress... Thanks to everyone's help, I have a simple app displaying a 
USGS Topo sheet using GeoTIFF.

Now, I'm trying to layer on information from the GRANIT website:

http://www.granit.unh.edu/

ftp://ftp.granit.sr.unh.edu/pub/GRANIT_Data/Vector_Data/Transportation_Networks/d-railroads/rrnh.e00

These files are "e00" files. Based on some other mailing list entries,  I used 
the ogr2ogr command to convert the file into Shapefiles:

ogr2ogr -f "ESRI Shapefile" temp ../data/rrnh.e00

then copied the files into my data directory, and added a new layer to the map 
file:

LAYER
  NAME "rail"
  TYPE LINE
  STATUS on
  DATA nhrail
  CLASS
NAME "Railroads"
SYMBOL 0
COLOR 168 0 184
  END
END

Something is wrong since the line data does not get rendered. Maybe there is a 
better approach instead of converting to Shapefiles?



  
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users