[Xastir] Difficulty loading Geotif
Brian D Heaton
brian.heaton at janusresearch.com
Fri Mar 26 14:59:36 EST 2004
Curt,
Yep, I looked at the the DOI10s against a couple of DRGs late last
night. The biggest difference on first inspection is that the DRGs are
projected and the DOI10s are GCS.
I've got it locked down to the line in map_tif.c where we first fail
part of the process. Its line 1234 and I've pasted it below. I'm still
poking around to see exactly what is failing. I'm leaning towards
proper handling of the GCS_WGS84 projection definition.
Is there a different function call to convert GCS (non-projected)
rasters to the internal coordinate systems?
------------------------------------------------------
if ( GTIFProj4FromLatLong( &defn, 1, &xxx, &yyy ) )
------------------------------------------------------
THX/BDH
On Fri, 2004-03-26 at 14:15, Curt, WE7U wrote:
> On Thu, 25 Mar 2004, Brian D Heaton wrote:
>
> > XASTIR debugging output.....
> > ---------------------------------------------
> > Load maps start
> > Load maps Open map file
> > Found mapname: TESTING/test1.tif
> > Load maps -TESTING/test1.tif
> > draw_map: Searching for map driver
> > draw_map: Found map driver: tif: 3
> > draw_map: /usr/local/share/xastir/maps/TESTING/test1.tif
> > /usr/local/share/xastir/maps/TESTING/test1.tif
> > /usr/local/share/xastir/maps/TESTING/test1.fgd
> > /usr/local/share/xastir/maps/metadata/test1.fgd
> > Couldn't open '.fgd' file, assuming no map collar to chop
> > /usr/local/share/xastir/maps/TESTING/test1.tif
> > XTIFFOpen
> > GTIFNew
> > GTIFDirectoryInfo
> > GTIFGetDefn
> > GCS: 4326/WGS 84
> > Datum: 6326/World Geodetic System 1984
> > Ellipsoid: 7030/WGS 84 (6378137.00,6356752.31)
> > Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
> >
> > Corner Coordinates:
> > Upper Left ( 43.000, 35.000)
> > Lower Right ( 44.000, 35.000)
> > Lower Right ( 43.000, 34.500)
> > Lower Right ( 44.000, 34.500)
> > Datum: 6326/World Geodetic System 1984
> > ***** Skipping Datum Conversion
> > Width: 8406
> > Height: 5562
> > Photometric: 3
> > Rows Per Strip: 1
> > Bits Per Sample: 8
> > Samples Per Pixel: 1
> > Planar Config: 1
> > Loading geoTIFF map: /usr/local/share/xastir/maps/TESTING/test1.tif
> >
> > NW neat-line corner = 43.000000 35.000000
> > Problem in translating
> > NE neat-line corner = 44.000000 35.000000
> > Problem in translating
> > SW neat-line corner = 43.000000 34.500000
> > Problem in translating
> > SE neat-line corner = 44.000000 34.500000
> > Problem in translating
> > Crop points (pixels):
> > Top: 3 Bottom: 5561 Left: 3 Right: 8405
> >
> > Initial Bytes Per Row: 8406
> > xastir_left_x_increment: 0.000000 80280000 80280000 0.000000 3 3 5561 3
> > xastir_left_y_increment: 30000.000000 19980000 19800000 1.000000 -3 3 5561 3
> > xastir_right_x_increment: 0.000000 80640000 80640000 0.000000 -3 -3 5561 3
> > xastir_right_y_increment: 30000.000000 19980000 19800000 1.000000 -3 3 5561 3
> > Left x increments: 0.000000 0.000000
> > Left y increments: 30000.000000 1.000000
> > Right x increments: 0.000000 0.000000
> > Right y increments: 30000.000000 1.000000
> > xastir_top_y_increment: 0.000000 19800000 19800000 0.000000 3 3 8405 3
> > xastir_bottom_y_increment: 0.000000 19980000 19980000 0.000000 -3 -3 8405 3
> > PixelScale: 0.000119 0.000090 0.000000
> > Size: x 262144, y 262144
> > SkipRows: 556
> > 5563 rows read in
> > Load maps stop
> > ---------------------------------------------------
>
>
> Here's mine from the o48123h1.tif file, again for comparison:
>
> Load maps start
> draw_map: Searching for map driver
> draw_map: Found map driver: tif: 3
> draw_map: /usr/local/share/xastir/maps/World/USA/WA/Tif/o48123h1.tif
> /usr/local/share/xastir/maps/World/USA/WA/Tif/o48123h1.tif
> /usr/local/share/xastir/maps/World/USA/WA/Tif/o48123h1.fgd
> West Bounding: -123.125000
> East Bounding: -123.000000
> North Bounding: 49.000000
> South Bounding: 48.875000
> 48.875000 49.000000 -123.125000 -123.000000
> FGD: W:-123.125000 E:-123.000000 N:49.000000 S:48.875000
> Map is viewable
> XTIFFOpen
> GTIFNew
> GTIFDirectoryInfo
> GTIFGetDefn
> PCS = 26710 (NAD27 / UTM zone 10N)
> Projection = 16010 (UTM zone 10N)
> Projection Method: CT_TransverseMercator
> ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
> ProjNatOriginLongGeoKey: -123.000000 (123d 0' 0.00"W)
> ProjScaleAtNatOriginGeoKey: 0.999600
> ProjFalseEastingGeoKey: 500000.000000 m
> ProjFalseNorthingGeoKey: 0.000000 m
> GCS: 4267/NAD27
> Datum: 6267/North American Datum 1927
> Ellipsoid: 7008/Clarke 1866 (6378206.40,6356583.80)
> Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
> Projection Linear Units: 9001/metre (1.000000m)
> Datum: 6267/North American Datum 1927
> ***** Attempting Datum Conversions
> Width: 5030
> Height: 6729
> Photometric: 3
> Rows Per Strip: 1
> Bits Per Sample: 8
> Samples Per Pixel: 1
> Planar Config: 1
> Loading geoTIFF map:
> /usr/local/share/xastir/maps/World/USA/WA/Tif/o48123h1.tif
>
> NW neat-line corner = -123.125000 49.000000
> 490856.892,5427244.298
> X/Y Pixels: 485.613243, 351.925915
> X/Y Pixels: 486, 352
> NE neat-line corner = -123.000000 49.000000
> 500000.000,5427236.771
> X/Y Pixels: 4235.247518, 355.012843
> X/Y Pixels: 4235, 355
> SW neat-line corner = -123.125000 48.875000
> 490834.034,5413348.801
> X/Y Pixels: 476.239248, 6050.538346
> X/Y Pixels: 476, 6051
> SE neat-line corner = -123.000000 48.875000
> 500000.000,5413341.270
> X/Y Pixels: 4235.247518, 6053.627114
> X/Y Pixels: 4235, 6054
> Crop points (pixels):
> Top: 355 Bottom: 6051 Left: 479 Right: 4232
>
> Initial Bytes Per Row: 5030
> xastir_left_x_increment: 0.000527 20474538 20474535 -0.001757 479 489 6051 355
> xastir_left_y_increment: 7.904444 14805072 14760072 1.000000 6048 355 6051 355
> xastir_right_x_increment: 0.000000 20519538 20519538 0.000000 4232 4232 6051 355
> xastir_right_y_increment: 7.904444 14805072 14760072 1.000000 6051 358 6051 355
> Left x increments: 0.000527 -0.001757
> Left y increments: 7.904444 1.000000
> Right x increments: 0.000000 0.000000
> Right y increments: 7.904444 1.000000
> xastir_top_y_increment: 0.000000 14760072 14760072 0.000801 358 355 4232 479
> xastir_bottom_y_increment: 0.000000 14805072 14805072 0.000799 6051 6048 4232 479
> PixelScale: 2.438400 2.438400 0.000000
> Size: x 3, y 2
> SkipRows: 1
> 6052 rows read in
> Load maps stop
>
>
> --
> Curt, WE7U archer at eskimo dot com
> Arlington, WA, USA http://www.eskimo.com/~archer
> "Lotto: A tax on people who are bad at math." -- unknown
> "Windows: Microsoft's tax on computer illiterates." -- WE7U
> "The world DOES revolve around me: I picked the coordinate system!"
> _______________________________________________
> Xastir mailing list
> Xastir at xastir.org
> https://krypton.hscs.virginia.edu/mailman/listinfo/xastir
--
------------------------------------------------------------------
Brian D Heaton
Senior Network Engineer
Janus Research Group
(706) 791-8342
GPG Fingerprint: C99E 3E9C E23A 4E47 46F4 0A77 3A45 CB65 9E19 5B0A
More information about the Xastir
mailing list