Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-02-01 Thread Rahkonen Jukka
Hi,

EPSG:4326 bounds do work all the time but in order to see your image the 
epsg:4326 extents must be the same than the epsg:26719 extents of the image 
projected to epsg:4326.  

This may go wrong, I have not used cs2cs program much, but let's hope I did 
this correctly
So I tried to take the 26716 coordinates and reproject them to 4321 

d:\FWTools2.4.6>cs2cs +init=epsg:26716 +to +init=epsg:4326   #reproject from 
26716 -> to 4326
295125.999  4761593.049  #input coordinates 
 
89d30'45.961"W  42d58'52.872"N 0.000 #cs2cs output
308264.099 4778305.843   #input coordinates 
 
89d21'27.097"W  43d8'6.525"N 0.000   #cs2cs output

If this want OK it means that your image should be inside an epsg:4326 box 
defined by EXTENT -90 42 -89 44  
If it went wrong some more experienced cs2cs user will certainly help you 
further.

-Jukka-


> -Alkuperäinen viesti-
> Lähettäjä: mapserver-users-boun...@lists.osgeo.org 
> [mailto:mapserver-users-boun...@lists.osgeo.org] Puolesta Jon Miner
> Lähetetty: 1. helmikuuta 2010 14:46
> Vastaanottaja: mapserver-users@lists.osgeo.org
> Aihe: Re: [mapserver-users] New to mapserver;question about 
> USGS DRG/GeoTIFF files
> 
> 
> 
>  Thanks Jukka,
> 
> That did it.  I guess I need to better understand the 
> relationship between EXTENT and UNITS.  I thought that my 
> original bounds, in degrees, would have worked?
> 
> EXTENT -72 42 -71 44
> UNITS DD
> 
> - Jon
> 
> 
> 
> 
> 
> - Original Message 
> From: Rahkonen Jukka 
> To: Jon Miner ; mapserver-users@lists.osgeo.org
> Sent: Mon, February 1, 2010 12:28:09 AM
> Subject: Re: [mapserver-users] New to mapserver;question 
> about USGS  DRG/GeoTIFF files
> 
> Hi,
> 
> One possibility is that your image is not inside the spatial 
> window of your request.  Extents or the image in epsg:26719 
> system seem to be
> Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 
> 42d58'52.88"N)
> Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)
> 
> Extents of the MAP are in epsg:4326 system EXTENT -72 42 -71 
> 44. If the image is not inside this box you won't see it 
> which is natural.
> 
> Make a test by changing the MAP estents and projection 
> temporary to suit the image
> EXTENT 295125.999  4761593.049 308264.099 4778305.843
>   PROJECTION
>   "init=epsg:26719"
>   END
> 
> If you get the image now with shp2img you know that 
> everything is OK with the geotiff.
> 
> -Jukka Rahkonen-
> 
> Jon Miner wrote:
> 
> Hi Roger -
> 
> Thanks for the advice, but no change.  I did stumble across 
> the advice to try to run "shp2img" and even with debug set to 
> "5", get very little debugging information.
> 
> # shp2img -m test.map -all_debug 5 -o test.png
> GDAL: 
> GDALOpen(/opt/fgs/apps/gmap-demo-cvs_MS_VERSION_54/htdocs/./..
> /data/o43071a4.tif, this=0x9bd61b0) succeeds as GTiff.
> GDAL: GDALDeregister_GTiff() called.
> 
> The resulting test.png ends up being all white.
> 
> The actual TIFF file (10MB) can be found here:
> 
> http://trainguy.dyn.dhs.org/~jminer/o43071a4.tif
> 
> Thanks
> 
> 
> 
> 
> 
> From: Roger André 
> To: Jon Miner 
> Cc: mapserver-users@lists.osgeo.org
> Sent: Fri, January 29, 2010 4:16:08 PM
> Subject: Re: [mapserver-users] New to mapserver; question 
> about USGS  DRG/GeoTIFF files
> 
> Hi Jon,
> 
> You don't appear to have any projection info in your map 
> file.  Since you have "UNITS dd" in it, I assume you want to 
> display the map in lat/lon WGS84?  If so, your MAP file needs 
> an entry that looks like this:
> 
>   PROJECTION
> "init=epsg:4326"
>   END
> 
> and your usgs LAYER (the name should be quoted btw) should 
> have the UTM projection defined in it, like this:
> 
> LAYER
>   NAME "usgs"
>   PROJECTION  "init=epsg:26719" END
> 
> Best of luck,
> 
> Roger
> --
> 
> 
> 
> 
> 
> On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner  wrote:
> 
> 
> >
> >
> >>Hi -
> >
> >>I'm new to Mapserver, and so far impressed by it's 
> capabilities.  I've installed the FGS package via 
> fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played 
> with their demo program.
> >
> >>I'm having a problem with a USGS DRG file that is encoded 
> as a GeoTIFF file.  It doesn't display on the webpage, and 
> I'm no

Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-02-01 Thread Jon Miner


 Thanks Jukka,

