Re: [UMN_MAPSERVER-USERS] Image location inside its corner coordinates

Wed, 09 Aug 2006 12:55:11 -0700

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.
 
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?  Are there other LAYER statements in the file?
 
     - Ed
 

Ed 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 coordinates

Hi 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.0004163
My extent is indeed like you said it should be:
    EXTENT -48.0004159 -25.0004164 -46.5004159 -24.0004163

It was because of the line below that you state that the units should be Decimal Degrees?
    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



    - Ed

Ed 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 coordinates

I 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

Reply via email to