geotiff problem since update

[William Reid]

Frank Warmerdam wrote:

> William Reid wrote:
>
>> Since upgrading my rasters have left the map.  
>
>
> ...
>
>> The image itself is readable by gdal.
>>
>> Driver: GTiff/GeoTIFF
>> Size is 1600, 1200
>> Coordinate System is `'      (Is this a problem?)
>> GCP Projection = GEOGCS["unnamed",DATUM["WGS_1984",SPHEROID["WGS 
>> 84",6378137,298.257223563],TOWGS84[0,0,0,-0,-0,-0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]] 
>>
>> GCP[  0]: Id=, Info=
>>          (0,0) -> (-78.7842,38.1912,0)
>> GCP[  1]: Id=, Info=
>>          (0,1200) -> (-78.7824,37.7893,0)
>> GCP[  2]: Id=, Info=
>>          (1600,0) -> (-78.118,38.1912,0)
>> GCP[  3]: Id=, Info=
>>          (1600,1200) -> (-78.1198,37.7893,0)
>
>
> Bill,
>
> It was my understanding that MapServer has never supported GDAL
> images with only GCPs for control information.  So I don't know how
> this used to work.

I am trying to digest this as I am not super familiar with all the 
various control information.

>
> Is it possible you have a .wld file associated with this GeoTIFF?   It
> may be that the internal geotiff info is now taking precidence over the
> .wld.
>
Frank,
You are tireless.  Thank you for the response.

I have a TAB file with the Geotiff (sorry I didn't mention this).  If I 
don't let gdal see it then the output from gdalinfo looks like this.
Driver: GTiff/GeoTIFF
Size is 1600, 1200
Coordinate System is `'
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 1200.0)
Upper Right ( 1600.0,    0.0)
Lower Right ( 1600.0, 1200.0)
Center      (  800.0,  600.0)
Band 1 Block=1600x8 Type=Byte, ColorInterp=Red
Band 2 Block=1600x8 Type=Byte, ColorInterp=Green
Band 3 Block=1600x8 Type=Byte, ColorInterp=Blue

This is also the same output I get after using your gdaltranslate 
suggestion.  I may not have fully understood the PROFILE=BASELINE 
portion of the command.


The TAB file has this.  (originates from EDX software for doing RF 
predictions, so if I find the problem I would like to help clue them in 
as well)

!table
!version 300
!charset WindowsLatin1

Definition Table
  File "Carters2.tif"
  Type "RASTER"
  (-78.784164,38.19123) (0,0) Label "",
  (-78.782352,37.789336) (0,1200) Label "",
  (-78.117974,38.19123) (1600,0) Label "",
  (-78.119787,37.789336) (1600,1200) Label ""
  CoordSys Earth Projection 1, 104
  Units "degree"

Finally,
Earlier in the week I attempted to create a simple .wld file to fix the 
situation but it didn't seem to help (I no longer have the attempt but I 
recall it being pretty easy).

I am off to read about GCP.

-=Bill

> If that is the case, you can do an operation like:
>
>  gdal_translate -co PROFILE=BASELINE old.tif new.tif
>  mv new.tif old.tif
>
> to strip away the geotiff information from a geotiff file.
>
> Best regards,