That did it.  I guess I need to better understand the relationship between 
EXTENT and UNITS.  I thought that my original bounds, in degrees, would have 
worked?

EXTENT -72 42 -71 44
UNITS DD

- Jon





- Original Message 
From: Rahkonen Jukka 
To: Jon Miner ; mapserver-users@lists.osgeo.org
Sent: Mon, February 1, 2010 12:28:09 AM
Subject: Re: [mapserver-users] New to mapserver;question about USGS  
DRG/GeoTIFF files

Hi,

One possibility is that your image is not inside the spatial window of your 
request.  Extents or the image in epsg:26719 system seem to be
Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)

Extents of the MAP are in epsg:4326 system EXTENT -72 42 -71 44. If the image 
is not inside this box you won't see it which is natural.

Make a test by changing the MAP estents and projection temporary to suit the 
image
EXTENT 295125.999  4761593.049 308264.099, 4778305.843
  PROJECTION
  "init=epsg:26719"
  END

If you get the image now with shp2img you know that everything is OK with the 
geotiff.

-Jukka Rahkonen-

Jon Miner wrote:

Hi Roger -

Thanks for the advice, but no change.  I did stumble across the advice to try 
to run "shp2img" and even with debug set to "5", get very little debugging 
information.

# shp2img -m test.map -all_debug 5 -o test.png
GDAL: 
GDALOpen(/opt/fgs/apps/gmap-demo-cvs_MS_VERSION_54/htdocs/./../data/o43071a4.tif,
 this=0x9bd61b0) succeeds as GTiff.
GDAL: GDALDeregister_GTiff() called.

The resulting test.png ends up being all white.

The actual TIFF file (10MB) can be found here:

http://trainguy.dyn.dhs.org/~jminer/o43071a4.tif

Thanks





From: Roger André 
To: Jon Miner 
Cc: mapserver-users@lists.osgeo.org
Sent: Fri, January 29, 2010 4:16:08 PM
Subject: Re: [mapserver-users] New to mapserver; question about USGS  
DRG/GeoTIFF files

Hi Jon,

You don't appear to have any projection info in your map file.  Since you have 
"UNITS dd" in it, I assume you want to display the map in lat/lon WGS84?  If 
so, your MAP file needs an entry that looks like this:

  PROJECTION
"init=epsg:4326"
  END

and your usgs LAYER (the name should be quoted btw) should have the UTM 
projection defined in it, like this:

LAYER
  NAME "usgs"
  PROJECTION  "init=epsg:26719" END

Best of luck,

Roger
--





On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner  wrote:


