[GRASS-user] Import/conversion problem

Adam Dershowitz, Ph.D., P.E. adershowitz at exponent.com
Fri Feb 11 12:23:39 EST 2011


I am trying to import some data in GRASS and have run into a problem.  I hope someone can give me a little insight to help out.
I downloaded some aerial imagery from here:
http://www.michigan.gov/dnr/0,1607,7-153-10371_14546-30211--,00.html

I want to import it into a UTM zone 17 project.  So I did the following:
gdalwarp -t_srs '+proj=utm +zone=17 +datum=WGS84' input.sid output_utm.sid

gdapwarp responds with an error, but then continues with the conversion:
ERROR 6: Failed to initialize PROJ.4 with `+proj=omerc +lat_0=0 +lonc=0 +alpha=0 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs '.
lat_0 = 0 or 90 or alpha = 90
0...10...20...30...40...50...60...70...80...90...100 - done.

If I then try to import the image, into my project, it does work, but places it WAY off from where it should be (millions of meters!).  I guess I should not be surprised, since the error says that it is trying to use a lat/long of 0.  But I don't understand why it is not working correctly.  

So, clearly something is wrong with how I am doing the conversion.  Am I missing a flag or something?  It looks to me like gdalinfo is doing a correct read, but somehow gdalwarp is missing something on the input.
Any help would be greatly appreciated.  

If I do gdalinfo on the input file, here is what I get:

Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: oakgrove_sw.sid
       oakgrove_sw.sdw
       oakgrove_sw.sid.aux.xml
Size is 6415, 8061
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE",
    GEOGCS["NAD83",
        DATUM["North_American_Datum_1983",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6269"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4269"]],
    PROJECTION["Hotine_Oblique_Mercator"],
    PARAMETER["latitude_of_center",0],
    PARAMETER["longitude_of_center",0],
    PARAMETER["azimuth",0],
    PARAMETER["rectified_grid_angle",90],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (663144.500000000000000,240035.500000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  IMAGE__INPUT_NAME=K:\blk18\oakgrove_sw.tif
  IMAGE__INPUT_FILE_SIZE=155199578.000000
  GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
  GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
  GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
  GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
  GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
  GEOTIFF_NUM__2048__GeographicTypeGeoKey=4269
  GEOTIFF_CHAR__GeographicTypeGeoKey=GCS_NAD83
  GEOTIFF_NUM__2060__GeogAzimuthUnitsGeoKey=9102
  GEOTIFF_CHAR__GeogAzimuthUnitsGeoKey=Angular_Degree
  GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32767
  GEOTIFF_CHAR__ProjectedCSTypeGeoKey=User-Defined
  GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
  GEOTIFF_NUM__3074__ProjectionGeoKey=32767
  GEOTIFF_CHAR__ProjectionGeoKey=User-Defined
  GEOTIFF_NUM__3075__ProjCoordTransGeoKey=3
  GEOTIFF_CHAR__ProjCoordTransGeoKey=CT_ObliqueMercator
  GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001
  GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter
  GEOTIFF_NUM__3088__ProjCenterLongGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  GEOTIFF_NUM__3089__ProjCenterLatGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  GEOTIFF_NUM__3090__ProjCenterEastingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  GEOTIFF_NUM__3091__ProjCenterNorthingGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  GEOTIFF_NUM__3093__ProjScaleAtCenterGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  GEOTIFF_NUM__3094__ProjAzimuthAngleGeoKey=0.999600,45.309167,2546731.496000,-4354009.816000,337.255560,-86.000000
  IMAGE__Z_RESOLUTION=0.000000
  GEO__ModelTypeGeoKey=1
  GEO__ProjectedCSTypeGeoKey=32767
  GEO__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Oblique Mercator (Hotine)|IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Oblique Mercator (Hotine)
Units = meters
GeoTIFF Units = meters|
  IMAGE__Z_ORIGIN=0.000000
  IMAGE__COMPRESSION_VERSION=1,6,1
  IMAGE__TARGET_COMPRESSION_RATIO=20.000000
  IMAGE__COMPRESSION_NLEV=4
  IMAGE__COMPRESSION_WEIGHT=4.000000
  IMAGE__COMPRESSION_GAMMA=2.000000
  IMAGE__COMPRESSION_BLOCK_SIZE=512
  IMAGE__CREATION_DATE=Wed Mar 10 17:11:37 2004

  VERSION=MG2
  IMAGE__BITS_PER_SAMPLE=8
  IMAGE__WIDTH=6415
  IMAGE__HEIGHT=8061
  IMAGE__COLOR_SCHEME=0
  IMAGE__DATA_TYPE=0
  IMAGE__XY_ORIGIN=663145.000000,240035.000000
  IMAGE__X_RESOLUTION=1.000000
  IMAGE__Y_RESOLUTION=1.000000
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  663144.500,  240035.500) 
Lower Left  (  663144.500,  231974.500) 
Upper Right (  669559.500,  240035.500) 
Lower Right (  669559.500,  231974.500) 
Center      (  666352.000,  236005.000) 
Band 1 Block=1024x128 Type=Byte, ColorInterp=Red
  Minimum=0.000, Maximum=255.000, Mean=97.677, StdDev=49.163
  Overviews: 3208x4031, 1604x2016, 802x1008, 401x504
Band 2 Block=1024x128 Type=Byte, ColorInterp=Green
  Minimum=0.000, Maximum=255.000, Mean=97.555, StdDev=49.371
  Overviews: 3208x4031, 1604x2016, 802x1008, 401x504
Band 3 Block=1024x128 Type=Byte, ColorInterp=Blue
  Minimum=0.000, Maximum=255.000, Mean=97.402, StdDev=49.024
  Overviews: 3208x4031, 1604x2016, 802x1008, 401x504



--Adam





More information about the grass-user mailing list