[gdal-dev] Coordinate Transformation Question

[Joel Odom]

Hi, Frank.  Sometimes I have to admit to idiocy.  Being new to OGR, I
didn't realize that the projection code was dependent on the proj
library.  Maybe I should have looked at my output where it said,
"ERROR 6: Unable to load PROJ.4 library (proj.dll), creation of
OGRCoordinateTransformation failed."

At any rate, it's working beautifully now that I've got the dll
installed.  Again, thanks for your help.





On Tue, Dec 16, 2008 at 3:32 PM, Frank Warmerdam <warmerdam at pobox.com> wrote:
> Joel Odom wrote:
>>
>> Thanks for checking it out, Frank.  Here is my full projection string:
>>
>> PROJCS[" Projection Name = Lambert Conformal Conic Units = meters
>> GeoTIFF Units =
>> meters",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",38.66666666666667],PARAMETER["standard_parallel_2",33.33333333333336],PARAMETER["latitude_of_origin",34.10833333333333],PARAMETER["central_meridian",-84.50000000000003],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
>
> Joel,
>
> This works fine for me with the testepsg utility:
>
> WKT[test.wkt] =
> PROJCS[" Projection Name = Lambert Conformal Conic Units = meters GeoTIFF
> Units =meters",
>    GEOGCS["NAD83",
>        DATUM["North_American_Datum_1983",
>            SPHEROID["GRS1980",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",38.66666666666667],
>    PARAMETER["standard_parallel_2",33.33333333333336],
>    PARAMETER["latitude_of_origin",34.10833333333333],
>    PARAMETER["central_meridian",-84.50000000000003],
>    PARAMETER["false_easting",0],
>    PARAMETER["false_northing",0],
>    UNIT["metre",1,
>        AUTHORITY["EPSG","9001"]]]
>
> warmerda at gdal64[75]% testepsg -t test.wkt WGS84 100000 200000
> OGRCT: Source: +proj=lcc +lat_1=38.66666666666667 +lat_2=33.33333333333336
> +lat_0=34.10833333333333 +lon_0=-84.50000000000003 +x_0=0 +y_0=0
> +ellps=GRS80 +datum=NAD83 +units=m +no_defs
> OGRCT: Target: +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> (100000.000000,200000.000000,0.000000) -> (-83.390980,35.907560,0.000000)
>
> Actually, I neglected to ask before if you were confident you had
> PROJ.4 properly installed and findable by OGRCoordinateTransformation.
> If coordinate transformation *never* works then it is likely the PROJ
> DLL/sharedlib is not being found at runtime.
>
> Best regards,
> --
> ---------------------------------------+--------------------------------------
> I set the clouds in motion - turn up   | Frank Warmerdam,
> warmerdam at pobox.com
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush    | Geospatial Programmer for Rent
>
>



-- 
http://giscoder.blogspot.com/