[mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
Rahkonen Jukka
Jukka.Rahkonen at mmmtike.fi
Mon Feb 1 05:07:58 PST 2010
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-bounces at lists.osgeo.org
> [mailto:mapserver-users-bounces at lists.osgeo.org] Puolesta Jon Miner
> Lähetetty: 1. helmikuuta 2010 14:46
> Vastaanottaja: mapserver-users at 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 <Jukka.Rahkonen at mmmtike.fi>
> To: Jon Miner <cp_n18 at yahoo.com>; mapserver-users at 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é <randre at gmail.com>
> To: Jon Miner <cp_n18 at yahoo.com>
> Cc: mapserver-users at 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 <cp_n18 at yahoo.com> 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",500000],
> >> PARAMETER["false_northing",0],
> >> UNIT["metre",1,
> >> AUTHORITY["EPSG","9001"]],
> >> AUTHORITY["EPSG","26719"]]
> >>Origin = (295125.999493999988772,4778305.842659999616444)
> >>Pixel Size = (2.438400000000000,-2.438400000000000)
> >>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 at lists.osgeo.org
> >http://lists.osgeo.org/mailman/listinfo/mapserver-users
> >
>
>
>
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/mapserver-users
>
More information about the MapServer-users
mailing list