Come determinare se un tiff è georeferenziato o no?

40

Ho ricevuto un file tiff da qualcuno, ma non sono sicuro che sia georeferenziato. C'è un modo semplice per controllare questo?

geotiff-tiff  georeferencing 
— Jeff Storey
fonte
3
Che software stai usando? Se si utilizza ArcMap, caricarlo in un frame di dati, fare clic con il tasto destro del mouse su di esso e selezionare le proprietà, è possibile trovare informazioni sui file da lì. Sono sicuro che altri software hanno procedure simili.
grazie, abbastanza facile. Non sapevo se esistesse uno strumento da riga di comando che avrei potuto eseguire (magari usando gdal o qualcosa del genere).
— Jeff Storey,

Risposte:

42

Sembra che gdalinfo ti dirà:

Una figura non georeferenziata esportata dal layout ArcMap:

C:\Temp>gdalinfo figure1.tif
Driver: GTiff/GeoTIFF
Files: figure1.tif
Size is 244, 210
Coordinate System is `'
Metadata:
  TIFFTAG_XRESOLUTION=96
  TIFFTAG_YRESOLUTION=96
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  210.0)
Upper Right (  244.0,    0.0)
Lower Right (  244.0,  210.0)
Center      (  122.0,  105.0)
Band 1 Block=244x11 Type=Byte, ColorInterp=Red
Band 2 Block=244x11 Type=Byte, ColorInterp=Green
Band 3 Block=244x11 Type=Byte, ColorInterp=Blue

Una foto aerea:

C:\Temp>gdalinfo Adona_1.tif
Driver: GTiff/GeoTIFF
Files: Adona_1.tif
       Adona_1.tif.ovr
       Adona_1.tfw
       Adona_1.tif.aux.xml
Size is 5426, 5522
Coordinate System is:
PROJCS["NAD83 / Arkansas North",
    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["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",36.23333333333333],
    PARAMETER["standard_parallel_2",34.93333333333333],
    PARAMETER["latitude_of_origin",34.33333333333334],
    PARAMETER["central_meridian",-92],
    PARAMETER["false_easting",1312333.333333333],
    PARAMETER["false_northing",0],
    UNIT["US survey foot",0.3048006096012192,
        AUTHORITY["EPSG","9003"]],
    AUTHORITY["EPSG","26951"]]
Origin = (1045317.554567784100000,261147.746197238540000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
  TIFFTAG_SOFTWARE=Adobe Photoshop CS4 Windows
  TIFFTAG_DATETIME=2010:08:04 07:39:30
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 1045317.555,  261147.746) ( 92d53'31.49"W, 35d 2'50.94"N)
Lower Left  ( 1045317.555,  258386.746) ( 92d53'31.19"W, 35d 2'23.64"N)
Upper Right ( 1048030.555,  261147.746) ( 92d52'58.86"W, 35d 2'51.19"N)
Lower Right ( 1048030.555,  258386.746) ( 92d52'58.56"W, 35d 2'23.88"N)
Center      ( 1046674.055,  259767.246) ( 92d53'15.03"W, 35d 2'37.41"N)
Band 1 Block=5426x1 Type=Byte, ColorInterp=Red
  Overviews: 2713x2761, 1357x1381, 679x691, 340x346, 170x173
Band 2 Block=5426x1 Type=Byte, ColorInterp=Green
  Overviews: 2713x2761, 1357x1381, 679x691, 340x346, 170x173
Band 3 Block=5426x1 Type=Byte, ColorInterp=Blue
  Overviews: 2713x2761, 1357x1381, 679x691, 340x346, 170x173
— Chad Cooper
fonte
16

Aggiungendo alla risposta del Ciad: esiste un file .TFW che accompagna il TIF? In tal caso, l'immagine ha informazioni georef che non sono incorporate nell'intestazione. Non so se GDALINFO lo prenderà o no.

Sfortunatamente un TFW ti fornisce solo coordinate e dimensioni delle celle. Non specifica a quale CRS / proiezione / datum si fa riferimento all'immagine, né quali unità vengono utilizzate per la proiezione.

Buon punto, Kevin.
— Chad Cooper,
non c'è tfw, quindi speravo fosse incorporato nell'intestazione. Farò alcuni test per vedere quali informazioni gdal fornisce sul file. Grazie.
— Jeff Storey,
Non direi che sia una buona indicazione; ArcMap 10 crea GeoTIFF con georeferenziazione integrata e file tfw ridondanti.
— Sarà l'
7

Nel mondo ESRI, se l'immagine è georeferenziata, ma non rettificata, gdalinfo non restituirà le informazioni di proiezione

Driver: GTiff/GeoTIFF
Files: COVE_COGEO.tif
       COVE_COGEO.aux
       Cove_COGEO.rrd
       COVE_COGEO.tif.aux.xml
Size is 6618, 8176
Coordinate System is `'
Metadata:
  TIFFTAG_IMAGEDESCRIPTION=
  TIFFTAG_MAXSAMPLEVALUE=1
  TIFFTAG_MINSAMPLEVALUE=0
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=300
  TIFFTAG_YRESOLUTION=300
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 8176.0)
Upper Right ( 6618.0,    0.0)
Lower Right ( 6618.0, 8176.0)
Center      ( 3309.0, 4088.0)
Band 1 Block=6618x512 Type=Byte, ColorInterp=Red
  Overviews: 1655x2044, 828x1022, 414x511, 207x256, 104x128, 52x64
