On 8/9/06, Ed McNierney <[EMAIL PROTECTED]> wrote:
My original motivation to adjust the ows_extent was to treat this problem but following your advice e changing the units to decimal degrees my zoom problem was solved, thanks a lot!
I did not had a PROJECTION element other than the one in my LAYER element. Now based on the lines bellow outputted by gdalinfo:
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
END
But I'm not sure if that's correct.
Here goes my hole mapfile (based on the mapserver demo Itasca mapfile):
No, I'm just trying with one at the moment for the sake of simplicity.
In fact I want to display a DEM (Digital Elevation Model) image through WCS (Web Coverage Service) in order to open it in the OpenJUMP application ( The use of OpenJUMP is a mandatory requirement). For that I want first to test the availability of WCS of this DEM in the MapServer Demo application.
That's my current strategy to achieve this objective, any suggestions would be most appreciated :) ( but I know, this can be a subject for a hole new post to this list).
Bruno -I would suggest not mixing problems at this point - let's not worry about ows_ metadata settings until you get the simple image display to work properly. Please leave out the METADATA statements until the image displays properly; then you'll be able to add the METADATA back in without being confused about what settings are helping and what settings are causing trouble.
My original motivation to adjust the ows_extent was to treat this problem but following your advice e changing the units to decimal degrees my zoom problem was solved, thanks a lot!
Yes, I suggested UNITS DD because your output and input image units are decimal degrees.
You should have EXTENT and UNITS statements at the top level in your map file,
and a PROJECTION block there, too (with the same settings as below, as you seem to want to display the image in its native unprojected coordinates). Is that the case? Can you post the whole map file?
I did not had a PROJECTION element other than the one in my LAYER element. Now based on the lines bellow outputted by gdalinfo:
GEOGCS["WGS 84",
DATUM["WGS_1984",
I've created the following PROJECTION element:DATUM["WGS_1984",
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
END
But I'm not sure if that's correct.
Here goes my hole mapfile (based on the mapserver demo Itasca mapfile):
MAP
NAME ITASCA
STATUS ON
SIZE 600 600
EXTENT -48.0004159 -25.0004164 -46.5004159 -24.0004163
UNITS DD
SHAPEPATH "data"
IMAGECOLOR 255 255 255
TEMPLATEPATTERN "itasca"
IMAGETYPE PNG
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
END
#
# Start of web interface definition (including WMS enabling metadata)
#
WEB
HEADER templates/header.html
TEMPLATE "set in index.html"
FOOTER templates/footer.html
MINSCALE 1000
MAXSCALE 1550000
IMAGEPATH "set in index.html"
IMAGEURL "set in index.html"
METADATA
WMS_TITLE "UMN MapServer Itasca Demo"
WMS_ABSTRACT "This is a UMN MapServer application for Itasca County located in north central Minnesota."
WMS_ACCESSCONSTRAINTS "none"
# change this value to match your setup
WMS_ONLINERESOURCE "http://10.10.5.171/cgi-bin/mapserv.exe "
WMS_SRS "EPSG:26915"
END
END
#
# Start of reference map
#
REFERENCE
IMAGE graphics/reference.png
EXTENT 324152.222 7384541.141 336175.222 7392179.141
SIZE 120 120
STATUS ON
MINBOXSIZE 5
MAXBOXSIZE 100
COLOR 255 0 0
OUTLINECOLOR 0 0 0
MARKERSIZE 8
MARKER 'star'
END
OUTPUTFORMAT
NAME GEOTIFFINT16
DRIVER "GDAL/GTiff"
MIMETYPE "image/tiff"
IMAGEMODE "INT16"
EXTENSION "tif"
END
LAYER
NAME "MyLayer"
STATUS OFF
TYPE RASTER
DUMP TRUE
DATA "SG-23-V-A.tif"
PROJECTION
"init=epsg:4326"
END
END
END
NAME ITASCA
STATUS ON
SIZE 600 600
EXTENT -48.0004159 -25.0004164 -46.5004159 -24.0004163
UNITS DD
SHAPEPATH "data"
IMAGECOLOR 255 255 255
TEMPLATEPATTERN "itasca"
IMAGETYPE PNG
PROJECTION
"proj=latlong"
"ellps=WGS84"
"datum=WGS84"
END
#
# Start of web interface definition (including WMS enabling metadata)
#
WEB
HEADER templates/header.html
TEMPLATE "set in index.html"
FOOTER templates/footer.html
MINSCALE 1000
MAXSCALE 1550000
IMAGEPATH "set in index.html"
IMAGEURL "set in index.html"
METADATA
WMS_TITLE "UMN MapServer Itasca Demo"
WMS_ABSTRACT "This is a UMN MapServer application for Itasca County located in north central Minnesota."
WMS_ACCESSCONSTRAINTS "none"
# change this value to match your setup
WMS_ONLINERESOURCE "http://10.10.5.171/cgi-bin/mapserv.exe "
WMS_SRS "EPSG:26915"
END
END
#
# Start of reference map
#
REFERENCE
IMAGE graphics/reference.png
EXTENT 324152.222 7384541.141 336175.222 7392179.141
SIZE 120 120
STATUS ON
MINBOXSIZE 5
MAXBOXSIZE 100
COLOR 255 0 0
OUTLINECOLOR 0 0 0
MARKERSIZE 8
MARKER 'star'
END
OUTPUTFORMAT
NAME GEOTIFFINT16
DRIVER "GDAL/GTiff"
MIMETYPE "image/tiff"
IMAGEMODE "INT16"
EXTENSION "tif"
END
LAYER
NAME "MyLayer"
STATUS OFF
TYPE RASTER
DUMP TRUE
DATA "SG-23-V-A.tif"
PROJECTION
"init=epsg:4326"
END
END
END
Are there other LAYER statements in the file?
No, I'm just trying with one at the moment for the sake of simplicity.
In fact I want to display a DEM (Digital Elevation Model) image through WCS (Web Coverage Service) in order to open it in the OpenJUMP application ( The use of OpenJUMP is a mandatory requirement). For that I want first to test the availability of WCS of this DEM in the MapServer Demo application.
That's my current strategy to achieve this objective, any suggestions would be most appreciated :) ( but I know, this can be a subject for a hole new post to this list).
- EdEd McNierney
President and Chief Mapmaker
TopoZone.com / Maps a la carte, Inc.
73 Princeton Street, Suite 305
North Chelmsford, MA 01863
Phone: +1 (978) 251-4242
Fax: +1 (978) 251-1396
[EMAIL PROTECTED]
From: Bruno Patini Furtado [mailto:[EMAIL PROTECTED]]
Sent: Wednesday, August 09, 2006 1:47 PM
To: Ed McNierney
Cc: [email protected]
Subject: Re: [UMN_MAPSERVER-USERS] Image location inside its corner coordinatesHi Ed,
On 8/9/06, Ed McNierney <[EMAIL PROTECTED]> wrote:Bruno -Could you post the EXTENT statement you're using? Two things to check - UNITS DD should be there (should not affect this problem, however) and your EXTENT needs to be in the order MINX MINY MAXX MAXY, so it should be:EXTENT -48.0004159 -25.0004164 -46.5004159 -24.0004163My extent is indeed like you said it should be:It was because of the line below that you state that the units should be Decimal Degrees?
EXTENT -48.0004159 -25.0004164 -46.5004159 -24.0004163
UNIT["degree",0.0174532925199433],
Your image is 1.5 degrees wide (X dimension) and 1 degree tall (Y dimension) so your output image size should also be that shape (600 pixels wide by 400 pixels high, for example) or you will see the map shrunk to fit the shape of the image.
Your example shows the right proportion as gdalinfo tells the image is [1800, 1200] but I would like to set the layer metadata ows_extent with the actual image coordinates, not size, in order to show the image a bit better in the MapServer demo. Below goes my mapinfo layer element as I'm trying to configure it.
LAYER
NAME "my DEM Image"
STATUS OFF
TYPE RASTER
DUMP TRUE
DATA "SG-23-V-A.tif"
PROJECTION
"init=epsg:4326"
END
METADATA
wcs_label "Some label"
ows_extent "- 48.000416 -24.000416 -48.001816 -24.001816"
wcs_resolution "0.00083 -0.00083"
ows_srs "EPSG:4326"
wcs_formats "GEOTIFFINT16"
wcs_nativeformat "geotiff"
END
END
- EdEd McNierney
President and Chief Mapmaker
TopoZone.com / Maps a la carte, Inc.
73 Princeton Street, Suite 305
North Chelmsford, MA 01863
Phone: +1 (978) 251-4242
Fax: +1 (978) 251-1396
[EMAIL PROTECTED]
From: UMN MapServer Users List [mailto:[email protected] ] On Behalf Of Bruno Patini Furtado
Sent: Wednesday, August 09, 2006 10:59 AM
To: [email protected]
Subject: [UMN_MAPSERVER-USERS] Image location inside its corner coordinatesI have a DEM image of which below I've copied the gdalinfo output.
Driver: GTiff/GeoTIFF
Size is 1800, 1200
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-48.000416,-24.000416)
Pixel Size = (0.00083333,-0.00083333 )
Metadata:
AREA_OR_POINT=Area
TIFFTAG_SOFTWARE=IMAGINE TIFF Support
Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: etif.c $ $Revision: 1.10.1.9 $ $Date: 2003/11/13 19:25:49EST $
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
Corner Coordinates:
Upper Left ( -48.0004159, -24.0004163) ( 48d 0'1.50"W, 24d 0'1.50"S)
Lower Left ( -48.0004159, -25.0004164 ) ( 48d 0'1.50"W, 25d 0'1.50"S)
Upper Right ( -46.5004159, -24.0004163) ( 46d30'1.50"W, 24d 0'1.50"S)
Lower Right ( -46.5004159, -25.0004164) ( 46d30'1.50"W, 25d 0'1.50"S)
Center ( -47.2504159, -24.5004163) ( 47d15'1.50"W, 24d30'1.50"S)
Band 1 Block=1800x2 Type=UInt16, ColorInterp=Gray
Overviews: 450x300, 225x150, 113x75, 57x38
I could take from it the EXTENT to use in my Mapfile using the Lower Left and Upper Right coordinates. But the problem is that the image is to small inside this coordinates and I can't find the image location inside then by this output. Is there some other utility to get this coordinates?
--
"Minds are like parachutes, they work best when open."
Bruno Patini Furtado
Software Developer
webpage: http://bpfurtado.net
software development blog: http://bpfurtado.livejournal.com
--
"Minds are like parachutes, they work best when open."
Bruno Patini Furtado
Software Developer
webpage: http://bpfurtado.net
software development blog: http://bpfurtado.livejournal.com
--
"Minds are like parachutes, they work best when open."
Bruno Patini Furtado
Software Developer
webpage: http://bpfurtado.net
software development blog: http://bpfurtado.livejournal.com
