Craig,
thanks, but that did not help. The code that really bugs me is this one:
$res=$self*$other #
$res->sethdr($self->hdr_copy);
say "keys self: ",%{$self->hdr};
say "keys other: ",%{$other->hdr};
say "keys res: ",%{$res->hdr};
When I use PDL Transform to generate $other, I get this:
keys self:
typusphaseorientationTraselfImage1max3.13852469201402xHASH(0x52257c8)
channelHASH(0x4c1db68)SetHASH(0x51e19c8)ascconvHASH(0x5213078)
yHASH(0x52130c0)slzroxndims7echoHASH(0x5176930)tHASH(0x5203338)linksHASH(0x4bc0e20)parentroi1min-3.14159265358979active1dimnamesARRAY(0x4bcff00)updating1dicomHASH(0x4c9bb88)peyzHASH(0x4bca8b8)3d1
keys other:
yHASH(0x4b6dd80)dimnamesARRAY(0x4bcfa98)ndims3xHASH(0x4ae3df8)zHASH(0x4ca63f0)
keys res:
yHASH(0x491b420)dimnamesARRAY(0x51ed588)ndims3xHASH(0x51b25d8)zHASH(0x4994820)
You see, the header is copied from $other, not from $self, but only if
$other was the result of map(t_linear ...), see below.
It works as expected otherwise.
That's why I wanted a "clean" piddle.
Ingo
On 03/17/2014 08:38 PM, Craig DeForest wrote:
> Ah. It's putting a scientific coordinate system in there in FITS
> form. You want to use the "pix=>1" option to map to work in pixel
> coordinates. Otherwise the output image is autoscaled to fit most of
> the output into the output grid, and the mapping between output
> scientific coordinates and output pixel coordinates is stored in WCS
> FITS-header form.
>
> Try:
>
> $u = $ts->map($x,{pix=>1});
>
>
> On Mar 17, 2014, at 1:35 PM, Ingo Schmid <[email protected]
> <mailto:[email protected]>> wrote:
>
>> $a=zeroes(10,10)
>> $ts=t_linear(s=>3)
>> $u=$ts->map($x)
>> p %{$u->gethdr}
>> CRVAL2 13.5 CDELT1 3 CRPIX1 5.5 SIMPLE T CRVAL1 13.5 HISTORY Header
>> written by PDL::Transform::Cartography::map NAXIS 2 CRPIX2 5.5
>> CUNIT1 NAXIS1 10 CDELT2 3 CTYPE1 CUNIT2 CTYPE2 NAXIS2 10
>>
>> yes, they are piddles, but you see?
>>
>> My real transform is like that (simplified, so it may have some typos):
>> $maskb=$maska->copy->map(t_linear(s=>1,d=>3),[$sizex,$sizey,$sizez],{m=>'g',
>> ir=>$a,},);
>> $maskb->sethdr({});
>>
>> Ingo
>>
>> On 03/17/2014 08:22 PM, Craig DeForest wrote:
>>> What method are you using to map? The map() and apply() methods
>>> both return PDLs...
>>>
>>> On Mar 17, 2014, at 1:21 PM, Ingo Schmid <[email protected]
>>> <mailto:[email protected]>> wrote:
>>>
>>>> Hi,
>>>>
>>>> I use PDL::Transform to map piddles from one coordinate system to
>>>> another. My problem is, that after the transform, the result is a
>>>> PDL::Transform object, which seems to interfere with my own module.
>>>>
>>>> How do I get a naked piddle back after a map? I tried
>>>> $res->sethdr({})
>>>> which clears the header alright, but that doesn't cut it for me
>>>> somehow. I need fullaccess to the piddle header, my guess is, that
>>>> after the hard reset of the header, something bad happens.
>>>>
>>>> Ingo
>>>
>>
>
_______________________________________________
Perldl mailing list
[email protected]
http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