Band 2 Block=6618x512 Type=Byte, ColorInterp=Green
  Overviews: 1655x2044, 828x1022, 414x511, 207x256, 104x128, 52x64
Band 3 Block=6618x512 Type=Byte, ColorInterp=Blue
  Overviews: 1655x2044, 828x1022, 414x511, 207x256, 104x128, 52x64

Se guardi nel tif.aux.xml troverai le informazioni sulla proiezione: 

<?xml version="1.0" encoding="UTF-8"?>-
<GeodataXform xmlns:typens="http://www.esri.com/schemas/ArcGIS/9.2" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:type="typens:PolynomialXform">
<PolynomialOrder>1</PolynomialOrder>-<SpatialReference xsi:type="typens:ProjectedCoordinateSystem">  <WKT>PROJCS["NAD_1983_UTM_Zone_15N",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-93.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]</WKT>
      <XOrigin>-5682839.5944962</XOrigin>
      <YOrigin>-10997760.6594148</YOrigin>
      <XYScale>409501512.793437</XYScale>
      <ZOrigin>-100000</ZOrigin>
      <ZScale>10000</ZScale>
      <MOrigin>-100000</MOrigin>
      <MScale>10000</MScale>
      <XYTolerance>0.001</XYTolerance>
      <ZTolerance>0.001</ZTolerance>
      <MTolerance>0.001</MTolerance>
      <HighPrecision>true</HighPrecision>
    </SpatialReference>-<SourceGCPs xsi:type="typens:ArrayOfDouble">
    <Double>1.57734892030567</Double>
    <Double>25.8359122977818</Double>
    <Double>20.4046417720032</Double>
    <Double>25.9084153582944</Double>
    <Double>20.4733053142072</Double>
    <Double>3.14623602371285</Double>
    <Double>1.61802731006583</Double>
    <Double>3.1060187781668</Double>
    </SourceGCPs>-<TargetGCPs xsi:type="typens:ArrayOfDouble">
    <Double>362267.039129157</Double>
    <Double>3818630.52804079</Double>
    <Double>373747.691530917</Double>
    <Double>3818468.79669514</Double>
    <Double>373558.508499628</Double>
    <Double>3804605.04409216</Double>
    <Double>362064.5523863</Double>
    <Double>3804768.11127776</Double>
    </TargetGCPs>
    </GeodataXform>

Tuttavia, una volta corretto gdalinfo funzionerà e i dati non saranno più nel tif.aux.xml, quindi potresti dover cercare in entrambi i posti.

— JasonT
fonte
3

C'è anche listgeo , parte di libgeotiff .

Inoltre, se stai usando Windows, c'è una comoda applicazione GUI per esso e geotifcp (anche parte di libgeotiff) che ti consente di scaricare e ripristinare facilmente i metadati su un geotiff - utile se vuoi modificare il tuo geotiff in modo non spaziale software di modifica delle immagini abilitato o per l'applicazione di metadati spaziali ai risultati di una classificazione, ecc.

— FOD
fonte
2

Come determinare se un tiff è georeferenziato o no?

Scarica ftp://ftp.remotesensing.org/geotiff/libgeotiff/listgeo_GUI.zip

In listgeo_GUI.zip:

Strumenti GeoTIFF in forma di interfaccia grafica

(1). "listgeo.exe" - programma di utilità per il dumping dei metadati di un file GeoTIFF. (2). "geotifcp.exe" - programma di utilità per l'applicazione di metadati a un file TIFF, rendendolo un file GeoTIFF.

Sono stati distribuiti liberamente con la popolare libreria libgeotiff.

— Paganel
fonte
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.