Fwd: Re: [Proj] PROJ 4.5.0 beta (and Krovak!)

Frank Warmerdam warmerdam at pobox.com
Wed Oct 18 13:04:32 EDT 2006


Gerald I. Evenden wrote:
 >Frank Warmerdam wrote:
>> warmerda at gdal2200[1926]% ./proj +proj=krovak +ellps=bessel
>> 10 45
>> -1164802.15     -1546506.69
>>
>> I don't really know any test points for Krovak or even what options
>> are required in the two implementations.
>>
>> For now, I'll stick with the existing Krovak code for the 4.5.0 release,
>> but if there are things we can do to make the two implementations
>> compatible with the aim of a seamless transition to use of libproj4
>> in the future, I'd be all for that!

Gerald I. Evenden wrote:
> I inserted the 10 45 number into an old script:
> 
> #!/bin/sh
> # This is a test script to check the correct execution
> # of the Krovak projection
> cat <<EOF >temp
> 35d42'35.2147 48d7'46.2973
> 10 45
> EOF
> echo "forward projection\n"
> ../lproj +proj=kocc +ellps=bessel +czech temp \
>      +lon_0=42d30 +lat_0=49d30 \
>      +lat_t=59d42\'42.69689 +lat_1=78d30 +k_0=.9999
> ../lproj +proj=kocc +ellps=bessel +czech -o temp \
>      +lon_0=42d30 +lat_0=49d30 \
>      +lat_t=59d42\'42.69689 +lat_1=78d30 +k_0=.9999 \
> 
> | ../lproj -I +proj=kocc +ellps=bessel +czech -i \
> 
>      +lon_0=42d30 +lat_0=49d30 \
>      +lat_t=59d42\'42.69689 +lat_1=78d30 +k_0=.9999 -W6
> #
> echo "above should agree with input of\n"
> cat temp
> rm temp
> 
> which produces the results:
> 
> forward projection\n
> -16643462.09	2422822.35
> -15043707.84	7571503.40
> 35d42'35.214700"E	48d07'46.297300"N
> 10d00'00.000000"E	45d00'00.000000"N
> above should agree with input of\n
> 35d42'35.2147 48d7'46.2973
> 10 45

Gerald,

OK, this leaves me with the impression that the PJ_krovak.c that was
contributed to me includes built in defaults that are compatible with
common use in the Czech Republic while, perhaps, yours (kocc) uses standard
defaults (ie. lon_0=0, lat_0=0, lat_t=0, lat_1=0?)

Do you have a suggestion on how we can align the two implementations so they
work the same from the commandline?

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    | President OSGeo, http://osgeo.org



More information about the Proj mailing list