[Proj] rHEALPix new parameters
spethm at landcareresearch.co.nz
Thu Nov 3 14:21:47 EST 2011
I am in the process of switching my test scripts to use the cs2cs
instead of proj. I am having trouble running the inverse tests.
This command works for proj and produces the inverse projection of HEALPix.
proj +proj=healpix -f '%.'$PERCISION'f' -I +lon_0=$1 +a=$2 +R_A
+ellps=WGS84 -E <<EOF
Can you help me convert this to cs2cs?
On 02/11/11 13:13, Frank Warmerdam wrote:
> On Tue, Nov 1, 2011 at 3:17 PM, Michael Speth
> <spethm at landcareresearch.co.nz> wrote:
>> An example best explains my question.
>> Command line
>> proj +proj=healpix -f '%.'$PERCISION'f' -I +lon_0=$1 +a=$2 +R_A +ellps=WGS84
>> $3 $4
>> ==== HEALPix Spheriod Inverse Test ====
>>  0 0 0.0000 0.0000
>>  0 0 45.0000 0.0000
>>  0 0 -45.0000 0.0000
>>  2 -2 * *
>>  2 -2 * *
>>  1 2 * *
>>  2 1 * *
>> So in the output, notice that the last 4 test result in an error (noted by
>> the *). The output I am looking for is "invalid x or y" for the last 4 test
>> results. Does this make sense?
> It does not seem to have been typical in the past to issue
> errors during point transformation - instead HUGE_VAL was
> returned and that was it. So the "proj" command has no code
> to report errors other than by printing "*" after a point transformation.
> cs2cs on the other hand does - partly motifivated by my issuing of
> errors from the datum shifting code that comes into play in cs2cs
> and not in proj. I would suggest you try with cs2cs.
> Best regards,
More information about the Proj