<div dir="ltr">Antii,<div><br></div><div>Thank you for the testing. It looks like Arc/GIS works while FME and Intergraph don't.</div><div><br></div><div>That file had both StdParallel1 and ScaleAtNatOrigin defined.</div>
<div><br></div><div>I've uploaded two new files. This one has only StdParallel1 defined with no ScaleAtNatOrigin.</div><div><a href="https://drive.google.com/file/d/0BybuTedE9CLxRnNlNUl2Zm5JRzg/edit?usp=sharing">https://drive.google.com/file/d/0BybuTedE9CLxRnNlNUl2Zm5JRzg/edit?usp=sharing</a><br>
</div><div><br></div><div>This file has the lat_ts stored in NatOriginLat and no ScaleAtNatOrigin defined.</div><div><a href="https://drive.google.com/file/d/0BybuTedE9CLxMzltYlBFWU51RzA/edit?usp=sharing">https://drive.google.com/file/d/0BybuTedE9CLxMzltYlBFWU51RzA/edit?usp=sharing</a><br>
</div><div><br></div><div>I think the latter might have better compatibility. Yet it seems to clearly be the incorrect way to do it.</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Nov 18, 2013 at 4:54 AM, Antti Castrén <span dir="ltr"><<a href="mailto:antti.castren@iki.fi" target="_blank">antti.castren@iki.fi</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Trent,<br>
<br>
The chart opens in ArcMap (ArcGIS 10.0) well, and it is in right<br>
location (Seattle).<br>
<br>
Relevant properties of the file as seen by ArcGIS:<br>
Cell Size: 240.003787, 240.003787<br>
Extent<br>
Top:4143530.20898<br>
Left: -9278434.53415<br>
Right: -9165872.75807<br>
Bottom: 3984167.69444<br>
Spatial Reference: Global Mercator<br>
Linear Unit: Meter<br>
false_easting: 0<br>
false_northing: 0<br>
central_meridian: 0<br>
standard_parallel_1: 47.667<br>
Datum: D_North_American_1983<br>
<br>
<br>
I opened the file in FME Data Inspector (FME Desktop 2013 SP1) also.<br>
It did not read the latitude of true scale from the<br>
ProjStdParallel1GeoKey as you can see from the following parameters.<br>
Therefore it places the chart to NE Georgia/NW South Carolina as if<br>
the latitude of true scale was the Equator.<br>
<br>
<FME Parameters><br>
Coordinate System Parameters<br>
CS_NAME: _FME_0<br>
DESC_NM: Global Mercator<br>
DT_NAME: NAD83<br>
PROJ: MRCAT<br>
UNIT: METER<br>
<br>
Datum Parameters<br>
DESC_NM: NAD 1983, Alaska, Canada,Continental US,Mexico,Central America<br>
ELLIPSOID: GRS1980<br>
SOURCE: US Defense Mapping Agency, TR-8350.2-B, August 1993<br>
USE: NAD83<br>
<br>
Ellipsoid Parameters<br>
DESC_NM: Geodetic Reference System of 1980<br>
E_RAD: 6378137<br>
P_RAD: 6356752.3141403478<br>
SOURCE: Stem, L.E., Jan 1989, State Plane Coordinate System of 1983<br>
<br>
OGC WKT Description<br>
PROJCS["Global Mercator",<br>
GEOGCS["NAD83",<br>
DATUM["North_American_Datum_1983",<br>
SPHEROID["Geodetic Reference System of<br>
1980",6378137,298.2572221008916,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6269"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4269"]],<br>
PROJECTION["Mercator_1SP"],<br>
PARAMETER["central_meridian",0],<br>
PARAMETER["scale_factor",1],<br>
<div class="im"> PARAMETER["false_easting",0],<br>
PARAMETER["false_northing",0],<br>
</div> UNIT["METER",1]]<br>
<br>
Esri WKT Description<br>
PROJCS["Global_Mercator",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["Geodetic_Reference_System_of_1980",6378137,298.2572221008916]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Mercator"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],PARAMETER["standard_parallel_1",0.0]]<br>
</FME Parameters><br>
<br>
Intergraph's Geomedia Professional recognises the file as GeoTIFF, but<br>
the chart appears to be in the same false place as with FME.<br>
<br>
If you want, I could update the file's georeferencing by FME and<br>
Geomedia so that location is correct. Then I could send listgeo<br>
information or even put the actual files available if desired.<br>
<br>
The goal is, of course, to get various software vendors to use the<br>
GeoTIFF-information the same way. Hopefully the common way would also<br>
be logically correct, but I guess that's too much to ask.<br>
<br>
Best regards,<br>
<br>
Antti<br>
<br>
<br>
Here is the listgeo info from the original file for those who are interested:<br>
<br>
Geotiff_Information:<br>
Version: 1<br>
Key_Revision: 1.0<br>
Tagged_Information:<br>
ModelTiepointTag (2,3):<br>
0 0 0<br>
-9278434.53 4143530.21 0<br>
ModelPixelScaleTag (1,3):<br>
240.003787 240.003787 0<br>
End_Of_Tags.<br>
Keyed_Information:<br>
GTModelTypeGeoKey (Short,1): ModelTypeProjected<br>
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea<br>
GTCitationGeoKey (Ascii,16): "Global Mercator"<br>
GeographicTypeGeoKey (Short,1): GCS_NAD83<br>
GeogCitationGeoKey (Ascii,6): "NAD83"<br>
GeogAngularUnitsGeoKey (Short,1): Angular_Degree<br>
GeogSemiMajorAxisGeoKey (Double,1): 6378137<br>
GeogInvFlatteningGeoKey (Double,1): 298.257222<br>
Unknown-2062 (Double,3): 0 0 0<br>
ProjectedCSTypeGeoKey (Short,1): User-Defined<br>
ProjectionGeoKey (Short,1): User-Defined<br>
ProjCoordTransGeoKey (Short,1): CT_Mercator<br>
ProjLinearUnitsGeoKey (Short,1): Linear_Meter<br>
ProjStdParallel1GeoKey (Double,1): 47.667<br>
ProjNatOriginLongGeoKey (Double,1): 0<br>
ProjNatOriginLatGeoKey (Double,1): 0<br>
ProjFalseEastingGeoKey (Double,1): 0<br>
ProjFalseNorthingGeoKey (Double,1): 0<br>
ProjScaleAtNatOriginGeoKey (Double,1): 1<br>
End_Of_Keys.<br>
End_Of_Geotiff.<br>
Projection Method: CT_Mercator<br>
ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)<br>
ProjNatOriginLongGeoKey: 0.000000 ( 0d 0' 0.00"E)<br>
ProjScaleAtNatOriginGeoKey: 1.000000<br>
ProjFalseEastingGeoKey: 0.000000 m<br>
ProjFalseNorthingGeoKey: 0.000000 m<br>
GCS: 4269/NAD83<br>
Datum: 6269/North American Datum 1983<br>
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)<br>
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)<br>
Projection Linear Units: 9001/metre (1.000000m)<br>
<div class="im">Corner Coordinates:<br>
Upper Left (-9278434.534, <a href="tel:4143530.209" value="+14143530209">4143530.209</a>)<br>
</div>Lower Left (-9278434.534, 3984167.694)<br>
Upper Right <a href="tel:%28-9165872.758" value="+19165872758">(-9165872.758</a>, <a href="tel:4143530.209" value="+14143530209">4143530.209</a>)<br>
Lower Right <a href="tel:%28-9165872.758" value="+19165872758">(-9165872.758</a>, 3984167.694)<br>
Center (-9222153.646, 4063848.952)<br>
<br>
<br>
2013/11/14 Trent Piepho <<a href="mailto:tpiepho@gmail.com">tpiepho@gmail.com</a>>:<br>
<div class="HOEnZb"><div class="h5">> It looks like arcgis does not support Mercator_1SP in geotiff files.<br>
> It just ignores the scale factor and uses 1.0. Which means in GDAL<br>
> and libgeotiff as they are now, there is no way to export Mercator<br>
> projections that don't have a scale of 1.0 to arcgis.<br>
><br>
> ArcGIS does support Mercator_2SP, but instead of placing the latitude<br>
> of true scale in ProjStdParallel1GeoKey, it uses<br>
> ProjNatOriginLatGeoKey. Which seems to be plainly wrong, yet that's<br>
> what it does.<br>
><br>
> I wonder if ESRI used libgeotiff? With unpatched libgeotiff, the only<br>
> way to get a latitude from a Mercator geotiff it to use one of the<br>
> origin latitude GeoKeys. It's not correct, but it does allow one to<br>
> stick the number in the file and get it back out without having to<br>
> modify or understand any libgeotiff code.<br>
><br>
> It's possible that ArcGIS supports getting the latitude of true scale<br>
> from ProjStdParallel1GeoKey, even it it doesn't store it there.<br>
> Antii, could you test this file with ArcGIS?<br>
><br>
> <a href="https://drive.google.com/file/d/0BybuTedE9CLxak9iakZRTC00cHM/edit?usp=sharing" target="_blank">https://drive.google.com/file/d/0BybuTedE9CLxak9iakZRTC00cHM/edit?usp=sharing</a><br>
><br>
> Pick file->download to get the geotiff file.<br>
><br>
> It should have projection information like this:<br>
> PROJECTION["Mercator_2SP"],<br>
> PARAMETER["standard_parallel_1",47.667],<br>
> PARAMETER["central_meridian",0],<br>
> PARAMETER["false_easting",0],<br>
> PARAMETER["false_northing",0],<br>
> UNIT["metre",1,<br>
> AUTHORITY["EPSG","9001"]]]<br>
> Origin = (-9278434.534151820465922,4143530.208984711207449)<br>
> Pixel Size = (240.003786961970206,-240.003786961970206)<br>
> Corner Coordinates:<br>
> Upper Left (-9278434.534, 4143530.209) (123d32'26.22"W, 48d23'56.59"N)<br>
> Lower Left (-9278434.534, 3984167.694) (123d32'26.22"W, 46d57'58.64"N)<br>
> Upper Right (-9165872.758, 4143530.209) (122d 2'30.76"W, 48d23'56.59"N)<br>
> Lower Right (-9165872.758, 3984167.694) (122d 2'30.76"W, 46d57'58.64"N)<br>
> Center (-9222153.646, 4063848.952) (122d47'28.49"W, 47d41'15.43"N)<br>
</div></div></blockquote></div><br></div>