comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado
Home » Public Forums » archive » Re: getting GeoTiff Information
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend  		Switch to threaded view of this topic Create a new topic Submit Reply
Re: getting GeoTiff Information [message #72190] Mon, 16 August 2010 13:01
Steve McAlister is currently offline  Steve McAlister
Messages: 3
Registered: August 2010
Junior Member
On Aug 16, 11:22 am, Paulo Penteado <pp.pente...@gmail.com> wrote:
> On Aug 16, 2:04 pm, Steve McAlister <s.mcalist...@gmail.com> wrote:
>
>
>
>
>
>> geotiff keyword from read_tiff will give me projection information in
>> structure of GeoTiff like:
>
>>    MODELPIXELSCALETAG
>>                    DOUBLE    Array[3]
>>    MODELTIEPOINTTAG
>>                    DOUBLE    Array[6, 1]
>>    GTMODELTYPEGEOKEY
>>                    INT              1
>>    GTRASTERTYPEGEOKEY
>>                    INT              1
>>    GEOGLINEARUNITSGEOKEY
>>                    INT           9001
>>    GEOGANGULARUNITSGEOKEY
>>                    INT           9102
>>    PROJECTEDCSTYPEGEOKEY
>>                    INT          26713
>
>> But, I want something more sensible for users to read like:
>
>>                   PCS: 'NAD83 / Massachusetts Mainland'
>>            Projection: 'SPCS83 Massachusetts Mainland zone (m)'
>>                MapSys: 'STATE_PLANE_83'
>>                  Zone: 2001
>>          CTProjection: 'CT_LambertConfConic_2SP'
>>              ProjParm: [7x1 double]
>>            ProjParmId: {7x1 cell}
>>                   GCS: 'NAD83'
>>                 Datum: 'North American Datum 1983'
>>             Ellipsoid: 'GRS 1980'
>
> Yes, the structure uses the numeric codes. Maybe David's geocoord
> (http://www.dfanning.com/documents/programs.html#GEOCOORD) can be
> useful in translating them into names.
>
> Also, see
>
> http://www.dfanning.com/map_tips/tiffoverlay.html

Thank you Paul! David's code was very helpful.

Report message to a moderator

Re: getting GeoTiff Information [message #72195 is a reply to message #72190] Mon, 16 August 2010 10:22 Go to previous message
penteado is currently offline  penteado
Messages: 866
Registered: February 2018
Senior Member
Administrator
On Aug 16, 2:04 pm, Steve McAlister <s.mcalist...@gmail.com> wrote:
> geotiff keyword from read_tiff will give me projection information in
> structure of GeoTiff like:
>
>    MODELPIXELSCALETAG
>                    DOUBLE    Array[3]
>    MODELTIEPOINTTAG
>                    DOUBLE    Array[6, 1]
>    GTMODELTYPEGEOKEY
>                    INT              1
>    GTRASTERTYPEGEOKEY
>                    INT              1
>    GEOGLINEARUNITSGEOKEY
>                    INT           9001
>    GEOGANGULARUNITSGEOKEY
>                    INT           9102
>    PROJECTEDCSTYPEGEOKEY
>                    INT          26713
>
> But, I want something more sensible for users to read like:
>
>                   PCS: 'NAD83 / Massachusetts Mainland'
>            Projection: 'SPCS83 Massachusetts Mainland zone (m)'
>                MapSys: 'STATE_PLANE_83'
>                  Zone: 2001
>          CTProjection: 'CT_LambertConfConic_2SP'
>              ProjParm: [7x1 double]
>            ProjParmId: {7x1 cell}
>                   GCS: 'NAD83'
>                 Datum: 'North American Datum 1983'
>             Ellipsoid: 'GRS 1980'

Yes, the structure uses the numeric codes. Maybe David's geocoord
(http://www.dfanning.com/documents/programs.html#GEOCOORD) can be
useful in translating them into names.

Also, see

http://www.dfanning.com/map_tips/tiffoverlay.html

Report message to a moderator

Re: getting GeoTiff Information [message #72196 is a reply to message #72195] Mon, 16 August 2010 10:04 Go to previous message
Steve McAlister is currently offline  Steve McAlister
Messages: 3
Registered: August 2010
Junior Member
On Aug 16, 10:47 am, Paulo Penteado <pp.pente...@gmail.com> wrote:
> On Aug 16, 1:41 pm, Steve McAlister <s.mcalist...@gmail.com> wrote:
>
>> Hello all!
>
>> I am aware of the 'geotiffinfo' command in Matlab.http://www.mathworks.de/access/helpdesk/help/toolbox/ map/ref/geotiffi...
>
>> Do you know of equivalent way of converting geotiff structure to a
>> readable format using IDL only (without ENVI)?
>
>> Thanks a lot for help!
>
>> Steve
>
> Sounds like the geotiff keyword to read_tiff().

Hi Paulo,
geotiff keyword from read_tiff will give me projection information in
structure of GeoTiff like:

MODELPIXELSCALETAG
DOUBLE Array[3]
MODELTIEPOINTTAG
DOUBLE Array[6, 1]
GTMODELTYPEGEOKEY
INT 1
GTRASTERTYPEGEOKEY
INT 1
GEOGLINEARUNITSGEOKEY
INT 9001
GEOGANGULARUNITSGEOKEY
INT 9102
PROJECTEDCSTYPEGEOKEY
INT 26713

But, I want something more sensible for users to read like:

PCS: 'NAD83 / Massachusetts Mainland'
Projection: 'SPCS83 Massachusetts Mainland zone (m)'
MapSys: 'STATE_PLANE_83'
Zone: 2001
CTProjection: 'CT_LambertConfConic_2SP'
ProjParm: [7x1 double]
ProjParmId: {7x1 cell}
GCS: 'NAD83'
Datum: 'North American Datum 1983'
Ellipsoid: 'GRS 1980'

Report message to a moderator

Re: getting GeoTiff Information [message #72197 is a reply to message #72196] Mon, 16 August 2010 09:47 Go to previous message
penteado is currently offline  penteado
Messages: 866
Registered: February 2018
Senior Member
Administrator
On Aug 16, 1:41 pm, Steve McAlister <s.mcalist...@gmail.com> wrote:
> Hello all!
>
> I am aware of the 'geotiffinfo' command in Matlab.http://www.mathworks.de/access/helpdesk/help/toolbox/ map/ref/geotiffi...
>
> Do you know of equivalent way of converting geotiff structure to a
> readable format using IDL only (without ENVI)?
>
> Thanks a lot for help!
>
> Steve

Sounds like the geotiff keyword to read_tiff().

Report message to a moderator

  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Combing structures
Next Topic: How do I plot two images with !P.Multi, while keep the aspect ration reflecting the image size

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 11:44:33 PDT 2025

Total time taken to generate the page: 0.00693 seconds
.:: Contact :: Home ::.

Powered by: FUDforum 3.0.7.
Copyright ©2001-2013 FUDforum Bulletin Board Software