>
>
>>Hi -
>
>>I'm new to Mapserver, and so far impressed by it's capabilities.  I've 
>>installed the FGS package via 
>>fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played with their demo 
>>program.
>
>>I'm having a problem with a USGS DRG file that is encoded as a GeoTIFF file.  
>>It doesn't display on the webpage, and I'm not sure where to start debugging 
>>this.  I've checked the file with gdalinfo:
>
>>Driver: GTiff/GeoTIFF
>>Files: o43071a4.tif
>>Size is 5388, 6854
>>Coordinate System is:
>>PROJCS["NAD27 / UTM zone 19N",
>>GEOGCS["NAD27",
>>DATUM["North_American_Datum_1927",
>>SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
>>AUTHORITY["EPSG","7008"]],
>>AUTHORITY["EPSG","6267"]],
>>PRIMEM["Greenwich",0],
>>UNIT["degree",0.0174532925199433],
>>AUTHORITY["EPSG","4267"]],
>>PROJECTION["Transverse_Mercator"],
>>PARAMETER["latitude_of_origin",0],
>>PARAMETER["central_meridian",-69],
>>PARAMETER["scale_factor",0.9996],
>>PARAMETER["false_easting",50],
>>PARAMETER["false_northing",0],
>>UNIT["metre",1,
>>AUTHORITY["EPSG","9001"]],
>>AUTHORITY["EPSG","26719"]]
>>Origin = (295125.99949388772,4778305.84265616444)
>>Pixel Size = (2.4384000,-2.4384000)
>>Metadata:
>>  AREA_OR_POINT=Area
>>  TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of Manchester North. 
>> Product:522644
>>  TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0
>>  TIFFTAG_DATETIME=1997:01:10 07:16:28
>>  TIFFTAG_XRESOLUTION=72
>>  TIFFTAG_YRESOLUTION=72
>>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>>Image Structure Metadata:
>>  COMPRESSION=PACKBITS
>>  INTERLEAVE=BAND
>>Corner Coordinates:
>>Upper Left  (  295125.999, 4778305.843) ( 71d31'7.66"W, 43d 7'54.16"

Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-01-31 Thread Rahkonen Jukka
Hi,

Another thing, your mapfile has:
 MINSCALEDENOM 24000
 MAXSCALEDENOM 24000

I means that the image is shown only at a scale 1:24000, exactly.  Comment out 
these lines with #.  The demo application is obviously asking maps always in 
this scale but it is difficult for you do the same with direct requests. And 
learn to send request with a browser by using CGI calls or WMS, then you will 
have full control on all the parameters.

-Jukka-


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


Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-01-31 Thread Rahkonen Jukka
Hi,

One possibility is that your image is not inside the spatial window of your 
request.  Extents or the image in epsg:26719 system seem to be
Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)

Extents of the MAP are in epsg:4326 system EXTENT -72 42 -71 44. If the image 
is not inside this box you won't see it which is natural.

Make a test by changing the MAP estents and projection temporary to suit the 
image
EXTENT 295125.999  4761593.049 308264.099, 4778305.843
  PROJECTION
  "init=epsg:26719"
  END

If you get the image now with shp2img you know that everything is OK with the 
geotiff.

-Jukka Rahkonen-

Jon Miner wrote:

Hi Roger -

Thanks for the advice, but no change.  I did stumble across the advice to try 
to run "shp2img" and even with debug set to "5", get very little debugging 
information.

# shp2img -m test.map -all_debug 5 -o test.png
GDAL: 
GDALOpen(/opt/fgs/apps/gmap-demo-cvs_MS_VERSION_54/htdocs/./../data/o43071a4.tif,
 this=0x9bd61b0) succeeds as GTiff.
GDAL: GDALDeregister_GTiff() called.

The resulting test.png ends up being all white.

The actual TIFF file (10MB) can be found here:

http://trainguy.dyn.dhs.org/~jminer/o43071a4.tif

Thanks





From: Roger André 
To: Jon Miner 
Cc: mapserver-users@lists.osgeo.org
Sent: Fri, January 29, 2010 4:16:08 PM
Subject: Re: [mapserver-users] New to mapserver; question about USGS  
DRG/GeoTIFF files

Hi Jon,

You don't appear to have any projection info in your map file.  Since you have 
"UNITS dd" in it, I assume you want to display the map in lat/lon WGS84?  If 
so, your MAP file needs an entry that looks like this:

  PROJECTION
