Steve, This resolved the issue. I ran `gdalwarp /tmp/image.tiff /tmp/image_inverted.tiff` and it inverted the y axis coordinates again.
Chris On Sun, Apr 17, 2011 at 9:09 AM, Christopher Dedels <cded...@gmail.com>wrote: > I'll look into the N-S coordinates. I am using gdal_grid to create a grid > from a point layer and then using gdal_contour to create a line layer from > the grid. The input point layer and output line layer seem to have the > right N-S coordinates and display on my map just fine. > > Thanks, > > Chris > > > > On Sun, Apr 17, 2011 at 12:11 AM, Stephen Woodbridge < > wood...@swoodbridge.com> wrote: > >> Sorry, you had posted the gdalinfo, I was rushing out the door and missed >> that. But looking at it the image extents looks strange, it seems that the >> north-south coordinates are inverted. >> >> >> Upper Left (-123.8600311, 38.0116386) (123d51'36.11"W, 38d 0'41.90"N) >> Lower Left (-123.8600311, 53.7055626) (123d51'36.11"W, 53d42'20.03"N) >> Upper Right ( -97.1885986, 38.0116386) ( 97d11'18.96"W, 38d 0'41.90"N) >> Lower Right ( -97.1885986, 53.7055626) ( 97d11'18.96"W, 53d42'20.03"N) >> >> lower left should be the xmin, ymin and upper right should be the xmax, >> ymax values and so that looks weird. >> >> Otherwise I'm not seeing a problem. You might need to specify an image >> mode that you want mapserver to render you image in: >> >> IMAGETYPE jpeg >> >> >> -Steve W >> >> On 4/16/2011 5:52 PM, Christopher Dedels wrote: >> >>> gdalinfo output for my raster file is: >>> >>> Driver: GTiff/GeoTIFF >>> Files: /tmp/image.tiff >>> Size is 800, 600 >>> Coordinate System is: >>> GEOGCS["WGS 84", >>> DATUM["WGS_1984", >>> SPHEROID["WGS 84",6378137,298.257223563, >>> AUTHORITY["EPSG","7030"]], >>> AUTHORITY["EPSG","6326"]], >>> PRIMEM["Greenwich",0], >>> UNIT["degree",0.0174532925199433], >>> AUTHORITY["EPSG","4326"]] >>> Origin = (-123.860031127930000,38.011638641357401) >>> Pixel Size = (0.033339290618897,0.026156539916992) >>> Metadata: >>> AREA_OR_POINT=Area >>> Image Structure Metadata: >>> INTERLEAVE=BAND >>> Corner Coordinates: >>> Upper Left (-123.8600311, 38.0116386) (123d51'36.11"W, 38d 0'41.90"N) >>> Lower Left (-123.8600311, 53.7055626) (123d51'36.11"W, 53d42'20.03"N) >>> Upper Right ( -97.1885986, 38.0116386) ( 97d11'18.96"W, 38d 0'41.90"N) >>> Lower Right ( -97.1885986, 53.7055626) ( 97d11'18.96"W, 53d42'20.03"N) >>> Center (-110.5243149, 45.8586006) (110d31'27.53"W, 45d51'30.96"N) >>> Band 1 Block=800x2 Type=Float32, ColorInterp=Gray >>> >>> Regards, >>> >>> Chris >>> >>> On Sat, Apr 16, 2011 at 5:47 PM, Stephen Woodbridge >>> <wood...@swoodbridge.com <mailto:wood...@swoodbridge.com>> wrote: >>> >>> What does: >>> >>> gdalinfo /tmp/image.tif >>> >>> report? >>> >>> -Steve W >>> >>> >>> >>> On 4/16/2011 4:52 PM, Christopher Dedels wrote: >>> >>> Hmmm...missed that one. But changing the units to dd still >>> doesn't have >>> any effect. >>> >>> I am using mapserver 4.8.1. DEBUG is not an option in this >>> version. >>> mapserv -v produces the following output: >>> >>> MapServer version 4.8.1 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG >>> OUTPUT=WBMP >>> OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER >>> SUPPORTS=GEOS INPUT=EPPL7 INPUT=POSTGIS INPUT=GDAL INPUT=SHAPEFILE >>> >>> My request uri is: >>> >>> http://localhost/cgi-bin/mapserv?MAP=/tmp/mapNl7O7H.map&LAYERS=ALL&MODE=MAP >>> < >>> http://localhost/cgi-bin/mapserv?MAP=/tmp/mapNl7O7H.map&LAYERS=ALL&MODE=MAP >>> > >>> < >>> http://localhost/cgi-bin/mapserv?MAP=/tmp/mapNl7O7H.map&LAYERS=ALL&MODE=MAP >>> < >>> http://localhost/cgi-bin/mapserv?MAP=/tmp/mapNl7O7H.map&LAYERS=ALL&MODE=MAP >>> >> >>> >>> gdalinfo output for my raster file is: >>> >>> Driver: GTiff/GeoTIFF >>> Files: /tmp/image.tiff >>> Size is 800, 600 >>> Coordinate System is: >>> GEOGCS["WGS 84", >>> DATUM["WGS_1984", >>> SPHEROID["WGS 84",6378137,298.257223563, >>> AUTHORITY["EPSG","7030"]], >>> AUTHORITY["EPSG","6326"]], >>> PRIMEM["Greenwich",0], >>> UNIT["degree",0.0174532925199433], >>> AUTHORITY["EPSG","4326"]] >>> Origin = (-123.860031127930000,38.011638641357401) >>> Pixel Size = (0.033339290618897,0.026156539916992) >>> Metadata: >>> AREA_OR_POINT=Area >>> Image Structure Metadata: >>> INTERLEAVE=BAND >>> Corner Coordinates: >>> Upper Left (-123.8600311, 38.0116386) (123d51'36.11"W, 38d >>> 0'41.90"N) >>> Lower Left (-123.8600311, 53.7055626) (123d51'36.11"W, >>> 53d42'20.03"N) >>> Upper Right ( -97.1885986, 38.0116386) ( 97d11'18.96"W, 38d >>> 0'41.90"N) >>> Lower Right ( -97.1885986, 53.7055626) ( 97d11'18.96"W, >>> 53d42'20.03"N) >>> Center (-110.5243149, 45.8586006) (110d31'27.53"W, >>> 45d51'30.96"N) >>> Band 1 Block=800x2 Type=Float32, ColorInterp=Gray >>> >>> Regards, >>> >>> Chris >>> >>> >>> On Sat, Apr 16, 2011 at 4:34 PM, Stephen Woodbridge >>> <wood...@swoodbridge.com <mailto:wood...@swoodbridge.com> >>> <mailto:wood...@swoodbridge.com >>> <mailto:wood...@swoodbridge.com>>> wrote: >>> >>> On 4/16/2011 4:33 PM, Christopher Dedels wrote: >>> >>> Thanks for the suggestion. I have tried it and it >>> appears to >>> have no >>> effect. >>> >>> Chris >>> >>> >>> On Sat, Apr 16, 2011 at 4:26 PM, Stephen Woodbridge >>> <wood...@swoodbridge.com <mailto:wood...@swoodbridge.com> >>> <mailto:wood...@swoodbridge.com <mailto:wood...@swoodbridge.com>> >>> <mailto:wood...@swoodbridge.com <mailto:wood...@swoodbridge.com> >>> >>> <mailto:wood...@swoodbridge.com >>> <mailto:wood...@swoodbridge.com>>>> wrote: >>> >>> Try changing the layer status to default, ie: >>> >>> STATUS DEFAULT >>> >>> -Steve W >>> >>> >>> On 4/16/2011 2:22 PM, Christopher Dedels wrote: >>> >>> I am having some problems adding a raster layer >>> to my >>> map. The >>> image >>> produced by mapserver is blank. My mapfile is >>> as follows: >>> >>> MAP >>> SIZE 800 600 >>> UNITS meters >>> >>> >>> This is a problem, Your UNITs are meters, but your extents and >>> projection are DD! >>> >>> >>> EXTENT -123.8600311 38.0116386 >>> -97.1885986 >>> 53.7055626 >>> >>> PROJECTION >>> "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" >>> END >>> >>> LAYER >>> NAME "grid" >>> TYPE RASTER >>> DATA "/tmp/image.tiff" >>> STATUS ON >>> PROJECTION >>> "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" >>> END >>> END >>> END >>> >>> I have also noticed that changing the DATA path to >>> a >>> file which >>> does not >>> exist in the filesystem (or is not a >>> GDAL-supported file >>> format) >>> has no >>> effect on the rendered output and does not >>> produce an >>> error message. >>> >>> Can somebody help me where I am going wrong? >>> >>> Thanks >>> >>> >>> >>> _______________________________________________ >>> mapserver-users mailing list >>> mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org> >>> <mailto:mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org>> >>> <mailto:mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org> >>> <mailto:mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org>>> >>> >>> http://lists.osgeo.org/mailman/listinfo/mapserver-users >>> >>> >>> >>> >>> >>> _______________________________________________ >>> mapserver-users mailing list >>> mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org> >>> <mailto:mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org>> >>> http://lists.osgeo.org/mailman/listinfo/mapserver-users >>> >>> >>> _______________________________________________ >>> mapserver-users mailing list >>> mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org> >>> <mailto:mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org>> >>> >>> http://lists.osgeo.org/mailman/listinfo/mapserver-users >>> >>> >>> >>> >>> _______________________________________________ >>> mapserver-users mailing list >>> mapserver-users@lists.osgeo.org >>> <mailto:mapserver-users@lists.osgeo.org> >>> http://lists.osgeo.org/mailman/listinfo/mapserver-users >>> >>> >>> _______________________________________________ >>> mapserver-users mailing list >>> mapserver-users@lists.osgeo.org <mailto: >>> 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 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