"init=epsg:4326"
  END

and your usgs LAYER (the name should be quoted btw) should have the UTM 
projection defined in it, like this:

LAYER
  NAME "usgs"
  PROJECTION  "init=epsg:26719" END

Best of luck,

Roger
--





On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner  wrote:


>
>
>>Hi -
>
>>I'm new to Mapserver, and so far impressed by it's capabilities.  I've 
>>installed the FGS package via 
>>fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played with their demo 
>>program.
>
>>I'm having a problem with a USGS DRG file that is encoded as a GeoTIFF file.  
>>It doesn't display on the webpage, and I'm not sure where to start debugging 
>>this.  I've checked the file with gdalinfo:
>
>>Driver: GTiff/GeoTIFF
>>Files: o43071a4.tif
>>Size is 5388, 6854
>>Coordinate System is:
>>PROJCS["NAD27 / UTM zone 19N",
>>GEOGCS["NAD27",
>>DATUM["North_American_Datum_1927",
>>SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
>>AUTHORITY["EPSG","7008"]],
>>AUTHORITY["EPSG","6267"]],
>>PRIMEM["Greenwich",0],
>>UNIT["degree",0.0174532925199433],
>>AUTHORITY["EPSG","4267"]],
>>PROJECTION["Transverse_Mercator"],
>>PARAMETER["latitude_of_origin",0],
>>PARAMETER["central_meridian",-69],
>>PARAMETER["scale_factor",0.9996],
>>PARAMETER["false_easting",50],
>>PARAMETER["false_northing",0],
>>UNIT["metre",1,
>>AUTHORITY["EPSG","9001"]],
>>AUTHORITY["EPSG","26719"]]
>>Origin = (295125.99949388772,4778305.84265616444)
>>Pixel Size = (2.4384000,-2.4384000)
>>Metadata:
>>  AREA_OR_POINT=Area
>>  TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of Manchester North. 
>> Product:522644
>>  TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0
>>  TIFFTAG_DATETIME=1997:01:10 07:16:28
>>  TIFFTAG_XRESOLUTION=72
>>  TIFFTAG_YRESOLUTION=72
>>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>>Image Structure Metadata:
>>  COMPRESSION=PACKBITS
>>  INTERLEAVE=BAND
>>Corner Coordinates:
>>Upper Left  (  295125.999, 4778305.843) ( 71d31'7.66"W, 43d 7'54.16"N)
>>Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
>>Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)
>>Lower Right (  308264.099, 4761593.049) ( 71d21'5.95"W, 42d59'5.20"N)
>>Center  (  301695.049, 4769949.446) ( 71d26'6.44"W, 43d 3'29.80"N)
>>Band 1 Block=5388x1 Type=Byte, ColorInt

Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-01-31 Thread Jon Miner
Hi Roger -

Thanks for the advice, but no change.  I did stumble across the advice to try 
to run "shp2img" and even with debug set to "5", get very little debugging 
information.

# shp2img -m test.map -all_debug 5 -o test.png
GDAL: 
GDALOpen(/opt/fgs/apps/gmap-demo-cvs_MS_VERSION_54/htdocs/./../data/o43071a4.tif,
 this=0x9bd61b0) succeeds as GTiff.
GDAL: GDALDeregister_GTiff() called.

The resulting test.png ends up being all white.

The actual TIFF file (10MB) can be found here:

http://trainguy.dyn.dhs.org/~jminer/o43071a4.tif

Thanks





From: Roger André 
To: Jon Miner 
Cc: mapserver-users@lists.osgeo.org
Sent: Fri, January 29, 2010 4:16:08 PM
Subject: Re: [mapserver-users] New to mapserver; question about USGS  
DRG/GeoTIFF files

Hi Jon,

You don't appear to have any projection info in your map file.  Since you have 
"UNITS dd" in it, I assume you want to display the map in lat/lon WGS84?  If 
so, your MAP file needs an entry that looks like this:

  PROJECTION
"init=epsg:4326"
  END

and your usgs LAYER (the name should be quoted btw) should have the UTM 
projection defined in it, like this:

LAYER
  NAME "usgs"
  PROJECTION  "init=epsg:26719" END

Best of luck,

Roger
--





On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner  wrote:


>
>
>>Hi -
>
>>I'm new to Mapserver, and so far impressed by it's capabilities.  I've 
>>installed the FGS package via 
>>fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played with their demo 
>>program.
>
>>I'm having a problem with a USGS DRG file that is encoded as a GeoTIFF file.  
>>It doesn't display on the webpage, and I'm not sure where to start debugging 
>>this.  I've checked the file with gdalinfo:
>
>>Driver: GTiff/GeoTIFF
>>Files: o43071a4.tif
>>Size is 5388, 6854
>>Coordinate System is:
>>PROJCS["NAD27 / UTM zone 19N",
>>GEOGCS["NAD27",
>>DATUM["North_American_Datum_1927",
>>SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
>>AUTHORITY["EPSG","7008"]],
>>AUTHORITY["EPSG","6267"]],
>>PRIMEM["Greenwich",0],
>>UNIT["degree",0.0174532925199433],
>>AUTHORITY["EPSG","4267"]],
>>PROJECTION["Transverse_Mercator"],
>>PARAMETER["latitude_of_origin",0],
>>PARAMETER["central_meridian",-69],
>>PARAMETER["scale_factor",0.9996],
>>PARAMETER["false_easting",50],
>>PARAMETER["false_northing",0],
>>UNIT["metre",1,
>>AUTHORITY["EPSG","9001"]],
>>AUTHORITY["EPSG","26719"]]
>>Origin = (295125.99949388772,4778305.84265616444)
>>Pixel Size = (2.4384000,-2.4384000)
>>Metadata:
>>  AREA_OR_POINT=Area
>>  TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of Manchester North. 
>> Product:522644
>>  TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0
>>  TIFFTAG_DATETIME=1997:01:10 07:16:28
>>  TIFFTAG_XRESOLUTION=72
>>  TIFFTAG_YRESOLUTION=72
>>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
>>Image Structure Metadata:
>>  COMPRESSION=PACKBITS
>>  INTERLEAVE=BAND
>>Corner Coordinates:
>>Upper Left  (  295125.999, 4778305.843) ( 71d31'7.66"W, 43d 7'54.16"N)
>>Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
>>Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)
>>Lower Right (  308264.099, 4761593.049) ( 71d21'5.95"W, 42d59'5.20"N)
>>Center  (  301695.049, 4769949.446) ( 71d26'6.44"W, 43d 3'29.80"N)
>>Band 1 Block=5388x1 Type=Byte, ColorInterp=Palette
>>  Color Table (RGB with 256 entries)
>>0: 0,0,0,255
>>1: 255,255,255,255
>>2: 0,151,164,255
>>3: 203,0,23,255
>>4: 131,66,37,255
>>5: 201,234,157,255
>>6: 137,51,128,255
>>7: 255,234,0,255
>>8: 167,226,226,255
>>9: 255,184,184,255
>>   10: 218,179,214,255
>>   11: 209,209,209,255
>>   12: 207,164,142,255
>>(etc...)
>
>>And my map file looks like this:
>
>>MAP
>>NAME DEMO
>>STATUS ON
>>SIZE 400 300
>>SYMBOLSET ../etc/symbols.sym
>>EXTENT -72 42 -71 44
>>UNITS dd
>>SCALEDENOM 24000
>>SHAPEPATH "../data"
>>IMAGECOLOR 255 255 255
>>FONTSET ../etc/fonts.txt
>>CONFIG "MS_ERRORFILE" "/tmp/ms_error.txt"
>>#
>># Start of web interface definition
>>#
>>WEB
>>  MINSCALEDENOM 24000
>>  MAXSCALEDENOM 24000
>>#
>># On Windows systems, /tmp and /opt/fgs/tmp/ms_tmp/ should be created at the 
>>root
>># of the drive where the .MAP file resides.
>>#
>>  IMAGEPATH "/opt/fgs/tmp/ms_tmp/"
>>  IMAGEURL "/ms_tmp/"
>>  LOG "/tmp/gmap.log"
>>END
>
>>LAYER
>>  NAME usgs
>>  METADATA
>>"DESCRIPTION" "USGS GeoTIFF"
>>  END
>>  TYPE RASTER
>>  STATUS ON
>>  DATA o43071a4.tif
>>END
>
>>END # Map File
>
>
>
>>___
>>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] New to mapserver; question about USGS DRG/GeoTIFF files

2010-01-29 Thread Roger André
Hi Jon,

You don't appear to have any projection info in your map file.  Since you
have "UNITS dd" in it, I assume you want to display the map in lat/lon
WGS84?  If so, your MAP file needs an entry that looks like this:

  PROJECTION
"init=epsg:4326"
  END

and your usgs LAYER (the name should be quoted btw) should have the UTM
projection defined in it, like this:

LAYER
  NAME "usgs"
  PROJECTION  "init=epsg:26719" END

Best of luck,

Roger
--




On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner  wrote:

>
>
>
> Hi -
>
> I'm new to Mapserver, and so far impressed by it's capabilities.  I've
> installed the FGS package via
> fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played with their demo
> program.
>
> I'm having a problem with a USGS DRG file that is encoded as a GeoTIFF
> file.  It doesn't display on the webpage, and I'm not sure where to start
> debugging this.  I've checked the file with gdalinfo:
>
> Driver: GTiff/GeoTIFF
> Files: o43071a4.tif
> Size is 5388, 6854
> Coordinate System is:
> PROJCS["NAD27 / UTM zone 19N",
>GEOGCS["NAD27",
>DATUM["North_American_Datum_1927",
>SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
>AUTHORITY["EPSG","7008"]],
>AUTHORITY["EPSG","6267"]],
>PRIMEM["Greenwich",0],
>UNIT["degree",0.0174532925199433],
>AUTHORITY["EPSG","4267"]],
>PROJECTION["Transverse_Mercator"],
>PARAMETER["latitude_of_origin",0],
>PARAMETER["central_meridian",-69],
>PARAMETER["scale_factor",0.9996],
>PARAMETER["false_easting",50],
>PARAMETER["false_northing",0],
>UNIT["metre",1,
>AUTHORITY["EPSG","9001"]],
>AUTHORITY["EPSG","26719"]]
> Origin = (295125.99949388772,4778305.84265616444)
> Pixel Size = (2.4384000,-2.4384000)
> Metadata:
>  AREA_OR_POINT=Area
>  TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of Manchester
> North. Product:522644
>  TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0
>  TIFFTAG_DATETIME=1997:01:10 07:16:28
>  TIFFTAG_XRESOLUTION=72
>  TIFFTAG_YRESOLUTION=72
>  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
> Image Structure Metadata:
>  COMPRESSION=PACKBITS
>  INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (  295125.999, 4778305.843) ( 71d31'7.66"W, 43d 7'54.16"N)
> Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
> Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)
> Lower Right (  308264.099, 4761593.049) ( 71d21'5.95"W, 42d59'5.20"N)
> Center  (  301695.049, 4769949.446) ( 71d26'6.44"W, 43d 3'29.80"N)
> Band 1 Block=5388x1 Type=Byte, ColorInterp=Palette
>  Color Table (RGB with 256 entries)
>0: 0,0,0,255
>1: 255,255,255,255
>2: 0,151,164,255
>3: 203,0,23,255
>4: 131,66,37,255
>5: 201,234,157,255
>6: 137,51,128,255
>7: 255,234,0,255
>8: 167,226,226,255
>9: 255,184,184,255
>   10: 218,179,214,255
>   11: 209,209,209,255
>   12: 207,164,142,255
> (etc...)
>
> And my map file looks like this:
>
> MAP
> NAME DEMO
> STATUS ON
> SIZE 400 300
> SYMBOLSET ../etc/symbols.sym
> EXTENT -72 42 -71 44
> UNITS dd
> SCALEDENOM 24000
> SHAPEPATH "../data"
> IMAGECOLOR 255 255 255
> FONTSET ../etc/fonts.txt
> CONFIG "MS_ERRORFILE" "/tmp/ms_error.txt"
> #
> # Start of web interface definition
> #
> WEB
>  MINSCALEDENOM 24000
>  MAXSCALEDENOM 24000
> #
> # On Windows systems, /tmp and /opt/fgs/tmp/ms_tmp/ should be created at
> the root
> # of the drive where the .MAP file resides.
> #
>  IMAGEPATH "/opt/fgs/tmp/ms_tmp/"
>  IMAGEURL "/ms_tmp/"
>  LOG "/tmp/gmap.log"
> END
>
> LAYER
>  NAME usgs
>  METADATA
>"DESCRIPTION" "USGS GeoTIFF"
>  END
>  TYPE RASTER
>  STATUS ON
>  DATA o43071a4.tif
> END
>
> END # Map File
>
>
>
> ___
> 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


[mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

2010-01-29 Thread Jon Miner



Hi -

I'm new to Mapserver, and so far impressed by it's capabilities.  I've 
installed the FGS package via fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin 
and played with their demo program.

I'm having a problem with a USGS DRG file that is encoded as a GeoTIFF file.  
It doesn't display on the webpage, and I'm not sure where to start debugging 
this.  I've checked the file with gdalinfo:

Driver: GTiff/GeoTIFF
Files: o43071a4.tif
Size is 5388, 6854
Coordinate System is:
PROJCS["NAD27 / UTM zone 19N",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-69],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",50],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","26719"]]
Origin = (295125.99949388772,4778305.84265616444)
Pixel Size = (2.4384000,-2.4384000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of Manchester North. 
Product:522644
  TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0
  TIFFTAG_DATETIME=1997:01:10 07:16:28
  TIFFTAG_XRESOLUTION=72
  TIFFTAG_YRESOLUTION=72
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  COMPRESSION=PACKBITS
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  295125.999, 4778305.843) ( 71d31'7.66"W, 43d 7'54.16"N)
Lower Left  (  295125.999, 4761593.049) ( 71d30'45.52"W, 42d58'52.88"N)
Upper Right (  308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N)
Lower Right (  308264.099, 4761593.049) ( 71d21'5.95"W, 42d59'5.20"N)
Center  (  301695.049, 4769949.446) ( 71d26'6.44"W, 43d 3'29.80"N)
Band 1 Block=5388x1 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 255,255,255,255
2: 0,151,164,255
3: 203,0,23,255
4: 131,66,37,255
5: 201,234,157,255
6: 137,51,128,255
7: 255,234,0,255
8: 167,226,226,255
9: 255,184,184,255
   10: 218,179,214,255
   11: 209,209,209,255
   12: 207,164,142,255
(etc...)

And my map file looks like this:

MAP
NAME DEMO
STATUS ON
SIZE 400 300
SYMBOLSET ../etc/symbols.sym
EXTENT -72 42 -71 44
UNITS dd
SCALEDENOM 24000
SHAPEPATH "../data"
IMAGECOLOR 255 255 255
FONTSET ../etc/fonts.txt
CONFIG "MS_ERRORFILE" "/tmp/ms_error.txt"
#
# Start of web interface definition
#
WEB
  MINSCALEDENOM 24000
  MAXSCALEDENOM 24000
#
# On Windows systems, /tmp and /opt/fgs/tmp/ms_tmp/ should be created at the 
root
# of the drive where the .MAP file resides.
#
  IMAGEPATH "/opt/fgs/tmp/ms_tmp/"
  IMAGEURL "/ms_tmp/"
  LOG "/tmp/gmap.log"
END

LAYER
  NAME usgs
  METADATA
"DESCRIPTION" "USGS GeoTIFF"
  END
  TYPE RASTER
  STATUS ON
  DATA o43071a4.tif
END

END # Map File


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